发明名称 一种可消除底单干扰的印鉴真伪识别方法
摘要 本发明公开了一种可消除底单干扰的印鉴真伪识别方法,该方法按照以下步骤具体实施:印鉴图像的预处理;预处理结果图像的修复;印鉴图像与预存标准印鉴图像的配准;膨胀处理;对消除底单干扰的图像进行滤波;将配准后的标准印鉴图像进行一次腐蚀处理;进一步消除残留干扰;对印鉴图像的真伪进行粗判别;对印鉴图像的真伪进行精细判别。本发明的真伪识别方法适用于实际盖印条件下的票据印鉴,印鉴图像的预处理环节主要是消除票据叠加在印鉴之上的底单表格、文字等干扰,采用了颜色提取,印鉴的签字叠加缺损修复,局部干扰滤波去除等步骤完成;印鉴图像的真伪鉴别,则通过定义左差及右差的方式,根据残差比的分布特性给出真伪鉴别的结果。
申请公布号 CN101894268A 申请公布日期 2010.11.24
申请号 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><mi>x</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>,</mo><mi>y</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,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><mi>x</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>,</mo><mi>y</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>则预处理后的印鉴图像[I(x,y)]<sub>m×n</sub>,为:<img file="FSA00000194908700013.GIF" wi="1583" he="140" />其中,P<sub>th</sub>为阈值,在20%-50%的范围内选择;步骤2,对上步预处理后的印鉴图像[I(x,y)]<sub>m×n</sub>进行修复;步骤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="FSA00000194908700021.GIF" wi="1595" he="141" />步骤5,对消除底单干扰的图像[I<sub>R</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="FSA00000194908700022.GIF" wi="1427" he="140" />式(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="FSA00000194908700023.GIF" wi="1650" he="145" />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="FSA00000194908700031.GIF" wi="220" he="55" />为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="FSA00000194908700032.GIF" wi="1729" he="140" /><img file="FSA00000194908700033.GIF" wi="1733" he="138" />记左差不为0的像素点个数为N<sub>L</sub>,右差不为0的像素点个数为N<sub>R</sub>;9.2)对左差图[e<sub>L</sub>(x,y)]<sub>m×n</sub>进行贴标签处理,统计所有连通域的大小,保留的面积大于阈值<img file="FSA00000194908700034.GIF" wi="75" he="60" />的连通域,这些连通域为疑似伪造残差,设这些连通域有K<sub>L</sub>个,记作<img file="FSA00000194908700035.GIF" wi="397" he="65" />其连通域的面积分别为<img file="FSA00000194908700036.GIF" wi="462" he="70" />9.3)对K<sub>L</sub>个疑似伪造残差区域<img file="FSA00000194908700037.GIF" wi="367" he="65" />分别进行三次膨胀,以获得其邻域范围,设为<img file="FSA00000194908700038.GIF" wi="388" he="65" />其邻域面积分别为<img file="FSA00000194908700039.GIF" wi="450" he="72" />9.4)分别统计右差图中相应的邻域<img file="FSA000001949087000310.GIF" wi="363" he="65" />范围内为1的像素点数,分别记作<img file="FSA000001949087000311.GIF" wi="404" he="65" />9.5)分别按照下面公式计算右残差率<img file="FSA000001949087000312.GIF" wi="439" he="68" />即:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mrow><msub><mi>R</mi><mi>k</mi></msub><mo>=</mo></mrow></msub><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><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>K</mi><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>9.6)如果<img file="FSA000001949087000314.GIF" wi="63" he="46" />大于阈值<img file="FSA000001949087000315.GIF" wi="254" he="47" />为20%-50%,则保留该左差连通域,否则,将其删除;9.7)统计所有保留下来的左差连通域的面积总和,记为<img file="FSA000001949087000316.GIF" wi="107" he="58" />计算左差的总残差率ρ<sub>L</sub>为:<img file="FSA00000194908700041.GIF" wi="1097" he="127" />9.8)对右差图[e<sub>R</sub>(x,y)]<sub>m×n</sub>进行贴标签处理,统计所有连通域的大小,保留的面积大于阈值<img file="FSA00000194908700042.GIF" wi="75" he="60" />的连通域,这些连通域为疑似伪造残差,设这些连通域有K<sub>R</sub>个,记作<img file="FSA00000194908700043.GIF" wi="405" he="63" />其连通域的面积分别为<img file="FSA00000194908700044.GIF" wi="462" he="68" />9.9)对K<sub>R</sub>个疑似伪造残差区域<img file="FSA00000194908700045.GIF" wi="375" he="63" />分别进行三次膨胀,以获得其邻域范围,设为<img file="FSA00000194908700046.GIF" wi="387" he="63" />其邻域面积分别为<img file="FSA00000194908700047.GIF" wi="458" he="63" />9.10)分别统计左差图中相应的邻域<img file="FSA00000194908700048.GIF" wi="374" he="58" />范围内为1的像素点数,分别记作<img file="FSA00000194908700049.GIF" wi="409" he="63" />9.11)分别按照下面公式计算左残差率<img file="FSA000001949087000410.GIF" wi="421" he="66" />即:<maths num="0005"><![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><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>K</mi><mi>R</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths>9.12)如果<img file="FSA000001949087000412.GIF" wi="62" he="46" />大于阈值<img file="FSA000001949087000413.GIF" wi="258" he="46" />为20%-60%,则保留该右差连通域,否则,将其删除;9.13)统计所有保留下来的右差连通域的面积总和,记为<img file="FSA000001949087000414.GIF" wi="111" he="58" />计算右差的总残差率ρ<sub>R</sub>:<img file="FSA000001949087000415.GIF" wi="997" he="124" />9.14)计算总残差率ρ:ρ=(ρ<sub>L</sub>+ρ<sub>R</sub>)/2,如果ρ大于阈值ρ<sub>Th</sub>,ρ<sub>Th</sub>为1%-5%,则判断为假印鉴,否则判断为真印鉴。
地址 710048 陕西省西安市金花南路5号