专利名称：指纹识别方法
技术领域：
本发明涉及一种生物识别技术，具体涉及一种指纹识别方法。
技术背景由于人的指纹千差万别，几乎没有重复，因此指纹识别是身份识别的一种 重要方法。指纹识别技术是指对采集到的指纹图像信息进行计算机图像处理和 模式识别的一项综合性技术。指纹识别的大致过程是首先通过指纹采集设备获取指纹图像，然后对该 图像进行预处理，以方便提取指纹特征，特征提取后，再通过模式匹配的方法 判别两者的一致性，若特征值相似，则可判定为同一人的指纹，反之，则判定 两指纹不是来自同一枚手指。现有指纹识别方法一般仅考虑指纹的细节点局部特征，普遍存在识别率低、 识别速度慢等问题。发明内容本发明的目的在于克服现有技术的上述缺陷，提供一种识别率高，识别速 度快的指纹识别方法，通过分析细节点的邻域关系和脊线曲率构建稳定性更好 的特征矢量，然后利用匹配的细节点对总数和匹配分值之和构建多级判别体系。本发明用到的一些定义含义如下指纹的细节点是指指纹脊线上的端点或者分歧点；脊线间的平均距离X:是指细化后脊线间的平均宽度，用象素的个数表示， 对于大小为300X300，分辨率为500DPI的指纹图像，脊线平均距离为10个象 素，也即A =10;细节特征的方向根据特征的类型不同，其方向的计算方法也有差异。对 于端点，它的方向可由以下方法计算而得从端点出发，沿脊线跟踪入个象素， 用一元线性回归求出这些点的拟合直线，该直线的方向就是所求特征端点的方 向。而对于分歧点特征，它的方向可由从分歧点出发的两条同方向脊线的方向 角的平均值确定，而脊线的方向角可通过上面介绍的一元回归线性拟合的方法 求得。特征点间的距离d:是指细节特征点间的欧氏距离。假定两个特征点的位置坐标分别为(xl， yl) ， (x2， y2)，则距离d为特征点的结构数Nm:是指以细节特征为中心，R为半径的圆或正方形内细 节特征的个数。下标m可用e或b取代，分别表示端点和分歧点的个数。例如， 可以采用正方形模板，半径为3入/2。特征点垂直对比度Varo:该参数的计算需要借助于没有进行二值化的灰度 图像。在16X16的灰度图像块中，选取中心象素作为待测点。根据预处理确定 的块方向，过待测点作该方向的垂线，沿垂线在待测点两边各取4个象素点， 则这8个象素的灰度方差就被认为是待测点的垂直对比度，为方便起见，这一 块中所有象素的垂直对比度都等于中心点的垂直对比度，那么如果细节特征所 处的图像块已知，则其垂直对比度也就可以测得。特征点平行对比度Varp:该参数的计算类似于垂直对比度，仍需要借助于 没有进行二值化的灰度图像。在16X16的灰度图像块中，选取中心象素作为待 测点。根据预处理确定的块方向，沿该方向在待测点两边各取4个象素点，则这 8个象素的灰度方差就被认为是待测点的平行对比度，为方便起见，这一块中所 有象素的平行对比度都等于中心点的平行对比度，那么如果细节特征所处的图 像块已知，则其平行对比度也就可以测得。以上定义在论文"指纹细节特征提取与剪枝"(王崇文，李见为等，光电 工程，2002， No. 4， pp67-70)中已有详细定义。在该论文中提出的指纹图像的 表示方法为将指纹图像表示为一个二维矩阵，每一个象素就是矩阵的一个元素， 取值为0 255，矩阵的维度就是图像的宽和高。在本发明中，继续采用这些定 义和图像的表示方法。本发明提出了一种指纹识别方法，由指纹特征提取和特征匹配两个过程组成所述指纹特征提取的过程为针对于预处理的图像进行细节特征初选；对 初选后的细节特征进行剪枝；对剪枝后的细节点特征进行可靠性验证；最后根 据细节点的邻域关系和脊线曲率特征构建特征矢量；所述特征匹配的过程为首先利用脊线的曲率特征进行匹配参考点对的识 别；然后利用极坐标对图像进行校准；接着利用自适应弹性匹配方法进行细节 点配对和匹配分值计算；最后利用匹配的点对数和匹配分值，采用多级判别体 系进行指纹匹配。一、指纹特征提取 所述指纹特征提取的具体过程为1、 指纹细节点初选初步选出所有的细节点，记录下这些细节特征的类型 和位置坐标，以便于下一步的剪枝处理。2、 细节点剪枝对于细化后的点线指纹图，所述细节特征剪枝的伪特征包括因为脊线断 裂产生的两个错误端点；因为脊线的错误连接和断开产生的伪分歧点和伪端点； 在对非平滑脊线进行细化后产生的短枝；因为脊线间的错误连接产生伪分歧点。 可以针对每一类细节伪特征分别进行剪枝。删除了伪细节特征后，可以进一步根据指纹图像的拓扑性质来对存留的细 节特征进行可靠性检测。由指纹图像可以看出，分歧点侧面的脊线与形成分歧 点的脊线分支平行，可以基于这个特性对分歧点的可靠性进行检测。3、 构建特征矢量本发明构建特征矢量的方法为在本发明中，为了减少指纹形变对特征匹 配的影响，定义了一种新的细节特征矢量，它总共包含24个元素，表示方法如 下FC = {jc， ；/， cy， var， /,c^y，厂，c/1， "1， sl， /1， c/2， "2， s2， /2， d3, "3, s3, /3， t/4， "4， s4， /4}各元素的含义为mO是当前待研究的细节特征点，即中心细节点，它的方 向为o0 (该方向值可以在预处理阶段获得)。以m0为坐标原点，以o0为x轴的正 方向，构建坐标系。ml， m2和m3是分别位于第一、第二和第三象限的细节特征 点，这些特征点应该是在各象限所有细节点中距离mO最近的。oi是各细节点mi 的方向，x和y是中心细节点mO在原来坐标系中位置坐标，"是包含中心细节点 的局部脊线方向，即图中的o0。考虑到特征提取过程中误差，可以用参数var和f 取代类型参数，var是在原始指纹图像中以mO为中心的局部窗口内各象素的灰度 方差，f是在原始指纹图像中包含mO的局部窗口的频率，它代表局部窗口中脊线 的频率。这两个参数都可以在指纹预处理阶段获得。如果仔细观察指纹图像，就会发现大多数的脊线在有限长度内可以当作是 圆弧的一部分，基于这一点，本发明使用了3个参数来表示与细节特征点相连的 脊线，它们分别是cx， cy和r。把与细节特征点相连的脊线看作是一段圆弧， 利用圆弧上的三个点就可以得到所在圆的中心坐标和半径。考虑到弹性形变的 影响，利用最小二乘法去拟合脊线，最后用中心的平均坐标cx 、 cy和半径的平 均值r来表示与细节特征点相连的脊线。di,si,ai和ii(i-l，2,3，4)是中心细节特征点的邻域关系特征，它们描述了中心细节点(m0)与周围细节点(mi)之间的关系。di分别是mi到mO的距离，ai 是mO方向角所在直线与mi方向角所在直线的夹角，si是mO与mi连线所在直线分 别和mO方向角所在直线的夹角。il， i2， i3和i4分别是mi特征点在所有细节点集 合中的索引。考虑到在某些情况下无法为某个象限确定合适的细节特征，就用 一个特殊的极小值来标识d4， s4， a4和i4，如一l。注意到在特征矢量FC中，各分量特征都来自于相对量或统计量，因此该特 征矢量具有旋转不变和平移不变等特性。另外，对于任意的细节点，它们的特 征矢量FC的长度都是一样的，这为我们计算矢量之间的相似性带来很大的方便。二、指纹特征匹配所述特征匹配的过程为首先利用脊线的曲率特征进行匹配参考点对的识 别；然后利用极坐标对图像进行校准；接着利用自适应弹性匹配方法进行细节 点配对和匹配分值计算；最后利用匹配的点对数和匹配分值，采用多级判别体 系进行指纹匹配。具体处理流程如图l右半部分所示。1.识别匹配参考点对如果用穷举法对输入指纹和模板指纹的每个点对的 相似性进行计算，显然是非常耗时的，本发明提出了一种基于脊线的曲率特征 来识别存在对应关系细节点对的新方法。首先，需要对已经获得的两个点集的 矢量特征按照曲率大小的顺序排序，如果输入指纹的第i个细节点和模板指纹的第j个特征点存在对应关系，那么与这两个特征点相连的脊线的曲率 和^也 应该是基本相同的，并且该曲率是不随图像旋转或平移所变化的。于是有，如| - |< ;，则可以利用公式0，W/zen^e 计算两个细节特征的相似度。如果S(i， j)>Ta，则可以认为对应点对已经找到。否则如果计算了所有曲率相同的点对之间的相似度，没有哪个点对的相似度大于给定的阈值Ta，则可以把具有最大相似度的点对作为后续图像校准的中心。如果点对的最大相似度小于给定的阈值Tb，则可以直接判断两个图像不匹配。相似度计算公式中T是预先定义的阈值，"ql是两个特征矢量的欧氏距 离。相似度S(i,j)指明了两个细节特征点的匹配可能性，值越大匹配的可能性越 高。可以看出0《S(i,j)《1，对于S(i,j)-l说明两个细节特征完全匹配，而对于S(i, j) = 0则说明两个细节点完全不相同。2.图像校准图像校准的目的在于消除两幅指纹图像的几何形变。利用平 移变换参数和旋转变换参数可以修正图像间的刚性形变。对于非线性形变，因 为其主要出现在图像中心区域，其表现形式为非线性的向外辐射。因此，采用 极坐标表示方法是一种更有效的解决方法。首先把输入指纹图像和模板指纹图 像原来笛卡儿坐标系下的细节点转换到极坐标下，极坐标的原点为上一步骤中 确定的参考对应点对；然后把新的特征矢量集合按照极角的升序进行排列。进行图像校准的步骤如下1 )把输入指纹图像和模板指纹图像原来笛卡儿坐标系下的细节点转换到极 坐标下，极坐标的原点为上一步确定的对应参考点对，转换公式如下所示。P。<formula>formula see original document page 9</formula>乂 乂 式中"，x，w'"是细节点在原来坐标系下的取值，"，^，^f是图像校准中心在原来坐标系的值，(pri，ei,ei)T是转换后细节点极坐标的值，其中pri代表极径， ei代表极角，6i代表细节点i的方向与校准中心细节点方向的差值。这样一来，就 可以把原来的特征矢量FC改成一种更为简洁的形式FO(pri， ei， & var, f, r)，其中 var，俘B r的含义和FC中定义一样。2)把新的特征矢量集合按照极角的升序进行排列，结果如下 尸。={[<， ， w, ，//>,].. [ ,《,v《/ V:]}K[，M,《vaif,/iV,Q]…[/^,々f,varf，/"V^3. 细节点配对与匹配计算由于存在非线性变形，这就要求匹配算法要有 一定的弹性。通常，在匹配过程中，会利用一个边界盒来进行弹性匹配，即在 模板指纹的细节点上放一个边界盒，然后检查对应的输入指纹位置，看是否有 细节点落在边界盒内。4. 指纹匹配计算完所有存在对应关系的细节点对的匹配分值后，则可以 得到匹配的细节点对总数Nm和匹配分值之和Mm，其中^"^I^("力。下面将介绍多级判别体系首先，如果^^12或《/她W，似)^'6，其中N和M分别是指纹模板和输 入图像细节点的数量，则可以直接判断两幅指纹图像来自同一个手指； 否则，如果7^〈6，则可以直接判断两幅指纹图像不匹配；如果U^气Sii并且M^^i，则也可以判断两幅指纹图像匹配成功； 如果"A^"并且U2，则也可以判断两幅指纹图像匹配成功；其他情况都可以判断两幅图像不匹配。其中T1和T2都是预先设定的阈值，可根据一定数量的指纹样本统计而得。 本发明的优点在于(1)由于构建的指纹特征矢量不仅包含细节点本身的 位置、方向、类型信息，还包含细节点周围的邻域关系信息和脊线曲率信息， 从而保证了它的平移、旋转不变性，也更能反映指纹纹线的全局信息；(2)基于这种新的特征矢量，在指纹匹配时一方面利用脊线曲率特征进行图像校准， 降低了校准的复杂度，另一方面利用指纹的邻域关系特征进行匹配分值的计算， 利用匹配的细节点对总数和匹配分值之和构建多级判别体系，增加了匹配的鲁棒性；(3)细节点的伪特征剪枝以及可靠性验证，保证了细节点提取的可信度。 本发明的目的、特征和优点将结合实施例并参照附图进一步说明。
图l为本发明的指纹识别流程框图。图2为指纹细节点初选使用的3X3模板。图3为常见伪特征结构示意图。图4为删除因图像模糊带来的伪特征流程图。图5为删除因脊线断裂造成的伪特征流程图。图6为删除短枝和桥形伪特征流程图。图7为分歧点可靠性检验流程图。图8为分歧点可靠性检验示意图。图9为端点可靠性检验流程图。图10为端点可靠性检测示意图。图ll为特征矢量示意图。
具体实施方式
现结合附图对本发明做进一步说明。本发明所述的指纹识别方法由指纹特征提取和特征匹配两个步骤组成，下 面分别加以叙述。一、指纹特征提取所述特征提取的过程为针对于预处理的图像进行细节特征初选；对初选后的细节特征进行剪枝；对剪枝后的细节点特征进行可靠性验证；最后根据细 节点的邻域关系和脊线曲率特征构建特征矢量；具体处理流程如图l左半部分所 示。现对每一处理流程分述如下(一)指纹细节点初选本实施例采用论文"指纹细节特征提取与剪枝" (王崇文，李见为等，光电工程，2002， No. 4， pp67-70)中提出的指纹细节点 初选方法。使用了一个3X3的小窗口在细化后的指纹点线图上确定细节特征点 的位置。假定脊线上的点用"1"表示，背景用"0"表示，如图2所示，M是待 检测点，XI ......X8为M的八邻域点，呈逆时针排列。如果M是端点，那么必须满足下面公式(1):如果M是分歧点，则满足下面的公式(2):<formula>formula see original document page 11</formula> (2)用3X3的模板逐点对细化后的指纹图像的脊线进行检测，就可初步选出所 有的细节特征，记录下这些细节特征的类型和位置坐标，以便于下一步的剪枝 处理。(二)细节点剪枝对于细化后的点线指纹图，主要有如图3所示的几种 伪特征结构。图3 (a)是因为脊线断裂产生了两个错误端点；而3 (b)禾B3 (c) 是因为脊线的错误连接和断开产生了伪分歧点和伪端点；3 (d)是在对非平滑 脊线进行细化后产生的短枝；3 (e)是因为脊线间的错误连接产生伪分歧点。 本实施例采用论文"指纹细节特征提取与剪枝"(王崇文，李见为等，光电工 程，2002， No. 4， pp67-70)中的细节点剪枝算法，针对每一类细节伪特征分别 进行剪枝。1消除图像的边缘效应在图像的边缘，由于指纹传感器本身的原因以及图像预处理过程的不完整 性， 一般有大量的伪特征。可采用以下办法消除它们对初选出的特征，计算 它们与边界的距离，当距离小于预先设定的阈值Td时(本实施例中，针对256X 256的指纹图像，Td取25)，则认为该特征点不可靠，从指纹特征中删除该点的 记录。2删除图像的模糊部分带来的伪特征点图像的模糊部分，在细化之后，会出现大量的断点和分叉，这些对于指纹 识别来说，都是些不可靠的信息。因此要去除这部分的伪特征点。去除这些伪 特征点可以利用两个标准，首先是在单位区域内的细节密度标准，接着又考虑 到在图像清晰可靠的区域，特征的垂直对比度应比较大，而平行对比度则比较 小。所以，去除步骤如下算法一 (如图4所示)(1) 计算各特征点的结构数Ne和Nb;(2) 如果(Ne>Te)，其中Te为预先设定的端点的个数(本实施例中取5)， 执行步骤(6);(3) 如果(Nb〉Tb),其中Tb为预先设定的分歧点的个数(本实施例中取3)， 执行步骤(6);(4) 如果(Ne+Nb>Ts)，其中Ts为预先设定的结构数(本实施例中取6)， 执行步骤(6);(5) 如果((Varo—Varp) /logB<Tc)，其中B为灰度图像块的均值，Tc为 预先设定的对比度阈值(Tc的值可以用统计方法获得。本实施例利用一定数量 的指纹样本图像，用人工方法精确标定包含细节特征的图像块，计算各图像块 的(Varo—Varp) /1ogB，求得均值M，以0.75M作为Tc的备选值)，执行步骤(6)， 否则执行步骤(7);(6) 删除该特征点，算法结束；(7) 保留该特征点。3、 删除因脊线断裂造成的伪特征点如图3(a)、 (b)、 (c)所示的由于指纹脊线断裂和误连会产生一些伪细节特征，对于这些伪特征，主要考虑细节点间的距离和方向关系，故采用以下的方 法对其进行处理算法二 (如图5所示)(1) 确定一特征点ml，计算其方向6 1;(2) 如果存在另一特征点m2，满足(d (ml, m2) <Tdl)，其中Tdl为预 设的距离阈值(本实施例中Tdl设置为单倍指纹宽度)，删除特征点ml和m2;(3) 如果存在另一特征点m2，满足(d(ml， m2)〈Td2且9 1—92"±冗)， 其中0 2为特征点m2的方向角，Td2为预设的距离阈值(本实施例中Td2设置为 1.5倍指纹宽度)，则删除特征点ml和m2。4、 删除短枝和桥形伪特征点如图3 (d)和(e)是由短枝和桥形造成的伪特征点。仔细观察这两种结构，短枝分歧点的两个分支一般在一条直线上，而第三条分支的方向一般与另两条 有明显差异；对于桥形分歧点， 一般是两条分支共线，而第三条与之近似垂直。 本发明在算法中用到了三个预设值9bsl、 9bs2、 9bs3， 9bsl表示两条分支垂 直(本实施例中设为75度)，e bs2和9 bs3都表示两分支共线(本实施例中分别 设置为15度和25度)，下面介绍处理的算法过程 算法三(如图6所示)(1) 对于每一个分歧点，搜索其各条分支，如果各条分支的长度均大于入，则执行以下步骤，否则执行步骤(2);(a) 用一元线性回归的方法估算各分支的方向；(b) 如果有两条分支共线即夹角小于9 bs2且第三条分支垂直于共线的 两分支，沿垂直的分支搜索，搜索长度为5入/6;如果找到另一个特征点——分 歧点或端点)，则找到短枝或桥形，删除这两个伪细节；(c) 如果有两条分支共线即夹角小于e bs2且第三条分支不垂直于共线 的两分支，沿第三条分支搜索，搜索长度为3入/2;如果找到另一个端点，则找 到短枝，删除这两个伪细节；(2) 如果只有两条分支长度大于A ，则执行以下步骤，否则执行步骤(3);(a) 用一元线性回归的方法估算这两条分支的方向；(b) 如果这两条分支共线即夹角小于0bs3，沿第三条分支搜索，如 果找到端点或分歧点，则为短枝或短桥形，删除这两个伪细节；(3) 如果最多只有一条分支长度大于入，则删除该分歧点和长度不满入分 支上的端点。删除了伪细节特征后，可以进一步根据指纹图像的拓扑性质来对存留的细 节特征进行可靠性检测。由指纹图像可以看出，分歧点侧面的脊线与形成分歧 点的脊线分支平行，基于这个特性，本实施例用论文"指纹细节特征提取与剪 枝"(王崇文，李见为等，光电工程，2002， No, 4， pp67-70)中提到的可靠性 验证方法对分歧点的可靠性进行检测算法四(如图7所示)(1) 对于每一个分歧点的每一条分支，用3入/2个点来进行直线拟合，估算 出各分支的方向角；(2) 沿各分支跟踪，如果在A长度内出现另一细节特征，则该分歧点不可 靠，算法结束；(3) 沿各分支跟踪至A/2处，在垂直于各分支的方向上，检测是否在3入/2 的范围内存在临近脊线(如图8(a)所示)；(4) 如果有一条以上脊线未找到，则该分歧点不可靠，算法结束；(5) 如果都找到了临近脊线，则标记该分歧点为一般可靠；(6) 定义一个矩形ABCD，其中AC的长度为3入/2， AB的长度为4入(如图 8(b)所示)；(7) 从临近脊线与矩形的左交点P0和P1出发，沿临近脊线搜索；(8) 如果发现一个端点，则该分歧点不可靠，算法结朿；(9) 如果存在右交点P2和P3，则该分歧点高度可靠。对于端点特征，由指纹图像可知，脊线端点两侧的脊线逐步靠近，以保证 指纹脊线间的距离保持基本恒定。基于此特性，我们用以下算法来验证端点的可靠性算法五(如图9所示)(1) 用入个点来进行直线拟合，估算端点所在脊线的方向角；(2) 如果端点所在脊线的长度小于A，则该端点不可靠，算法结束；(3) 从端点出发，沿脊线跟踪至入/2处，在垂直于脊线的方向上，检测是否在士3A/2的范围内存在临近脊线；(4) 如果找不到两根脊线，则该端点不可靠，算法结束；(5) 定义一个矩形ABCD，其中AB的长度为3入/2, AD的长度为3入(如图 IO所示)；(6) 从临近脊线与矩形的左交点P0和P1出发，沿临近脊线搜索；(7) 如果一条脊线与AC或BD相交且另一条临近脊线与DC相交，则标记该 端点一般可靠；(8) 如果在搜索过程中发现另一特征点，则该端点不可靠，算法结束；(9) 如果两条临近脊线都与DC相交，则检测(a) 如果两临近脊线的夹角小于预设阈值9w (本实施例设为10度) 也即两脊线近似平行，则该端点高度可靠；否则该端点一般可靠，算法结束；(b) 如果在P2和P3之间检测到另一脊线，则该端点不可靠，算法结束；(三)构建特征矢量在本发明中，为了减少指纹形变对特征匹配的影响，定义了一种新的细节特征矢量，它总共包含24个元素，表示方法如下 = {x, 乂 &>， var， /, cx, cy， r,sl, /1, d2, 。2， s2, /2, d3， o3, s3, /3, d4, o4,54, /4}各元素的含义如图6所示mO是当前待研究的细节特征点，即中心细节点， 它的方向为oO (该方向值可以在预处理阶段获得)。以mO为坐标原点，以oO为x 轴的正方向，构建坐标系(如图虚线所示)，可以看到ml， m2和m3是分别位于 第一、第二和第三象限的细节特征点，这些特征点应该是在各象限所有细节点中距离m0最近的。oi是各细节点mi的方向，x和y是中心细节点mO在原来坐标系 中位置坐标，"是包含中心细节点的局部脊线方向，即图中的o0。考虑到特征提 取过程中误差，可以用参数var和f取代了类型参数，var是在原始指纹图像中以 mO为中心的局部窗口内各象素的灰度方差，f是在原始指纹图像中包含mO的局部 窗口的频率，它代表局部窗口中脊线的频率。这两个参数都可以在指纹预处理 阶段获得。如果仔细观察指纹图像，就会发现大多数的脊线在有限长度内可以当作是 圆弧的一部分，基于这一点，本发明使用了3个参数来表示与细节特征点相连的 脊线，它们分别是cx， cy和r。把与细节特征点相连的脊线看作是一段圆弧， 利用圆弧上的三个点就可以得到所在圆的中心坐标和半径。考虑到弹性形变的 影响，利用最小二乘法去拟合脊线，最后用中心的平均坐标cx 、 cy和半径的平 均值r来表示与细节特征点相连的脊线。di， si， ai and ii (i-l，2,3，4)是中心细节特征点的邻域关系特征，它们描述了中 心细节点(m0)与周围细节点(mi)之间的关系。di分别是mi到mO的距离，ai 是mO方向角所在直线与mi方向角所在直线的夹角，si是mO与mi连线所在直线分 别和mO方向角所在直线的夹角。il， i2， i3和i4分别是mi特征点在所有细节点集 合中的索引。考虑到在某些情况下无法为某个象限确定合适的细节特征，例如， 在图6中，在第四象限找不到合适的细节点，那么就用一个特殊的极小值来标识 d4， s4， a4和i4，如一1。二、指纹特征匹配所述特征匹配的过程为首先利用脊线的曲率特征进行匹配参考点对的识 别；然后利用极坐标对图像进行校准；接着利用自适应弹性匹配方法进行细节 点配对和匹配分值计算；最后利用匹配的点对数和匹配分值，采用多级判别体 系进行指纹匹配。具体处理流程如图l右半部分所示。现对每一处理流程分述如 下1.识别匹配参考点对首先，需要对己经获得的两个点集的矢量特征按照 曲率大小的顺序排序，如果输入指纹的第i的细节点和模板指纹的第j个特征点 存在对应关系，那么与这两个特征点相连的脊线的曲率 和^也应该是基本 相同的，并且该曲率是不随图像旋转或平移所变化的。于是有，如果卜p—乂1<7;(在本实施例中，对于分辨率为500DPI的指纹图像，其z;取25)，则可以利用<formula>formula see original document page 16</formula>公式0，W/zmW^ 计算两个细节特征的相似度。如果S(i,j)>Ta (在本实施例中，Ta的取值为0.75)，则可以认为对应点对已经找到。否则 如果计算了所有曲率相同的点对之间的相似度，没有哪个点对的相似度大于给 定的阈值Ta，则可以把具有最大相似度的点对作为后续图像校准的中心。如果 点对的最大相似度小于给定的阈值Tb (在本实施例中，Tb的取值为0.25)，则可以直接判断两个图像不匹配。相似度计算公式中T是预先定义的阈值(T的值可以用统计方法获得。本实 施例利用一定数量的指纹样本图像，用人工方法精确标定具有完全对应关系的 指纹细节，计算各对应点对的特征矢量的欧式距离，求得均值M和方差S，以 M+2 S作为T的备选值)，^^l是两个特征矢量的欧氏距离。2.图像校准图像校准的目的在于消除两幅指纹图像的几何形变。利用平 移变换参数和旋转变换参数可以修正图像间的刚性形变。对于非线性形变，因 为其主要出现在图像中心区域，其表现形式为非线性的向外辐射。因此，采用 极坐标表示方法是一种更有效的解决方法。进行图像校准的步骤如下1)把输入指纹图像和模板指纹图像原来笛卡儿坐标系下的细节点转换到极 坐标下，极坐标的原点为上一步确定的对应参考点对，转换公式如下所示。<formula>formula see original document page 16</formula>式中"，^"')7是细节点在原来坐标系下的取值，"，夂，^)7是图像校准中心在原来坐标系的值，(pri,ei，ei)T是转换后细节点极坐标的值，其中pri代表极径， ei代表极角，ej代表细节点i的方向与校准中心细节点方向的差值。这样一来，就 可以把原来的特征矢量FC改成一种更为简洁的形式FC'气pri, eb & var, f， r)，其中 讶口 r的含义和FC中定义一样。2)把新的特征矢量集合按照极角的升序进行排列，结果如下<formula>formula see original document page 16</formula>3.细节点配对与匹配计算由于存在非线性变形，这就要求匹配算法要有 一定的弹性。通常，在匹配过程中，会利用一个边界盒来进行弹性匹配，即在 模板指纹的细节点上放一个边界盒，然后检查对应的输入指纹位置，看是否有 细节点落在边界盒内。本发明所采用的细节点配对算法利用了自适应弹性匹配 方法，即边界盒的范围是随着细节点距离匹配中心的距离的远近而不断调整的。 判断两个细节点存在对应关系的条件如下两个细节点的极径差值小于给定阈值Apr (本实施例取单倍纹线宽度入)； 两个细节点的极角差值小于给定阈值Ae (本实施例取15度)。假设Pi模板指纹上一个细节点，Qj是输入指纹上与Pi存在对应关系的细节点，则它们直接匹配分值S (i,j)可以用公式 计算。如果匹配分值S(i,j)大于给定的阈值Tm(在本实施例中，Ta的取值为0.70)，则可以认为这两个细节点匹配，在记录它们匹配分值的同时，把已经匹配的细节点总数Nm加l;否则，只记录两个细节点匹配的分值。4.指纹匹配计算完所有存在对应关系的细节点对的匹配分值后，则可以 得到匹配的细节点对总数Nm和匹配分值之和Mm，其中^=2>(^')。下面将介绍多级判别体系首先，如果《^2或A^/M一7V，M)^X6，其中N和M分别是指纹模板和输 入图像细节点的数量，则可以直接判断两幅指纹图像来自同一个手指；否则，如果^<6，则可以直接判断两幅指纹图像不匹配； 如果H^i^S11并且12Tl，则也可以判断两幅指纹图像匹配成功；如果^1^9并且^"2，则也可以判断两幅指纹图像匹配成功； 其他情况都可以判断两幅图像不匹配。其中Tl和T2都是预先设定的阈值，并且满足n〈r2。其中T1和T2都是预先设定的阈值，可根据一定数量的指纹样本统计而得， 在本实施例中分别为11和12.5 。
权利要求
1. 一种指纹识别方法，其特征在于由指纹特征提取和特征匹配两个步骤组成所述指纹特征提取的过程为针对于预处理的图像进行细节特征初选；对初选后的细节特征进行剪枝；对剪枝后的细节点特征进行可靠性验证；最后根据细节点的邻域关系和脊线曲率特征构建特征矢量；所述特征匹配的过程为首先利用脊线的曲率特征进行匹配参考点对的识别；然后利用极坐标对图像进行校准；接着利用自适应弹性匹配方法进行细节点配对和匹配分值计算；最后利用匹配的点对数和匹配分值，采用多级判别体系进行指纹匹配。
2、 根据权利要求1所述的指纹识别方法，其特征在于，所述细节特征剪枝 的伪特征包括因为脊线断裂产生的两个错误端点；因为脊线的错误连接和断 开产生的伪分歧点和伪端点；在对非平滑脊线进行细化后产生的短枝；因为脊 线间的错误连接产生的伪分歧点。
3、 根据权利要求1所述的指纹识别方法，其特征在于，所述构建特征矢量 的方法为定义一种新的细节特征矢量，它总共包含24个元素，表示方法如下 FC = {x, j, w, var， /， cx， c少,r，"1, /1, d2, "2,512, /2, d3, a3, s3, /3, d4, a4, s4, /4}其中x和y是参考细节点在原来坐标系中位置坐标，w是包含参考细节点 的局部脊线方向，参数var是在原始指纹图像中以参考细节点为中心的局部窗口 内各象素的灰度方差，f是在原始指纹图像中包含参考细节点的局部窗口的频 率，它代表局部窗口中脊线的频率；cx 、 cy和r表示与细节特征点相连的脊线 所拟合圆弧的中心位置坐标和半径；di， si, ai and ii (i-l,2,3,4)是参考细节特征点 的邻域关系特征，它们描述了参考细节点(m0)与周围细节点(mi)之间的关 系；di分别是mi到mO的距离，ai是m0方向角所在直线与mi方向角所在直线 的夹角，si是mO与mi连线所在直线分别和mO方向角所在直线的夹角；il， i2， i3和i4分别是mi特征点在所有细节点集合中的索引。
4、 根据权利要求1所述的指纹识别方法，其特征在于，所述匹配参考点对 的识别方法为首先，需要对已经获得的两个点集的矢量特征按照曲率大小的顺序排序，如 果输入指纹的第i个细节点和模板指纹的第j个特征点存在对应关系，那么如果s(/，；)=1 —夂1<7^，则可以利用公式r一 FC, — FC,G，W&nv^ 计算两个细节 特征的相似度，如果S(i，j)〉Ta，则可以认为对应点对已经找到；否则如果计算了所有曲率相同的点对之间的相似度，没有哪个点对的相似度大于给定的阈值Ta，则可以把具有最大相似度的点对作为后续图像校准的中心；如果点对的最大相似度小于给定的阔值Tb，则可以直接判断两个图像不匹配；相似度计算公式中T是预先定义的阈值，T的值可以用统计方法获得； IW—^。l是两个特征矢量的欧氏距离。
5、根据权利要求1所述的指纹识别方法，其特征在于，所述对图像进行校 准的方法为-1 )把输入指纹图像和模板指纹图像原来笛卡儿坐标系下的细节点转换到极 坐标下，极坐标的原点为已经确定的对应参考点对，转换公式如下所示>,、 、《乂tan—'(:、一 xr式中"，^"',是细节点在原来坐标系下的取值，"，》A"是图像校准中心在原来坐标系的值，(pri,ei,ei)T是转换后细节点极坐标的值，其中pri代表极径， ei代表极角，ei代表细节点i的方向与校准中心细节点方向的差值；这样一来，就 可以把原来的特征矢量FC改成一种更为简洁的形式FO(pri, ei, & var, f， r)，其中 var， fS] r的含义和FC中定义一样；2)把新的特征矢量集合按照极角的升序进行排列，结果如下 《={[ //>,]...[力:,《， ，。}2。 = {[P。2, e卩，《e ， var,e, 一 ]…[;^ e, e e,《e, varf, / e, r e ]}。
6、根据权利要求1所述的指纹识别方法，其特征在于，所述指纹匹配的方 法为计算完所有存在对应关系的细节点对的匹配分值后，则可以得到匹配的细 节点对总数Nm和匹配分值之和Mm，其中M" =2>("');首先，如果^^12或^^"W，^0^'6，其中N和M分别是指纹模板和输入图像细节点的数量，则可以直接判断两幅指纹图像来自同一个手指；否则，如果气<6,则可以直接判断两幅指纹图像不匹配； 如果10SA^《11并且M^n，则也可以判断两幅指纹图像匹配成功；如果"A^^9并且M^T2，则也可以判断两幅指纹图像匹配成功； 其他情况都可以判断两幅图像不匹配；其中T1和T2都是预先设定的阈值，可根据一定数量的指纹样本统计而得。
全文摘要
本发明涉及一种指纹识别方法，由指纹特征提取和特征匹配两个步骤组成特征提取包括针对于预处理的指纹图像提取细节点特征，对初选的细节点进行剪枝，对剪枝后的细节点特征进行可靠性验证，以及进行特征矢量编码；特征匹配首先利用脊线曲率特征进行图像校准，降低了校准的复杂度，然后利用指纹的邻域关系特征进行匹配分值的计算和匹配点对的统计，根据匹配的点对数和匹配分值，利用多级判别体系来判断两枚指纹的相似度。本方法具有识别率高、识别速度快，可靠性强等优点。
文档编号G06K9/00GK101276411SQ200810106358
公开日2008年10月1日 申请日期2008年5月12日 优先权日2008年5月12日
发明者丁刚毅, 唐明湘, 王崇文 申请人:北京理工大学