一种基于mscr区域特征的图像复制检测方法
【专利摘要】本发明提供了一种基于MSCR区域特征的图像区域复制检测方法。首先对彩色图像提取MSCR特征，然后将这些特征区域归一化为圆形区域。接下来计算每个圆形特征区域的Zernike矩，作为该区域的特征向量。然后计算这些特征向量的欧式距离，找出候选的特征向量匹配对，及对应的特征区域对。通过这些特征区域对在图像当中的位置关系，估计区域复制过程当中的仿射变换矩阵。最后根据仿射矩阵来确定图像是否经过区域复制，并且定位复制区域的位置。本发明使用了一种新的彩色图像特征，并且优化了区域定位方法，具有很好的检测准确性和定位精确度。
【专利说明】
-种基于MSCR区域特征的图像复制检测方法
技术领域
[0001] 本发明设及图像数字取证领域，更具体地，设及一种基于MSCR区域特征的图像复 制检测方法。
【背景技术】
[0002] 随着数码相机和智能手机的普及，数字图像已经和日常生活密不可分。但由于图 像处理软件使用的低口槛和互联网的飞速传播，数字图像内容的真实性和图像的来源变得 无法确定。而在许多场合，例如司法、新闻出版等领域，我们又急需保证图像内容的真实性 和来源的可靠性，因此数字图像取证技术应运而生。
[0003] 图像复制检测是数字图像取证领域中的一个重要分支，它主要是为了检测一副图 像中是否存在区域复制行为，即图像中的一部分区域复制到本图像中的其它区域。图像区 域复制的目的是为了夸大或是隐藏图像中部分内容，复制的过程中往往会包含缩放、旋转、 加噪声等处理手段，使得复制的痕迹无法肉眼查别。因此，一个好的复制检测方法应该要能 考虑到运些干扰措施，并能准确定位复制区域。
[0004] 现有检测技术通常分为两种:基于图像块和基于图像特征。基于图像块的方法由 于计算时间过长，适应性不高等缺点逐渐不再被使用。目前主流的方法是基于图像特征的 方法。不同检测方法之间的差别主要是在图像特征的选取、特征匹配策略和复制区域定位 方面。
【发明内容】
[0005] 本发明提供一种基于MSCR区域特征的图像复制检测方法，能够有效对图像篡改中 区域复制行为进行检测和定位。
[0006] 为了达到上述技术效果，本发明的技术方案如下：
[0007] -种基于MSCR区域特征的图像复制检测方法，包括W下步骤：
[000引S1:提取MSCR图像特征:对彩色图像中的每个像素点，计算它与周围像素点之间的 颜色差，色差小于给定阔值加寸，运些点被聚为一类。重复运个过程，直到所有的点都完成计 算。每个类中的所有像素点构成一个MSCR区域；
[0009] S2:描述特征区域:将每个MSCR不规则区域采用构造放射不变的方法表示成楠圆， 然后将楠圆区域归一化为圆形区域。对每个圆形区域，计算其Zernike矩，生成一个12维特 征向量；
[0010] S3:匹配特征:对于每个特征区域，计算其与其它所有特征区域之间特征向量的欧 式距离。比较运些距离中的最小值do与次最小值山之间的比值do/di，如果比值小于0.5，则 认为距离为do的两个特征匹配；
[0011] S4:聚类并过滤特征区域:如果没有足够多匹配的特征对，则检测结束。否则，对于 所有匹配的一对特征，记录所有特征区域的中屯、点。对于所有的中屯、点，采用k-means聚类 算法，设定k = 2,将运些点聚为两类。删除掉无法被聚类的点。如果一类中少于3个点，则删 除所有点，检测结束；
[0012] S5:估计仿射矩阵:任意Ξ个不共线的特征点对，可计算得到一个仿射矩阵Τι。对 于每个矩阵Τι，计算所有特征点对之间的误差，如果误差值小于阔值β，则运个矩阵Τι获得一 票。最终得票数最多的矩阵即为所得仿射矩阵；
[0013] S6:定位复制区域:对于原始图像，使用所得的仿射矩阵进行坐标变换得到变换后 的图像。计算原始图像与变换后的图像之间对应位置的相似度，生成一副表示相似度的图 像。相似度的值在[0，1]之间，值越大表示相似度越高。如果相似度大于0.5,则认为该位置 为复制区域。将生成的相似度图像经过简单的滤波处理，生成最终检测结果图。
[0014] 本发明中，首次使用MSCR特征作为图像特征进行复制检测。此特征直接从彩色图 像中提取，实际应用中待检测的图像基本上都是彩色图像。此特征在数量上与主流的特征 点一致，且特征提取速度上也很快。
[0015] 进一步地，所述步骤S2中计算特征区域的ZernAe矩的过程如下：
[0016] 对于每个圆形的特征区域，其ZernAe矩可用如下公式计算
[0022] 公式（l)I(p，Θ)为圆形的特征区域的极坐标表示，P和Θ分别表示坐标的级半径和 角度，j为复数单位。公式(3)中S为求和变量。
[0023] 公式（1)、（2)、（3)中，η为ZernAe矩的阶数，系数m是满足(n-|m|)/2为非负整数的 任意整数。令公式（1)、（2)、（3)中n = 5,则每个特征区域可计算出的Zernike矩为12维的特 征向量。
[0024] 进一步地，所述步骤S5中得到仿射矩阵T的过程如下：
[0025] 一个仿射矩阵可W表示为一个3X3的矩阵
[0026]
州
[0027] 公式(4)中a~f都为待定系数。
[0028] 对于Ξ个不共线的点对，可通过下式计算出一个仿射矩阵
[0029]
巧）
[0030] 其中P。与Pd为一对点。
[0031 ]进一步地，所述步骤S6中计算相似度的过程如下：
[0032]对与原始图像中每个像素点，使用公式(6)进行坐标变换，生成一副变换后的图 像。式中τ为之前所得的仿射矩阵，（x〇，y〇，l)为原始坐标，（Xd，yd，l)为变换后的坐标。
[0033] (x〇,y〇,l)T=(xd,yd,l) (6)
[0034] 取原始图像和变换图像对应位置的5X5区域的像素值，生成一对25维的向量，计 算运对向量的相关系数，通过相关系数值来表示运个5X5区域的相似度。遍历原始图像和 变换图像所有对应位置的5X5区域，即可得到所有点的相似度。
[0035] 与现有技术相比，本发明技术方案的有益效果是：
[0036] 本发明利用彩色图像的特征，与传统的灰度图的特征相比具有更强的区分性。区 域特征在复制检测中即既具有图像块的优点一一能够覆盖更多的图像区域，又有图像特征 的优点一一计算速度快。且本发明采用的基于区域特征的描述子相比基于点特征的描述子 具有更强的表达能力和区分度。在对复制区域的定位上，采用了求取仿射矩阵的方法，能够 在像素级别定位复制区域，更加准确。
【附图说明】
[0037] 图1为本发明方法的流程图；
[0038] 图2为本发明方法的复制检测的实验效果图。
【具体实施方式】
[0039] 附图仅用于示例性说明，不能理解为对本专利的限制；
[0040] 为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品 的尺寸；
[0041] 对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可W理解 的。
[0042] 下面结合附图和实施例对本发明的技术方案做进一步的说明。
[0043] 实施例1
[0044] 如图1所示，一种基于MSCR区域特征的图像复制检测方法，包括W下步骤：
[0045] S1:提取MSCR图像特征:对彩色图像中的每个像素点，计算它与周围像素点之间的 颜色差，色差小于给定阔值加寸，运些点被聚为一类。重复运个过程，直到所有的点都完成计 算。每个类中的所有像素点构成一个MSCR区域；
[0046] S2:描述特征区域:将每个MSCR不规则区域采用构造放射不变的方法表示成楠圆， 然后将楠圆区域归一化为圆形区域。对每个圆形区域，计算其Zernike矩，生成一个12维特 征向量，本实施例中如图2(a)所示，图中楠圆区域即是提取出的MSCR区域；
[0047] S3:匹配特征:对于每个特征区域，计算其与其它所有特征区域之间特征向量的欧 式距离。比较运些距离中的最小值do与次最小值山之间的比值do/di，如果比值小于0.5，则 认为距离为do的两个特征匹配，本实施例中如图2(b)所示，所有匹配成功的区域对由直线 所连接；
[004引S4:聚类并过滤特征区域:如果没有足够多匹配的特征对(至少5对），则检测结束。 否则，对于所有匹配的一对特征，记录所有特征区域的中屯、点。对于所有的中屯、点，采用k- means聚类算法，设定k = 2,将运些点聚为两类。删除掉无法被聚类的点。如果一类中少于3 个点，则删除所有点，检测结束；
[0049] S5:估计仿射矩阵:任意Ξ个不共线的特征点对，可计算得到一个仿射矩阵Τι。对 于每个矩阵Τι，计算所有特征点对之间的误差，如果误差值小于阔值β，则运个矩阵Τι获得一 票。最终得票数最多的矩阵即为所得仿射矩阵；
[0050] S6:定位复制区域:对于原始图像，使用所得的仿射矩阵进行坐标变换得到变换后 的图像。计算原始图像与变换后的图像之间对应位置的相似度，生成一副表示相似度的图 像。相似度的值在[0，1]之间，值越大表示相似度越高。如果相似度大于0.5,则认为该位置 为复制区域。将生成的相似度图像经过简单的滤波处理，生成最终检测结果图。
[0051] 如图2(c)所示，该基于MSCR区域特征的图像复制检测方法的实验效果。对比实际 的复制区域图2(d)，从图2(c)中可W明显地看到非常精确的定位了复制区域。该方法利用 图像MSCR区域特征，对比传统的点特征具有覆盖面大，计算效率快的特点。且在定位方面采 用的仿射变换矩阵方法准确率高，在实际的检测过程中可W看出检测结果较为理想。
[0052] 相同或相似的标号对应相同或相似的部件；
[0053] 附图中描述位置关系的用于仅用于示例性说明，不能理解为对本专利的限制；
[0054] 显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对 本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可 W做出其它不同形式的变化或变动。运里无需也无法对所有的实施方式予W穷举。凡在本 发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求 的保护范围之内。
【主权项】
1. 一种基于MSCR图像区域特征的区域复制检测方法，其特征在于，包括以下步骤： S1:提取MSCR图像特征:对彩色图像中的每个像素点，计算它与周围像素点之间的颜色 差，色差小于给定阈值d时，这些点被聚为一类;重复这个过程，直到所有的点都完成计算； 每个类中的所有像素点构成一个MSCR区域； S2:描述特征区域:将每个MSCR不规则区域采用构造放射不变的方法表示成椭圆，然后 将椭圆区域归一化为圆形区域;对每个圆形区域，计算其Zernike矩，生成一个12维特征向 量； S3:匹配特征:对于每个特征区域，计算其与其它所有特征区域之间特征向量的欧式距 离；比较这些距离中的最小值do与次最小值di之间的比值do/di，如果比值小于0.5，则认为 距离为do的两个特征匹配； S4:聚类并过滤特征区域:如果没有达到设定的匹配的特征对，则检测结束;否则，对于 所有匹配的特征，记录所有特征区域的中心点;对于所有的中心点，采用k-means聚类算法， 设定k = 2,将这些点聚为两类;删除掉无法被聚类的点；如果一类中少于3个点，则删除所有 点，检测结束； S5:估计仿射矩阵:任意三个不共线的特征点对，可计算得到一个仿射矩阵;对于每个 矩阵h，计算所有特征点对之间的误差，如果误差值小于阈值β，则这个矩阵获得一票;最 终得票数最多的矩阵即为所得仿射矩阵； S6:定位复制区域:对于原始图像，使用所得的仿射矩阵进行坐标变换得到变换后的图 像;计算原始图像与变换后的图像之间对应位置的相似度，生成一副表示相似度的图像;相 似度的值在[〇，1 ]之间，值越大表示相似度越高;如果相似度大于〇. 5，则认为该位置为复制 区域;将生成的相似度图像经过简单的滤波处理，生成最终检测结果图。2. 根据权利要求1所述的基于MSCR区域特征的图像复制检测方法，其特征在于，所述步 骤S2中计算特征区域的Zern i ke矩的过程如下： 对于每个圆形的特征区域，其Zern i ke矩可用如下公式计算公式（1)Ι(Ρ，Θ)为圆形的特征区域的极坐标表示，p和Θ分别表示坐标的级半径和角度， j为复数单位。公式(3)中s为求和变量。 公式（1)、（2)、（3)中，η为Zernike矩的阶数，系数m是满足(n-|m|)/2为非负整数的任意 整数。令公式（1)、（2)、（3)中n = 5,则每个特征区域可计算出的Zernike矩为12维的特征向 量。3. 根据权利要求1所述的基于MSCR区域特征的图像复制检测方法，其特征在于，所述步 骤S5中得到仿射矩阵T的过程如下： 一个仿射矩阵表示为一个3 X 3的矩阵：公式(4)中a~f都为待定系数。 对于三个不共线的点对，可通过下式计算出一个仿射矩阵：其中P。与Pd为一对已知点的坐标。4.根据权利要求1所述的基于MSCR区域特征的图像复制检测方法，其特征在于，所述步 骤S6中计算相似度的过程如下： 对与原始图像中每个像素点，使用公式(6)进行坐标变换，生成一副变换后的图像;式 中T为之前所得的仿射矩阵，（X。，y。，1)为原始坐标，（xd，yd，1)为变换后的坐标； (x〇,y〇,l)T=(xd,yd,l) (6) 取原始图像和变换图像对应位置的5X5区域的像素值，生成一对25维的向量，计算这 对向量的相关系数，通过相关系数值来表示这个5X5区域的相似度;遍历原始图像和变换 图像所有对应位置的5 X 5区域，即可得到所有点的相似度。
【文档编号】G06T7/40GK105825504SQ201610142022
【公开日】2016年8月3日
【申请日】2016年3月11日
【发明人】卢伟, 杨帆
【申请人】中山大学