发明名称 自动指纹识别方法中的图像分割及指纹纹线距离提取技术
摘要 本发明公开了一种自动指纹识别方法中的图像分割及指纹纹线距离提取技术,针对现有指纹图像分割方法对指纹图像质量的适应性较低,对质量较差的指纹图像分割结果处理不准确的缺陷,引入了指纹图像灰度分布的概率密度模型,提出了基于马尔科夫链蒙特卡洛模拟对指纹图像进行分割处理的新方法。本方法即使对低质量的指纹图像也可以准确的实现对指纹图像的有效分割。本发明针对离散傅立叶谱分析方法以及现有统计方法求取的指纹纹线距离和真实的纹线距离偏差较大及对图像质量的适应性不高的问题,提出了通过二维信号采样定理将指纹图像的离散谱转换为连续谱,并在此基础上提取准确的连续纹线距离,提高整个指纹识别技术的准确率。
申请公布号 CN1737821A 申请公布日期 2006.02.22
申请号 CN200510090723.6 申请日期 2005.08.15
申请人 阜阳师范学院 发明人 詹小四;陈蕴;孙道德;陈超;王峰
分类号 G06K9/00(2006.01) 主分类号 G06K9/00(2006.01)
代理机构 安徽省阜阳市科颍专利事务所 代理人 徐宝泉
主权项 1、一种自动指纹识别方法中的指纹图像分割及指纹纹线距离提取方法,其特征是指纹图像分割方法为先构造边界曲线概率密度模型,再根据指纹图像所具有的内在特征,构造外背景区域概率密度函数和内指纹区域概率密度函数,然后形成边界曲线的Markov链,再采用蒙特卡洛模拟方法加速Markov链收敛到最优解,实现对指纹图像的快速分割处理,将有效指纹图像从背景区域中分割开来,其步骤是:①构造边界曲线概率密度模型:对一条闭合曲线Γ,Γ的外部环形区域在背景区的外背景概率密度函数记为Pout(Γ),同样,内指纹概率密度函数记为Pin(Γ),Γ的边界曲线概率密度函数记为PL(Γ),则有:PL(Γ)=Pin(Γ)Pout(Γ);②根据指纹图像所具有的内在特征,构造外背景区域概率密度函数和内指纹区域概率密度函数:A、脊线区域中像素的灰度值服从以脊线中心灰度为中心的正态分布,概率密度函数定义为:<math> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>i</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>|</mo> <mi>ridge</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msup> </mrow> </math> B、谷线区域和背景区域中像素的灰度值服从以谷线中心灰度为中心的正态分布,概率密度函数定义为:<math> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>i</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>|</mo> <mi>valley</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <mi>i</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>|</mo> <mi>back</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msup> </mrow> </math> C、对于准确地将有效指纹区域从背景中分割开的闭合曲线Γ,其外部区域一定在背景区,则定义外背景概率密度函数为:<math> <mrow> <mi>Pout</mi> <mrow> <mo>(</mo> <mi>&Gamma;</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Pi;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mfrac> <mn>1</mn> <mrow> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msup> </mrow> </math> 而其内部区域一定在有效指纹区,其上的像素既可能在脊线区域上,也可能在谷线区域上,而谷线的灰度分布跟背景灰度分布相同,脊线和谷线在指纹区域交替出现。因此,定义内指纹概率密度函数为:<math> <mrow> <mi>Pin</mi> <mrow> <mo>(</mo> <mi>&Gamma;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mo>|</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mi>k</mi> </mrow> <mi>N</mi> </mfrac> <mo>|</mo> <mo>)</mo> </mrow> <munderover> <mi>&Pi;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mfrac> <mn>1</mn> <mrow> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msup> </mrow> </math> ③构造指纹边界曲线的马尔科夫链:A、对当前状态下边界曲线上每个点xik做布朗运动,得到下一个状态下所有边界曲线上的点<math> <mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>k</mi> </msubsup> <mo>=</mo> <mi>B</mi> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>k</mi> </msubsup> <mo>)</mo> </mrow> <mo>;</mo> </mrow> </math> B、依顺序连接xi+1k,得到下一个状态下的边界曲线<math> <mrow> <msubsup> <mi>&Gamma;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>0</mn> </msubsup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>1</mn> </msubsup> <mo>,</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>k</mi> </msubsup> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>i</mi> </msub> </msubsup> <mo>}</mo> <mo>;</mo> </mrow> </math> C、对曲线进行整理,去掉一些重复的圈,从而得到Markov链上的下一个状态<math> <mrow> <msub> <mi>&Gamma;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mi>m</mi> <mrow> <mo>(</mo> <msubsup> <mi>&Gamma;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mo>;</mo> </mrow> </math> ④采用蒙特卡洛模拟方法加速收敛到最优解A、采用Metropolis-Hastings方法构造以下转移核:<math> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>&Gamma;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>|</mo> <msub> <mi>&Gamma;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>min</mi> <mo>{</mo> <mn>1</mn> <mo>,</mo> <mfrac> <mrow> <mi>PL</mi> <mrow> <mo>(</mo> <msub> <mi>&Gamma;</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>PL</mi> <mrow> <mo>(</mo> <msub> <mi>&Gamma;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>}</mo> </mrow> </math> 对马尔科夫链上的两个状态Γi,Γi+1,计算其转移核,再以概率P(Γi+1|Γi)决定是否转移;B、如果连续拒绝转移次数>50,则认为收敛到最优解,这时的边界曲线即为所求的最优分割曲线,否则,继续计算其转移概率,判断是否转移;指纹纹线距离提取方法是将连续谱分析技术引入到指纹纹线距离提取技术中,提出了一种连续谱分析的指纹纹线距离提取方法,求取准确的指纹纹线距离信息。
地址 236032安徽省阜阳市清河路741号师范学院计算机系詹小四