基于背景密度估计的印章图像提取方法
【专利摘要】一种基于背景密度估计的印章图像提取方法，首先对印章灰度图像用OTSU算法进行二值化，并对二值化后的图像进行去噪；其次，计算四种二值图积分图像，自适应地计算分形维数窗口尺寸，计算背景密度能量图，估计印章位置，依据真实印章大小提取印章图像，即完成算法目标；本发明方法利用密度估计的思想，创新性地对图像背景像素点密度进行估计，而不是对图像前景点密度进行估计，使印章作为低密度区域，非印章部分作为高密度区域，将待鉴定印章从简单背景中提取，作为印章配准的前处理步骤，大大提高了印章图像配准的准确率，使图像配准结果可以作为印章鉴定的依据。
【专利说明】基于背景密度估计的印章图像提取方法
【技术领域】
[0001] 本发明涉及印章图像提取【技术领域】，具体涉及一种基于背景密度估计的印章图像 提取方法。
【背景技术】
[0002] 如图1所示，印章是我国特有的历史文化产物，古代主要用作身份凭证和行驶职 权的工具。如今成为一种用作印于文件上表示鉴定或签署的文具。被制作出的每枚印章都 是独一无二的，所以在现代社会中，印章成为判定文书，证件，金融票据等真伪性的实用方 法，与签名具有同等的作用。
[0003]图像配准技术是图像处理和计算机视觉领域中的一个典型问题和技术难点，是将 同一对象的不同图像变换到相同坐标系的过程。它被广泛地应用在遥感图像、医学影像、三 维重构、机器人视觉等诸多领域中。其目的在于比较或融合针对同一对象在不同条件下获 取的图像，使两图中对应点一一对应起来，从而达到相似性度量或者差异检出等目的。
[0004] 分形维数是描述分形最主要的参量，简称分维。它描述了一个分形对空间填充程 度的统计量，分形维数没有统一的定义，其定义方法包括有计盒维数，豪斯多夫维数和分配 维数等。分形维数模型可以用来衡量物体表面粗糙度，也可以用来对纹理进行描述，在本发 明中将其作为一种图像密度估计方法。
[0005] 对于仿制印章真假的鉴别，可以先采用图像配准技术对待鉴定印章与真实印章图 像进行配准，然后再交由专业鉴别师进行鉴定。在实际情况中从文件，文档上截取的待鉴定 印章图像一般会带有一些规则背景如图2所示，这些背景对图像配准算法会造成很大的干 扰，使配准结果具有较大的偏差，使鉴别师不能利用配准结果进行鉴定。
【发明内容】
[0006] 为了解决上述现有技术存在的问题，本发明目的在于提供一种基于背景密度估计 的印章图像提取方法，本方法利用密度估计的思想，创新性地对图像背景像素点密度进行 估计，而不是对图像前景点密度进行估计，使印章作为低密度区域，非印章部分作为高密度 区域，将待鉴定印章从简单背景中提取，作为印章配准的前处理步骤，大大提高了印章图像 配准的准确率，使图像配准结果可以作为印章鉴定的依据。
[0007] 为达到以上目的，本发明采用以下技术方案：
[0008] -种基于背景密度估计的印章图像提取方法，作为印章图像配准的前处理方法， 其特征在于：通过对背景像素点的密度进行估算，将印章位置作为背景像素点的低密度区 域进行定位；对印章待配准图像密度的估算采用到了有关分形维数的计算方法，利用图像 计盒维数的近似值来衡量背景像素点的密度；对于密度能量图分形维数窗口半径大小的 选择，采用标准模板图像半长轴按一定比例放缩的大小；在计算计盒维数近似值时采用了 四种不同的积分图像，通过其递推关系，使密度能量图的计算能在〇(N)的时间复杂度内实 现；该方法主要步骤为：首先对印章灰度图像用0TSU算法进行二值化，并对二值化后的图 像进行去噪；其次，计算四种二值图积分图像，自适应地计算分形维数窗口尺寸，计算背景 密度能量图，估计印章位置，依据真实印章大小提取印章图像，即完成算法目标；
[0009]所述计算四种二值图积分图像的具体方法为：分别计算图像的四种积分图像；第 一种图像每个像素点的值DP_SINGLE[I] [J]定义为此像素点到原图左上角围成的矩形中 背景点的个数；第二种图像每个像素点的值DP_D0UBLE[I] [J]定义为此像素点到原图左上 角围城的矩阵中相邻排列的至少包含一个背景点的2*2大小的方块的个数；第三种图像每 个像素点的值DP_R0W2[I] [J]定义为此像素点到相同行第一个像素点间包括第一个像素 点相邻排列的至少包含一个背景点的1*2大小的方块的个数；第三种图像每个像素点的值 DP_C0L2[I] [J]定义为此像素点到相同列第一个像素点间包括第一个像素点相邻排列的至 少包含一个背景点的2*1大小的方块个数；利用动态规划递推的算法，四种积分图像都能 够在〇(N)的时间复杂度内求解；
[0010] 所述自适应地计算分形维数窗口尺寸的具体方法为：分形维数窗口尺寸是下文计 算分形维数估计图像密度的重要参数。其大小的选择主要与印章大小，印章纹路宽度，背景 纹理宽度有关；需要保证分形维数窗口大小大于印章纹路宽度和背景纹理宽度，但又远小 于印章大小。在一般情况下，印章尺寸都是远大于印章纹路宽度与背景纹理宽度的，本发明 采用简单的与印章尺寸相关的比例系数计算分形维数窗口尺寸，这里采用印章半长轴放缩 一定比例作为分形维数窗口大小；
[0011] 所述计算背景密度能量图的具体方法为：背景密度能量图每个像素点的取值为在 以此像素点为中心，具有分形维数窗口尺寸的方块上的背景点的计盒维数大小；所述计盒 维数大小采用方块内背景点的个数Nsingle与相邻排列的至少包含一个背景点的2*2小的方 块的个数Nd()Uble的比例的以2为底的对数来近似代替；对于方块内Nsingle与Nd。#的计算， 利用上述所述的四种积分图像能够直接在0(1)的时间复杂度内求解；最后，对图像作归一 化处理，作为密度能量图的结果；
[0012]所述估计印章位置的具体方法为：对背景密度能量图用0TSU二值化算法作二值 化，将高密度区域作为背景，低密度区域作为前景；计算所有前景点像素坐标的平均值，以 此作为对印章位置的估计。
[0013]本发明和现有技术相比，具有如下优点：
[0014] 1、本发明方法利用密度估计的思想，创新性地对图像背景像素点密度进行估计， 而不是对图像前景点密度进行估计，使印章作为低密度区域，非印章部分作为高密度区域， 将待鉴定印章从简单背景中提取，作为印章配准的前处理步骤，大大提高了印章图像配准 的准确率，创新性地对背景像素点密度进行估计，对于确定印章位置更加鲁棒。
[0015] 2、算法运算效率高。
[0016] 3、对待配准图像背景的形状，规格等依赖性小。
[0017]4、可以对不同的印章进行提取，对印章样式的依赖性小。
【专利附图】
【附图说明】
[0018] 图1印章图像，以及印章配准时的标准模板图像。
[0019] 图2实际场景下，带有局部背景的待配准印章图像。
[0020] 图3直接对带局部背景的印章图像配准，与使用本方法做前处理后再对带局部背 景的印章图像配准的结果对比。
[0021] 图4印章配准时，模板图像上和待配准图像上实际需要配准的区域。
[0022] 图5实际场景下，执行印章提取算法前的待配准印章图像。
[0023] 图6执行本方法后，提取的印章图像。
[0024] 图7对印章图像二值化前后的图像。
[0025] 图8对二值化图像去噪前后的图像。
[0026] 图9图像解释积分图DP_R0W2定义。
[0027] 图10图像解释积分图DP_C0L2定义。
[0028] 图11图像解释积分图DP_D0UBLE定义。
[0029] 图I2模板图像长轴。
[0030] 图13对待配准图像计算出的背景密度能量图，越亮的地方代表背景的密度越大， 越暗的地方代表背景的密度越小。印章处于背景密度小的地方。
[0031] 图14在待配准图像上估计的印章位置（红点），以及以模板图像半长轴为半径的 印章图像提取范围（红圈以内）。
【具体实施方式】
[0032] 以下结合附图及具体实施例对本发明作进一步的详细描述。
[0033] 本发明方法利用密度估计的思想，创新性地对图像背景像素点密度进行估计，而 不是对图像前景点密度进行估计，使印章作为低密度区域，非印章部分作为高密度区域，将 待鉴定印章从简单背景中提取，作为印章配准的前处理步骤，大大提高了印章图像配准的 准确率，使图像配准结果可以作为印章鉴定的依据。采用印章提取和不采用印章提取进行 配准的结果如图3所示。
[0034] 本发明方法作为印章配准的前处理算法，将待配准印章图像中的有效图像与印章 模板相匹配，可以有效提高印章图像配准的准确率，如图4所示。
[0035] 本发明一种基于背景密度估计的印章图像提取方法，首先对印章灰度图像用0TSU 算法进行二值化，并对二值化后的图像进行去噪；其次，计算四种二值图积分图像，自适应 地计算分形维数窗口尺寸，计算背景密度能量图，估计印章位置，依据真实印章大小提取印 章图像，即完成算法目标。如图5、图6所示。
[0036] 下面分步详细描述：
[0037] 1、对印章灰度图像用最大类间方差0TSU算法（0TSU)进行二值化，如图7所示， 又叫大津法。其基本思想是使用某一个阈值将灰度图像根据灰度大小，分成目标部分和背 景部分两类，在这两类的类内方差最小和类间方差最大的时候，得到的阈值是最优的二值 化阈值。
[0038] 枚举图像二值化阈值Threshold,使灰度值大于Threshold的点为前景点，灰度 值小于Threshold的点为背景点。设前景点所占比例为Ratio_ForeGround,背景点所占 比例为Rati〇_BackGr〇und，前景点平均灰度为U，背景点平均灰度为V，最大类间方差G= Ratio_ForeGround*Ratio_BackGround* (U-V)2,用使最大类间方差G值最大的Threshold 对图像进行二值化。
[0039] 2.对二值化后的图像去噪。使用形态学操作对图像去噪，如图8所示，利用图像的 开闭操作将图像上结构尺寸较小的噪音滤除。
[0040] 3.计算二值图的四种积分图像，通过递推算法，四张积分图的计算都能在0(N)的 时间复杂度内求解。假设二值图像上像素点取值为〇或1。第一种积分图像每个像素点的 值DP_SINGLE[I] [J]定义为此像素点到原图左上角围成的矩形中背景点的个数，其递推表 达式为
[0041]DP_SINGLE[I][J] =DP_SINGLE[I-1][J]+DP_SINGLE[I][J-1]-
[0042]DP_SINGLE[I-1][J-l]+Pixel_Value[I][J]；
[0043] DP_R0W2[I] [J]定义为此像素点到相同行第一个像素点间（包括第一个像素点） 相邻排列的至少包含一个背景点的1*2大小的方块的个数，如图9所示，其递推表达式为
[0044]DP_R0W2[i] [j] =DP_R0W2 [i] [j-2]+Pixel[i] [j]*Pixel[i] [j-l]；
[0045] DP_C0L2[I] [J]定义为此像素点到相同列第一个像素点间（包括第一个像素点） 相邻排列的至少包含一个背景点的2*1大小的方块个数，如图10所示，其递推表达式为
[0046]DP_C0L2[i][j] =DP_C0L2[i-2][j]+Pixel[i,j]*Pixel[i-1][j]；
[0047]DP_D0UBLE[I] [J]定义为此像素点到原图左上角围城的矩阵中相邻排列的至少包 含一个背景点的2*2大小的方块的个数，如图11所示，其递推表达式为
[0048]DP_D0UBLE[i] [j] =DP_D0UBLE[i-l] [j]+DP_D0UBLE[i] [j-l]-DP_D0UBLE[i-l] [j-l]+Pixel[i] [j]*Pixel[i] [j-l]*Pixel[i_l] [j]*Pixel[i_l] [j-l];
[0049] 4.自适应地计算分形维数窗口尺寸。分形维数窗口尺寸根据模板印章尺寸而定， 如图12所示。首先，计算模板图像二阶中心矩ii^ii2(lli%，设X。，y。为图像质心的X，y坐 标，其定义为
【权利要求】
1. 一种基于背景密度估计的印章图像提取方法，作为印章图像配准的前处理方法，其 特征在于：通过对背景像素点的密度进行估算，将印章位置作为背景像素点的低密度区域 进行定位；对印章待配准图像密度的估算采用到了有关分形维数的计算方法，利用图像计 盒维数的近似值来衡量背景像素点的密度；对于密度能量图分形维数窗口半径大小的选 择，采用标准模板图像半长轴按一定比例放缩的大小；在计算计盒维数近似值时采用了四 种不同的积分图像，通过其递推关系，使密度能量图的计算能在〇(N)的时间复杂度内实 现。
2. 根据权利要求1所述的印章图像提取方法，其特征在于：其主要步骤为：首先对印章 灰度图像用OTSU算法进行二值化，并对二值化后的图像进行去噪；其次，计算四种二值图 积分图像，自适应地计算分形维数窗口尺寸，计算背景密度能量图，估计印章位置，依据真 实印章大小提取印章图像，即完成算法目标； 所述计算四种二值图积分图像的具体方法为：分别计算图像的四种积分图像；第一种 图像每个像素点的值DP_SINGLE[I] [J]定义为此像素点到原图左上角围成的矩形中背景 点的个数；第二种图像每个像素点的值DP_DOUBLE[I] [J]定义为此像素点到原图左上角围 城的矩阵中相邻排列的至少包含一个背景点的2*2大小的方块的个数；第三种图像每个像 素点的值DP_R0W2[I] [J]定义为此像素点到相同行第一个像素点间包括第一个像素点相 邻排列的至少包含一个背景点的1*2大小的方块的个数；第三种图像每个像素点的值DP_ C0L2[I] [J]定义为此像素点到相同列第一个像素点间包括第一个像素点相邻排列的至少 包含一个背景点的2*1大小的方块个数；利用动态规划递推的算法，四种积分图像都能够 在〇(N)的时间复杂度内求解； 所述自适应地计算分形维数窗口尺寸的具体方法为：采用印章半长轴放缩预设比例作 为分形维数窗口大小； 所述计算背景密度能量图的具体方法为：背景密度能量图每个像素点的取值为在以此 像素点为中心，具有分形维数窗口尺寸的方块上的背景点的计盒维数大小；所述计盒维数 大小采用方块内背景点的个数Nsingle与相邻排列的至少包含一个背景点的2*2小的方块的 个数Ndwble的比例的以2为底的对数来近似代替；对于方块内Nsingle与Ndwble的计算，利用 上述所述的四种积分图像能够直接在0(1)的时间复杂度内求解；最后，对图像作归一化处 理，作为密度能量图的结果； 所述估计印章位置的具体方法为：对背景密度能量图用0TSU二值化算法作二值化，将 高密度区域作为背景，低密度区域作为前景；计算所有前景点像素坐标的平均值，以此作为 对印章位置的估计。
【文档编号】G06T7/00GK104408721SQ201410690319
【公开日】2015年3月11日 申请日期:2014年11月25日 优先权日:2014年11月25日
【发明者】宋永红, 龚晨, 张元林 申请人:西安交通大学