虹膜图像预处理

虹膜图像预处理在实际获取的图像中，虹膜只占一部分区域。图像中还包括瞳孔、眼睑、睫毛及眼睛和人脸的其它部分。所以，在对虹膜纹理特征进行提取之前，要把虹膜区域从图像中分割出来。其分割准确与否直接影响到后面的特征提取及匹配精度，即如果分割不准确，则虹膜识别效果就差，很可能造成误判，使识别算法失效。针对Canny[2]算法优良的边缘检测性能及虹膜边界的特点，本文采用了该算法，与Hough变换相结合，用于定位虹膜内外边界、分割上下眼睑，采用阀值法剔除睫毛。2.1虹膜内外边界的定位虹膜的内边界即内边缘与瞳孔相交接的部分，外边界即外边缘与巩膜相交接的部分，均近似为圆形，因此定位虹膜内外边界就是要确定这两个拟合圆的圆心和半径。2.1.1Canny算子和Hough变换的基本原理Canny边缘检测算子边缘是图像的基本特征，保留了原始图像中相当重要的信息，而又使得总的数据减小了很多。准确提取边缘既能减少图像处理的信息量，又可重点描述物体的形态特征。传统的边缘检测算子如Roberts>Sobel、Prewitt、Kirsch和Laplacian等算子都是局域窗口梯度算子，由于它们对噪声敏感，所以在处理实际图像中效果并不理想。1986年，Canny提出边缘检测算子应满足以下3个判断准则：信噪比准则；定位精度准则；单边缘响应准则，并由此推导出了最佳边缘检测算子--Canny算子。该算子具有定位精度高、单一边缘和检测效果好等优点。它主要分以下四步：(?平滑图像Canny边缘检测算子是高斯函数的一阶导数，是对信噪比与定位精度之乘积的最优化逼近算子。选择合适的高斯低通滤波器来平滑原始的虹膜图像，既可以取得更加柔和的平滑效果，又不至于破坏图像中的边沿，同时还可以消除原始图像的高频噪声。Canny算法首先用二维高斯函数的一阶导数，对图像进行平滑，设二维高斯函数为：

rv、 1 ，尤2+y2G(x,y)=-——exp(———) (2.1)271(72 2b2其梯度矢量为：(2.2)dG/dy把G的2个滤波卷积模板分解为2个一维的行列滤波器：(2.3)(2.4)6G X2 (2.3)(2.4)__=^exp(____)exp(-—)=Vx)/Z2(J)—=Zryexp(--一)exp(--一)=h(y)h(x)dy 2c>22b2i2其中，k为常数，。为高斯滤波器参数，它控制着平滑程度。。小的滤波器，虽然定位精度高，但信噪比低；。大的情况则相反，因此要根据需要适当地选取高斯滤波器参数。。。计算梯度的幅值和方向传统Canny算法采用2x2邻域一阶偏导的有限差分来计算平滑后数据阵列的梯度幅值和梯度方向。其中，x和)方向偏导数的2个阵列P[ij]和尸[/"]分别为：pL,j]—(/L,j+1]_/L,j]+/L+1,j+1]_/L+1,j1)^2 (2.5)XPL,j]=(/L\j]-/L+1,j]+/L,j+1]-/D+1,j+1])/2 (2.6)y像素的梯度幅值和梯度方向用直角坐标到极坐标的坐标转化公式来计算，用二阶范数来计算梯度幅值和方向，分别为：M",j]=JpQ,j]2+Pkj]2 (2.7)'x y0Q,j]=arctanGkj]/PLj]) (2.8)yfx@对梯度幅值进行非极大值抑制为了精确定位边缘，必须细化梯度幅值图像mL,j]中的屋脊带，只保留幅值局部变化最大的点，这一过程就是非极大值抑制。如图2-1。在非极大值抑制过程中，Canny算法使用3x3大小、包含8个方向的邻域对

梯度幅值阵列ML,/•］的所有像素沿梯度方向进行梯度幅值的插值。在每一个点上，邻域的中心像素m上j］与沿梯度方向的2个梯度幅值的插值结果进行比较，fkj］是像素邻域中心处沿着梯度方向的扇形区域，非极大值抑制在此区域进行。如果邻域中心点的幅值mkj］不比梯度方向上的2个插值结果大，则将mkj］对应的边缘标志位赋值为0，这一过程把M上j］宽屋脊带细化为一个像素宽，并且保留了屋脊的梯度幅值。非极大值抑制(Non-MaximaSuppression,NMS)过程的数学表示为：Nkj］=NMS(ML,j］,匚L,j］) (2.9)图2-1Canny算子非极大值抑制Q检测和连接边缘双阈值算法是对经过非极大值抑制的图像N上j］分别使用高、低2个阈值S,h和S7分割得到2个阈值边缘图像Tkj］和T上j］。由于图像Tkj］是由高阈值得l h l h到的，因此它应该不含有假边缘，但T上j］可能在轮廓上有间断。因此双阈值算l法要在TL,j］中把边缘连接成轮廓，当达到轮廓端点时，该算法就在由低阈值得h到的边缘图像T上j］的8邻域位置寻找可以连接到轮廓上的边缘，这样，利用递l归跟踪的方法不断地在T上j］中搜集边缘，直到将［kj］中所有的间隙都连接起来为止。Hough圆检测Hough变换是一种用于区域形状描述的方法，经典的Hough变换常常被用于直线段、圆和椭圆的检测，广义的Hough变换可以推广至任意形状。两种变换的基本思想是将图像的空间域变换到参数空间，用大多数边界点满足的某种参数形式来描述图像中的曲线区域边界。Hough变换可以用于检测图像解析曲线f3,a=0，其中X为解析曲线上的点，a为参数空间上的点。对于圆，设其半径为,、圆心为(a,b)，则圆方程：(X-a)2+(y-b)2=r2 (2.10)图像空间的圆对应着参数空间(a,b,r)中的一个点，一个给定点(x,七)约束了通过该点的一簇圆参数(a,b,r)。传统的Hough圆检测是将空间域的每个轮廓点带入参数方程，其计算结果对参数空间(a,b,r)中的量化点进行投票，若投票超过某一门限，则认为有足够多的图像点位于该参数点所决定的圆上。由于其需要逐点投票、记录，故花费的时间较长，占用的计算机存储量也较大。2.1.2本文采用的Canny算子如图2-2(a)所示的虹膜图像，很多边缘对于检测虹膜内外边界来说都是干扰边缘，如果在计算梯度时，对x、y两个方向均进行计算，则会产生过多过细的边缘，而这些边缘又是不需要的，而且会对提取有用的边缘起干扰作用。此时，如果有选择性的重点检测某一或某些方向的边缘，则既能节省时间又有利于排除干扰。针对虹膜图像，由于上下部分常常会被眼睑遮挡，如果将眼睑边界也过清晰地检测出来了，则对下一步的Hough圆变换检测外边界产生很大的干扰，特别是眼睑基本成弧形，与圆形接近。为削弱非虹膜边界的影响，可以在进行边缘检测时重点检测有用的边缘，在保留足够的虹膜边界信息的同时尽量削弱其它边界的影响。注意到上下眼睑边界仅限于水平方向，在计算梯度幅值时对水平方向偏导值赋以较大的权重，而对竖直方向偏导值赋以较小的权重，从而削弱水平方向上的边界。本文在定位虹膜外边界时，按照这一思想对Canny算子进行了改进，在计算梯度时只计算水平方向上的梯度，以避免过多地检测出眼睑边缘。即：X方向偏导数：

pL,ji=(iL,j+1]-1L,j]+1L+1,j+1]-1L+1,j]),2 (2.ii)r方向偏导数：(2.12)P\i,j]=0(2.12)y(b)竖直权重为0(a)原始图像缩放(b)竖直权重为0(a)原始图像缩放0.4(c)两方向均匀权重 (d)水平权重为0图2-2虹膜图像以及Canny边缘检测的图像其中滤波时对边缘处，为简化，未处理，造成有虚假边缘，定位时除去，并不影响Hough定位圆。在非极大值抑制时，选择比较邻域为半径1.5像素范围之内，这样的到的边缘图像相对较好，边缘大于等于一个像素，小于等于两个像素，考虑到Hough参数的量化，比选择半径为1的邻域范围造成的误差小，而与选择半径2的邻域边缘图相比，边缘得到了较好的细化。如图2-3。(a)1.5邻域半径抑制(b)2邻域半径抑制(c)1邻域半径抑制(a)1.5邻域半径抑制(b)2邻域半径抑制(c)1邻域半径抑制(d)1.5邻域半径抑制(e)2邻域半径抑制 (f)1(d)1.5邻域半径抑制图2-3不同邻域半径边缘检测图(3)此外实际处理过程进行了少许改变：Q为了减小处理的数据量，在定位虹膜内外边界时，首先对虹膜图像进行了缩放，而且缩放定位后再折换成原始图像的圆，原始图像大小为320x280。Q对于Canny检测处理后的边缘图像，直接一次将内外两圆定位出，在定位外圆时比较准确，但是定位内圆时由于受其外围边缘的影响，很容易造成定位失败，因此首先定位出外圆，然后在外圆的范围中定位内圆，大大减小了干扰。@由于内圆相对较小，定位外圆时缩放比例0.4，定位内圆时缩放比例为0.6。2.1.3本文采用的Hough变换圆检测算法直接使用Hough变换是将图像上的每个边缘点都代入方程，且半径搜索范围从O到图像平面所能容纳的上限，因此计算量非常大。本文对其进行了修改，并根据先验知识及所采集的虹膜图像的特点尽可能地减少搜索范围及参与Hough变换的点数以高效率。过程如下：对图像进行边缘检测得到边缘图像；根据先验知识，分别确定内外边界圆半径的范围，从而减小搜索半径；圆的参数方程可改写为：a=x-rcos0,b=y-rsin0 (2.13)将图像空间中的边缘点而不是原图像中的每个点逐一代入上式，求出参数(a,b)值。从图中可以看出，虹膜外边界左右两侧的部分所受干扰较少，边界质量较高，因此可以限制0的取值范围，相当于只统计左右两侧部分的边界点。由于虹膜边界圆的圆心不可能太靠近图片边缘(否则就没有将整个虹膜区域采集进来，虹膜信息可能太少而无法用于识别)，即圆心应位于图像中间区域的某个范围内，因此如果(a,b)位于这一范围内，则将相应的累加阵H(a,b)中的元素加1，否则加0;找出H(a,b)中元素的最大值，即是对应半径为r，圆心为(a,b)，且圆周上边点最多的圆，即为边界圆。2.1.4虹膜内外边界的定位本文先定位虹膜外边界，然后再在外边界范围以内定位内边界，从而能更快更准地定位内边界。虹膜外边界上下部分常常被眼睑部分地遮挡，且往往还有睫毛的干扰，所以对外边界的定位先采用本文所述的改进的Canny算子进行边缘检测，从而在保留足够多的边缘信息的情况下尽量少地检测出干扰边缘，再利用改进的Hough变换进行圆检测，确定外边界圆的圆心和半径。由于瞳孔位于虹膜区域以内，所以定位内边界时，将搜索范围限制在外边界范围以内，以减少搜索范围。定位结果如图2-4图2-4虹膜定位结果2.2眼睑的分割在虹膜识别过程中，虹膜区域往往会受到干扰，由于尽量考虑实用性，不能对被采集者提出苛刻要求进行配合，这样在被采集者不配合的情况下，获得的虹膜有效区域是很有限的，往往会是部分虹膜区域，其余部分被眼皮、睫毛等遮挡。从采集的图像数据库中也可以看出，多数图像中虹膜区域都是不完全的，因此虹膜区域不完全是普遍存在的。眼睑包括上眼睑和下眼睑，上下眼睑开度不同，裸

露出的虹膜区域面积也不同。如果这部分眼睑也作为有效虹膜区域带入虹膜特征提取和编码模块，那么得到的虹膜编码可信度将大大降低。检测眼睑比较常用的方法有抛物线检测眼睑和直线检测法。在许多文献中，都采用了抛物线检测眼睑。MathieuADAM等在文献中首先在保留虹膜边界的同时，平滑了虹膜纹理，然后进行抛物线检测；Tae-HongMin等提出在归一化后的虹膜图像中采用抛物线检测，这样不需要考虑0的影响；YoungKyoonJang等提出检测眼睑候选区的曲线累积值，如果大于某一阈值则进行抛物线检测，否则进行直线检测。2.2.1抛物线检测眼睑检测眼睑时，可以根据眼睑的形状采用抛物线拟合的方法，在多维参数空间求最优解。抛物线曲线的定位可通过广义Hough变换来确定抛物线的各个参数。抛物线曲线可表示为：(-(x-h)sin0+(y-k)cos0)2=a((x-h)cos0+(y-k)sin0) (2.14)其中a是用来控制抛物线的曲率，(h,k)是抛物线的顶点，0是抛物线相对于x轴的旋转角。该方法的缺点是运算时间长。另外较常用的是将眼睑简化为直线的方法。这样可以大大的降低运算量。2.2.2Radon变换法分割眼睑对眼睑边界进行细分，可以近似认为其由一系列的线段连接而成。因此可以采用Radon变换直线检测法定位眼睑。Radon变换可在任意维空间定义，下面给出在2维空间的定义式:(2.15)P为坐标原点(2.16)R(p,0)=jjf(x,y)5(p-xcos0—ysin0)dxdy(2.15)P为坐标原点(2.16)D式中，D为整个图像平面，f(x,y)为图像在点(x,y)的灰度，到直线的距离，0为距离与x轴的夹角，5为冲激函数，定义为:3(x)=［顽离散时为1)x=°|0 x丰0它使f(x,y)沿直线p=xcos0+ysin0进行积分。如图2-5。Radon变换可以理解为图像在p-0空间的投影，p-0空间中的每一点对应着图像空间中的一条直线，而Radon变换是图像像素点在每条直线上的积分，也可理解为图像逆时针旋转0角度后在水平轴上的投影。因此图像中每条直线会在P-0空间形成一个“亮点”，直线的检测转化为在p-0变换域对“亮点”的检测。图2-5Radon变换法检测眼睑的步骤如下：分别确定可能存在上、下眼睑的虹膜区域，以减少搜索范围；对可能存在眼睑的虹膜区域采用改进的Canny算子进行边缘提取。由于眼睑边缘呈水平方向，因此在计算梯度时只取竖直方向的梯度，以避免虹膜边界的干扰，即：x方向偏导数：PL,j]=0 (2.17)尤Y方向偏导数：P上j]=(小 +1，j]+1上j+i】tG+1，j+山2 5对提取到的边缘图像进行Radon变换直线检测，变换后会得到若十个“亮点”，如果“亮点”的最大值大于某个阈值T，则视其对应的直线为眼睑边界直线L，否则认为没有上(或下)眼睑遮挡；如果检测到了眼睑边界直线L，则作一条水平直线，此直线过虹膜外边界与L的交点，为简化处理，将此直线所在的虹膜区域标为干扰点。对检测到眼睑像素用灰度0代替。眼睑检测结果如图2-6。可以发现采用这种方法能比较有效地将上下眼睑分割出去但检测方法过于粗略，易将部分有用虹膜当作眼睑分割出去，因此损失了一定量的有用虹膜信息。2.3剔除睫毛睫毛是用来保护眼睛的，可以有效阻挡异物进入眼睛，但是在虹膜识别时，落入虹膜图像区域的睫毛影响了虹膜特征的提取，如果不加以消除，会造成特征表示精度下降，难以反映实际的虹膜信息。由于虹膜编码识别算法本身具有一定的抗干扰能力，当睫毛稀疏且很短时，对