发明名称 一种求取重叠指纹方向场的方法和装置
摘要 本发明公开了一种求取重叠指纹方向场的方法和装置,方法步骤如下:标定出指纹有效区域,分块图像,将分块图像进行离散傅立叶变换,滤除低频分量并搜索频域幅值的极值点对,根据极值点对求得分块图像的初始方向场;将重叠有效区域的分块图像的初始方向场进行矩阵建模,计算相关性矩阵,设置初始标注概率矩阵,然后依次迭代进行随机标注、计算响应、更新标注概率矩阵、判断收敛,最终得到两个指纹的分离独立方向场;将独立有效区域和分离独立方向场进行组合并选择其中一种组合方式得到分离后的指纹方向场。装置包括标定模块、初始方向场估计模块、独立方向场分离模块和组合选择模块。本发明具有重叠指纹分离质量高、脊线检测准确性好、应用范围广的优点。
申请公布号 CN103500334B 申请公布日期 2016.08.31
申请号 CN201310424595.9 申请日期 2013.09.17
申请人 中国人民解放军国防科学技术大学 发明人 陈芳林;刘亚东;张熠;周宗潭;胡德文
分类号 G06K9/34(2006.01)I;G06K9/00(2006.01)I 主分类号 G06K9/34(2006.01)I
代理机构 湖南兆弘专利事务所(普通合伙) 43008 代理人 赵洪;谭武艺
主权项 一种求取重叠指纹方向场的方法,其特征在于实施步骤如下:1)在原始图像上标定出指纹有效区域,所述指纹有效区域包括两个指纹的独立有效区域以及两个指纹之间的重叠有效区域;2)将原始图像进行分块得到分块图像,将每一个分块图像进行离散傅立叶变换得到频域图,滤除频域图中的低频分量并搜索频域幅值的极值点对,在独立有效区域的分块图像中搜索获得1个极值点对,在重叠有效区域的分块图像中搜索获得2个极值点对,分别计算所述极值点对所对应二维正弦波的相位角,根据分块图像中极值点对所对应二维正弦波的相位角求得各个分块图像的初始方向场;3)针对所有包含2个极值点对的重叠有效区域的分块图像,将分块图像的初始方向场进行矩阵建模,对所述矩阵模型对应的对象集合计算相关性矩阵,设置初始标注概率矩阵,然后依次迭代进行在两个指纹对应的标注中随机选择一个标注、计算对象集合中各个对象在随机选择标注后的响应、更新标注概率矩阵、判断标注概率矩阵的概率密度是否收敛,如果尚未收敛则返回继续迭代,否则迭代完成,根据在所述两个指纹对应的标注中选择的标注将所述对象集合中各个对象进行分离,从而得到所述重叠有效区域对应所述两个指纹的两个分离独立方向场;4)将所述两个指纹的独立有效区域的初始方向场和两个分离独立方向场进行组合得到两种组合方式,从两种组合方式中选择一种组合方式得到两个指纹分离后的指纹方向场;所述步骤3)的详细步骤如下:3.1)针对所有包含2个极值点对的重叠有效区域的分块图像,将对应极值点对的初始方向场进行矩阵建模;3.2)将所述矩阵模型中的每一个元素视为一个对象,针对所述矩阵模型对应的对象集合根据式(2)计算相关性矩阵;<img file="FDA0000952758970000011.GIF" wi="1397" he="311" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>s</mi><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mi>&delta;</mi><mrow><mo>(</mo><mo>|</mo><mi>O</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>O</mi><mrow><mo>(</mo><msup><mi>i</mi><mo>&prime;</mo></msup><mo>,</mo><msup><mi>j</mi><mo>&prime;</mo></msup><mo>,</mo><msup><mi>k</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow><mrow><mi>&pi;</mi><mo>/</mo><mn>2</mn></mrow></mfrac></mrow>]]></math><img file="FDA0000952758970000012.GIF" wi="677" he="126" /></maths>式(2)中,R<sub>IJ</sub>为对象集合的相关性矩阵;D为预设的控制边界宽度参数;s表示两个对象之间取相同标注的支持度因子,1‑s表示两个对象之间取不同标注的支持度因子;δ(·)为归一化函数,用于将输入值归一化到区间(‑π/2,π/2]内,当输入值小于或等于π/2时,直接将输入值作为输出值输出,否则用π减去输入值后作为输出值输出;O(i,j,k)和O(i′,j′,k′)为对象集合中被计算取相同标注的支持度因子的两个对象,i,j,k为对象O(i,j,k)在矩阵模型中的坐标,i′,j′,k′为对象O(i′,j′,k′)在矩阵模型中的坐标;3.3)针对所述对象集合设置标注概率矩阵,所述标注概率矩阵的每一个元素为由对应所述两个指纹的概率组成的转置矩阵,且所述标注概率矩阵的每一个元素的两个概率之和为1,通过均匀分布的方式初始化设置标注概率矩阵中所有元素的概率初始值;3.4)针对当前的标注概率矩阵,为所述对象集合中的所有对象在所述两个指纹对应的标注中随机选择一个标注;3.5)根据式(3)分别计算所述对象集合中各个对象在随机选择标注后的响应;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>/</mo><mi>N</mi><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mi>j</mi></munder><msub><mi>R</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>q</mi><mo>,</mo><msub><mi>s</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952758970000021.GIF" wi="1398" he="111" /></maths>式(3)中,β<sub>iq</sub>为所述对象集合中序号为i的对象在随机选择标注后的响应,N为所述对象集合中的对象数量,j为所述对象集合中与序号为i的对象计算关联关系另一个对象的序号,R<sub>ij</sub>(q,s<sub>j</sub>)为所述对象集合中序号为i的对象和对象集合中序号为j的对象之间的相关性矩阵,q为序号为i的对象在步骤3.4)中选择的标注,s<sub>j</sub>为所述对象集合中序号为j的对象的在步骤3.4)中选择的标注;3.6)根据式(4)重新计算更新标注概率矩阵的每一个元素;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>p</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>p</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&alpha;&beta;</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>p</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>p</mi><mrow><mi>i</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>p</mi><mrow><mi>i</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&alpha;&beta;</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>r</mi><mo>&NotEqual;</mo><mi>q</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952758970000022.GIF" wi="1398" he="150" /></maths>式(4)中,p<sub>iq</sub>(t+1)为对象集合中序号为i的对象在本次迭代时对应第一个指纹的概率,p<sub>iq</sub>(t)为对象集合中序号为i的对象在上一次迭代时对应第一个指纹的概率,q为第一个指纹的标注;p<sub>ir</sub>(t+1)为对象集合中序号为i的对象在本次迭代时对应第二个指纹的概率,p<sub>ir</sub>(t)为对象集合中序号为i的对象在上一次迭代时对应第二个指纹的概率,r为第二个指纹的标注;α为预设的权值参数,β<sub>iq</sub>为所述对象集合中序号为i的对象在随机选择标注后的响应;3.7)计算第一个指纹在本次迭代时的概率和上一次迭代时的概率之间的概率差值,计算第二个指纹在本次迭代时的概率和上一次迭代时的概率之间的概率差值,如果所述两个概率差值均小于预设阈值,则判定所述标注概率矩阵的概率密度收敛,跳转执行步骤3.8),否则判定所述标注概率矩阵的概率密度尚未收敛,继续跳转执行步骤3.4)进行下一轮迭代处理;3.8)根据在所述两个指纹对应的标注中选择的标注将所述对象集合中各个对象进行分离,从而得到所述重叠有效区域对应所述两个指纹的两个分离独立方向场;所述步骤4)的详细步骤如下:4.1)将重叠有效区域对应所述两个指纹的分离独立方向场进行膨胀运算处理;4.2)将所述两个指纹的独立有效区域的初始方向场和膨胀运算处理后的两个分离独立方向场进行组合得到两种组合方式,并根据式(5)计算所述两种组合方式的相似性因子;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&lsqb;</mo><mfrac><mn>1</mn><msub><mi>N</mi><mn>1</mn></msub></mfrac><munder><mo>&Sigma;</mo><mi>x</mi></munder><munder><mo>&Sigma;</mo><mi>y</mi></munder><mi>&delta;</mi><mrow><mo>(</mo><mo>|</mo><msubsup><mi>O</mi><mrow><mi>o</mi><mo>,</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><msub><mi>O</mi><mrow><mi>n</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>|</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>B</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><mn>1</mn><msub><mi>N</mi><mn>2</mn></msub></mfrac><munder><mo>&Sigma;</mo><mi>x</mi></munder><munder><mo>&Sigma;</mo><mi>y</mi></munder><mi>&delta;</mi><mrow><mo>(</mo><mo>|</mo><msubsup><mi>O</mi><mrow><mi>o</mi><mo>,</mo><mn>2</mn></mrow><mo>&prime;</mo></msubsup><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><msub><mi>O</mi><mrow><mi>n</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>|</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&lsqb;</mo><mfrac><mn>1</mn><msub><mi>N</mi><mn>1</mn></msub></mfrac><munder><mo>&Sigma;</mo><mi>x</mi></munder><munder><mo>&Sigma;</mo><mi>y</mi></munder><mi>&delta;</mi><mrow><mo>(</mo><mo>|</mo><msubsup><mi>O</mi><mrow><mi>o</mi><mo>,</mo><mn>2</mn></mrow><mo>&prime;</mo></msubsup><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><msub><mi>O</mi><mrow><mi>n</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>|</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>B</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><mn>1</mn><msub><mi>N</mi><mn>2</mn></msub></mfrac><munder><mo>&Sigma;</mo><mi>x</mi></munder><munder><mo>&Sigma;</mo><mi>y</mi></munder><mi>&delta;</mi><mrow><mo>(</mo><mo>|</mo><msubsup><mi>O</mi><mrow><mi>o</mi><mo>,</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><msub><mi>O</mi><mrow><mi>n</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>|</mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952758970000031.GIF" wi="1398" he="567" /></maths>式(5)中,c<sub>1</sub>为第一种组合方式的相似性因子,c<sub>2</sub>为第二种组合方式的相似性因子,O′<sub>o,1</sub>(x,y)为膨胀运算处理后的第一个分离独立方向场,O′<sub>o,2</sub>(x,y)为膨胀运算处理后的第二个分离独立方向场,O<sub>n,1</sub>(x,y)为第一个独立有效区域,O<sub>n,2</sub>(x,y)为第二个独立有效区域,B<sub>1</sub>(x,y)为第一种组合方式中重叠有效区域与第一个独立有效区域之间形成的重合边界,B<sub>2</sub>(x,y)为第二种组合方式中重叠有效区域与第二个独立有效区域之间形成的重合边界,N<sub>1</sub>为重合边界B<sub>1</sub>(x,y)的所有像素总和,N<sub>2</sub>为重合边界B<sub>2</sub>(x,y)的所有像素总和;x,y为像素点在图像中的坐标;4.3)比较两种组合方式的相似性因子,选择相似性因子数值较小的一种组合方式得到两个被分离的指纹方向场,分别通过均值滤波器滤波后输出。
地址 410073 湖南省长沙市砚瓦池正街47号中国人民解放军国防科学技术大学三院