虹膜识别算法研究及实现

摘要

在当今信息化时代，如何精确鉴定个人的真实身份、保护信息安全，已成为一个急待解决的关键性问题。传统的身份认证极易伪造和丢失，难以满足急速发展的社会需求，目前最为便捷、安全的解决方案就是生物体识别技术。它不仅简洁快速，而且安全、可靠、准确。同时更易于配合网络和安全、监控、管理系统整合，实现自动便捷管理。虹膜识别是一种基于人体生理特征的生物体特征识别技术，与人体纹理、掌纹、脸相、音频、步频、血液等特征识别相比，具有唯一、高稳定、识别率高、检测方便等优点，因此虹膜识别技术己成为当前身份鉴别研究的热门领域。

本论文详细阐述了虹膜识别技术的研究背景和现状、虹膜生理结构和虹膜识别系统的构成。较深入的研究了虹膜识别算法，实现了三个步骤，即虹膜数字图像预处理、特征码提取和模式匹配。本论文的研究工作主要集中在对虹膜数字图像预处理的归一化和模式匹配两方面的研究。

在虹膜数字图像预处理的归一化过程中，采用一种基于Dangman橡皮片的辐射线段的归一化方式，将环型区域改变成为矩形区域。这种方式采用若千条线段表示两个非同心圆周之间的区域，只要设定这些线段上的点数就可以起到很好的归一化成果，并且仅出一种分析模型。更加定位搜索的效率，节约时间，简化定位的过程。采取Hough变换算子对虹膜进行精定位。

在特征提取及编码匹配过程中，本论文首先将归一化后的2D虹膜数字图像转换为1D灰度信号，从而减小了运算量;然后运用1D Gabor小波对构造的1D 灰度信号进行分析，选取一定尺度上的小波变换结果进行量化，生成二进制的特征向量，从而提高了虹膜识别技术的效率;最后通过改进的Hamming距离移位匹配，实现了识别中的旋转不可逆性。

本论文通过使用中科院自动化研究所的虹膜数据库(CASIA-IRISV1)的虹膜数字图像进行实验。经过一定工作量的科学实验对这一算法进行了验证。实验结果：本论文选取的CASIA-a中心波长为20像素效果最好。最好的模板大小为

辐射分辨率22像素，角分辨率为240像素。这些参数产生了信息量为9600比特的生物体模板。为了修正不一致旋转，平移10个单位在CASIA-a中效果最好。FAR与FRR分别为0.004%与0.236%，这个识别率己经足够用于实际;而传统的海明距离FAR与FRR分别为0.012%与0.475%。

关键词:虹膜识别，模式识别，Hough变换，Gabor小波，Hamming距离

Abstract

Keyword: Iris 1Dentify, Mode 1Dentifies, Hough transformation, Gabor small wave, Hamming is apart from

目录

摘要 (1)

Abstract (3)

第1章绪论 (5)

1.1课题研究背景 (5)

1.2几种典型的生物体识别技术 (5)

1.3国内外研究现状分析 (7)

1.4本论文完成的主要工作 (9)

第2章虹膜识别原理 (10)

2.1虹膜特征 (10)

2.2虹膜识别技术流程 (12)

2.3几种典型的虹膜识别方式 (13)

2.3.1 微积分算子法 (13)

2.3.2 Wildes虹膜识别法 (14)

2.3.3 Boles虹膜识别法 (14)

2.3.4 中科院自动化所的最小二乘拟合法 (14)

2.4技术难点 (14)

2.5本章小结 (15)

第3章虹膜数字图像预处理 (16)

3.1虹膜边缘定位 (16)

3.1.1虹膜内边缘粗定位 (16)

3.1.2虹膜外边缘粗定位...................... 错误！未定义书签。

3.2睫毛与噪音检测 (20)

3.3虹膜数字图像归一化及增强 (23)

3.4实现结果 (26)

3.5本章小结 (27)

第4章虹膜特征码 (28)

4.1小波函数.................................... 错误！未定义书签。

4.1.1小波概念.............................. 错误！未定义书签。

4.1.2小波的性质............................ 错误！未定义书签。

4.1.3常用小波.............................. 错误！未定义书签。

4.2基于2D Gabor滤波器的编码方式............... 错误！未定义书签。

4.2.1 2D Gabor小波特性..................... 错误！未定义书签。

4.2.2 2D Gabor小波的特征提取............... 错误！未定义书签。

4.3局部过零检测提取算法 (33)

4.4 1D Log-Gabor滤波器虹膜特征提取 (34)

4.5本章小结 (36)

第5章模式匹配 (37)

5.1模式匹配概论 (37)

5.2加权欧式距离 (38)

5.3改进的Hamming距离 (39)

5.4实验效果 (43)

5.4.1虹膜模式唯一性 (43)

5.4.2个体识别 (45)

5.4.3结论 (51)

5.5本章小结 (53)

第6章总结与展望 (54)

6.1研究工作与成果 (55)

6.2展望 (55)

参考文献 (57)

致谢 (59)

第1章绪论

1.1课题研究背景

进入二十一世纪以来，随着科学技术不断进步发展，电子信息日益影响这人们的日常工作生活，并起到越来越重要的作用。网络、安全、金融、电子商务、机场、地铁等，都需要可靠而准确的身份鉴定。传统的身份鉴别方式主要基于身份表示物品和身份标示知识。身份标识物品如证件、钥匙、银行卡等。身份标识知识如使用者名、密码等。随着科技进步这些方式的缺陷越来越突出，身份识别物可能丢失，也可能被伪造，使用者名、密码可能被忘记，也可能被破解泄露。一旦他人获取这些身份标识物品和知识，系统无法区别真正的用有者和仿冒者。后来又发展同时带有身份识别知识的身份识别物品，如数字证书和智能卡等，当其效果仍然满足不了人们对信息安全的需求。

为克服传统身份鉴别方式，人们提出基于人体生物体特征的新身份识别技术。所谓生物体识别技术就是，通过电子信息与光学、声学、生物学等学科密切结合，利用人体固有的生理特性，（如人体纹理、脸型、虹膜等）和行为特征（如笔迹、音频、步态等）来进行个人身份的鉴定。它具不易遗忘、防伪性能好、不易伪造或被盗、随身“携带”和随时随地可用等优点。

1.2几种典型的生物体识别技术[1][2][3]

随着科技的进步，我们已经实现了多种生物体识别技术，如人体纹理识别、虹膜识别、人脸识别、签名识别、语音识别等，其中部分技术性高的识别手段目前还存在于实验室中。我们有理由相信科技的进步，将会有越来越多的生物体识别技术应用到实际环境中。

人体纹理识别

在目前实际应用的多种识别技术中，最为广泛的人体纹理识别。它是基于对比我们手掌及其手指、脚、脚趾内侧表面的皮肤凸凹不平产生的纹路所形成各种各样的图案识别的。由于这些纹路的存在增加了皮肤表面的物理摩擦力，使得人

们可以用手来抓起重物。同时包括指纹在内的这些皮肤的纹理在图案、断点和交叉点上都不尽相同，也就是说，这些纹理是仅有的。也就是这种唯一性的存在，使得我们可以通过把使用者个人同他的指纹建立起来的对应，将使用者的指纹和预先保存的指纹进行对比，验证使用者的真实身份。

这种依靠人体的独有纹理特征来进行身份验证的技术称为人体纹理识别技术，指纹识别是生物识别技术当前应用最广泛的一种。它比较适用于室内安全系统，首先它可以有充分的条件为使用者经行讲解和培训，能较好的适应系统运行环境。相对其他识别系统，指纹识别具有体积小、价格底的优点，键盘、鼠标、打卡机等设备可以轻易的将其集成在内，因此目前指纹识别广泛使用的在计算机系统安全访问系统和考勤系统等小型设备中。

●语音识别

语音识别就是通过把使用者的语言和语声转换成可进行处理的信息的识别技术。广泛应用在语音邮件、语声开锁、文字输入等方面。由于音频传感设备和人的语言发音可变性很大，使用的步骤比较复杂，在需要安静的场合下不合适使用，使得这种识别系统应用起来不太方便。

●签名识别

签名识别在应用中具有其他生物体识别所没有的优势，人们已经习惯将签名作为一种在交易中确认身份的方式，它的进一步的发展也不会让人们觉得有太大不同。实践证明，签名识别是相当准确的，因此签名很容易成为一种可以被接受的识别符。但与其他生物体识别产品相比，这类产品目前数量很少。

●人脸识别

目前人脸识别技术是比较热门的计算机技术研究领域，它是对我们人体本身的生物特征来鉴别区分人的。它特指通过对比分析使用者脸部视觉特征信息进行鉴别使用者身份的识别技术。由于近年通信技术的发展使得视频通讯变得简单方便，因此有人夸大了这种识别的相关功效。实际应用当中一般要求使用人多次采集图像，这是由于它易受到光线变化变化、采集或比对时脸部位置的影响。

当前最有效的是基于主动近红外图像的多光源人脸识别技术。它能有效克服光线变化对采集人脸图像时的影响，取得了相当不错的识别能力，系统整体性能具有很高的稳定性、精确性、识别速度。

●基因识别

基因识别是当前生物信息学的一个重要分支，利用生物学实验或计算机等识别人体细胞DNA序列上的生物学特征进行比对识别。识别的主要对象是蛋白质编码基因，当然也包括RNA基因和调控因子等具有一定生物学功能的因子。目前人类对于基因组研究就是建立在基因识别上的。

●虹膜识别

虹膜识别是一种与人的眼睛有关的生物识别。它使用普通的摄像元件，不需要使用人与机器发生密切接触。据有关科研人员的考证和相关标本证明，人眼睛的虹膜是人还在胚胎时就形成的。其虹膜个体由于胚胎发育和母亲子宫环境的影响，形成了虹膜个体独有特征，在人体的成长发育中，它会保持罕有的稳定性。由于虹膜的生理结构独特，同时还有这高度得仿伪性和不可更改的特点，使得虹膜识别比其它的生物体识别技术具有更好的防伪造能力。到目前为止，虹膜识别相对其它各种生物体特征识别技术的识别正确率很高。因为虹膜识别算法相对简单、有效，使得相关研究越来越多的吸引广泛学者与企业的高度关注。

1.3国内外研究现状分析

国外研究现状

人体虹膜最早用于身份识别的尝试是由法国人Alphonse Bertillon首先做到的，他在1885年对法国巴黎刑事监狱中通过虹膜对的犯人进行身份确认，当时他主要通过人眼对人的虹膜颜色和形状的进行观察。其后一直到了上世纪八九十年代，人体虹膜识别技术才取得了快速发展。期间在1987年，眼科专家Lenonard Flom和Aran Safir医生首次提出利用虹膜数字图像进行自动虹膜识别的概念[4]，但令人遗憾的是他们并没有开发出能用于实际的虹膜识别系统，仅在理论上对人体虹膜识别技术做出了研究。在1991年，一名叫Johnson的学者在美国洛斯阿拉莫斯国家实验室首次实现了一个真正可以自动识别人体虹膜的系统，这是有文献记载以来的虹膜识别技术在应用系统上的首次使用[5]。如今，虹膜己成为生物体特征识别领域的热门研究方向。人体虹膜识别技术目前在国外的研究较多，这方面主要有美国的IrDian&SenSar和LG等公司在研究、应用，而且Ir1Dian公司还在1999年位列美国1998年度发展最快的企业500强中。在从多虹膜识别算法中当以英国剑桥大学的学者J.Daugman和美国麻省理工大学的学者

R.P.Wilde s提出的相关算法最为常用。当前，大多数虹膜识别商业系统都是以这两种算法做基础来实际应用的，特别是J.Daugman研发的基于人体虹膜分析的身份识别系统在1994年通过了美国专利申请，1993年他提出基于2D Gabor小波的虹膜纹理编码算法,目前这种算法已经成为当今所有商用虹膜识别应用系统的根本。相关的虹膜产品目前已经小范围的试运行在金融系统中。虽然虹膜识别技术尚未得到大规模的实验验证，但是在国外一些国家己经开始在实验检测虹膜识别在生活中大规模应用的可能性。在美国得克萨斯州联合银行的三个营业部内Ir1Dian公司研制出的虹膜识别系统已经应用，使用者无需任何其他标识物就可以来银行办理相关业务。当使用者在银行设备上操作时，用于采集人体虹膜的摄像机会首先对使用者的眼部进行扫描，然后系统会将扫描到的使用者虹膜图像转化成数字信息的，同时比对数据库中存储的使用者虹膜资料，这样就可以验证出使用者的实际身份是否真实。

国内研究现状

目前在国内识别技术的实际应用上人体纹理识别技术相对比较成熟使用范围比较广，如打卡器，计算机开机保护等。我国虹膜识别的实际应用研究起步相对较晚，在上世纪末期，研究虹膜识别技术实际应用的工作才刚刚开始进行。1998年，上海交通大学开始跟踪研究虹膜识别技术的相关工作。华中科技大学有相关科研人员在2000年开始研究虹膜图像的相关性。虹膜识别的演示系统同年在中国科技大学得到了实现。本世纪初，一套较高效的虹膜识别方式在北京大学信息科学中心研究出来。中科院自动化所相关研究人员，采用不同的方式对人体虹膜经行识别研究，并取得了很好成果。在本世纪初，具有我国自主知识产权的人体虹膜识别原型系统，在中科院自动化研究所拥有与国际同步的自主产权的核心技术和算法的研究人体虹膜识别的国家某重点实验室成功开发出来，同时该实验室还利用研制的人体虹膜数字图像采集设备，建立了用于科学研究的虹膜数据库——CASIA人体虹膜数字图象数据库，并使其成为我国用于人体虹膜识别技术研究领域的公共数据库。同时中科院自动化研究所还成立了相关企业，目前该企业是世界上第二个拥有核心算法的企业。2007年11月，我国颁布实施《信息安全技术虹膜识别系统技术要求》(GB/T20979-2007)标准。

总的来说，虽然我国人体虹膜识别技术的理论性研究起步较早，但研究虹膜

识别技术在实际环境的应用相对国外还比较晚。

1.4本论文完成的主要工作

论文预期完成研究虹膜识别的最新理论，完善小波理论在生物体识别技术中的应用，在前人提出的有关虹膜识别算法的基础上，展开更加系统和深入的讨论，力求探索模式分类的一些新方式。主要工作包括:

1）虹膜识别技术的基本理论：包括虹膜的生理结构和该识别技术的技术流程。介绍虹膜识别方式。从图像采集和算法的通用性两方面说明了虹膜识别的技术难点。

2）虹膜预处理的各个环节的说明：首先是虹膜预处理。通过采用Canny算子对虹膜数字图像进行边缘检测。其次使用有关边缘信息的Hough变换提取虹膜。最后，利用Daugman算子实现虹膜边缘的精确定位。小波变换说Log-Gabor小波的虹膜编码方式。

3）利用Hamming距离进行匹配，以获取识别率。

第2章虹膜识别原理

2.1虹膜识别特征

虹膜识别技术是通过利用虹膜个体终身不变和单个差异大的特点来识别身份的。通过相应的算法与虹膜数字图像相结合，获得性当高的精准度，就算某一数据库内具有当前人类所有的虹膜个体图像，虹膜比对错误也不会出现。但是这项技术的无法录入问题已经成了它同其它识别技术抗衡的最大障碍，由于虹膜识别技术的高精度使其能够在众多识别技术中占有一席之地。如图（2-1）所示，简单的展示出人类眼球的结构，其中虹膜是围绕瞳孔的一层生理薄膜。虹膜是包裹在眼球上的彩色环状物，每一个虹膜都包含一个独—无二的基于像冠、水晶体、细丝、斑点、结构、凹点、射线、皱纹和条纹等特征的结构。眼睛的虹膜是由相当复杂的纤维组织构成，其细部结构在出生之前就以随机组合的方式决定下来了，虹膜识别技术将虹膜的可视特征转换成一个512个字节的Iris Code（虹膜代码），这个代码模板被存储下来以便后期识别所用，512个字节，对生物体识别模板来说是一个十分紧凑的模板，但它对从虹膜获得的信息量来说是十分巨大的。

从直径11mm的虹膜上，Dr. Daugman的微积分算法用3.3个字节的数据来代表每平方毫米的虹膜信息，这样，一个虹膜约有266个量化特征点，而一般的生物体识别技术只有13个到60个特征点。266个量化特征点的虹膜识别算法在众多虹膜识别技术资料中都有讲述，在算法和人类眼部特征允许的情况下， Dr. Daugman指出，通过他的算法可获得173个二进制自由度的独立特征点。这在生物体识别技术中，所获得特征点的数量是相当大的。

图2-1人类眼球构造

目前统计表明虹膜识别的错误率是各种生物体识别错误率最低的。虹膜识别比起其他识别技术的优势是:

1）生物体活性: 虹膜是人眼体的可见部分，在巩膜的保护下，具有极强的生物体活性。例如，瞳孔的大小随光线强弱变化；视物时有无意识调节过程；有每秒可达十余次的瞳孔缩放；在人体脑死亡、处在深度昏迷状态或眼球组织脱离人体时，虹膜组织会立即完全收缩，出现散瞳现象。这些生物体活性与人体生命现象同时存在，共生共息，所以想用图像、视频、尸体的虹膜代替活体虹膜数字图像都是不可能的，从而达到保证了生理组织的真实性。

2）非接触性:从一定距离即可获得虹膜数字图像,无需人接触系统,对人身没有侵犯，易被公众所接受；

3）唯一性: 每个虹膜所包含的生物信息都不相同，出现形态完全相同的虹膜组织的可能性远远低于其他组织。虹膜的纤维组织细节复杂而丰富，并且它的形成与胚胎发生阶段该组织局部的物理化学条件有关，具有极大的随机性，即便使用克隆技术也无法复制某个虹膜。同卵双胞胎的虹膜纹理信息不同，同一个人左右眼的虹膜纹理都不会相互认同；

4）稳定性: 虹膜在人的一生中都极其稳定，出生前（胎儿7个月时）已经形成，出生6—18个月后定型，此后终身不变。一般性疾病不会对虹膜组织造成损伤，不会因职业等因素造成磨损。

5）防伪性: 不可能在对人视觉无严重影响的情况下通过外科手术改变虹膜特征，更不可能将一个人的活体虹膜组织特征改变得与某个特定对象的虹膜组织特征相一致，用图像、视频、尸体的虹膜代替活体的虹膜数字图像都可被检验出来。

2.2虹膜识别技术流程

虹膜识别系统一般由：图像提取、图像预处理、特征码提取和特征码匹配四部分组成。

1）图像提取：通过虹膜采集系统取得含有丰富信息的可供识别的虹膜数字图像，并将得到的数字图像传输到系统内存储。

2）图像预处理：获取到的眼部图像包括了很多无用途信息。应先消除采集到的虹膜数字图像噪声斑点和光照对图像的影响；其次根据虹膜几何、灰度特征定位虹膜内外边缘和眼睑，将虹膜从原始的图像中分离出来。然后对其进行归一化处理以达到消除图像提取时的因旋转、尺度变化等所代带来的影响。

将归一化处理后的图像进行增强，得到预处理好的图像；最后将预处理好的图像供给虹膜特征提取模块提取相应的特征。

3）特征提取：通过一定的算法从分离出的图像提取出独特的特征点，并对其进行编码。虹膜的特征提取方式跟具体的识别系统相关，因为不同的系统可能采用不同的虹膜特征。现有的虹膜特征提取算法可以大致分为以下三类：①对局部纹理定性分析的相位编码法②纹理分析法③基于1D小波变换的过零点检测的方式。如何有效地描述丰富的虹膜纹理，并提取适当的特征来描述其细节，是虹膜识别算法的核心部分．

4）特征匹配：根据采集到的虹膜数字图像进行特征提取获得的特征编码与数据库事先存储的虹膜数字图像特征编码进行比较、验证，从而达到识别的目的。

针对不同的特征所设计的检验方式，在特征提取之后，虹膜特征被记录为特征向量。如何衡量特征向量之间距离或相似程度属于特征匹配问题，也是识别系统中的重要的组成部分。主要关心的指标为：准确率、匹配速度。

图2-2虹膜识别系统流程

2.3几种典型的虹膜定位

2.3.1 微积分算子法

John Daugman 是虹膜识别技术的创始人，他开发的虹膜识别系统是虹膜识别技术中相对比较成熟的。在实践中，这套系统被验证具有非常好的效果。该系统使用Hamming 距离（Hamming Distance ）对待识别的虹膜编码位与数据库中的虹膜编码进行匹配。简单来说就是一个圆形检测匹配器，查找随着半径变化的相应圆周上像素平均值变化最大值所对应的圆心坐标和半径三个参数，以此来确定瞳孔和虹膜的边缘以及虹膜与巩膜的边缘。

在定位时使用了从粗到精的策略最后达到单像素的精度，并估计出虹膜和瞳孔的中心及半径。一般瞳孔和虹膜是不同圆心的，它的半径是虹膜半径的0.1到0.8倍。因此瞳孔圆周的三个参数应与虹膜圆周的参数分开估算。

对于人眼图像灰度分布I(x,y)，首先粗略获得瞳孔的圆心区域，得到一个瞳孔圆心点集。然后依次以每个瞳孔圆心点(x 0,y 0)作为瞳孔圆心，沿半径r 方向搜寻圆形模板。通过微积分操作在半径为r 的圆周上的积分，然后计算在半径方向r 的梯度变化。对于每一个候选圆心，以半径r 为变量，搜索梯度变化最大处对应的半径。梯度变化最大处对应的候选圆心和半径即为所求瞳孔圆心和半径。在已知瞳孔圆心和半径的基础上，利用经验确定一个较小的虹膜圆心点集，并用与瞳孔半径搜寻相同的方式确定虹膜外边界的圆心和半径。

微积分操作定义如下：

ds r

y x I r r G y x r y x r ???*πσ2),()(max ),,(),,(0000 （2-1） 其中：*为卷积分；),(y x I 为采集到的虹膜图像；r 是搜索的圆半径；)(r G σ是一个尺度σ的高斯函数；2202)(21)(σσσπr r e r G --=，由于高斯函数的傅里叶变换仍然是高斯函数，所以)(r G σ同时可以在空域和频域中其平滑滤波的作用；ds r y x I y x r ?π2),(),,(00 是),(y x I 在以r 为半径，(x 0,y 0)为圆心的曲线积分。

由于微积分算子法的求解方式运算时间长，因此影响采取此算法的虹膜识别速度。

2.3.2 Wildes虹膜圆定位法

Wildes识别法就是将瞳孔的形状模拟成一个椭圆，计算输入图像和输出图像的均值和方差等统计信息，然后将两幅图像的协方差作为相关性系数，最后采用Fisher线性判断确定两幅图像是否匹配。

2.3.3 Boles虹膜定位法

Boles识别法就是通过1D Gabor对沿虹膜中心同心圆的一条采样曲线进行过零点描述，利用两个自定义的不相似度函数计算两个虹膜纹理间的距离，从而得到识别结果。

2.3.4 中科院自动化所的最小二乘拟合虹膜定位法

利用Gabor滤波和Daubechies-4小波变换对特征进行提取。通过选择中心频率和相位角，获得24个Gabor滤波通道，对每个通道的滤波结果，提取平均值和方差作为代表该通道的特征，从而抽取48个特征，然后利用Daubechies-4小波分解为7个通道进行编码。

2.4技术难点

一、人体虹膜数字图像的采集

虹膜识别技术主要通过虹膜纹理信息来识别不同人的身份，因此得到清晰的虹膜数字图像是虹膜识别技术的基础。目前，图像采集方面存在的问题有以下两点:

1)图像采集系统对光线要求很较高。光线明、暗会影响人的瞳孔的大小，进而影响图像采集;不稳定的光照条件会改图像的灰度分布特征，从而影响虹膜定位算法精确度，进而影响识别结的精准度;识别过程中，虹膜区域产生的光斑会减少虹膜纹理的识别率，降低识别准确率。

2)图像采集系统不易使用。在图像采集过程中，首先，使用者往往要长时间地注视摄像镜头，同时不断调整姿势，以免虹膜数字图像偏出镜头;其次，由于摄像机镜头的焦距限制，需要使用者的眼睛与镜头间的距离处于指定的范围内，否则将导致图像模糊、影响虹膜识别系统识别。这些条件无疑会增加虹膜识别设备的使用难度，并引起人的反感。当然，可以通过很多技术手段克服上述困难，如:增加图像质量评估系统、摄像机镜头自动调焦等软、硬件辅助系统，但

这无疑会使整个系统的成本增加，同时增强使用难度。

二、现有虹膜识别算法不通用性

目前各虹膜图像采集装置工作原理、软硬件条件存在明显差异，所以还没有一个标准化的、通用性强的虹膜数字图像库提供研究者使用。很多科研人员为方便进行算法研究都研制了自己的虹膜数字图像采集系统，并利用采集系统所采集的图像库进行算法研究。由于不同的图像采集系统采到的图像的灰度分布特征存在一定差异，而虹膜的定位和特征提取算法又是建立在特定的灰度分布特征基础上，因此，很难在不同类型的虹膜数字图像上应用同一虹膜识别算法。

2.5本章小结

本章简单介绍了人体虹膜的生理结构，描述了虹膜识别系统实现的技术流程，详细介绍了几种成熟的虹膜识别方式，最后讨论了虹膜识别技术的难点。

第3章虹膜数字图像预处理

虹膜识别的第一阶段就是把虹膜从图像中提取出来。需要解决的问题是不但要定位虹膜区域，还要分离并去除人为现象。分割阶段对于虹膜识别系统的成功非常重要，因为错误表达的虹膜模式数据会破坏生成的生物体模板，最终导致低识别率。

成功的分割取决于人眼图像的质量。本论文使用的虹膜数字图像来自中国科学院自动化所（Chinese Academy of Sciences Institution of Automation,简称CASIA）CASIA-IrisV1虹膜数据库，因为使用红外光来照明所以没有镜面反射。

3.1虹膜边缘定位

虹膜定位的目的是找出人体眼睛图像中瞳孔与虹膜，虹膜与巩膜之间的两个边界，即虹膜的内边界和外边界，把虹膜从其他无用信息中提取出来。本论文针对虹膜特点，采用不同的定位方式定位虹膜内、外边缘。

3.1.1虹膜内边缘粗定位

瞳孔粗定位的目的是获取粗略的瞳孔圆心和半径，为精定位提供搜索范围。其具体步骤如下：

l ）首先分析虹膜图像的灰度直方图，如图3-1所示。

根据虹膜图像的灰度分布特点，可以得出这样的结论：第一个波峰对应瞳孔区域，第二个波峰对应虹膜区域，第三个波峰对应巩膜区域。因此，找到第一个和第二个波峰之间的波谷作为二值化阈值，从虹膜原图像中分割出瞳孔，如图3-2所示

图3-1 虹膜图像的直方图图3-2分割的瞳孔图像

二值化阈值分割公式如下：

???≤>=T

y x f T y x f y x G ),(,0),(,1),( （3-1） 其中，T 为分割阈值，),(y x f 为图像灰度值。

2 ）采用水平和垂直方向灰度投影法粗定位瞳孔圆心),(00y x 和半径0r 。

∑=x

y x I x I ),()( (3-2)

∑=y

y x I x I ),()( (3-3)

如上公式分别为水平和垂直方向的灰度投影和，其中),(y x I 是原虹膜图像灰度值矩阵。灰度投影累计和最小点),(00y x 对应瞳孔圆心

(a)水平方向投影 (b)垂直方向投影图

3-3 瞳孔的灰度投影

在二值化后的瞳孔图像上，选取行30±y 。对选取的每行，分别向左和向右搜索第一个灰度值为0 的点x 1和x 2 ,每行对应的半径r 为点x 1和x 2平均值。

最后将所有的半径取平均值，得到粗定位的瞳孔半径r 0 。

3 ）己经得到粗定位瞳孔圆周的圆心和半径，是对整幅虹膜图像粗定位的结果，易受到睫毛，眼皮等噪声的影响。因此，可以选取我们感兴趣的瞳孔区域，

如图3-4所示，此瞳孔图像包含了瞳孔，部分虹膜，以及少许睫毛等噪声。

图3-4选取的瞳孔区域矩形图

分离瞳孔后，采用相同的灰度投影法再次定位瞳孔圆心),(11y x 及半径r 1 。

这时分割的瞳孔因为受睫毛和眼睑等影响很小准确率非常高。

3.1.2虹膜内边缘精定位

瞳孔的圆心和半径大致确定后，就可以对内边缘进行精定位。本论文采用索贝尔（Sobel ） 边缘检测算子结合Hough 投票机制精定位虹膜内边缘。

3.1.2.1索贝尔（Sobel ） 边缘检测算子

边缘是分界两块灰度值相对一致区域的届线，是梯度大或灰度函数不连续的区域。图像灰度变化可以用梯度变化反映，所以可以用局部微分技术得到边缘检测算子。

由灰度分布特点可知，虹膜的边缘轮廓属于典型阶跃性边缘。对阶跃性边缘点来说，图像灰度在点的两边变化规律应是该点灰度变化的一阶导数的极值点，也就是二阶导数过零点，阶跃性边缘点使用一阶微分子算法比较有效。另外，虹膜与巩膜之间的灰度差异要比瞳孔与虹膜之间的灰度差异小，虹膜内边缘非常容易检测出来。在一阶微分算子中，索贝尔（Sobel ）算子比较简单，且计算速度快，因此采用索贝尔（Sobel ）算子检测瞳孔边缘。

对图3-4采用索贝尔（Sobel ）算子进行边缘检测，生成如图3-5所示检测图像，检测的边缘信息基本上都是有效的瞳孔边缘点，噪声点非常少。

图3-5索贝尔（Sobel ）边缘检测图像

3.1.2.2Hough 圆检测投票竞选策略

标准Hough 变换（S TH ）是一种用于提取边缘特征的简便而有效的方式，它能够提取直线、圆、椭圆、二次曲线甚至是任意形状的边缘。它所实现的是一种从图像空间到参数空间的映射关系。STH 的基本思想是将图像空间内具有一定关系的像素进行聚类，寻找能把这些像素用某一解析形式联系起来的参数空间累积对应点。

Hough 变换检测圆的原理：假设{}

1,2,3,...n j y ,x j j =为圆周上的点，则圆周上的点),(y x 满足在参数),,(r βα系下的方程：

222)()(r y x =-+-βα (3-4)

其中r 表示半径，),(βα表示圆心坐标。上式表示一个锥面，对于图像中任意一个点，均有参数空间中的一个3D 锥面与之对应。对于同一个圆周上的点集，这些3D 锥面构成圆锥面族，并且相交与3D 空间中的一点),,(000r βα，则这个圆周的圆心坐标为),(00βα，半径为r 0 。

Hough 变换在参数空间不超过2D 的情况下，有着理想的效果。但是，伴随着参数空间增大，计算量便会大幅上升，同时会耗费掉巨大的存储空间。为提高虹膜定位速度，本论文采用Chiu [25]提出的Hough 圆检测投票策略方式。

根据内边缘粗定位结果，确定候选圆心集及其半径搜索范围。以)3,3(11±±y x 作为圆心候选集，以15r 1±作为瞳孔半径进行搜索，并建立相应的3D 数组（大小为7 * 7 * 31 ) ，初始化为0 。对图2 一中的每个边缘点，计算时每个圆心候选点的距离r x ，在数组中对应r x -1 , r x , r x +1位置加1 ，最后将数组中累加值最大的元素对应的候选圆心及半径，作为瞳孔圆心),(p p y x 和半径r p 。与标准Hough 变换相比，这种Hough 圆检测投票策略可以有效提高速度。

3.1.3虹膜外边缘精定位

虹膜内边缘精定位后，可以根据此结果定位外边缘。目前主要的两种外边缘

定位方式是活动圆检测法和Hough 圆检测法。由于虹膜外边缘较模糊，且虹膜纹理非常丰富，用活动圆检测法很难准确定位虹膜外边缘。因此，最普遍的外边缘定位方式是边缘检测结合Hough 圆投票机制。

由于虹膜与巩膜间的灰度差异较小，因此一般的边缘检测算子很难准确的检测到虹膜外边界。Canny算法是目前公认最好的边缘检测算法，可以最大限度地消除噪声的影响，又不至于过多的丢失真正的边缘信息。理论上，Canny 算法通过调整参数可以检测到所有的边界，因此也是目前最常用的检测生物边缘特征的算法。

Canny 边缘检测算子的流程如图(3-6)所示，其基本步骤为：

①使用2D高斯滤波模板对采集到的图像用滤波消除噪声；

②通过微分算子来计算图像灰度两个方向的偏导数，并求出梯度的大小和

方向；

③对梯度幅值进行非极大值抑制: 首先将边缘的梯度方向大致地分为水

平、竖直、π/4、3π/4方向四种，然后通过利用各个方向不同的邻近像

素进行比较，决定出来局部极大值。若某个像素的灰度值与其梯度方向

上前后两个像素的灰度值相比不是最大的，那么该像素值置为0 ，即不

是边缘；

④用双阂值算法检测和连接边缘。使用累计直方图计算两个阂值。凡是大

于高闭值的一定是边缘；凡是小于低阂值的一定不是边缘：如果检测结

果大于低闭值但又小于高闭值，那就要看这个像素的邻接像素中有没有

超过高闭值的边缘像素，如果有，它就是边缘，否则就不是边缘。

图3-6Canny算子的基本流程图

3.2干扰检测

虹膜图像中中包含了眼睑、睫毛等噪声。如不消除这些干扰，而作为真实信息看待，势必会错误的表示虹膜特征，导致错误率增大。因此，虹膜边界定位后，