发明名称 一种可消除底单干扰的印鉴真伪识别方法
摘要 本发明公开了一种可消除底单干扰的印鉴真伪识别方法,该方法按照以下步骤具体实施:印鉴图像的预处理;预处理结果图像的修复;印鉴图像与预存标准印鉴图像的配准;膨胀处理;对消除底单干扰的图像进行滤波;将配准后的标准印鉴图像进行一次腐蚀处理;进一步消除残留干扰;对印鉴图像的真伪进行粗判别;对印鉴图像的真伪进行精细判别。本发明的真伪识别方法适用于实际盖印条件下的票据印鉴,印鉴图像的预处理环节主要是消除票据叠加在印鉴之上的底单表格、文字等干扰,采用了颜色提取,印鉴的签字叠加缺损修复,局部干扰滤波去除等步骤完成;印鉴图像的真伪鉴别,则通过定义左差及右差的方式,根据残差比的分布特性给出真伪鉴别的结果。
申请公布号 CN101894268B 申请公布日期 2012.05.23
申请号 CN201010229733.4 申请日期 2010.07.16
申请人 西安理工大学 发明人 朱虹;王翔;邓颖娜;王栋;邢楠
分类号 G06K9/62(2006.01)I;G06K9/36(2006.01)I;G06T7/00(2006.01)I;G06T5/00(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 罗笛
主权项 1.一种可消除底单干扰的印鉴真伪识别方法,其特征在于,该方法按照以下步骤具体实施:步骤1,印鉴图像的预处理:假设盖印的印泥颜色为红色,设输入的大小为m×n的印鉴票据图像的红、绿、蓝三颜色通道值分别为[R(x,y)]<sub>m×n</sub>,[G(x,y)]<sub>m×n</sub>,[B(x,y)]<sub>m×n</sub>,则其灰度值为[F(x,y)]<sub>m×n</sub>:F(x,y)=0.299·R(x,y)+0.587·G(x,y)+0.114·B(x,y),x=1,2,...,m,y=1,2,...,n  (1)对灰度化后的印章图像[F(x,y)]<sub>m×n</sub>采用最大类间类内方差比法进行二值化处理,消除纸张背景影响,得到二值图像[F<sub>b</sub>(x,y)]<sub>m×n</sub>,即:<maths num="0001"><![CDATA[<math><mrow><msub><mi>F</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>Th</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&lt;</mo><mi>Th</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>x=1,2,...,m,y=1,2,...,n            (2)其中,Th为阈值,采用类间类内方差比法确定,之后计算每个像素点的红色占比:<maths num="0002"><![CDATA[<math><mrow><msub><mi>P</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>x=1,2,...,m,y=1,2,...,n    (3)则预处理后的印鉴图像[I(x,y)]<sub>m×n</sub>为:<img file="FSB00000656209900013.GIF" wi="683" he="158" />x=1,2,...,m,y=1,2,...,n            (4)其中,P<sub>th</sub>为阈值,在20%-50%的范围内选择;步骤2,对上步预处理后的印鉴图像[I(x,y)]<sub>m×n</sub>进行修复,具体按照以下步骤进行修复:首先,提取签名信息,对[I(x,y)]<sub>m×n</sub>和[F<sub>b</sub>(x,y)]<sub>m×n</sub>进行异或运算得[S(x,y)]<sub>m×n</sub>,即:<img file="FSB00000656209900021.GIF" wi="1182" he="140" />x=1,2,...,m,y=1,2,...,n    (5)根据手写体签名一般是采用黑色或者是蓝色笔书写而成,所以,根据原图的灰度分布消除[S(x,y)]<sub>m×n</sub>中的噪声,得到签名信息[S<sub>b</sub>(x,y)]<sub>m×n</sub>即:<img file="FSB00000656209900022.GIF" wi="1042" he="139" />x=1,2,...,m,y=1,2,...,n    (6)其中,Th<sub>0</sub>为灰度图像中较小的阈值,以判断黑色或蓝色签名的灰度值,之后,用[S<sub>b</sub>(x,y)]<sub>m×n</sub>来修复因叠加导致的印鉴预处理图像[I(x,y)]<sub>m×n</sub>的缺损,采用以下的方法来进行修复:找到所有[I(x,y)]<sub>m×n</sub>为0,但[S<sub>b</sub>(x,y)]<sub>m×n</sub>为1的点,设为{(x<sub>1</sub>,y<sub>1</sub>),...,(x<sub>K</sub>,y<sub>K</sub>)},分别在以点(x<sub>i</sub>,y<sub>i</sub>),i=1,2,...,K为中心,覆盖一个N×N的模板,N为奇数,统计该模板中为1像素个数,记作N<sub>b</sub>,如果N<sub>b</sub>≥(N×N)/2,则该点置I(x<sub>i</sub>,y<sub>i</sub>),i=1,2,...,K为1,否则,不做处理;步骤3,对上步修复后的印鉴图像[I(x,y)]<sub>m×n</sub>与预存标准印鉴图像[I<sub>s</sub>(x,y)]<sub>m×n</sub>进行配准:对预存标准印鉴图像[I<sub>s</sub>(x,y)]<sub>m×n</sub>进行平移与旋转,找到与[I(x,y)]<sub>m×n</sub>相比残差最小的位置,即完成两者的配准,配准后的标准印鉴图像仍为[I<sub>s</sub>(x,y)]<sub>m×n</sub>;步骤4,将[I<sub>s</sub>(x,y)]<sub>m×n</sub>进行膨胀处理:设膨胀后的标准印鉴图像为[I<sub>sd</sub>(x,y)]<sub>m×n</sub>,用其消除印鉴图像[I(x,y)]<sub>m×n</sub>的底单干扰,特别是当底单的表格与文字均为红色时的干扰,消除底单干扰的图像[I<sub>R</sub>(x,y)]<sub>m×n</sub>计算公式为:<img file="FSB00000656209900023.GIF" wi="791" he="141" />x=1,2,...,m,y=1,2,...,n     (7)步骤5,对消除底单干扰的图像[I<sub>R</sub>(x,y)]<sub>m×n</sub>进行滤波,消除残留的局部干扰,消除残留的局部干扰的具体操作如下:5.1)设计长度为N<sub>m</sub>的十字形模板,N<sub>m</sub>为奇数,设模板的中心为(x<sub>0</sub>,y<sub>0</sub>),则模板为:Tem={(x<sub>0</sub>-N<sub>m</sub>/2,y<sub>0</sub>),...,(x<sub>0</sub>,y<sub>0</sub>),...(x<sub>0</sub>+N<sub>m</sub>/2,y<sub>0</sub>),(x<sub>0</sub>,y<sub>0</sub>-N<sub>m</sub>/2),...,(x<sub>0</sub>,y<sub>0</sub>-1),(x<sub>0</sub>,y<sub>0</sub>+1),...,(x<sub>0</sub>,y<sub>0</sub>+N<sub>m</sub>/2)}5.2)将模板Tem依次作用于[I<sub>R</sub>(x,y)]<sub>m×n</sub>上所有值为1的点,模板的中心点(x<sub>0</sub>,y<sub>0</sub>)依次为待处理像素的位置,计算[I<sub>R</sub>(x,y)]<sub>m×n</sub>上模板范围内的所有为1的像素点的个数,如果该个数超过(N<sub>m</sub>+1)/2,则该点修正为背景点,由此,消除局部的残留干扰,经过处理后的印鉴图像记为[I<sub>Rs</sub>(x,y)]<sub>m×n</sub>;步骤6,将由步骤3得到的配准后的标准印鉴图像[I<sub>s</sub>(x,y)]<sub>m×n</sub>进行腐蚀处理,结构元素采用步骤5中的十字形模板,得到的结果记作[I<sub>se</sub>(x,y)]<sub>m×n</sub>;步骤7,将[I<sub>se</sub>(x,y)]<sub>m×n</sub>与[I<sub>Rs</sub>(x,y)]<sub>m×n</sub>进行“或”操作,以在保持输入印鉴图像的基本特征的同时,进一步消除残留干扰,结果为[I<sub>Rf</sub>(x,y)]<sub>m×n</sub>,即:<img file="FSB00000656209900031.GIF" wi="1414" he="143" />式(8)中x=1,2,...,m,y=1,2,...,n;步骤8,对印鉴图像的真伪进行粗判别:8.1)将步骤7得到的印鉴图像[I<sub>Rf</sub>(x,y)]<sub>m×n</sub>与配准后的标准印鉴[I<sub>s</sub>(x,y)]<sub>m×n</sub>做差异检测,得到残差图为[e(x,y)]<sub>m×n</sub>,即:<img file="FSB00000656209900032.GIF" wi="759" he="157" />x=1,2,...,m,y=1,2,...,n    (9)8.2)对[e(x,y)]<sub>m×n</sub>进行贴标签处理,统计所有连通域的大小,保留的面积大于阈值S<sub>Th</sub>的连通域,统计保留下来的所有连通域中的像素个数,记作N<sub>k</sub>;统计标准印鉴图像[I<sub>s</sub>(x,y)]<sub>m×n</sub>的目标像素点数,记作N<sub>s</sub>,残差率ρ<sub>e</sub>计算公式为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>e</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><mi>k</mi></msub><msub><mi>N</mi><mi>s</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>8.3)当残差率ρ<sub>e</sub>大于阈值<img file="FSB00000656209900042.GIF" wi="208" he="46" />为5%-10%,则表明该印鉴为伪造精度不高的假印鉴,如果ρ<sub>e</sub>不大于阈值ρ<sub>Th</sub>,则进入精细判别环节;步骤9,对印鉴图像的真伪进行精细判别:9.1)将步骤7得到的印鉴图像[I<sub>Rf</sub>(x,y)]<sub>m×n</sub>与配准后的标准印鉴[I<sub>s</sub>(x,y)]<sub>m×n</sub>,分别进行左差检测和右差检测,左差图[e<sub>L</sub>(x,y)]<sub>m×n</sub>和右差图[e<sub>R</sub>(x,y)]<sub>m×n</sub>的计算公式如下:<img file="FSB00000656209900043.GIF" wi="933" he="136" />x=1,2,...,m,y=1,2,...,n    (11)<img file="FSB00000656209900044.GIF" wi="900" he="172" />x=1,2,...,m,y=1,2,...,n    (12)记左差不为0的像素点个数为N<sub>L</sub>,右差不为0的像素点个数为N<sub>R</sub>;9.2)对左差图[e<sub>L</sub>(x,y)]<sub>m×n</sub>进行贴标签处理,统计所有连通域的大小,保留的面积大于阈值<img file="FSB00000656209900045.GIF" wi="74" he="61" />的连通域,这些连通域为疑似伪造残差,设这些连通域有K<sub>L</sub>个,记作<img file="FSB00000656209900046.GIF" wi="394" he="65" />其连通域的面积分别为<img file="FSB00000656209900047.GIF" wi="458" he="69" />9.3)对K<sub>L</sub>个疑似伪造残差区域<img file="FSB00000656209900048.GIF" wi="365" he="64" />分别进行三次膨胀,以获得其邻域范围,设为<img file="FSB00000656209900049.GIF" wi="386" he="70" />其邻域面积分别为<img file="FSB000006562099000410.GIF" wi="446" he="68" />9.4)分别统计右差图中相应的邻域<img file="FSB000006562099000411.GIF" wi="360" he="64" />范围内为1的像素点数,分别记作<img file="FSB000006562099000412.GIF" wi="400" he="64" />9.5)分别按照下面公式计算右残差率<img file="FSB000006562099000413.GIF" wi="434" he="71" />即:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&rho;</mi><msub><mi>R</mi><mi>k</mi></msub></msub><mo>=</mo><mfrac><msub><mi>N</mi><msub><mi>R</mi><mi>k</mi></msub></msub><msub><mi>N</mi><msub><mi>&Omega;</mi><msub><mi>L</mi><mi>k</mi></msub></msub></msub></mfrac><mo>,</mo></mrow></math>]]></maths>k=1,2,...,K<sub>L</sub>                     (13)9.6)如果<img file="FSB000006562099000415.GIF" wi="62" he="46" />大于阈值<img file="FSB000006562099000416.GIF" wi="254" he="58" />为20%-50%,则保留该左差连通域,否则,将其删除;9.7)统计所有保留下来的左差连通域的面积总和,记为<img file="FSB00000656209900051.GIF" wi="108" he="58" />计算左差的总残差率ρ<sub>L</sub>为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>L</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><msub><mi>&Omega;</mi><mi>L</mi></msub></msub><msub><mi>N</mi><mi>L</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>9.8)对右差图[e<sub>R</sub>(x,y)]<sub>m×n</sub>进行贴标签处理,统计所有连通域的大小,保留的面积大于阈值<img file="FSB00000656209900053.GIF" wi="74" he="60" />的连通域,这些连通域为疑似伪造残差,设这些连通域有K<sub>R</sub>个,记作<img file="FSB00000656209900054.GIF" wi="401" he="62" />其连通域的面积分别为<img file="FSB00000656209900055.GIF" wi="457" he="69" />9.9)对K<sub>R</sub>个疑似伪造残差区域<img file="FSB00000656209900056.GIF" wi="373" he="63" />分别进行三次膨胀,以获得其邻域范围,设为<img file="FSB00000656209900057.GIF" wi="386" he="64" />其邻域面积分别为<img file="FSB00000656209900058.GIF" wi="454" he="62" />9.10)分别统计左差图中相应的邻域<img file="FSB00000656209900059.GIF" wi="372" he="57" />范围内为1的像素点数,分别记作<img file="FSB000006562099000510.GIF" wi="405" he="89" />9.11)分别按照下面公式计算左残差率<img file="FSB000006562099000511.GIF" wi="417" he="65" />即:<maths num="0006"><![CDATA[<math><mrow><msub><mi>&rho;</mi><msub><mi>L</mi><mi>k</mi></msub></msub><mo>=</mo><mfrac><msub><mi>N</mi><msub><mi>L</mi><mi>k</mi></msub></msub><msub><mi>N</mi><msub><mi>&Omega;</mi><msub><mi>R</mi><mi>k</mi></msub></msub></msub></mfrac><mo>,</mo></mrow></math>]]></maths>k=1,2,...,K<sub>R</sub>                        (15)9.12)如果<img file="FSB000006562099000513.GIF" wi="60" he="48" />大于阈值<img file="FSB000006562099000514.GIF" wi="255" he="49" />为20%-60%,则保留该右差连通域,否则,将其删除;9.13)统计所有保留下来的右差连通域的面积总和,记为<img file="FSB000006562099000515.GIF" wi="110" he="57" />计算右差的总残差率ρ<sub>R</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>R</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><msub><mi>&Omega;</mi><mi>R</mi></msub></msub><msub><mi>N</mi><mi>R</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow></math>]]></maths>9.14)计算总残差率ρ:ρ=(ρ<sub>L</sub>+ρ<sub>R</sub>)/2,如果ρ大于阈值ρ<sub>Th</sub>,ρ<sub>Th</sub>为1%-5%,则判断为假印鉴,否则判断为真印鉴。
地址 710048 陕西省西安市金花南路5号