一种基于多域耦合的图像人工模糊检测方法与流程

本发明涉及数字图像信息技术领域，具体涉及一种基于多域耦合的图像人工模糊检测方法。

背景技术：

随着网络和智能设备的普及，多媒体信息的传输变得无处不在。数字图像作为传输信息其中重要的一种载体，更是成为重中之重。但如今，数字图像内容的篡改变得越来越容易，人们不在需要专业的技术都可以几乎不留痕迹的对图像进行编辑修改，这为数字图像的真实性和安全性带来了挑战。因此对数字图像取证的研究对维护网络信息的安全具有重要的意义。

图像的人工模糊篡改检测作为数字图像取证技术的一个重要分支，其目的是为了检测数字图像中是否存在人工模糊的痕迹并定位出篡改的区域。在图像中人工模糊的目的通常是为了掩饰原始图像中一些敏感的信息内容，或者作为后处理的润饰功能掩饰图像中拼接等篡改的痕迹。通过这种精心巧妙的模糊设计，在图像中几乎很难通过人眼观察看到图像中篡改的痕迹。即使在图像中仍留有一些可以辨别的痕迹，也需要有力的方法去证明篡改的存在。因此有效的数字图像人工模糊检测可以在信息安全保障体系中可以发挥非常重要的作用。图像人工模糊检测技术被广泛的应用在司法取证、新闻纪实、媒体娱乐等领域中。

现有的图像人工模糊检测技术主要分为两种：一种是基于图像模糊的边缘检测，人工模糊区域和原始区域的交界边缘作为图像中特殊的存在不同于图像中正常的纹理边缘，通过滤波等方式可以强化提取到这个边缘，从而实现人工模糊的定位；另一种是基于图像像素的统计特征，对图像的每一个像素点进行特征的求取，并利用机器学习的方法对统计特征进行分类识别，这类方法之间的区别主要在于图像变换域和统计特征选取。

技术实现要素：

本发明的目的在于提供一种基于多域耦合的图像人工模糊检测方法，能够有效地检测定位出图像中是否存在人工模糊的操作，具有检测效率快、定位准确率高、不需要人为判断的优点。

为了达到上述技术目的，本发明的技术方案如下：

一种基于多域耦合的图像人工模糊检测方法，包括以下步骤：

s1：灰度转化：对于输入的原始测试图像首先转换成单通道的灰度图像i；

s2：二次模糊：对步骤s1中的灰度图像i使用高斯模糊进行全局模糊操作得到二次模糊后的图像ib，用于之后在各个域上i和ib相似度的比较，用公式表示为ib＝i*g

其中，g是高斯模糊核，*是卷积操作；

s3：提取dct域内特征：对于灰度图像i上每一个像素点p，找到它在二次模糊后的图像ib上的对应点p′，分别截取它们n*n邻域的图像块b和bb，并进行dct操作，得到dct域内特征b＝dct(b)，bb＝dct(bb)；

s4：dct域内相似度比较：利用相对差异计算步骤s3得到的b和bb的相似度，得到图像在dct域内的人工模糊度估计，有

其中，b(i,j)和bb(i,j)分别表示在坐标(i,j)上对应的dct系数值；tf是在dct上得到计算相似度的距离，用来表示dct域内的人工模糊程度；

s5：对于步骤s4中比值当分母接近或者为0时会出现异常，将比值利用sigmoid映射到[0,0.5]的范围内，具体的有，

其中，r表示比值和

s6：获得候选的可疑人工模糊区域：利用步骤s4求得图像中每个像素点对应的人工模糊度估计值tf，即可得到图像dct域内的人工模糊度估计映射图；

而后使用阈值t1对该映射图进行分割，利用图像形态学的开闭操作进行去燥和填充孔洞，即可在二值图中提取出候选的可疑人工模糊区域其中rb是候选的可疑人工模糊区域的集合；

s7：在步骤s6中得到的每一个可疑候选区域在二次模糊后的图像ib中找到其对应的区域其中rg是二次模糊后的图像ib中对应的集合，在空域内利用空域纹理的描述子包括异常灰度级统计、平滑度以及信息熵作为特征f计算和的相似度s，有

其中fb和fg分别是和计算特征f所得到的结果，然后以t2为阈值进行最终的判断得到人工模糊区域的定位结果，有

其中是求出的最终人工模糊区域的集合，表示对求出的s值。

进一步，所述步骤s7中特征f计算的处理过程如下：

s701：计算候选区域的图像梯度：对于每一个候选区域上的像素点计算它的梯度幅值，有

其中[.]表示取整操作；

s702：异常灰度级统计：利用灰度共生矩阵对提取特征，有

其中m为利用灰度共生矩阵提取到的矩阵，δ(·)为脉冲函数，用来对像素对进行计数，写为

式子中分别当di＝0，dj＝1和di＝1，dj＝0时可以得到m在水平和垂直方向上相邻梯度幅值的统计，分别用mh和mv来表示；

将统计中较大幅值的统计对视为异常的点对，异常元素表示为则对异常点对的统计有,

s703：平滑度估计：首先求得梯度值的方差σ，

其中n是中点的个数，是的均值，l是梯度幅值的变化范围；根据σ求得候选区域的平滑度，用f2表示有

s704：纹理信息熵估计：候选图像区域梯度的纹理信息熵可以表示为

其中是所出现的概率，通过计算直方图的频数求得，

其中指的是所对应直方图的统计数，因为f1、f2以及f3取值都正相关与s，采用如下式子将三种描述子结合起来，

得到候选区域在图像空域内的特征值。

与现有技术相比，本发明技术方案的有益效果是：

本发明联合了dct域和图像空域的提取的特征，能够更好的对图像中的细节信息进行表征，对人工模糊区域更加的敏感。首先是在dct内采用相对差异的方式和二次模糊后的图像进行相似度的比较，能够抓住图像高频分量在人工模糊和原始图像中的不同；在空域中利用三种纹理的描述子对图像的候选人工模糊区域进行筛选甄别，有效利用了人工模糊区域在空域上的一些内容的先验知识，进一步提高了定位的准确率。相比于传统的基于特征的检测算法，速度更快，实用性更强，能够准确的给出人工模糊区域的轮廓。

本发明首次联合dct域和空域对图像的人工模糊进行检测定位，能够充分利用两个域的优势，在和二次模糊的相似度比较中得到准确的检测定位结果，避免了人工按照经验来判断的操作，具有很好的检测效率和定位的准确率。

附图说明

图1为本发明方法流程图

图2为本发明实施例中待处理的人工模糊篡改后的图像

图3为本发明实施例中待处理篡改图像的掩膜

图4为本发明实施例中得到的dct域内的人工模糊度映射图

图5为本发明实施例定位结果图

具体实施方式

下面结合附图对本发明作进一步的说明。

为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；

对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。

实施例1

如图1所示，一种基于多域耦合的图像人工模糊检测方法，包括以下步骤：

s1：灰度转化：对于输入的原始测试图像首先转换成单通道的灰度图像i，转化公式为，

i＝0.299*r+0.587*g+0.114*b(1)

其中，r，g，b分别为图像在彩色三通道上的像素值。

s2：二次模糊操作：对步骤s1中的灰度图像i使用高斯模糊进行全局模糊操作得到二次模糊后的图像ib用于之后在各个域上i和ib相似度的比较，用公式表示为

ib＝i*g(2)

其中*是卷积操作，g是高斯模糊核，本实施例中g的参数模糊核大小为n＝25，标准差σ＝1。

s3：提取dct域内特征：对于i上像素点p，它在dct域上表现出的人工模糊程度由它所在的图像块的人工模糊程度所决定。截取出它n*n邻域的图像块b和该图像块在二次模糊图像ib上所对应的图像块bb，分别进行dct操作，b＝dct(b)，bb＝dct(bb)。本实施例中参数n＝25。

s4：dct域内相似度比较：为了突出表现高频分量的作用，利用相对差异计算s3得到的b和bb的相似度，得到像素点p在dct域内的人工模糊度估计值，有

其中，b(i,j)和bb(i,j)分别表示在坐标(i,j)上对应的dct系数值，tf是在dct上得到计算相似度的距离，可以用来表示像素点p在dct域内的人工模糊程度。

s5：对于步骤s4中比值当分母接近或者为0时会出现异常，在计算公式(3)中的比值时，利用sigmoid函数将比值映射到[0,0.5]的范围内，具体的有，

其中，r表示比值和则在公式(3)中求解的两个比值中最大值转化为映射后数值的最大值。

s6：获得候选的可疑人工模糊区域：通过步骤s3、s4以及s5可以得到像素点p在dct域内的人工模糊程度的距离表示tf。对图像中所有的像素点进行该操作即可得到测试图像在dct域内的人工模糊度估计映射图m，而后使用阈值t1对该映射图进行分割，

其中t′f代表tf二值分割后的结果，标记为‘1’所组成的区域就是当前检测出的人工模糊的区域。而后利用图像形态学的开闭操作进行去燥和填充孔洞，即可在二值图中提取出标记为‘1’候选的可疑人工模糊区域其中rb是候选的可疑人工模糊区域的集合。本实施例中参数t1＝0.3989。

s7：图像空域中候选区域的筛选：在步骤s6中得到的每一个可疑候选区域在二次模糊后的图像ib中找到其对应的区域其中rg是二次模糊后的图像ib中对应的集合。在空域内利用空域纹理的描述子包括异常灰度级统计、平滑度以及信息熵作为特征f计算和的相似度s，有

其中fb和fg分别是和计算特征f所得到的结果。然后以t2为阈值进行最终的判断得到人工模糊区域的定位结果，有

其中是求出的最终人工模糊区域的集合，表示对求出的s值。在二值的篡改定位图中将所对应的区域标记为‘1’，其他区域标记为‘0’，即可得到最终的检测定位图。本实施例中参数t2＝1.333。

进一步地，所述步骤s7中特征f计算的处理过程如下：

s701：计算给定区域的图像梯度：对于区域上每一个像素点计算它的梯度幅值，有

其中[.]表示取整操作。

s702：异常灰度级统计：利用灰度共生矩阵对提取特征，有

其中δ(·)为脉冲函数，用来对像素对进行计数，可写为

式子中分别当di＝0，dj＝1和di＝1，dj＝0时可以得到m在水平和垂直方向上相邻梯度幅值的统计，分别用mh和mv来表示。将统计中较大幅值的统计对视为异常的点对，m中的异常元素可以表示为则对异常点对的统计有,

本实施例中参数τ＝3。

s703：给定区域整体平滑度估计：对于给定图像区域内的所有点，首先求得梯度值的标准差σ，

其中n是中点的个数由给定图像区域的大小决定，是的均值，l是梯度幅值的变化范围。根据σ求得候选区域的平滑度，用f2表示有

s704：纹理信息熵估计：给定图像区域梯度的纹理信息熵可以表示为

其中是所出现的概率，可以通过计算直方图的频数求得，

其中指的是所对应直方图的统计数。直方图的统计间隔为1。因为f1、f2以及f3取值都正相关与s，采用如下式子将三种描述子结合起来，

得到候选区域在图像空域内的特征值。

如图2至5所示，是基于多域耦合人工模糊检测方法的实验效果。图2为人工模糊后的待检测的篡改试图像，图3为该篡改图像的掩膜图，图4是在dct域上得到人工模糊度估计的映射图，图5是最终的篡改定位的结果，从图中可以明显地看出定位效果比较理想。

相同或相似的标号对应相同或相似的部件；

附图中描述位置关系的用于仅用于示例性说明，不能理解为对本专利的限制；

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。