一种图像篡改盲取证方法

摘要

一种基于Sobel边缘检测与图像块亮度特征的图像篡改盲取证方法，其特征在于，包括以下步骤：对待检测的可疑图像转换为灰度图像，对灰度图像I(i,j)进行卷积处理，获得该图像的梯度图像G(i,j)，对G(i,j)进行阙值分割得到梯度图像G(i,j)，对梯度图像G(i,j)进行二值化处理得到二值化图像W(i,j)，对二值化图像W(i,j)进行1像素化，通过判断两子图像集的相似度，本发明通过对图像集进行分块处理，通过亮度均值排序能够有效提高取证算法效率，而且通过比较图像块的相似亮度值，进一步提高了准确率，而且进一步提高了图像取证的效率；本发明解决了由于大比例的缩放的克隆篡改行为而导致不能准确检测篡改图片的问题，通过归一化处理后的图像块亮度值，从而进一步提高了不同亮度值的克隆图像的检测结果。

权利要求

1.一种图像篡改盲取证方法，其特征在于，包括以下步骤：

S1)、对待检测的可疑图像进行灰度判决，如果为非灰度图像，则转换为灰度图像；

S2)、通过垂直和水平方向的两个3×3的滤波器分别与待检测灰度图像I(i,j)进行卷积处理，获得该图像的梯度图像G(i,j)；对梯度图像G(i,j)进行阈值分割，其计算表达式为： 其中，k1取值为4，m(G(i,j)为梯度图像G(i,j)所有像素点灰度的平均值；以遍历图像中的像素点，灰度大于阈值的置为白点，灰度小于阈值的则置为黑点的方式对梯度图像G(i,j)进行二值化处理，得到二值化图像W(i,j)；

S3)、利用非极大值抑制算法对二值化图像W(i,j)进行1像素化，对于二值化图像W(i,j)内的任意一点p(i,j)采用从左至右、从上至下的顺序遍历二值化图像W(i,j)，二值化图像W(i,j)中的任意一点p(i,j)竖直边缘的梯度值大于该点左边和右边的点的梯度值或水平边缘的梯度值大于该点上端和下端的点的梯度值，则设为白点，反之则设为黑点，从而得到非极大值抑制后的二值化图像P(i,j)；

S4)、按边缘闭合区域分割图像P(i,j)，形成图像集A＝(A1,A2,A3,…Ai)；

S5)、将待检测灰度图像I(i,j)进行分块处理，采用从左至右，从上到下的顺序对每块图像进行标号，每一个图像块标记为Mi,j；

S6)、将所有图像块Mi,j与图像集A＝(A1,A2,A3,…Ai)的子图像集进行相应置换，其中，子图像集Ai＝(Mi,1,Mi,2,…Mi,j)；

S6)、计算得到每个子图像集Ai的亮度均值Ci，并按照亮度均值Ci从大到小的顺序对图像集A的子图像集进行排序；

S7)、对每个子图像集Ai的每个图像块Mi,j的亮度值Ci,j进行归一化处理，每个图像块Mi,j归一化处理后的亮度值为Di,j，其计算式为： 其中， 为子图像集Ai的最大亮度值， 为子图像集Ai的最小亮度值；

S8)、将每个图像块Mi,j归一化处理后的亮度值Di,j按照从左至右、从上至下的顺序排列；

S9)、循环选取排序后的亮度均值Ci相近的子图像集Aq和Ap，计算子图像集Aq和Ap的相似度Sp,q，其表达式为： 假设子图像集Aq包含的图像块数量小于子图像集Ap包含的图像块数量，其中，nq为子图像集Aq的图像块的数量，qi为子图像集Aq的第i个图像块归一化后的亮度值， 为子图像集Aq包含图像块的数量nq和子图像集Ap图像块的数量np的比值，pi×r为子图像集Ap的第[i*r]个图像块归一化后的亮度值；

S10)、如果子图像集Aq和Ap的相似度Sp,q≤δ，则子图像集Aq和Ap不相似，按照步骤S9)继续处理；如果子图像集Aq和Ap的相似度Sp,q＞δ，则子图像集Aq和Ap相似，将相似的图像集Aq和Ap标记为克隆图像，输出取证结果。

2.根据权利要求1所述的图像篡改盲取证方法，其特征在于：步骤S10)中δ取值0.1。

说明书

一种图像篡改盲取证方法

技术领域

[0001] 本发明涉及一种信息安全与信息处理技术领域，尤其是一种图像篡改盲取证方法。

背景技术

[0002] 图像被动取证技术，也称为图像盲取证技术，是一个近年来快速发展的研究领域。它仅通过图像本身的特征鉴别其真伪，而不需要任何的数字签名或数字水印等先验知识，也不需要往相机设备中添加入水印生成程序，因此有着应用范围较广的优点；同时，由于不需要在图像中嵌入任何信息，因此也不会带来图像质量的降低。自然图像在生成时会受到成像设备内部软硬件或自然图像本身统计规律等因素的影响，输出的图像会存在某些自然图像独有的一致性和连续性固有特征，而篡改行为往往会打破这些图像固有特征的一致性和连续性。

[0003] 图像盲取证技术通过曝光这些被破坏的自然图像特征即可鉴别出图像的真伪，甚至定位出篡改的位置，逼真的篡改图像在视觉上虽然让人们不易察觉，但是使用各种相应技术来检测图像受到的篡改类型，甚至可以定位出篡改的具体位置。

[0004] 在数字图像篡改中最常用到的手段之一是图像克隆，它通过复制拷贝同一幅图像中的某些区域以清除或者添加一些重要目标。由于在同一幅图像中克隆区域间的亮度、色彩度都较为相近，进行复制粘贴不会引起明显的视觉变化，因此比较难于准确检测出来。近十年来，大部分数字图像取证研究都针对克隆类型的篡改。

[0005] Fridrich等人首先提出了一种通过对图像进行图像块分割的方法来进行图像克隆盲取证的技术，并被作为了一类经典算法。他们首先对图像进行分块操作；然后提取每个图像块的离散余弦变换(Discrete Cosine Transform，简称DCT)特征系数，并组成特征向量；最后通过对特征向量进行排序找出相似的图像块；

[0006] AC.Popescu和H.Farid随后提出了一种跟Fridrich等人类似的算法，其主要改进在于每个分割后的图像块的特征提取方面，他们使用主成份分析法(PCA)代替了DCT，以得到一个维数更小的特征系数，从而提高特征向量匹配的效率。

[0007] 但上述方法无法检测一些经过旋转操作的克隆物体，这也影响了该方法的推广应用前景。

[0008] 随后的基于图像块匹配的取证技术形成了两个研究分支，一部分学者着重提高对图像几何变换攻击的鲁棒性，而另一部分文献则针对效率的改进进行研究。

[0009] 在提高克隆检测算法的鲁棒性方面，B.Mahdian和S.Saic首先研究数字图像在经过模糊卷积后的特征不变性问题，从数学角度建立了一个模糊不变性等式：B(f)＝B(D(f))，其中D表示模糊，f表示图像，模糊不变性问题就是找到一个函数B，能够满足以上等式，通过数学推导和实验分析，他们最终确定了一个24维的模糊不变特征，该特征可以较好的减少模糊，噪点或镜像变换等后期操作带来的影响。他们对图像块进行模糊不变量运算，根据模糊不变特征的相似性鉴别出克隆的图像块；

[0010] Wu等提出了一种能抵御90°和270°的普通旋转攻击的克隆图像取证算法，他们首先将图像转换为YUV色彩模式并提取Y分量得到灰度图像；然后对图像分块并计算每个图像块的平均亮度值；最后根据每个图像块的平均亮度与其相邻8个图像块的平均亮度的之间的关系来判断图像块是否被克隆。

[0011] Ryu等人则通过对图像块进行傅立叶-梅林变化(FMT)及一维对数投影的方法来得到一个较鲁棒的特征向量，以此来提高克隆图像块检测的鲁棒性。

[0012] 但是，上述基于图像块的克隆盲取证方法会产生大量的图像块，特别是在检测高分辨率图像时，图像块的数量是非常巨大的，增加了取证算法的运算时间。

[0013] Khan等人首先对检测图像进行离散小波变换(Discrete Wavelet Transform，简称DWT)，然后提取低频子频带信息再进行克隆检测。类似的方法还有Li等人的DWT-SVD方法，Ghorbani等人的DWT-DCT(QCD)方法，Huang等人在Fridrich等人算法的基础上提出了一种改进的基于离散余弦变换的方法，他们通过引入一个阈值p，来控制图像块特征向量的维数，上述方法平均能减少p％的匹配时间，但是随着p值的增加，越来越多的特征维被舍弃，算法的准确率也随之降低。

发明内容

[0014] 针对现有技术的不足，本发明提供一种基于Sobel边缘检测与图像块亮度特征的图像篡改盲取证方法。

[0015] 本发明的技术方案为：一种图像篡改盲取证方法，其特征在于，包括以下步骤：

[0016] S1)、对待检测的可疑图像进行灰度判决，如果为非灰度图像，则转换为灰度图像，其表达式为：Gray＝0.299R+0.587G+0.114B；

[0017] S2)、通过垂直和水平方向的两个3×3的滤波器分别与待检测灰度图像I(i,j)进行卷积处理，获得该图像的梯度图像G(i,j)；对梯度图像G(i,j)进行阈值分割，其计算表达式为： 其中，k1取值为4，m(G(i,j)为梯度图像G(i,j)所有像素点灰度的平均值；以遍历图像中的像素点，灰度大于阈值的置为白点，灰度小于阈值的则置为黑点的方式对梯度图像G(i,j)进行二值化处理，得到二值化图像W(i,j)；

[0018] S3)、利用非极大值抑制算法对二值化图像W(i,j)进行1像素化，对于二值化图像W(i,j)内的任意一点p(i,j)采用从左至右、从上至下的顺序遍历二值化图像W(i,j)，二值化图像W(i,j)中的任意一点p(i,j)竖直边缘的梯度值大于该点左边和右边的点的梯度值或水平边缘的梯度值大于该点上端和下端的点的梯度值，则设为白点，反之则设为黑点，从而得到非极大值抑制后的二值图像P(i,j)；

[0019] S4)、按边缘闭合区域分割图像P(i,j)，形成图像集A＝(A1,A2,A3,…Ai)；

[0020] S5)、将待检测灰度图像I(i,j)进行分块处理，采用从左至右，从上到下的顺序对每块图像进行标号，每一个图像块标记为Mi,j；

[0021] S6)、将所有图像块Mi,j与图像集A＝(A1,A2,A3,…Ai)的子图像集进行相应置换，其中，子图像集Ai＝(Mi,1,Mi,2,…Mi,j)；

[0022] S6)、计算得到每个子图像集Ai的亮度均值Ci，并按照亮度均值Ci从大到小的顺序对图像集A的子图像集进行排序；

[0023] S7)、对每个子图像集Ai的每个图像块Mi,j的亮度值Ci,j进行归一化处理，每个图像块Mi,j归一化处理后的亮度值为Di,j，其计算式为： 其中，为子图像集Ai的最大亮度值， 为子图像集Ai的最小亮度值；

[0024] S8)、将每个图像块Mi,j归一化处理后的亮度值Di,j按照从左至右、从上至下的顺序排列；

[0025] S9)、循环选取排序后的亮度均值Ci相近的子图像集Aq和Ap，计算子图像集Aq和Ap的相似度Sp,q，其表达式为： 假设子图像集Aq包含的图像块数量小于子图像集Ap包含的图像块数量，其中，nq为子图像集Aq的图像块的数量，qi为子图像集Aq的第i个图像块归一化后的亮度值， 为子图像集Aq包含图像块的数量nq和子图像集Ap图像块的数量np的比值；

[0026] S10)、如果子图像集Aq和Ap的相似度Sp,q≤δ，则子图像集Aq和Ap不相似，按照步骤S9)继续处理；如果子图像集Aq和Ap的相似度Sp,q＞δ，则子图像集Aq和Ap相似，将相似的图像集Aq和Ap标记为克隆图像，输出取证结果。

[0027] 上述技术方案中，步骤S10)中δ取值0.1。

[0028] 本发明的有益效果为：通过对图像集进行分块处理，通过亮度均值排序能够有效提高取证算法效率，而且通过比较图像块的相似亮度值，进一步提高了准确率，而且进一步提高了图像取证的效率；本发明解决了由于大比例的缩放的克隆篡改行为而导致不能准确检测篡改图片的问题，通过归一化处理后的图像块亮度值，从而进一步提高了不同亮度值的克隆图像的检测结果。

附图说明

[0029] 图1为本发明的流程示意图。

具体实施方式

[0030] 下面结合附图对本发明的具体实施方式作进一步说明：

[0031] 如图1所示，一种图像篡改盲取证方法，其特征在于，包括以下步骤：

[0032] S1)、对待检测的可疑图像进行灰度判决，如果为非灰度图像，则转换为灰度图像，其表达式为：Gray=0.299R+0.587G+0.114B；

[0033] S2)、通过垂直和水平方向的两个3×3的滤波器分别与待检测灰度图像I(i,j)进行卷积处理，获得该图像的梯度图像G(i,j)；对梯度图像G(i,j)进行阈值分割，其计算表达式为： 其中，k1取值为4，m(G(i,j)为梯度图像G(i,j)所有像点灰度的平均值；以遍历图像中的像素点，灰度大于阈值的置为白点，灰度小于阈值的则置为黑点的方式对梯度图像G(i,j)进行二值化处理，得到二值化图像W(i,j)；

[0034] S3)、利用非极大值抑制算法对二值化图像W(i,j)进行1像素化，对于二值化图像W(i,j)内的任意一点p(i,j)采用从左至右、从上至下的顺序遍历二值化图像W(i,j)，二值化图像W(i,j)中的任意一点p(i,j)竖直边缘的梯度值大于该点左边和右边的点的梯度值或水平边缘的梯度值大于该点上端和下端的点的梯度值，则设为白点，反之则设为黑点，从而得到非极大值抑制后的二值图像P(i,j)；

[0035] S4)、按边缘闭合区域分割图像P(i,j)，形成图像集A＝(A1,A2,A3,…Ai)；

[0036] S5)、将待检测灰度图像I(i,j)进行分块处理，采用从左至右，从上到下的顺序对每块图像进行标号，每一个图像块标记为Mi,j；

[0037] S6)、将所有图像块Mi,j与图像集A＝(A1,A2,A3,…Ai)的子图像集进行相应置换，其中，子图像集Ai＝(Mi,1,Mi,2,…Mi,j)；

[0038] S6)、计算得到每个子图像集Ai的亮度均值Ci，并按照亮度均值Ci从大到小的顺序对图像集A的子图像集进行排序；

[0039] S7)、对每个子图像集Ai的每个图像块Mi,j的亮度值Ci,j进行归一化处理，每个图像块Mi,j归一化处理后的亮度值为Di,j，其计算式为： 其中，为子图像集Ai的最大亮度值， 为子图像集Ai的最小亮度值；

[0040] S8)、将每个图像块Mi,j归一化处理后的亮度值Di,j按照从左至右、从上至下的顺序排列；

[0041] S9)、循环选取排序后的亮度均值Ci相近的子图像集Aq和Ap，计算子图像集Aq和Ap的相似度Sp,q，其表达式为： 假设子图像集Aq包含的图像块数量小于子图像集Ap包含的图像块数量，其中，nq为子图像集Aq的图像块的数量，qi为子图像集Aq的第i个图像块归一化后的亮度值， 为子图像集Aq包含图像块的数量nq和子图像集Ap图像块的数量np的比值；

[0042] S10)、如果子图像集Aq和Ap的相似度Sp,q≤δ，则子图像集Aq和Ap不相似，按照步骤S9)继续处理；如果子图像集Aq和Ap的相似度Sp,q＞δ，则子图像集Aq和Ap相似，将相似的图像集Aq和Ap标记为克隆图像，输出取证结果。

[0043] 上述技术方案中，步骤S10)中δ取值0.1。

[0044] 上述实施例和说明书中描述的只是说明本发明的原理和最佳实施例，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。