0 引言 基于人眼虹膜信息特征的生物识别技术具有独特性、稳定性、防伪性、非接触性、高精确性等优点,在安全验证系统、信用卡验证、医学、档案管理、视频会议、人机交互系统等方面具有巨大的应用前景。Daugman[1]和Wildes[2]等人提出的识别方法由于对噪声敏感、要求硬件捕获高质量的虹膜图像而难以在实际中广泛应用。本文提出了一种基于修正的精确定位和一维多尺度小波变换下提取归一化虹膜图像特征及基于SVM的WED与DSIM虹膜分类识别ISD算法,实验结果表明该算法具有较高的识别效率和精度。 1 虹膜的精定位与修正 1.1 精定位 虹膜与巩膜、瞳孔的边界均近似为圆形。瞳孔内部灰度分布均匀且与虹膜的接界处存在极大灰度梯度;巩膜与虹膜的接界也存在较大的梯度,由此可以检测出内外圆形边界。由于虹膜具有良好的环状特性,所以可采用如下定位算子对虹膜内边缘进行精定位 :
(1-1) 其中 ,是高斯函数,起到了平滑的作用,可以去除图象中部分噪声干扰。虹膜内边界检测出来后,通过 及 的值限定合适的范围,采用同样的方法可搜索虹膜的外边界。由于虹膜外边缘一致性较差(主要是指虹膜外边缘偏离圆形,以及它与巩膜之间的过渡带较宽),所以可能使某些虹膜上存在的非常明显的向心性的环状纹理。为防止式(1-1)的算法掉入陷阱而失去作用, 用面积分式(1-2)代替式(1-1)中的线积分。
(1-2) 另外,在人眼正常凝视状态下,上下眼睑可能会遮盖一部分虹膜,进行虹膜外边缘精定位时,对其积分域必须进行限定,积分域一般限定在 到 及 到 之间的扇形区域内。 1.2定位修正 虹膜定位时基于圆模板搜索得到的灰度梯度值为在该半径下的模板圆与虹膜内边缘或是外边缘相符的最大圆弧处的灰度梯度值, 这就带来了不可避免的误差。为使模板圆的圆心尽量靠近虹膜边缘的几何中心,就应尽可能多的采集到真实的虹膜边缘点。虹膜与巩膜之间边界对应于灰度曲线上的跳变边沿,采用如下的灰度差分累加和的方法检测。它是一种简便、快捷、稳健的波形边沿检测方法。

图 1-1 园模板 设灰度曲线上的点表示为 , 。由灰度曲线,经过差分运算,便可求得灰度差分曲线( )。差分运算如下:

在灰度差分曲线上,m个相邻点的累加和为:
(1-3) 边沿的确定可以采用阀值的判断法。设定一个阀值,对于右边沿,差分的累加和大于阀值,即 ;对于左边沿,差分累加和的负值大于阀值,即- 。边沿确定后,取 为边界点。由于人的虹膜与巩膜边界上有一个一定宽度的模糊带,因此在灰度曲线上所反映的跳变是有不同的坡度的。为了更准确更省时寻找虹膜边界点,可以采用如下方法:设定一个合适的阈值,以精定位确定的边缘点为中心,设对应的 ,则
(1-4) 在灰度差分曲线上,2p个点的差分累加和为: (1-5) 然后,变化k的取值,k=k+q , ,检测到差分累加和大于阈值时就把该 值列为候选点,在所有候选点中,值最大者选为最终的边界点。 2 虹膜图像归一化 虹膜的内外圆形边界在许多情景下不是同心圆,同时由于瞳孔的牵引使得虹膜伸缩变换,因此必须把虹膜区域经过坐标变换并归一化展开成极坐标下的带状信号以便后续研究。 如图:外圆半径 , 圆心 O(Ox,Oy), 内圆半径 , 圆心I(Ix,Iy)

图 2-1虹膜的归一化 设 上一点P(X,Y),归一化 ,由图有:
(2-1) 由此原图中的点对P(X,Y)一一映射到极坐标中的点对P(r,θ)。由于此坐标变换实质上是对不同的θ方向采用不同的量化间距 ,因此对内外圆的伸缩具有不变特性,同时具有平移和尺度不变特性。由于虹膜纹理集中分布于内环边缘, 实验中取瞳孔附近环带的70%, 量化级Nr=64,Nθ=256,θ=90º为归一化时θ轴的原点。 3 特征提取与分类识别 3.1 基于Haar小波基的特征提取 我们对归一化的虹膜图像 在角向作N次一维小波变换( ),以变换后的图像和原图像的差值 作为待识别图像, 保留了虹膜的纹理结构特征并且减少了光照和噪声的影响。对待识别目标 提取其角向相邻两行的均值与方差 以及尺度为16*16的计盒维数 作为分类的二组有效特征量向量 。 标志着图像表面的粗超程度和自相似构造规律,表明了虹膜的纹理信息。
(3-1) 其中膨胀函数 和尺度因子S之间的关系为:
(3-2) 函数 在尺度S和 处的一维小波变换可写成下列卷积:
(3-3) 实验中取正交归一的Haar小波基,对应的滤波器系数: (3-4)

图 3-1 Harr一维小波变换 3.2 特征训练 根据实验对某类的 个训练样本特征向量 (  )进行如下训练以获取该类特征。
(3-5) 3.3 分类器的设计 本算法采用基于支持向量机SVM的相异性测度DSIM 和倒数加权欧式WED距离的ISD分类。设有模式 ,则
(3-6)
(3-7) 拒识判定准则:设模板数为m,待分样本到第i个模板的WED距离di; dj和 dk为最小的两个距离。定义拒识系数 。设 为拒识阈值,若 则表明该样本归属清楚,直接送入WED距离分类器;否则拒识,对于获得较小WED的几个模板计算DSIM值并利用上述准则判定,如仍拒识表明该样本应送入SVM分类器分类。. 支持向量机SVM: 对于线性可分的集合S,如果存在 对于 由条件 定义超平面 ,使得分类间隔 最大,转化为二次规划问题,定义Lagrange函数( 为与样本对应的Lagrange因子):
(3-8) 用核函数 代替样本内积 由极值求法可得到最优化分类函数:
(3-9) 设测试集H,确定前两级分类器的阈值 ,自动确定参数 ,则ISD算法如下: 1) 若 ,则算法结束,输出结果.否则取样本 ; 2) 计算样本特征向量V1并求其到模板的WED距离,调用拒识判定准则; 3) 将未被拒识的样本送WED分类器分类,对于拒识样本计算V2特征并求其到模板的DSIM距离, 调用拒识判定准则: 4) 将未被拒识样本送DSIM分类器分类,否则送往SVM分类器分类. 5) ,转第1步. 4 实验与性能分析 4.1 虹膜预处理

图 4-1虹膜的定位 归一化及Harr一维小波变换 4.2 分类测试 实验中提取不同个体在两个不同阶段共423个64*256的虹膜图像样本,第二阶段样本经训练后行成89个模板,类内测试185次,内间测试1408次,性能如表1.从表中可知随着 的增大,识别率相应提高,但由于前两级拒识的样本增加SVM分类的时间代价也提高,综合考虑 = =0.4时具有较佳的性能.

图 4-2 单纯WED与 DSIM分类示意图 表1 不同 下ISD与单纯WED和DSIM分类器性能比较

4.3 结论 从实验结果可知本文提出的算法很好的实现了虹膜的定位与分类识别,误识率低,具有比较高的识别精度和效率,其性能优于单纯采用SVM,DSIM或WED距离度量的方法
|