虹膜识别特征提取及鉴别

1、-PAGE . z摘 要随着信息社会的快速开展，对平安的需求也日益增长。虹膜识别技术作为一种身份识别，以其很高的可靠性得到人们的重视。虹膜识别系统核心一般由图像采集、虹膜定位、归一化、特征提取及编码和训练识别五局部构成。本文介绍了目前虹膜识别的现状，简单阐述了一些经典的虹膜识别算法和技术，并完成识别系统。在虹膜的定位阶段，首先对图像进展缩放，在不影响后续处理的情况下减小了处理的数据量，然后采用梯度加权的Canny算法进展边缘检测，再对边缘图像，采用圆Hough定位方法，分别定位了虹膜的外边界。接着采用Radon变换检测直线的方法分割上下眼睑，阈值法除去睫毛干扰。同时也研究了一些文献中分割眼睑和

2、睫毛的方法。归一化阶段，采用了文献中普遍使用的Rubber-Sheet模型，将虹膜归一化为大小的矩形，以利于特征比对。在虹膜的特征提取及编码阶段，基于信号处理中的空间/频域技术，采用一维Log Gabor滤波器提取虹膜的纹理信息，对滤波结果的实部和虚局部别进展相位量化和编码，同时也对噪声进展处理，获得相应的掩码。训练识别阶段，采用海明距离度量虹膜之间的相似度，选取最小距离分类器和具有最小错误率的分类阈值形成组合的分类决策规则。整个识别系统主要在中科院V3.0虹膜数据库上进展了测试。关键词：虹膜识别；虹膜定位；圆Hough变换；Log Gabor小波-. zABSTRACTWith the ra

3、pid growth of information technology, the demands of information security areever-growing. As the technology of identification, iris recognition, for its high reliability, gets great attention. Iris recognition system consists of image capturing, iris location, iris normalization, feature e*traction

4、 and coding and decision training. In this dissertation, the situation of iris recognition is presented. Some practical algorithms and technique are briefly introduced. A system of iris recognition is fulfilled.In iris location stage, image is zoomed to reduce the data volume with less influence on

5、post processing. Afterword, the algorithm of Canny edge detection, with weighted gradient, is adopted. In the edge image, Circular Hough transform is applied to locate the inside and out boundary of iris. Then, linear Radon transform is put to use to detect the eyelids. Eyelashes are eliminated by t

6、hreshold. In normalization stage, rubber-sheet model, in general use, is usedt-o unwrap iris image into a rectangle of the same size, for the parison of characteristics.In feature e*traction and coding stage. 1D Log Gabor filters are used to filter the iris te*ture features in the space-frequency do

7、main. Then, the real part and the image part is quantized by phase encoding respectively. Besides, the processing of noise to generates mask.In decision training stage, the Hamming distance is employed to measure the similarity of two iris. The minimum distance classifier bining with the threshold w

8、ith minimum error ratio are used as the decision rules. The e*periments implemented on CASIA iris database V3.0 show that the system performs well.Key Words: Iris Recognition; Iris Location; Circular Hough Transform; Log Gabor Wavelets-. z目录 TOC o 1-3 h z u HYPERLINK l _Toc263860304第一章绪论1HYPERLINK l

9、 _Toc2638603051.1 生物特征识别1HYPERLINK l _Toc263860306非生物特征身份鉴别方法1HYPERLINK l _Toc2638603071.1.2 生物特征识别技术的提出2HYPERLINK l _Toc2638603081.1.3 研究生物特征识别技术的意义3HYPERLINK l _Toc2638603091.2 虹膜识别研究3HYPERLINK l _Toc2638603101.2.1 虹膜识别的机理3HYPERLINK l _Toc2638603111.2.2 虹膜识别的开展与现状6HYPERLINK l _Toc2638603121.2.3 虹膜

10、库9HYPERLINK l _Toc2638603131.3 本文研究的容9HYPERLINK l _Toc263860314第二章虹膜图像预处理11HYPERLINK l _Toc2638603152.1 虹膜外边界的定位11HYPERLINK l _Toc2638603162.1.1 Canny算子和Hough变换的根本原理11HYPERLINK l _Toc2638603172.1.2 本文采用的Canny算子14HYPERLINK l _Toc2638603182.1.3 本文采用的Hough变换圆检测算法16HYPERLINK l _Toc2638603192.1.4 虹膜外边界的定

11、位17HYPERLINK l _Toc2638603202.2 眼睑的分割18HYPERLINK l _Toc2638603212.2.1 抛物线检测眼睑19HYPERLINK l _Toc2638603222.2.2 Radon变换法分割眼睑19HYPERLINK l _Toc2638603232.3 剔除睫毛21HYPERLINK l _Toc2638603242.3.1 阈值法剔除睫毛21HYPERLINK l _Toc2638603252.4 小结22HYPERLINK l _Toc263860326第三章虹膜图像归一化及特征提取23HYPERLINK l _Toc2638603273

12、.1 虹膜纹理归一化23HYPERLINK l _Toc2638603283.1.1 平移23HYPERLINK l _Toc2638603293.1.2 旋转24HYPERLINK l _Toc2638603303.1.3 伸缩24HYPERLINK l _Toc2638603313.1.4 分辨率24HYPERLINK l _Toc2638603323.2 虹膜的Rubber-Sheet模型25HYPERLINK l _Toc2638603333.2.1 虹膜边界表示26HYPERLINK l _Toc2638603343.2.2 虹膜区域的表示27HYPERLINK l _Toc2638

13、603353.2.3 虹膜区域规化28HYPERLINK l _Toc2638603363.2.4 灰度级插值29HYPERLINK l _Toc2638603373.2.5 具体实现31HYPERLINK l _Toc2638603383.3 小结31HYPERLINK l _Toc263860339第四章虹膜特征编码32HYPERLINK l _Toc2638603404.1 图像纹理分析的根本理论32HYPERLINK l _Toc2638603414.2 虹膜特征编码方案33HYPERLINK l _Toc2638603424.3 基于Log Gabor小波的特征编码35HYPERLI

14、NK l _Toc2638603434.3.1 Log Gabor函数性能分析35HYPERLINK l _Toc2638603444.3.2 Log Gabor小波的构造39HYPERLINK l _Toc2638603454.3.3 特征编码40HYPERLINK l _Toc2638603464.4 小结42HYPERLINK l _Toc263860347第五章虹膜模式匹配43HYPERLINK l _Toc2638603485.1 识别判决43HYPERLINK l _Toc2638603495.2 虹膜模式匹配44HYPERLINK l _Toc2638603505.2.1 海明距

15、离匹配44HYPERLINK l _Toc2638603515.2.2 最小距离分类器46HYPERLINK l _Toc2638603525.3 统计分析49HYPERLINK l _Toc2638603535.4 小结52HYPERLINK l _Toc263860354第六章总结53HYPERLINK l _Toc2638603556.1论文容总结与展望53HYPERLINK l _Toc2638603566.2系统演示54HYPERLINK l _Toc263860357参考文献61HYPERLINK l _Toc263860358致63-. z第一章 绪论1.1生物特征识别生物特征识

16、别有时也被称为生物测定技术。生物特征识别是指利用人体所固有的生理特征或行为特征，进展个人身份识别的技术19,21,22。近年来在维护国家平安、航空平安、金融平安、社会平安、网络平安等应用领域，身份识别和认证变得越来越重要，对身份的有效认证要求更准确、更平安、更实用的鉴别方法，生物技术的开展和进步为身份鉴别提供了新的方法和手段，基于生物特征的身份识别方法成了近年身份识别领域研究的热点。非生物特征身份鉴别方法个人身份鉴别可以分为：认证和识别。认证是指验证用户是否是他所声明的身份，识别指的是确定用户的身份。目前，身份鉴别大多采用、钥匙、密码、用户名等等，通过验证这些标识身份的身外之物来识别个人身份，

17、这些身份鉴别方法是把身份识别问题转化为鉴别一些标识个人身份的事物。这些均为身外之物，而不是生物特征，这类方法统称为非生物特征身份鉴别方法。这些身份鉴别方法确实在一定程度上提供了简单有效的身份鉴别，给人们的生活带来了方便，但是随着信息平安性和身份鉴别可靠性要求的提高，这些非生物特征的身份鉴别方法很难满足要求。采用钥匙、证件的门禁系统和采用用户名和密码的网上交易，其平安性容易受到攻击。对于日益增加的计算机网络用户，登陆不同的，可能需要注册不同的用户名，相当多的用户为了便于记忆而设定的用户名和密码容易被猜出和破译，而非生物特征的身份鉴别方法很难识别真正的拥有者和取得身份标识物品的冒充者，所以一旦他人

18、获得了这些身份标识事物，则就拥有了相应的权利，造成的后果会很严重的。防伪性和防欺骗性差是非生物特征身份鉴别方法平安性低的另一个缺点。犯罪分子伪造证件的手段越来越高明，这使目前广泛使用的依靠证件、个人识别码、口令或钥匙等来确认个人身份的技术平安性降低。生物特征识别技术的提出二十一世纪是网络化、信息化时代，随之而来的一大特征就是身份的数字化和隐性化，如何准确鉴别一个人的身份、保护信息的平安是当今信息化时代必须解决的问题。随着网络的日益普及，社会信息化程度越来越高，对平安、可靠的身份识别技术需求变得越来越迫切。人们希望有一种更平安、更可靠、携带使用更方便且不会被遗忘的事物来表示个人身份，显然个体本身

19、的生物特征将是一种可行的、比拟理想的选择。生物特征识别技术利用人类的生理或者行为特征进展身份识别和认证，认证的是人身之物，而不是身外之物。人们可能会遗忘或丧失标识他们身份的证件物品或用户名及密码，但是人们绝不会遗忘或者丧失他们的生物特征(如人脸、指纹、虹膜、掌纹等)，另外，个人的生物特征也不会被分享，所以生物识别系统很难被欺骗或欺骗本钱很高。目前，生物特征识别采用的生理特征主要有：指纹、掌纹、眼睛(视网膜和虹膜)、人体气味、脸型、皮肤毛孔、手腕/手的血管纹理和DNA等，是先天具有的，采用的行为特征主要包括：签名、声音、行走的步态、击打键盘的力度等，是后天习惯养成的。基于生物特征的识别利用计算机

20、技术很容易实现身份自动识别，它们的根本工作原理一样，如图1-1所示。首先是从独立个体采集生物样本，这些样本可以是虹膜图像、指纹图像、人脸图像，声音的数字化描述，步态时序图像等；接着是进展预处理，主要进展特征区域定位或者去噪处理；然后进展特征提取，并将提取的特征与数据库存储的身份特征进展比对，最后输出比对结果，做出身份判断。在基于生物特征的身份认证领域，身份信息全部是以数字形式存储于数据库或者智能IC卡中，鉴别身份时，能够对持有者合法性进展验证。图1-1 生物识别系统原理理想的生物特征识别系统应满足：(1)所有人都拥有这一生物特征，并且不同人的生物特征是可以区分的；(2)生物特征的采集不随采集的

21、条件而不同；(3)系统能够区分冒充者。近年来，随着计算机技术和信息处理与识别技术的不断进步，生物识别技术得到了迅猛的开展，并逐渐被群众所认可。研究生物特征识别技术的意义利用人类个体生理和行为特征进展个人身份识别己经取得了许多可喜成果。目前，国外许多高新技术公司用眼睛虹膜、指纹、面貌特征等取代人们手中的信用卡或密码，并且已经开场在机场、银行和各种电子器具上进展了实际应用。1.2虹膜识别研究虹膜识别的机理1.虹膜构造虹膜识别是最具潜力的生物识别方法之一，是识别率高、非接触、仿欺骗性好的识别方法。虹膜是人眼受保护的部组织，位于角膜和水样液后面、晶状体前面，如图1-2。虹膜是外部可见的，环绕瞳孔的有色

22、圆环，是一个肌肉组织，虹膜直径约12mm，厚约0.5mm，根部最薄。虹膜外表上下不平坦，有皱缓和凹陷，凹陷又称隐窝。由于虹膜血管分布不匀，使虹膜外表出现许多的放射形条纹。这其中包含的许多互相交织的类似于斑点、细丝、冠状、条纹、隐窝等的细微特征，就构成了我们所说的虹膜纹理信息，虹膜识别就是利用虹膜组织上这些丰富的纹理信息，作为重要的身份识别特征。图1-2 虹膜构造2.虹膜具有的优良特性虹膜模式具有作为身份特征的许多优良特性：(1)普遍性虹膜是每个人都具有的。(2)唯一性每个人错综复杂的虹膜是独一无二的，虹膜特征的唯一性与它的形成过程有关。虽然虹膜颜色由遗传决定，但是虹膜的纹理特征，就像指纹一样，

23、是随机形成的。John Daugman博士和剑桥大学同事使用Gabor算法3，将人眼虹膜区域的纹理转换成二进制代码，在公共场所对虹膜识别系统作了大量的虹膜图像比拟试验，超过两百万的眼睛图像分析，提供了平安系统的数学支持在伦敦的希思罗机场测试，结果说明眼睛可以提供个人身份信息，而且虹膜识别几乎是零误差样本。 John Daugman估计，两幅虹膜图像得到的代码完全一样的几率几乎为零，没有形状完全一样的虹膜，即使是双胞胎或者同一个人的左、右眼得到的虹膜代码都是毫不相关的，这正是虹膜识别的引人之处。(3)稳定性发育生物学家通过大量观察发现，当虹膜发育完全以后，它在人的一生中是稳定不变的，因而具有稳定

24、性。(4)非入侵检测和指纹识别不同，虹膜识别不需要物理接触，可采用非接触进展虹膜图像采集。(5)可承受程度较好虹膜识别以其认证准确度高、速度快、平安性高，被用户所承受。在识别过程中，用户不会有任何不舒服和不安的感觉，只需要在设备前停留片刻，无需为排长队等候而感到厌烦。3.虹膜识别机理生物特征识别通过捕获模板样本，然后采用数学函数把样本转化成生物学模板，而且这模板应该能够提供标准化、有效的、高度有区别的特征表示，这样可以客观地和其它模板进展比拟进而确定身份。虹膜识别系统的构成及其研究容大致包括以下五个方面，如图1-3所示。(1)虹膜图像采集：虹膜图像采集直接影响到后续的特征提取和识别的效果。在大

25、多数实际情况中，由于受到光照、噪音、饰物、遮挡等因素的干扰，采集到高质量、可承受的虹膜图像是一个比拟困难的问题。(2)虹膜检测与定位：从各种虹膜图像(包括不同的光照、背景、大小、分辨率等)中检测出虹膜的存在，并确定其准确位置及有效区域的过程。(3)虹膜纹理归一化：虹膜外圆的半径以及圆心的坐标后，将虹膜纹理区域归一化为统一大小和分辨率，减少虹膜在图像中的位置、整体尺度、瞳孔缩放尺度对识别产生的影响。(4)特征提取与编码：目的在于提取虹膜纹理中可用于识别的有效信息，并构造虹膜纹理特征编码。(5)匹配与分类：将虹膜纹理特征编码与数据库中的虹膜纹理特征编码进展匹配，进展识别。图1-3 虹膜识别系统的构

26、成虹膜识别的开展与现状用虹膜进展身份识别的设想最早出现于19世纪80年代，近20年，该项技术有了飞跃的开展。1885年在巴黎的监狱中曾利用虹膜的构造和颜色区分同一监狱中的不同犯人，这是最早利用虹膜进展的身份识别，真正的自动虹膜识别系统则是上世纪末才出现。1987年，眼科专家Aran Safir和Leonard Florm首次提出了利用虹膜图像进展自动身份识别的概念。发育生物学家通过大量观察发现，当虹膜发育完全后，它在人的一生中是稳定不变的，因而具有稳定性，可用来作为身份特征。影响力较大的虹膜识别系统主要有Daugman系统、Wildes系统、Boles系统和中科院虹膜系统等3,1,6, 11,

27、13,16。(1)Daugman系统目前国际上很多虹膜识别产品都使用了英国剑桥大学的Daugman博士提出的虹膜识别算法。算法中利用积分微分算子(Integro-differential operator)检测虹膜的外圆边界：(1.1)其中代表虹膜图像在处的灰度值：表示卷积；是标准差为的高斯算子，起平滑滤波的作用；是虹膜外边缘的参数(半径及圆心)。算子在以圆心，半径为的圆周上，对像素灰度值做积分并把它归一化，再求差分的极大值，从而得到圆的参数。然后，将虹膜区域视为各向同性的弹性体进展归一化，这就是Rubber-Sheet模型，其结果是将环形的虹膜纹理区域均匀拉伸到统一大小的矩形区域中。另外，D

28、augman设计了二维Gabor滤波器对虹膜纹理进展特征提取，滤波器表示如下：(1.2)滤波后提取相位信息，并将所得到的相位信息量化为二值的虹膜编码，共计2048比特。利用两幅虹膜图像所得到的二值编码间的归一化海明距离作为相似形度量对虹膜进展比对识别。另外此系统也采用了一种比拟简单、直观的方法解决了虹膜纹理的旋转性问题。原虹膜纹理的旋转对应着归一化后的左右平移，所以Daugman采用将虹膜编码左右平移较少的位数后分别比对，找出其中最小的海明距离作为两个虹膜之间的相似性度量。(2)Wildes系统Wildes系统与Daugman系统进展比照，采用了不同的识别过程。首先，利用扩展的Hough变换进

29、展虹膜外圆的检测。在克制虹膜的平移、缩放和旋转问题上，采用了图像配准的方法，而没有将其归一化到统一的矩形区域。Wildes认为仅采用2048比特的编码来表示虹膜的纹理特征可能包含的信息量过少，所以其系统中的数据量较大。他提出利用各向同向的高斯-拉普拉斯滤波器：(1.3)在不同分辨率下对图像进展滤波之后，逐次进展1/2采样，得到不同尺度的数据构成4层金字塔构造，利用这些数据作为虹膜识别的特征。最后使用Fisher分类器进展分类。(3)Boles系统Boles利用小波变换的过零点和两个连续过零点之间的小波变换的积分平均值来表示虹膜特征。在对虹膜纹理图像编码前，先沿着以虹膜中心为圆心的同心圆对虹膜图

30、像采样，把两维的虹膜图像变为一维的信号，然后利用特定的小波函数对它进展变换。这里特定的小波函数定义为*一光滑函数的两阶导数，即定义小波母函数为：(1.4)其中为*一光滑函数。根据小波变换的定义：(1.5)的小波变换正比于经过函数光滑化的的二阶导数。小波变换的零穿插对应于的变形点，即函数曲线剧烈变化的局部。在虹膜识别的应用中，表示虹膜图像样本，则二元点序列可以作为虹膜特征的编码。随后，通过其自定义的相异度函数完成分类。该算法能够抵抗光照变化，但只在很小规模的数据库上进展过测试。(4)中科院虹膜系统中国科学院的马力、谭铁牛等在他们的系统中不但提出了利用傅立叶变换来对虹膜图像进展质量评估，并且在Da

31、ugman利用Gabor滤波器进展特征提取的根底上，将Gabor滤波器的调制函数进展改良，使其具有圆对称的特性，表达式如下：(1.6)(1.7)然后将滤波后的图像分块，提取每块的均值和方差形成虹膜的特征向量。最后利用Fisher线形判据降低特征向量的维数，根据最近邻分类器进展特征匹配。与其它两种经典方法Daugman和Boles的方法相比，此方法也能取得较好的识别效果。虹膜库随着虹膜识别算法研究的开展，广阔研究者对用于测试算法性能的虹膜图像库的需求也越来越高。在此领域上应用最广泛的虹膜图像库是中国科学院自动化研究所公开的CASIA虹膜图像库17。另外还有样本质量较复杂的Ubiris数据库；只包

32、括清晰虹膜纹理的UPOL虹膜图像库；其它的还有MMU虹膜图像库、Bath虹膜图像库、ICE虹膜图像库以及WVU虹膜图像库等等。各个图像库都各自的特点。1.3本文研究的容虹膜识别包括虹膜图像采集、虹膜检测与定位、虹膜纹理归一化、特征提取与编码和匹配与分类。本文将完成一套虹膜识别算法实现，并在此根底上完成虹膜识别系统。论文分六章，容概述如下：(1)第一章：绪论。全面概述了生物识别技术研究背景及其意义，介绍了虹膜识别机理、开展现状和比拟有影响力的虹膜库。(2)第二章：虹膜图像预处理。采用Canny算子和Hough变换对虹膜的外边界进展定位，进展干扰检测，Radon直线检测除去眼睑和阈值法除去睫毛干扰

33、。(3)第三章：虹膜图像归一化及特征提取。针对虹膜区域图像的大小、分辨率受拍摄距离的影响，研究虹膜图像的归一化，采用Rubber-Sheet模型完成规化。(4)第四章：虹膜特征编码。简要介绍了几个重要系统所采用的编码方案，讨论比拟并实现了Log Gabor小波的特征编码。(5)第五章：虹膜模式匹配。简单分析了虹膜模式匹配的方法，计算虹膜编码的海明距离，基于最小距离分类器确定虹膜判决阈值。(6)第六章：总结。对全文进展系统总结展望，并给出系统使用流程。-. z第二章 虹膜图像预处理在实际获取的图像中，虹膜只占一局部区域。图像中还包括瞳孔、眼睑、睫毛及眼睛和人脸的其它局部。所以，在对虹膜纹理特征进

34、展提取之前，要把虹膜区域从图像中分割出来。其分割准确与否直接影响到后面的特征提取及匹配精度，即如果分割不准确，则虹膜识别效果就差，很可能造成误判，使识别算法失效。针对Canny2算法优良的边缘检测性能及虹膜边界的特点，本文采用了该算法，与Hough变换相结合，用于定位虹膜外边界、分割上下眼睑，采用阀值法剔除睫毛。2.1虹膜外边界的定位虹膜的边界即边缘与瞳孔相交接的局部，外边界即外边缘与巩膜相交接的局部，均近似为圆形，因此定位虹膜外边界就是要确定这两个拟合圆的圆心和半径。Canny算子和Hough变换的根本原理(1)Canny边缘检测算子边缘是图像的根本特征，保存了原始图像中相当重要的信息，而又

35、使得总的数据减小了很多。准确提取边缘既能减少图像处理的信息量，又可重点描述物体的形态特征。传统的边缘检测算子如Roberts、Sobel、Prewitt、Kirsch和Laplacian等算子都是局域窗口梯度算子，由于它们对噪声敏感，所以在处理实际图像中效果并不理想。1986年，Canny提出边缘检测算子应满足以下3个判断准则：信噪比准则；定位精度准则；单边缘响应准则，并由此推导出了最正确边缘检测算子-Canny算子。该算子具有定位精度高、单一边缘和检测效果好等优点。它主要分以下四步：eq oac(,1)平滑图像Canny边缘检测算子是高斯函数的一阶导数，是对信噪比与定位精度之乘积的最优化逼近

36、算子。选择适宜的高斯低通滤波器来平滑原始的虹膜图像，既可以取得更加柔和的平滑效果，又不至于破坏图像中的边沿，同时还可以消除原始图像的高频噪声。Canny算法首先用二维高斯函数的一阶导数，对图像进展平滑，设二维高斯函数为：(2.1)其梯度矢量为：(2.2)把的2个滤波卷积模板分解为2个一维的行列滤波器：(2.3)(2.4)其中，为常数，为高斯滤波器参数，它控制着平滑程度。小的滤波器，虽然定位精度高，但信噪比低；大的情况则相反，因此要根据需要适当地选取高斯滤波器参数。eq oac(,2)计算梯度的幅值和方向传统Canny算法采用邻域一阶偏导的有限差分来计算平滑后数据阵列的梯度幅值和梯度方向。其中，

37、和方向偏导数的2个阵列和分别为：(2.5)(2.6)像素的梯度幅值和梯度方向用直角坐标到极坐标的坐标转化公式来计算，用二阶数来计算梯度幅值和方向，分别为：(2.7)(2.8)eq oac(,3)对梯度幅值进展非极大值抑制为了准确定位边缘，必须细化梯度幅值图像中的屋脊带，只保存幅值局部变化最大的点，这一过程就是非极大值抑制。如图2-1。在非极大值抑制过程中，Canny算法使用大小、包含8个方向的邻域对梯度幅值阵列的所有像素沿梯度方向进展梯度幅值的插值。在每一个点上，邻域的中心像素与沿梯度方向的2个梯度幅值的插值结果进展比拟，是像素邻域中心处沿着梯度方向的扇形区域，非极大值抑制在此区域进展。如果邻

38、域中心点的幅值不比梯度方向上的2个插值结果大，则将对应的边缘标志位赋值为0，这一过程把宽屋脊带细化为一个像素宽，并且保存了屋脊的梯度幅值。非极大值抑制(Non-Ma*ima Suppression,NMS)过程的数学表示为：(2.9)图2-1 Canny算子非极大值抑制eq oac(,4) 检测和连接边缘双阈值算法是对经过非极大值抑制的图像分别使用高、低2个阈值和分割得到2个阈值边缘图像和。由于图像是由高阈值得到的，因此它应该不含有假边缘，但可能在轮廓上有连续。因此双阈值算法要在中把边缘连接成轮廓，当到达轮廓端点时，该算法就在由低阈值得到的边缘图像的8邻域位置寻找可以连接到轮廓上的边缘，这样，

39、利用递归跟踪的方法不断地在中搜集边缘，直到将中所有的间隙都连接起来为止。(2)Hough圆检测Hough变换是一种用于区域形状描述的方法，经典的Hough变换常常被用于直线段、圆和椭圆的检测，广义的Hough变换可以推广至任意形状。两种变换的根本思想是将图像的空间域变换到参数空间，用大多数边界点满足的*种参数形式来描述图像中的曲线区域边界。Hough变换可以用于检测图像解析曲线，其中为解析曲线上的点，为参数空间上的点。对于圆，设其半径为、圆心为，则圆方程：(2.10)图像空间的圆对应着参数空间中的一个点，一个给定点约束了通过该点的一簇圆参数。传统的Hough圆检测是将空间域的每个轮廓点带入参数

40、方程，其计算结果对参数空间中的量化点进展投票，假设投票超过*一门限，则认为有足够多的图像点位于该参数点所决定的圆上。由于其需要逐点投票、记录，故花费的时间较长，占用的计算机存储量也较大。本文采用的Canny算子如图2-2(a)所示的虹膜图像，很多边缘对于检测虹膜外边界来说都是干扰边缘，如果在计算梯度时，对、两个方向均进展计算，则会产生过多过细的边缘，而这些边缘又是不需要的，而且会对提取有用的边缘起干扰作用。此时，如果有选择性的重点检测*一或*些方向的边缘，则既能节省时间又有利于排除干扰。针对虹膜图像，由于上下局部常常会被眼睑遮挡，如果将眼睑边界也过清晰地检测出来了，则对下一步的Hough圆变换

41、检测外边界产生很大的干扰，特别是眼睑根本成弧形，与圆形接近。为削弱非虹膜边界的影响，可以在进展边缘检测时重点检测有用的边缘，在保存足够的虹膜边界信息的同时尽量削弱其它边界的影响。注意到上下眼睑边界仅限于水平方向，在计算梯度幅值时对水平方向偏导值赋以较大的权重，而对竖直方向偏导值赋以较小的权重，从而削弱水平方向上的边界。(1)本文在定位虹膜外边界时，按照这一思想对Canny算子进展了改良，在计算梯度时只计算水平方向上的梯度，以防止过多地检测出眼睑边缘。即：方向偏导数：(2.11)方向偏导数：(2.12)(a)原始图像缩放0.4(b)竖直权重为0(c)两方向均匀权重(d)水平权重为0图2-2虹膜图

42、像以及Canny边缘检测的图像其中滤波时对边缘处，为简化，未处理，造成有虚假边缘，定位时除去，并不影响Hough定位圆。(2)在非极大值抑制时，选择比拟邻域为半径1.5像素围之，这样的到的边缘图像相对较好，边缘大于等于一个像素，小于等于两个像素，考虑到Hough参数的量化，比选择半径为1的邻域围造成的误差小，而与选择半径2的邻域边缘图相比，边缘得到了较好的细化。如图2-3。(a)1.5邻域半径抑制(b)2邻域半径抑制(c)1邻域半径抑制(d)1.5邻域半径抑制(e)2邻域半径抑制(f)1邻域半径抑制图2-3不同邻域半径边缘检测图(3)此外实际处理过程进展了少许改变：eq oac(,1)为了减小

43、处理的数据量，在定位虹膜外边界时，首先对虹膜图像进展了缩放，而且缩放定位后再折换成原始图像的圆，原始图像大小为。eq oac(,2)对于Canny检测处理后的边缘图像，直接一次将外两圆定位出，在定位外圆时比拟准确，但是定位圆时由于受其外围边缘的影响，很容易造成定位失败，因此首先定位出外圆，然后在外圆的围中定位圆，大大减小了干扰。eq oac(,3)由于圆相对较小，定位外圆时缩放比例0.4，定位圆时缩放比例为0.6。本文采用的Hough变换圆检测算法直接使用Hough变换是将图像上的每个边缘点都代入方程，且半径搜索围从O到图像平面所能容纳的上限，因此计算量非常大。本文对其进展了修改，并根据先验知

44、识及所采集的虹膜图像的特点尽可能地减少搜索围及参与Hough变换的点数以高效率。过程如下：(1)对图像进展边缘检测得到边缘图像；(2)根据先验知识，分别确定外边界圆半径的围，从而减小搜索半径；(3)圆的参数方程可改写为：(2.13)将图像空间中的边缘点而不是原图像中的每个点逐一代入上式，求出参数值。从图中可以看出，虹膜外边界左右两侧的局部所受干扰较少，边界质量较高，因此可以限制的取值围，相当于只统计左右两侧局部的边界点。由于虹膜边界圆的圆心不可能太靠近图片边缘(否则就没有将整个虹膜区域采集进来，虹膜信息可能太少而无法用于识别)，即圆心应位于图像中间区域的*个围，因此如果位于这一围，则将相应的累

45、加阵中的元素加1，否则加0；(4)找出中元素的最大值，即是对应半径为，圆心为，且圆周上边点最多的圆，即为边界圆。虹膜外边界的定位本文先定位虹膜外边界，然后再在外边界围以定位边界，从而能更快更准地定位边界。虹膜外边界上下局部常常被眼睑局部地遮挡，且往往还有睫毛的干扰，所以对外边界的定位先采用本文所述的改良的Canny算子进展边缘检测，从而在保存足够多的边缘信息的情况下尽量少地检测出干扰边缘，再利用改良的Hough变换进展圆检测，确定外边界圆的圆心和半径。由于瞳孔位于虹膜区域以，所以定位边界时，将搜索围限制在外边界围以，以减少搜索围。定位结果如图2-4图2-4虹膜定位结果2.2眼睑的分割在虹膜识别

46、过程中，虹膜区域往往会受到干扰，由于尽量考虑实用性，不能对被采集者提出苛刻要求进展配合，这样在被采集者不配合的情况下，获得的虹膜有效区域是很有限的，往往会是局部虹膜区域，其余局部被眼皮、睫毛等遮挡。从采集的图像数据库中也可以看出，多数图像中虹膜区域都是不完全的，因此虹膜区域不完全是普遍存在的。眼睑包括上眼睑和下眼睑，上下眼睑开度不同，裸露出的虹膜区域面积也不同。如果这局部眼睑也作为有效虹膜区域带入虹膜特征提取和编码模块，则得到的虹膜编码可信度将大大降低。检测眼睑比拟常用的方法有抛物线检测眼睑和直线检测法。在许多文献中，都采用了抛物线检测眼睑。Mathieu ADAM等在文献中首先在保存虹膜边界

47、的同时，平滑了虹膜纹理，然后进展抛物线检测；Tae-Hong Min等提出在归一化后的虹膜图像中采用抛物线检测，这样不需要考虑的影响；Young Kyoon Jang等提出检测眼睑候选区的曲线累积值，如果大于*一阈值则进展抛物线检测，否则进展直线检测。抛物线检测眼睑检测眼睑时，可以根据眼睑的形状采用抛物线拟合的方法，在多维参数空间求最优解。抛物线曲线的定位可通过广义Hough变换来确定抛物线的各个参数。抛物线曲线可表示为：(2.14)其中是用来控制抛物线的曲率，是抛物线的顶点，是抛物线相对于轴的旋转角。该方法的缺点是运算时间长。另外较常用的是将眼睑简化为直线的方法。这样可以大大的降低运算量。R

48、adon变换法分割眼睑对眼睑边界进展细分，可以近似认为其由一系列的线段连接而成。因此可以采用Radon变换直线检测法定位眼睑。Radon变换可在任意维空间定义，下面给出在2维空间的定义式：(2.15)式中，为整个图像平面，为图像在点的灰度，为坐标原点到直线的距离，为距离与轴的夹角，为冲激函数，定义为：(2.16)它使沿直线进展积分。如图2-5。Radon变换可以理解为图像在空间的投影，空间中的每一点对应着图像空间中的一条直线，而Radon变换是图像像素点在每条直线上的积分，也可理解为图像逆时针旋转角度后在水平轴上的投影。因此图像中每条直线会在空间形成一个亮点，直线的检测转化为在变换域对亮点的检

49、测。图2-5Radon变换法检测眼睑的步骤如下：(1)分别确定可能存在上、下眼睑的虹膜区域，以减少搜索围；(2)对可能存在眼睑的虹膜区域采用改良的Canny算子进展边缘提取。由于眼睑边缘呈水平方向，因此在计算梯度时只取竖直方向的梯度，以防止虹膜边界的干扰，即：方向偏导数：(2.17)方向偏导数：(2.18)(3)对提取到的边缘图像进展Radon变换直线检测，变换后会得到假设干个亮点，如果亮点的最大值大于*个阈值，则视其对应的直线为眼睑边界直线，否则认为没有上(或下)眼睑遮挡；(4)如果检测到了眼睑边界直线，则作一条水平直线，此直线过虹膜外边界与的交点，为简化处理，将此直线所在的虹膜区域标为干扰

50、点。(5)对检测到眼睑像素用灰度0代替。眼睑检测结果如图2-6。可以发现采用这种方法能比拟有效地将上下眼睑分割出去但检测方法过于粗略，易将局部有用虹膜当作眼睑分割出去，因此损失了一定量的有用虹膜信息。2.3剔除睫毛睫毛是用来保护眼睛的，可以有效阻挡异物进入眼睛，但是在虹膜识别时，落入虹膜图像区域的睫毛影响了虹膜特征的提取，如果不加以消除，会造成特征表示精度下降，难以反映实际的虹膜信息。由于虹膜编码识别算法本身具有一定的抗干扰能力，当睫毛稀疏且很短时，对虹膜识别的影响还不是很大。但是如果睫毛的数量较多、比拟长，在进展虹膜特征编码时很容易被误认为是虹膜的纹理特征，得到错误的虹膜编码。文献中有一些睫

51、毛剔除的方法。John Dangman基于虹膜灰度直方图4，提出睫毛检测的统计模型；Zhao feng He 等人也提出了一种基于统计模型的睫毛剔除法；D Zhang等首先用Sobel算子处理睫毛所在区域，然后判决是否为睫毛遮挡，假设是，则沿边沿方向进展1D中值滤波15；Byung Jun Kang等将睫毛区域分成两种类型-多条睫毛重叠区域和睫毛可分区域，然后分别处理。本论文中采用了简单的阈值法剔除睫毛。阈值法剔除睫毛睫毛通常具有灰度值较低的特点，和虹膜区域其它局部的灰度相比，睫毛局部灰度值更低，因此从灰度上根本可以将睫毛与虹膜区分开，本文采用了灰度阈值法别离睫毛。而且睫毛具有多而分散、形状有

52、弯有曲和存在穿插现象的特点，因此采用阈值法检测是比拟适宜的。设用于别离睫毛的阈值为，在可能存在睫毛的虹膜区域，如果像素的灰度值于，则认为是睫毛，标记为干扰点，记作1：(2.19)图2-6是局部虹膜外圆定位、眼睑分割及睫毛剔除的最终结果。图2-6局部虹膜外圆定位、眼睑分割及睫毛剔除结果2.4小结本章完成虹膜图像的预处理，采用Canny边缘检测和Hough变换对虹膜外边界进展了定位，采用Radon变换进展直线检测，检测上下眼睑，采用阈值法除去睫毛干扰，并给出了处理结果。-. z第三章 虹膜图像归一化及特征提取3.1虹膜纹理归一化由于拍摄环境等条件的改变，摄取的每虹膜图片的位置与大小都可能会发生变化

53、，所引起的虹膜的平移、缩放和旋转都会对虹膜识别造成影响，而人的瞳孔也会因外界环境而缩放。为了消除这些影响，将每幅原始图像调整到一样的尺寸和对应位置，实现准确的匹配，就必须对虹膜进展归一化。因为虹膜外边界都近似为圆形，因此可以用几何的方法归一化。通常的虹膜识别算法都假定外界光照不发生变化或者外界环境可控制，所以瞳孔的缩放可忽略不计，或简化为线性形变。在John Daugman系统中，提出了弹性片模型，认为虹膜区域可以任意拉伸或者压缩，总可以通过拉伸(压缩)将圆环区域变为矩形区域即极坐标表示的矩形区域。在Wildes系统中，通过采集图像和数据库图像的缩放、旋转变换进展特征区域配准，然后注册图像。在

54、Boles系统中，为了使两幅虹膜图像可以比拟，将一幅虹膜图像作为参考，将另一幅虹膜图像缩放为同样半径大小，这一规化只在匹配时进展，而不是为了将来的匹配而存储同一大小的虹膜特征模板。无论哪种情况，只要进展模式比对，必须要把不同分辨率的虹膜区域转化为一样分辨率的虹膜区域。平移虹膜平移是指虹膜位置在实际捕捉到的图像中的位置是不固定的，总是存在相对的位置变化。分析造成虹膜区域在图像中位置相对平移的原因，主要是拍摄时采集设备位置与人眼的位置存在光轴不重合，导致了图像采集中虹膜位置发生偏离图像中心位置。平移变化只能改变虹膜在图像中的位置，对于这一变化总可以通过坐标的平移将虹膜的中心平移到图像的中心位置，其

55、特征不会发生变化，因此说虹膜特征是平移不变的。旋转虹膜旋转是指采集设备和虹膜所在的两个平面发生了相对的角度旋转变化，导致采集的不同虹膜图像发生了相对的旋转变化。造成虹膜图像发生旋转的原因是：(1)头部倾斜(或采集设备倾斜)；(2)眼球在眼眶的转动。伸缩虹膜的伸缩变化是指图像中虹膜面积大小的变化。发生伸缩变化的原因有：(1)拍摄时相机距离人眼的距离发生变化，造成虹膜大小不同；(2)瞳孔对光的反响，造成边界的伸缩变化，造成虹膜大小改变。这两种改变虹膜大小的情况在实际中是普遍存在的，拍摄距离的变化会造成虹膜整体大小(外边界和边界)的同时变化。在拍摄距离一定时，虹膜区域的大小变化主要来源于瞳孔的变化，

56、外边界是不变的。瞳孔的作用是控制进入瞳孔的光通量，当光线强时，瞳孔缩小，对虹膜区域产生拉伸作用，虹膜区域面积增大，当光线变暗时，瞳孔增大，对虹膜区域产生压缩作用，虹膜区域面积减小。由于虹膜区域近似为弹性模板，所以无论是拍摄距离导致的图像大小变化还是瞳孔大小使得虹膜发生拉伸变化，虹膜区域上任一点的位置都是相对不变的，都可用外边界的线性组合表示。分辨率虹膜图像存在大小变化使得虹膜区域的分辨率不一致，即使来自一样虹膜的图像，在不同条件下得到的虹膜区域分辨率也可能是不同的25，分辨率不一致使得虹膜的特征点无法进展有效比对。只有将不同大小的虹膜图像转化为一样大小的特征区域，才能实现不同大小虹膜图像的有效

57、比对。关于分辨率的选择没有一个统一的标准，在文献中采用的分辨率有 ，等，不同分辨率反映了不同频率的虹膜细节信息。对于分辨率不同的图像进展比对首先要统一分辨率，由于高分辨图像中同时包含着低分辨信息，但是低分辨图像无法获得高分辨信息，所以本章用极坐标对虹膜区域进展规化，规化区域的大小根据大尺度特征不变的思想，利用虹膜纹理的发散状分布以及邻近点的相关性，将虹膜区域统一用低分辨率表示，从而实现不同大小虹膜的特征点数统一。3.2虹膜的Rubber-Sheet模型虹膜图像归一化方法中比拟常用的Daugman的Rubber-Sheet模型。该算法把虹膜区域表示成一个具有各向同性和柔韧性的弹性模型。利用该方法

58、经虹膜圆环区域拉伸成指定大小的矩形区域。其算法原理的映射公式如下：(3.1)(3.2)其中、分别是虹膜边缘(也即瞳孔边缘)和虹膜外边缘的坐标。归一化过程可以用图3-1表示：图3-1归一化过程在采集得到的虹膜图像中，虹膜外圆心不重合的情况是很常见的，其原因可分为以下两种：(1)由眼睛的生理构造决定的。眼科学研究说明，人眼的虹膜与瞳孔不完全是圆形，而是椭圆形或接近椭圆的不规则多边形，所以将虹膜外边界当作圆形边界进展定位时，会出现外圆心偏差。大量观察发现，由生理原因造成的虹膜外圆心偏差一般比拟小。(2)采集过程中，由于眼球转动造成圆心偏差。眼球是生物活体，在拍摄过程中不可能始终保持固定不动。由于摄像

59、头微距拍摄时会出现桶形失真，即当眼球转动一个微小角度时，虹膜图像就会产生很大的畸变，从而导致外圆心出现偏差。由于通常情况下，虹膜、外边界的圆心是不同心的，二者的圆心有一定的偏离(不完美虹膜的定位分割)，偏离使得虹膜区域成了不同宽度的环状区域。在将其转换成矩形区域的工程中，必须对这一偏离进展处理。根据下面的推到过程就可以将虹膜的、外边界都用瞳孔中心为极点的极坐标展开。虹膜边界表示图3-2虹膜边界关系如图3-2所示虹膜边界关系。设为瞳孔圆心，坐标为；为虹膜外边界的圆心，坐标为。、外圆心间的距离为：(3.3)圆心连线与水平方向的夹角为：(3.4)虹膜外边界的半径为：，用瞳孔中心表示时半径为，则有：(

60、3.5)将此式展开并化简的到关于的方程：(3.6)求其解，并选择整数半径为：(3.7)此式中，对一幅虹膜图像，圆心连线距离和其与水平方向的夹角及虹膜两边界的半径均为常数，这样将虹膜的外边界表示成角度的函数，也是用瞳孔中心表示虹膜外边界的表达形式。由此虹膜的、外边界都可以用瞳孔中心为极点进展表示。虹膜区域的表示虹膜区域是介于瞳孔边界和虹膜外边界之间的局部，瞳孔的边界点的位置和虹膜的外边界点的位置均可由瞳孔的中心表示，将虹膜区域看作弹性模型，整个虹膜区域也就可以采用虹膜的、外边界坐标的加权来表示，权值一样的特征点组成了一个圆周，如图3-3。图3-3调整圆心虹膜区域上任一点的位置都可以看作虹膜、外边