虹膜图像处理与识别技术研究

1

訇出

虹膜图像处理与识别技术研究St y rs i age pr ud on ii m oce i g ss n and r ecognii ton

齐霁，孙光民Ql iSUN Gu n . n . J a g mi (北京工业大学电子信息与控制工程学院，北京 1 0 2 0 2) 0摘要：虹膜作为重要的身份鉴别特征，具有唯一性、稳定性、可采集性、非侵犯性等优点。据统计

目前虹膜识别的错误率在各种生物特征识别技术中是最低的。本文主要是通过对虹膜图像的预处理，纹理特征提取，编码及匹配达到对虹膜所有者的身份进行验证的目的。虹膜图像的预处理采用基于灰度差以及基于H u h o S变换的边界提取方法实现；然后用Ga o滤波器对虹膜 br图像进行纹理分析，提取图像的平均绝对偏差作为特征向量；最后用加权欧氏距离 ( D wE )对虹膜图像进行身份验证。此外，本文还使用B神经网络的方法对虹膜图像的特征进行分类 P识别。实验结果表明，这两种识别方法均达到了不错的识别效果。

关键词：虹膜识别；C n y a n算子；H u h o g变换；G b r a o"滤波中图分类号：T 3 1 P 9文献标识码：A 文章编号：1 0 - 14 21 ) ( ) 0 4 - 5 9 0 ( 0 1 7上一 08 0 0 3Doi . 9 9 J is . 0 -0 4.0 1. (I) 1 :1 3 6/ . n 1 9 1 2 1 7 - . 0 s 0 3 5

0引言生物特征识别技术是指计算机利用人体所固有的生理或行为特征来进行个人身份鉴定。常用的生物特征包括：指纹、掌纹、虹膜、脸像、声音、笔迹等 (如表 l所示 )。表 1各种生物识别方法安全等级比较方法编码方式误接受率安全等级应用场所

考勤；金融和证券；电子护照和第二代身份证； 医疗保险和养老保险的发放：公安和司法；电子商务；信息安全；特殊行业等，相对其他生物识别技术具有广阔的应用前景。因此，国内外学者对其高度重视，并进行深入而广泛的研究。 虹膜识别系统由虹膜图像预处理、特征提取

及编码和匹配识别等部分构成。 虹膜图像的预处理包括对虹膜内外边界的定位，对提取出的虹膜部分进行直方图均衡及归一化操作等。而

虹膜图像纹理的特征提取本文采用了多尺度二维Ga o滤波器对图像进行滤波，再 br

虹膜识别虹膜结构指纹识别指纹纹理掌型只别手掌面积，长度、厚度面部识别脸嚣轮廓、形状和眼、 鼻分布

l 200/l0 0 01 10/ 00 l 0/70 l 0/10

高中低低

需求安全高的场所普通设施安全需求低的设施安全需求低的设施

签名识别字型、书写顺序、用笔力度声音识5声音特征； I I

l 0/10I3/0

低低

安全需求低的设施电信服务

对输出的卷积结果计算平均绝对偏差值作为特征值，最后采用两种方法对特征值进行分类已达到

虹膜识别技术作为人体生物特征识别技术的尖端代表，具有的普遍性、唯一性、稳定性、可靠性、非侵犯性、可采集性等特点，错误接受率和错误拒绝率较低如表2示，因而具有无可比拟所的优势和很高的研究价值。表2虹膜与指纹识别方法性能比较性能虹膜识别指纹识别单指精度不高，十指指纹识别可以提高精度

识别图像的目的。一种方法是计算标准图像和待测图像特征值间的加权欧式距离 ( E;另一种 W D)方法是将标准图像的特征值输入B神经网络进行 P学习和训练，再用训练好的网络对待测图像的特征值进行分类，并输出判断结果。

1虹膜图像的预处理虹膜是瞳孔之外巩膜之内的圆环部分，包含了丰富的纹理信息。为了减少眼睑、睫毛、图像拍摄角度以及拍摄距离等因素对虹膜识别效果的影响，须对图像进行预处理。本文将预处理分为虹膜定位、虹膜图像归一化和虹膜图像增强三个

准确性最准确

稳定性

一岁半后终身不变终身不变体力劳动磨损指纹，老年指纹于涩难用可以手术移植手指可受污染

抗欺诈性 最强，无痕迹残留 有痕述残留可以伪造指纹印模

抗改变性最强抗污染性 不受污染

信号采集较难采集

易于采集

虹膜识别技术重点应用于领域包括：门禁和

步骤。虹膜与巩膜、瞳孔的边界均近似为圆形，

收稿日期：2 1- 2 2 01 0- 3作者简介：齐霁 (9 7一),男，北京人，硕士，

研究方向为图像与模式识别、通信与信号处理等。 18 【8第3卷 4】 3第7期 21— ( 0 1 7上)

l匐似生可将虹膜的内外边缘近似看作两个不同心的圆。 因此定位的过程就是找出内外两个圆的半径和圆心的过程。而归一化的目的是将每幅图像调整到果，如英国剑桥大学的J h a g n授提出了 o n D u ma教一

种利用多尺度二维Ga o滤波器提取虹膜纹理的 br

识别算法，是目前应用最广泛的算法，在很多场这合均取得了较高的识别效果。而本论文也采用了 D u ma教授的多尺度二维Ga o滤波算法对图像 ag n br滤波从而提取出特征信息。 虹膜图像是基于已提取的特征向量进行识别

相同的尺寸和位置，以消除拍摄角度和距离的不同而带来的影响。增强是为了消除光照不均的影响，这里对图像作直方图均衡使虹膜的纹理信息更加清晰。

1 . 1虹膜图像处理和识别系统虹膜识别系统由光源、摄像头、图像分析系

的，不同的图像特征提取和编码方式通常有不同的比较方法，要根据具体的特征选择和编码算

统 (算机 )三部分构成，对光源、摄像头、计计算机均没有特殊的要求。图像为 8 i2 6灰度 bt5级图。光源提供照明，使虹膜可以清晰地成像。光源应稳定，因为随入射光线强度的变化，瞳孔会产生收缩或扩张，牵动虹膜变化，增大匹配时的误差。稳定光源，即稳定了虹膜，可以提高匹配的精度。

法自行设计匹配算法。但也有多种通用的判别方法，如利用汉明距离 ( mmigDi ac )、加权 Ha n s n e t欧几里得距离 ( ihe u l enD s ne Weg tdE ci a i a c )、相 d t

关系数、Ha s of距离、不同虹膜间的归一化码 ud r问距等对待测图像和参考图像的相应信息对比， 从而达到鉴定身份的目的。本文中就应用加权欧几里得距离 ( E W D)对提取的特征向量进行了识别，取得了令人满意的鉴定效果。此外，神经网络 ( ua Ne r )由于其强大的分类和自适应 Ne rl t k wo学习能力而被广泛用于模

式识别领域，本论文也尝试了应用B神经网络对虹膜进行了识别。 P 本文采用的眼部图像取自中国科学院自动化研究所下属的虹膜图像数据库 ( CASA),均为 I8 i 2×2 0 bt 0 8大小的B位图如图2 3 MP所示。

图像分析系统对拍摄下的虹膜图像进行预处理，包括对虹膜内外边界的定位，提取出虹膜部分进行直方图均衡及归一化操作等。之后从处理后的虹膜图像中提取需要的纹理特征，包括像冠、水晶体、细丝、斑点、结构、凹点、射线、 皱纹和条纹等。识别阶段即将待识别的图像的特征向量与预存的参考图像的特征向量进行匹配的过程，是典型的模式识别问题。通过比较提取的

特征向量之间的相似性，进而确定图像是否来自同一对象，以达到鉴定身份的目的。

霉图 2 CAS A的虹膜图像样本 I

该数据库中有 1 8不同人的虹膜图像，每人 0个图 1虹膜图像处理和识别系统框图

虹膜图像识别问题的关键在于特征的选取和算法选择，特征是否有足够代表性，其包含的信息量是否够多直接影响到了识别效果，进而对最后的身份鉴别起决定性的影响。目前，国内外有许多虹膜图像的识别算法均能取得不错的识别效

在不同光照和角度的照片7,分为两组，分别张用于测试和检验。本文采用的两种虹膜识别方法均用到了两组图片，WE D识别中作为标准图片和待测图片，B神经网络中用作训练样本和测试样 P本。这样提高了系统识别的复杂性，使最终结果

更加具有普遍性，更具有说服力，也使得算法更第3卷 3第7期 2 1— ( )[0 0 1 7上 41

l加严谨，客观。

甸化 的减数行 ( )与差值最小的被减数行 ( )列列 的序号 (即坐标 )的差值即为瞳孔直径，再除以

12虹膜图像的内边界定位 .眼睛的外观图由巩膜、虹膜、瞳孔几部分构成。虹膜在角膜和水样液的后方，晶状体的前

2为内边界的半径；然后用差值最大的减数行即 ( )的序号与半径作差，得到的序号 (列即坐标 )值为圆心的纵 ( )横坐标，这

样便得到了列内边界的圆心坐标。最后依照圆心坐标及半径在原始图像上画出内边界如图3示。所

方，是瞳孔之外巩膜之内的圆环部分。它包含了最丰富的纹理信息。虹膜与巩膜、瞳孔的边界均

近似为圆形。因此，可将虹膜的内外边缘近似看作两个不同心的圆。定位的过程就是找出内外两个圆的半径和圆心的过程。 本文暨采用粗定位与精定位相结合的定位方法：先用灰度闽值分割法对图像进行粗定位，得到

1虹膜图像的外边界定位 . 3虹膜与巩膜的边界灰度差值不十分明显，用

灰度阈值法求外边界的效果十分不理想，必须另想他法。H u h换是理想的替代方法，其具体 o g变

内边界参数，然后利用瞳孔的圆心及半径来缩小搜索的范围，再检测图像边缘，对其使用 H u h o g变换提取外边界参数，这种方法定位虹膜的内外边界。 由于虹膜与瞳孔的灰度差值很大，而且虹膜

方法是：首先在图像直角坐标系中定义一个圆心和半径均可变化的圆。将参数适当离散化后，枚举出图像中所有可能的圆，统计每个圆经过的外

边界点数目，数目最多的就是与外边界最吻合的圆。但 Ho g变换需要遍历整幅图像，运算量巨 uh大，而且虹膜图像细节信息大，边缘特征复杂， 定位外边界的效果不十分理想，因此需要先验信息缩小遍历范围，降低运算量，提高定位准确率。本文先用Can算子对图像进行外边界边缘检 ny测，然后根据已知的内边界参数划定外边界参数的范围，最后再结合这些先验信息用H u h o g变换定位出虹膜外边界如图4示。所

和瞳孔本身内部的灰度值比较均匀，因此本文采用灰度阈值法分离瞳孔和虹膜，具体方法是先得

到整个眼部图像的灰度直方图，再对图像作二值化处理分离出瞳孔，即选择一个阈值，将图像转化为黑白二值。变换函数形式如式() 1,其中T阈为

值，如果图像中某个像素灰度值大于该阈值，则将该像素灰度值设定为2 5 5,否则置为0。.

r 0

f Y≤T ( ) x,

…

g,=15 f,≥ (Y 5 ( y T x ) 2 x )

■

()边缘检测 a ()定位出外边界 b

图4虹膜外边界的定位

()原始图像 a

()灰度直方图 b

14虹膜图像的归一化和增强 .由于光照和拍摄距离及位置等条件的不同， 即使同一个人的虹膜图像也会出现平移，缩放或

●()分离出的瞳孔 c ()定位出内边界 d

旋转，进而对识别效果造成影响，得到的虹膜图

像不尽相同。而多尺度二维Ga o滤波器对图像的 br旋转等因素十分敏感。为了尽可能消除这个因素的影响，需要对虹膜图像作归一化和增强处理。 归一化的目的就是将每幅原始图像调整到相同的尺寸和位置，以消除平移，放缩和旋转的影响。

图 3虹膜内边界的定位

此外为了便于纹理特征的提取，还需对图像作增之后对瞳孔图像的二值灰度矩阵各行 ( )列 灰度值依次求相邻行 ( )的灰度差，差值最大列[ 0第3卷 51 3第7期 2 1- ( ) 01 7上

强处理，而增强是为了消除光照不均的影响，使

得到的纹理信息更加清晰，突出细节信息，便于

务l下一阶段的特征提取。

匐化 率、4 b角度分别对每个子图以不同频率、不同/角度进行滤波，然后计算滤波后得到的子图信息的平均绝对偏差，并将这些绝对偏差值作为特征值，排列形成统一大小的特征向量用于识别。设计思路如图6示。所

虹膜的内外边界都为圆形，为此可以利用几何方法。虹膜拉伸数学模型由于比较接近虹膜真实变化因而是个较好的选择，利用该模型可将虹膜图像的直角坐标转化为极坐标形式，并展开成矩形。对虹膜图像进行归一化处理后，再用直

方图均衡处理，以突出虹膜的特征部分，便于观察、检测。 具体方法是利用拉伸模型将虹膜图像从直角坐标系转化到[,1问的极坐标系并展开为矩形， 0】

再作直方图均衡，以增强纹理特征如图5示。所

嗣露()归一化的矩形展开图 a

()直方图均衡化的矩形展开图 b

图5虹膜图像的归一化和增强

2虹膜图像的特征提取人眼的虹膜从形状到纹理结构具有许多可供区分、鉴定的特征，为了提取

这些特征信息，需对虹膜图像进行分解，对从中提取的信息编码。图6虹膜图像纹理特征的提取

现在最常用的特征提取算法为利用Ga o函数对图 br像滤波，提取相位信息编码。二维G br a o函数是一维Ga o函数的二维拓展 br

选取这 1个子图不但能降低眼睑、睫毛对识 0别效果的影响，还可以保证有足够丰富的纹理信息被提取来；而Ga o滤波器在空域和频域具 br

有良好的联合定位能力，够提供空间频率、位置能和方向的最大分辨率，因而非常适合纹理特征分析。本文选择滤波器的频率为 2、8 l周/,、4、 6度 每个频率分别从 0。、4。、9。、 1 5处对所 5 0 3。

形式，引入高斯函数作为其窗函数。在窗形区域内，该函数与人眼视觉细胞的二维接受域波形极其相似，适合分析图像不同频率的细节信息，而通过参数调整，可以使Ga o函数具有不同的空间 br

有子图滤波，每幅子图有 1/结果图，再计算每 6b幅结果图的平均绝对偏差值，输出共 1 0特征 6个

局部化特征，如方向选择性、频率选择性和二次相位关联性等。本文提取虹膜特征的具体方法是先将虹膜环形图像等分为 1个子图，由于上眼皮处的两区域 2

值，这样能够在最大程度上保证体现图像的纹理特征。而图像的平均绝对偏差定义如式 ( )其 2,中N是图像的像素数，r是图像的均值，fxy是 n (,)在 (,)的值： xy点

经常被眼睑、睫毛阻挡，这里除去不用，只选取 1个子图，用多尺度二维Ga o滤波器从4频 0 br个 V= ( x,) m1 y一 )第3卷 3第7期

() 22 1— ( I 1 0 1 7上) 5】

、訇 I生3虹膜图像的匹配与识别在提取了虹膜图像的特征信息后，要对其进输出误差最小、经过非线性转换的信息。 利用神经网络识别虹膜的具体方法是先输入标准图像的特征值让神经网络训练网络各层的权值和闽值，当输出与理想输出的误差达到某一设

行分类，以达到识别图像的目的。本文提出了两种匹配方法对虹膜图像进行识别，一

种是计算图

像特征值之间的加权欧式距离 ( D) WE的方法，另一

定值时，认为网络收敛已训练完毕，之后输入待测图像的特征值令神经网络根据训练结果自行分类，从而达到识别图像的目的。

种是利用学习了足够特征信息的 B神经网络对 P

图像进行识别的方法。

31基于WE . D的虹膜图像识别具体方法是先将参考图片的特征向量编译

4试验结果和结论在WE D法中采用 2人的虹膜图像，每人有不 O

成. t件保存，再对测试图片提取特征向量并分 ma文别载入，用算法计算待测图像特征值与每幅标准图像特征值间的WE D值，将结果一一对比，通过算法得到最 d WE￣ D值的序号和数值。若待测虹膜特征值的W E D在k最小，并且此WE时 D值小于设定的阈值，则认为待测虹膜是第 k虹膜，从而个识别出虹膜，达到身份鉴定的目的。加权欧式距

同拍摄角度和拍摄距离的图像3幅共 6幅图像进行 0测试。从中各选取 1幅共 2幅作为标准图像，将每 O人剩余的2共 4幅图像作为待测图像提取特征幅 0值，并一一与2幅标准图像的特征值计算WE O D距离，取每次对比中最d WE￣ D距离为该算法的判断结果，如表 1示。所 在B神经网络中采用 1人的虹膜图像，每人 P 0

离 ( D)定义如式 () WE的 3,其中Ai表第i代个子图的加权系数，B N和N分别代表子图的数目和从每个

有不同拍摄角度和拍摄距离的图像7共7幅图像幅 O参与测试。每人用其中的4幅共 4幅图像提取特征 0

子图提取的特征总数。f' f和分别代表未知虹膜 i i 和第k个虹膜中第i个子图的第- i个特征分量。WE ()∑A∑( (一 l) D k=/ fi fl kj i‘ _ j]】√ il= j=l

值作为训练样本输入到神经网络中进行训练，再用其余的3共3幅图像输入到训练好的神经网络幅 0中进行分类识别，结果如表3示。所表3两种方法对虹膜图像的识别结果样本数正确识别的正确率错误识别的错误率

————■

(1 3、

32基于神经网络的虹膜图像识别 .人工神

经网络是指利用工程技术手段模拟人脑神经网络的结构和功能的一宗技术。它的目的是使机器有人脑那样的感知、学习和推理能力。人工神经网络模型具有与人脑相似的高度并行WE D法

( 幅)4 0

图像数 ( )幅3" 7

( )%9 . 25

图像数 ( )幅3

( )%7. 5

B神经网络 P

3 0

2 3

7 . 67

7

2 . 33

结合试验结果得出结论：

性，良好的容错性和联想记忆功能、自适应和自 学习能力、强大的分类能力、鲁棒性，可以实现输入到输出的非线性映射，因此各类神经网络在图像识别领域得到越来越多的应用。 本文使用的B神经网络，是一种采用反向传 P播算法的多层前馈网络，能够逼近任意的非线性映射关系，而且有很好的泛化能力。BP络模网型处理信息的基本原理是：输入信号通过中间节

试验结果表明，用两种方法对虹膜图像进行

识别的效果都比较理想，WE D法的计算结果清楚明了，便于比较，准确率高，但判断耗费的时间

较B神经网络长；神经网络需要训练特征，对样 P本数量的需求比W E D法大，但判断时间短，输出结果简单直观，便于分析。若改进神经网络的输入样本特征类型，有望提高神经网络识别虹膜图像的准确率；而WE D法受图像的拍摄角度和拍摄距离的影响较大，若要减弱这种影响，提高识别的准确率，需要对提取特征的算法作进一步的改善。

点作用于输出节点，经过非线性变换产生输出信号。网络训练的每个样本包括输入向量和期望输出量，通过分别调整输入层与隐层间、隐层与输出层间的参数 (权值和阈值 )使误差沿梯度方向下降，经过反复学习训练，确定与最小误差相对

5结束语本文介绍了虹膜图像处理和识别系统，对各阶段从图像预处理到最终的虹膜识别用到的方法进行了研究与探讨。在图像预处理方面，选用灰度阈值和Ho g变换的方法分别对虹膜内外边界 uh

应的网络参数，训练即告停止。此时经过训练的神经网络既能对类似样本的输入信息，自行处理[2第

3卷 51 3第7期 21- ( ) 0 7上 1

、 l定位，并实现预期功能；通过算法实现了对虹膜图像归一化以及直方图均衡的处理以突出其纹理特征。在特征提取方面，选取了业内著名的 Da g n授的多尺度二维 Ga o滤波器法对 u ma教 br虹膜图像滤波，编程实现了提取滤波后图像的信息；选择滤波后图像的平均绝对偏差作为特征

訇化 习、训练样本特征值，进而对虹膜图像识别，鉴别效果达到预期。

参考文献：【]庞尚珍，跃荣， 1雷冯雪，艳 .膜识别技术研究[1川理杜玲虹 J四 .工学院学报，0 8 1 ) 59 . 2 0, (:—7 2 59 []黄惠芳， 2胡广书.虹膜识别算法的研究及实现[ _外与激 J红】光工程，0 2 5:0—0 . 2 0,1 ) 44 9 3( 4

值，利用算法将其编译成特征向量保存，为识别阶段做准备；在特征识别方面选用了两种：:号法 用加权欧几里得距离( D)虹膜图像的特征向 WE对量识别，鉴别效果十分理想；使用B神经网络学 P

[]孟浩，翠平 .识别算法的研究【 .尔滨工程大学学 3徐虹膜 J哈 J报，0 6 73: 0—0/5 . 2 0, ( )4 0 4 3 2 2 4

[】史鹏举 .光民.膜识别系统的研究[ _算机应用研究， 4孙虹 J计】2 0,3增 ) 0— 1 . 0 42 (: 9 5 5 1

.

.

{是I每赢I{&‘竞‘是I{品‘霸 .|出I翕‘ .‘{舞‘{&‘{高●岛‘意{离I翕‘赢●出{惫‘禹●毒‘ .翕

【上接第4 页】

点，根据这些期望路径点在低层关节位置伺服中实现关节的位置伺服控制，同时还设置了一定的系统保护功能，如当前实际位置值与规划点的误

分利用了DS运算的高速性，满足机械手控制的实 P时需求。

差过大、电机的转速过大或出现错误时将关掉伺服。图6关节位置伺服控制流程。为表2相关参数定义参数参数定义电机输出 E,

4结束语所开发的禽蛋吸运机械手具有以下特点： 1 )采用“业P工 C+ DS P运动控制器”分布式二级控制结构

,增加了系统的开放性、快速性及可靠

参数 K PK f

参数定义比例增益积分增益

性；2 )采用真空吸盘装置作为末端执行器，具有良好的缓冲性能，满足禽蛋移动的特殊性要求； 3 )更换机械手末端执行器，可以实现其他作业任务，具有良好的通用性；4 )控制软件采用 C十+编程，通过调用底层函数对硬件进行直接操作， 可视化环境提供良好的人机交互性能。经实践证

第 1个采样时刻的位置 3偏差第 n个采样时刻的目标位置和实际位置

K D

微分增益

Ⅲ 』

Y E积个样的 K 速前馈 第n采时刻误差 F度分值增益

明，很好地完成了禽蛋品质检测的分级吸运任务。

当前期望速度

BI AS

电机偏差

定时中断服务程序

参考文献：【】马香峰， .业机器人的操作机设计[ .京：金工 1等工 M】北冶业出版社，9 6 19 .

[]蔡自兴 .器人学[ .京：华大学出版社，0 3 2机 M】北清 20. []张兴国，海黎 . UCM一i型工业机器人结构及运动 3徐 F AN 6B学分析【 _通大学学报 ( J南】自然科学版 ) 0 91: 1., 0 ()—2 2 9

【]张兴国 .保压缩机装配机器人的运动学分析[】 4环 J.大南通学学报 (科学版 ) 0 41:—43 .自然， 0 () 23, 2 3 8[]S 5 MC(国)中有限公司 .代实用气动技术[ .京：械现 M]北机工业出版社，0 8 . 2 0, 8

【】王天然，奎 .业机器人控制系统的开放体系结构【 . 6曲道工 J J机器人，0 22 () 5— 6 . 2 0, 3: 62 1 4 2

图6关节位置伺服控制流程

[]深圳摩信科技有限公司. T 0 0列控制器使用手册 7 MC 80系[】圳：圳摩信科技有限公司，0 1 z.深深 20.

采用“业P工 C+DS运动控制器”的分布式 P控制结构，结合“ I速度前馈控制”算法，充 P D+

【】严晓照，兴国.量式 PD制在温控系统中的应用[】 8张增 I控 J.南通大学学报(自然科学版 ) 0 64:—1, 0 () 85 . 2 4

第3卷 3

第

7期

2 1- ( )[ 3 0 7上 5 1 1

本文来源：https://www.bwwdw.com/article/kni4.html