CN104933721A - 基于颜色滤波阵列特性的拼接图像篡改检测方法 - Google Patents

Abstract

本发明公开的基于颜色滤波阵列特性的拼接图像篡改检测方法，包括以下步骤：第1步，将待检测图像分成多个图像块的预处理；第2步，估算原始图像模式；第3步，利用边缘检测算子进行篡改定位检测。本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法利用颜色滤波阵列插值所引入的图像像素间的周期性相关模式的变化或差异性特点，进行拼接图像篡改检测，其不仅能够检测出图像是否被拼接篡改，而且能够检测被篡改区域的位置；在篡改定位阶段由于引进了Canny算子，使算法具有较高的篡改定位精度，即可以精确地定位出被篡改区域的边缘，并有效地拟制了虚假边缘；对内容保持的图像处理操作如JPEG压缩、不同类型的滤波、加噪处理等，具有较好的鲁棒性。

Description

基于颜色滤波阵列特性的拼接图像篡改检测方法

技术领域

本发明属于图像处理技术领域，具体涉及一种基于颜色滤波阵列特性的拼接图像篡改检测方法。

背景技术

在数字成像技术日新月异的发展过程中，数码照片被应用在我们的生活中的各个方面。然而，各种各样图像处理软件的广泛应用，可以方便地对图像进行一些处理操作，如局部修改、拼接、润饰等计算机处理，使得篡改图像无处不在，造成了数字图像的内容真实性变得不再可靠，无法作为一些法律案件、新闻传媒、科研成果、医疗诊断以及金融事件的强有力的证据。因此，如何检测数字图像内容的真实性已成为近年来法律界和信息产业界所面临的一个重要的热点问题和迫切需要解决的难点问题。展开对数字图像内容真实性的研究，对维护互联网的公共信任秩序、维护法律公正、新闻诚信、科学诚信等，具有十分重要的意义。

图像拼接是一种最普遍的图像篡改技术，是指把不同图像的部分内容拼接在一起生成合成图像，以伪造不存在的场景。拼接后的图像往往进行了一些后处理，如模糊、添加噪声、JPEG压缩，旋转/缩放等几何操作，以制造以假乱真的效果，使得人眼根本无法辨别真伪，机器识别也变得更加困难。

对于数码相机获取的全彩色图像，颜色滤波阵列(Color Filter Array，简称CFA)的运用为拼接图像的检测提供了理论基础：即CFA插值操作使图像相邻像素间具有相关性，拼接操作会破坏或者改变这种相关性模式。因此，可以通过在图像中检测这种相关模式的改变来追踪拼接伪造的痕迹。

首次将CFA插值所引入的图像相邻像素之间的周期性应用于数字图像篡改检测的方法出现在Popescu和Farid的文献中，作者首先估计了CFA插值模型的系数及插值后验概率图，并对后验概率图进行二维离散傅里叶变换，实现了空域到频域的转换，最后通过观察峰值的分布是否具有周期性实现篡改检测，该方法能够检测图像是否经历了拼接篡改，但不能检测被拼接的区域，而且对JPEG压缩不具有鲁棒性。除此之外，Dirik和Memon基于CFA的结构特征也提出了两种篡改检测方法：第一种，由于不同模式结构的CFA，通过插值得到的像素的残留误差不同，由此就可以判断待检测图像所使用的CFA模式结构，进而实现篡改检测与定位；第二种，给定一种相同模式结构的CFA，计算与之对应的由传感器直接获得的像素和由CFA插值得到的像素位置处的噪声强度比，最终实现篡改检测定位。这两种方法的不足之处也在于对JPEG压缩不具有鲁棒性。

通过大量调研我们发现，现有的基于CFA插值模式的图像拼接检测方法仍存在许多缺点，主要体现在两个方面：一是一些算法只能检测出图像是否经过了拼接操作，但无法确定被伪造区域的位置；二是一些算法虽然可以确定被伪造区域的位置，但对于JPEG压缩的鲁棒性较差，而JPEG是一种常用的图像压缩格式，目前使用的很多图像都是JPEG格式。因此，现存方法远远不能够满足图像取证的实际需求，发明篡改检测率高，篡改定位准确并且鲁棒的取证方法迫在眉睫。

发明内容

本发明的目的在于提供一种基于颜色滤波阵列特性的拼接图像篡改检测方法，解决了现有技术中不能精确定位被拼接的图像区域以及算法不具有鲁棒性的问题，其能够准确定位出拼接伪造的数字图像区域，并对于JPEG压缩、添加噪声、滤波、伽马校正等内容保持的图像处理操作具有鲁棒性。

本发明所采用的技术方案是：基于颜色滤波阵列特性的拼接图像篡改检测方法，包括以下步骤：

第1步，将待检测图像分成多个图像块的预处理；

第2步，估算原始图像模式；

第3步，利用边缘检测算子进行篡改定位检测。

本发明的特点还在于，

第1步，将待检测图像分成多个图像块的预处理：

将待测图像按像素点分为M×N大小的矩阵I，采用CFA差值模型将待检测图像的绿色分量记为I_CFA，将I_CFA划分为不重叠的64×64的图像块，即得到M×N/64²个图像块，用表示第k块：

$I_{CFA}^k=\left(\begin{array}{cccc}{I}_{CFA}^{11}& I_{CFA}^{12}& ...& I_{CFA}^{1s}\\ I_{CFA}^{21}& I_{CFA}^{22}& ...& I_{CFA}^{2s}\\ ...& ...& ...& ...\\ I_{CFA}^{s1}& I_{CFA}^{s2}& ...& I_{CFA}^{ss}\end{array}\right),(k=0,1,2,...,M\×N/({64}^2-1),s=64)---\left(1\right);$

第2步，估算原始图像模式时将I_CFA的像素分为M₁和M₂两类，其中M₁表示通过插值得到的像素值，M₂表示通过传感器直接获得的像素值，I_CFA(m,n)表示插值点(m,n)处的像素值，具体步骤如下：

第2.1步，对每一个图像块中插值点(m,n)处的像素值建立线性插值模型：

$I_{CFA}^k(m,n)=\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}^k(m+i,n+j)+r(m,n)---\left(2\right)$

其中，参数 $\stackrel{\‾}{v}=\left\{v(i,j)\right|-N_0\≤i,j\≤N_0\left\}\right(N_0\∈N^{+},v(0,0)=0),$ 参数r(m,n)是服从均值为0、方差为σ²正态分布的残余误差；

第2.2步，对参数进行初始化，令N₀＝1，方差σ＝2，属于M₂的条件概率为P₀＝1/256，对每一个图像块利用EM算法估算出其插值系数，记为计算所有的平均值，记为m_v：

$m_{\stackrel{\→}{v}}=\frac{{64}^2}{M\×N}\underset{k=1}{\overset{M\×N/{64}^2}{\Σ}}{\stackrel{\→}{v}}^k---\left(3\right);$

第2.3步，利用构造最终的插值系数矩阵，记为H：

$H=\left[\begin{array}{ccc}{m}_{\stackrel{\‾}{v}}\left(1\right)& m_{\stackrel{\‾}{v}}\left(2\right)& m_{\stackrel{\‾}{v}}\left(3\right)\\ m_{\stackrel{\‾}{v}}\left(4\right)& 0& m_{\stackrel{\‾}{v}}\left(5\right)\\ m_{\stackrel{\‾}{v}}\left(6\right)& m_{\stackrel{\‾}{v}}\left(7\right)& m_{\stackrel{\‾}{v}}\left(8\right)\end{array}\right]---\left(4\right);$

第2.4步，记绿色分量I_CFA插值点(m,n)的邻域矩阵为

$I_{CFA}^{m,n}=\left[\begin{array}{ccc}{I}_{CFA}\left(m-1,n-1\right)& I_{CFA}\left(m-1,n\right)& I_{CFA}\left(m-1,n+1\right)\\ I_{CFA}\left(m,n-1\right)& I_{CFA}\left(m,n\right)& I_{CFA}\left(m,n+1\right)\\ I_{CFA}\left(m+1,n-1\right)& I_{CFA}\left(m+1,n\right)& I_{CFA}\left(m+1,n+1\right)\end{array}\right]---\left(5\right);$

第2.5步，利用最终的插值系数矩阵H和差值点(m,n)邻域矩阵得到原始图像模式I'_CFA内的像素值I'_CFA(m,n)：

${I^{\′}}_{CFA}(m,n)=\mathrm{\Σ}\mathrm{\Σ}H\·I_{CFA}^{m,n}---\left(6\right);$

第3步利用边缘检测算子进行篡改定位检测具体步骤如下：

第3.1步，定义新矩阵I_C，其元素为I_CFA与I'_CFA的对应元素差的平方：

$I_C=\left[\begin{array}{cccc}{\left(I_{CFA}^{11}-I_{CFA}^{\′11}\right)}^2& {\left(I_{CFA}^{12}-I_{CFA}^{\′12}\right)}^2& ...& {\left(I_{CFA}^{1N}-I_{CFA}^{\′1N}\right)}^2\\ {\left(I_{CFA}^{21}-I_{CFA}^{\′21}\right)}^2& {\left(I_{CFA}^{22}-I_{CFA}^{\′22}\right)}^2& ...& {\left(I_{CFA}^{2N}-I_{CFA}^{\′2N}\right)}^2\\ ...& ...& ...& ...\\ {\left(I_{CFA}^{M1}-I_{CFA}^{\′M1}\right)}^2& {\left(I_{CFA}^{M2}-I_{CFA}^{\′M2}\right)}^2& ...& {\left(I_{CFA}^{MN}-I_{CFA}^{\′MN}\right)}^2\end{array}\right]---\left(7\right);$

第3.2步，对I_C进行二值化处理得到I'_C，然后利用Canny边缘检测算子对I'_C进行边缘检测，得到初步篡改定位结果I_L：

I_L＝E(I'_C,'canny')       (8)。

第3步还包括：

第3.3步，将初步篡改定位结果I_L使用形态学闭运算进行处理，得到最终的篡改定位结果I_Lend：

I_Lend＝imclose(I_L,SE)       (9)，

其中，SE是结构元素。

本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法利用颜色滤波阵列插值所引入的图像像素间的周期性相关模式的变化或差异性特点，进行拼接图像篡改检测，解决了现有技术中不能精确定位被拼接的图像区域以及算法不具有鲁棒性的问题，并具有以下有益效果：

(1)不仅能够检测出图像是否被拼接篡改，而且能够检测被篡改区域的位置；

(2)在篡改定位阶段由于引进了Canny算子，使算法具有较高的篡改定位精度，即可以精确地定位出被篡改区域的边缘，并有效地拟制了虚假边缘；

(3)对内容保持的图像处理操作如不同质量因子的JPEG压缩、不同类型的滤波、加噪处理等，具有较好的鲁棒性。

附图说明

图1a是本发明的一个实施例的原始测试图像；

图1b是图1a中拼接了其他图像部分内容生成的拼接篡改图像；

图1c是对图1b的检测结果图像；

图2a是本发明的另一实施例的原始测试图像；

图2b是图2a中拼接了其他图像部分内容生成的拼接篡改图像；

图2c是对图2b的检测结果图像；

图3b是在图3a中拼接了其他图像部分内容生成拼接篡改图像后再进行JPEG(QF＝80)压缩后的图像；

图3c是对图3b的检测结果图像；

图4a是本发明的另一实施例的原始测试图像；

图4b是在图4a中拼接了其他图像部分内容生成拼接篡改图像后再进行JPEG(QF＝60)压缩后的图像；

图4c是对图4b的检测结果图像；

图5a是本发明的另一实施例的原始测试图像；

图5b是在图5a中拼接了其他图像部分内容生成拼接篡改图像后再进行JPEG(QF＝40)压缩后的图像；

图5c是对图5b的检测结果图像；

图6a是本发明的另一实施例的原始测试图像；

图6b是在图6a中拼接了其他图像部分内容生成拼接篡改图像后再进行median(3×3)滤波后的图像；

图6c是对图6b的检测结果图像；

图7a是本发明的另一实施例的原始测试图像；

图7b是在图7a中拼接了其他图像部分内容生成拼接篡改图像后再进行wiener(3×3)滤波后的图像；

图7c是对图7b的检测结果图像；

图8a是本发明的另一实施例的原始测试图像；

图8b是在图8a中拼接了其他图像部分内容生成拼接篡改图像后再加入椒盐噪声(噪声因子为0.0006)后的图像；

图8c是对图8b的检测结果图像；

图9a是本发明的另一实施例的原始测试图像；

图9b是在图9a中拼接了其他图像部分内容生成拼接篡改图像后再加入椒盐噪声(噪声因子为0.001)后的图像；

图9c是对图9b的检测结果图像；

图10a是本发明的另一实施例的原始测试图像；

图10b是在图10a中拼接了其他图像部分内容生成拼接篡改图像后再进行伽马校正(校正因子为0.8)后的图像；

图10c是对图10b的检测结果图像。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

本发明提供的基于颜色滤波阵列特性的拼接图像篡改检测方法，包括以下步骤：

第1步，将待检测图像分成多个图像块的预处理：

将待测图像按像素点分为M×N大小的矩阵I，采用CFA差值模型将待检测图像的绿色分量记为I_CFA，将I_CFA划分为不重叠的64×64的图像块，即得到M×N/64²个图像块，用表示第k块：

$I_{CFA}^k=\left(\begin{array}{cccc}{I}_{CFA}^{11}& I_{CFA}^{12}& ...& I_{CFA}^{1s}\\ I_{CFA}^{21}& I_{CFA}^{22}& ...& I_{CFA}^{2s}\\ ...& ...& ...& ...\\ I_{CFA}^{s1}& I_{CFA}^{s2}& ...& I_{CFA}^{ss}\end{array}\right),(k=0,1,2,...,M\×N/({64}^2-1),s=64)---\left(1\right)$

第2步，估算原始图像模式：

将I_CFA的像素分为M₁和M₂两类，其中M₁表示通过插值得到的像素值，M₂表示通过传感器直接获得的像素值，I_CFA(m,n)表示插值点(m,n)处的像素值，具体步骤如下：

第2.1步，对每一个图像块中插值点(m,n)处的像素值建立线性插值模型：

$I_{CFA}^k(m,n)=\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}^k(m+i,n+j)+r(m,n)---\left(2\right)$

其中，参数 $\stackrel{\‾}{v}=\left\{v(i,j)\right|-N_0\≤i,j\≤N_0\left\}\right(N_0\∈N^{+},v(0,0)=0),$ 参数r(m,n)是服从均值为0、方差为σ²正态分布的残余误差；

第2.2步，对参数进行初始化，令N₀＝1，即与其相邻的8个像素值相关，方差σ＝2，属于M₂的条件概率为P₀＝1/256，对每一个图像块利用EM算法估算出其插值系数，记为具体地利用EM算法估算插值系数的步骤如下：

由于上述模型的系数和残余误差的方差σ²，一般用极大似然估计来估计，为了解决极大似然估计的迭代问题，使用期望最大化(简称EM)算法求得。该算法以两步迭代为过程，最终收敛为目的，分为E步和M步，E步估计插值点(m,n)属于M₁或M₂的概率，M步估计和σ²，进而估计出相邻像素间相关性的具体模式，

E步，已知插值点(m,n)处的像素值I_CFA(m,n)，由贝叶斯法则可以得到I_CFA(m,n)属于M₁的后验概率表示如下：

$\mathrm{Pr}\{I_{CFA}(m,n)\∈M_1|I_{CFA}(m,n)\}=\frac{\mathrm{Pr}\left\{I_{CFA}(m,n)\right|I_{CFA}(m,n)\∈M_1\left\}\mathrm{Pr}\right\{I_{CFA}(m,n)\∈M_1\}}{{\mathrm{\Σ}}_{i=1}^2\mathrm{Pr}\left\{I_{CFA}(m,n)\right|I_{CFA}(m,n)\∈M_i\left\}\mathrm{Pr}\right\{I_{CFA}(m,n)\∈M_i\}}---\left(10\right)$

这里假设先验概率Pr{I_CFA(m,n)∈M₁}和Pr{I_CFA(m,n)∈M₂}是常数并令初始值为1/2，I_CFA(m,n)属于M₂的条件概率P₀≡Pr{I_CFA(m,n)|I_CFA(m,n)∈M₂}服从均匀分布，即P₀等于I_CFA(m,n)可能取值范围的倒数，I_CFA(m,n)属于M₁的条件概率P(m,n)≡Pr{I_CFA(m,n)|I_CFA(m,n)∈M₁}表示如下：

$\mathrm{Pr}\left\{I_{CFA}\left(m,n\right)|I_{CFA}\left(m,n\right)\∈M_1\right\}=\frac{1}{\mathrm{\σ}\sqrt{2\mathrm{\π}}}\exp \[-\frac{1}{2{\mathrm{\σ}}^2}{\left(I_{CFA}\left(m,n\right)-\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v\left(i,j\right)I_{CFA}\left(m+i,n+j\right)\right)}^2\]---\left(11\right)$

其中，该步在估计模型系数时，第一次迭代的模型系数是随机选取的；

M步，通过对下面的二次误差函数进行最小化，使用加权最小二乘法重新估算出一组稳定的模型系数

$E\left(\stackrel{\→}{v}\right)=\underset{m,n}{\Σ}w\left(m,n\right){\left(I_{CFA}\left(m,n\right)-\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v\left(i,j\right)I_{CFA}\left(m+i,n+j\right)\right)}^2---\left(12\right)$

其中， $r(m,n)=|I_{CFA}(m,n)-\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}(m+i,n+j)|$ 代表差值点像素值的残余误差，w(m,n)≡Pr{I_CFA(m,n)∈M₁|I_CFA(m,n)}，即I_CFA(m,n)属于M₁的后验概率，

对中的一个元素求偏导，并设得如下两个线性方程：

$-2\underset{m,n}{\Σ}w(m,n)I_{CFA}(m+s,n+t)(I_{CFA}(m,n)-\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}(m+i,n+j))=0---\left(13\right)$

$\underset{m,n}{\Σ}w(m,n)I_{CFA}(m+s,n+t)\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}(m+i,n+j)=\underset{m,n}{\Σ}w(m,n)I_{CFA}(m+s,n+t)I_{CFA}(m,n)---\left(14\right)$

整理等式左边可得：

$\underset{i,j=-N_0}{\overset{N_0}{\Σ}}v\left(i,j\right)\left(\underset{m,n}{\Σ}w\left(m,n\right)I_{CFA}\left(m+s,n+t\right)I_{CFA}\left(m+i,n+j\right)\right)=\underset{m,n}{\Σ}w\left(m,n\right)I_{CFA}\left(m+s,n+t\right)I_{CFA}\left(m,n\right)---\left(15\right)$

对中所有的元素求偏导，就可以得到由一系列线性方程构成的方程组，对该方程组求解并带入初始化赋值即可重新得到一组系数。

为了得到稳定的系数，在E步和M步迭代过程中，对于第a次迭代，若则不稳定，令a＝a+1；否则，停止迭代，为最终求得的稳定的插值系数

为了使插值系数更稳定、更精确，因此计算所有的平均值，记为

$m_{\stackrel{\→}{v}}=\frac{{64}^2}{M\×N}\underset{k=1}{\overset{M\×N/{64}^2}{\Σ}}{\stackrel{\→}{v}}^k---\left(3\right);$

第2.3步，利用构造最终的插值系数矩阵，记为H：

$H=\left[\begin{array}{ccc}{m}_{\stackrel{\‾}{v}}\left(1\right)& m_{\stackrel{\‾}{v}}\left(2\right)& m_{\stackrel{\‾}{v}}\left(3\right)\\ m_{\stackrel{\‾}{v}}\left(4\right)& 0& m_{\stackrel{\‾}{v}}\left(5\right)\\ m_{\stackrel{\‾}{v}}\left(6\right)& m_{\stackrel{\‾}{v}}\left(7\right)& m_{\stackrel{\‾}{v}}\left(8\right)\end{array}\right]---\left(4\right);$

第2.4步，记绿色分量I_CFA插值点(m,n)的邻域矩阵为

$I_{CFA}^{m,n}=\left[\begin{array}{ccc}{I}_{CFA}\left(m-1,n-1\right)& I_{CFA}\left(m-1,n\right)& I_{CFA}\left(m-1,n+1\right)\\ I_{CFA}\left(m,n-1\right)& I_{CFA}\left(m,n\right)& I_{CFA}\left(m,n+1\right)\\ I_{CFA}\left(m+1,n-1\right)& I_{CFA}\left(m+1,n\right)& I_{CFA}\left(m+1,n+1\right)\end{array}\right]---\left(5\right);$

第2.5步，利用最终的插值系数矩阵H和差值点(m,n)邻域矩阵得到原始图像模式I'_CFA内的像素值I'_CFA(m,n)：

${I^{\′}}_{CFA}(m,n)=\mathrm{\Σ}\mathrm{\Σ}H\·I_{CFA}^{m,n}---\left(6\right);$

第3步，由于图像拼接会引入来自其它图像的区域，不同图像的CFA插值模式可能不尽相同，因此若测试图像为拼接图像，则其估算的原始图像模式I'_CFA中会存在不一致的区域。根据这一原理，结合I'_CFA和Canny算子检测拼接/合成图像的篡改区域，具体步骤如下：

第3.1步，定义新矩阵I_C，其元素为I_CFA与I'_CFA的对应元素差的平方：

$I_C=\left[\begin{array}{cccc}{\left(I_{CFA}^{11}-I_{CFA}^{\′11}\right)}^2& {\left(I_{CFA}^{12}-I_{CFA}^{\′12}\right)}^2& ...& {\left(I_{CFA}^{1N}-I_{CFA}^{\′1N}\right)}^2\\ {\left(I_{CFA}^{21}-I_{CFA}^{\′21}\right)}^2& {\left(I_{CFA}^{22}-I_{CFA}^{\′22}\right)}^2& ...& {\left(I_{CFA}^{2N}-I_{CFA}^{\′2N}\right)}^2\\ ...& ...& ...& ...\\ {\left(I_{CFA}^{M1}-I_{CFA}^{\′M1}\right)}^2& {\left(I_{CFA}^{M2}-I_{CFA}^{\′M2}\right)}^2& ...& {\left(I_{CFA}^{MN}-I_{CFA}^{\′MN}\right)}^2\end{array}\right]---\left(7\right);$

第3.2步，对I_C进行二值化处理得到I'_C，然后利用Canny边缘检测算子对I'_C进行边缘检测，得到初步篡改定位结果I_L：

I_L＝E(I'_C,'canny')      (8)；

第3.3步，将初步篡改定位结果I_L使用形态学闭运算进行处理，得到最终的篡改定位结果I_Lend：

I_Lend＝imclose(I_L,SE)       (9)，

其中，SE是结构元素。

本发明的实验验证过程及结果如下：

(1)篡改定位视觉效果

本实验的目的是测试本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法的准确性。实验所使用的图像选自国际通用的Columbia ImageSplicing Detection Evaluation Dataset[4](CISDED)图像数据库，用本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法对包含有不同大小拼接/合成区域的测试图像进行检测，实验步骤如下：

①图像预处理：提取待检测图像的绿色通道，对绿色通过图像分块，得到图像块

②估计图像模式：首先，对建立线性插值模型；然后，利用EM算法计算每个的一组模型系数计算所有的平均值并作为最终的插值系数；最后，通过对I_CFA进行双线性插值，估计得出I'_CFA；

③篡改定位：用I_CFA和I'_CFA建立矩阵I_C，接着用Canny算子对I_C进行边缘检测，定位出拼接区域，最终利用形态学处理定位结果。

本实验的目的是为了展示本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法的效果，即检测被拼接区域的位置的能力。实验中测试了大量大小不同的图像，图1a-图10c展示了实验结果，其中，用本发明的篡改定位方法检测出的拼接区域用二值图标出(注：原图是彩色的，很醒目，目前不醒目的原因是因为灰度图像造成的)。图1a为原始图像(来自CISDED)，图1b为图1a的拼接/合成篡改图像(来自CISDED)，其中的拼接区域是人眼视觉容易识别的，图1c为图1b的检测结果图像；图2b为图2a的拼接/合成篡改图像(其中，图2a和图2b均来自CISDED)，图2c分别为图2b的检测结果。

由实验结果可以看出，本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法对恶意篡改很敏感，而且能够精确地检测被拼接区域的位置。

(2)对常规图像处理操作的鲁棒性实验

常规图像处理操作是指内容保持的图像处理操作。本实验目的是检测本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法对内容保持的图像处理操作具有鲁棒性。

为此，我们选用CISDED数据库中的图像和部分自主获得的图像，选用的图像的特点是其拼接/合成篡改不易被肉眼察觉，需要利用定位算法定位出拼接区域。实验中对经历了不同的内容保持性图像处理操作的图像进行检测：

图3a是来自CISDED图像库的原始图像，图3b是在图3a中拼接了其它图像的部分内容生成拼接篡改图像，再进行JPEG(QF＝80)压缩图像，图3c是图3b的检测结果图像；

图4a是来自CISDED图像库的原始测试图像，图4b是在图4a中拼接了其它图像的部分内容生成拼接篡改图像，再进行JPEG(QF＝60)压缩生成的图像，图4c是图4b的检测结果图像；

图5a是自主获取的原始测试图像，图5b是在图5a中拼接了其它图像的部分内容生成拼接篡改图像，再进行JPEG(QF＝40)压缩生成的图像，图5c是图5b的检测结果图像；

图6a是来自CISDED图像库的原始测试图像，图6b是在图6a中拼接了其它图像的部分内容生成拼接篡改图像，再进行median(3×3)滤波生成的图像，图6c是图6b的检测结果图像；

图7a是自主获取的原始测试图像，图7b是在图7a中拼接了其它图像的部分内容生成拼接篡改图像，再进行wiener(3×3)滤波后的图像，图7c是图7b的检测结果图像；

图8a是来自CISDED图像库的原始测试图像，图8b是在图8a中拼接了其它图像的部分内容生成拼接篡改图像，再加入椒盐噪声(噪声因子为0.0006)生成的图像，图8c是图8b的检测结果图像；

图9a是自主获得的原始测试图像，图9b是在图9a中拼接了其它图像的部分内容生成拼接篡改图像，再加入椒盐噪声(噪声因子为0.001)生成的图像，图9c是图9b的测试结果图像；

图10a是来自CISDED图像库的原始测试图像，图10b是在图10a中拼接了其它图像的部分内容生成拼接篡改图像，再进行伽马校正(矫正因子为0.8)生成的图像，图10c是图10b的检测结果图像。

由实验结果可以看出，本发明的基于颜色滤波阵列特性的拼接图像篡改检测方法具有较好的鲁棒性。

Claims (5)

1.基于颜色滤波阵列特性的拼接图像篡改检测方法，其特征在于，包括以下步骤：

第1步，将待检测图像分成多个图像块的预处理；

第2步，估算原始图像模式；

第3步，利用边缘检测算子进行篡改定位检测。

2.如权利要求1所述的基于颜色滤波阵列特性的拼接图像篡改检测方法，其特征在于，所述第1步中将待检测图像分成多个图像块的预处理时，所述待测图像按像素点分为M×N大小的矩阵I，采用CFA差值模型将待检测图像的绿色分量记为I_CFA，将I_CFA划分为不重叠的64×64的图像块，即得到M×N/64²个图像块，用表示第k块：

$I_{CFA}^k=\left(\begin{array}{cccc}{I}_{CFA}^{11}& I_{CFA}^{12}& ...& I_{CFA}^{1s}\\ I_{CFA}^{21}& I_{CFA}^{22}& ...& I_{CFA}^{2s}\\ ...& ...& ...& ...\\ I_{CFA}^{s1}& I_{CFA}^{s2}& ...& I_{CFA}^{ss}\end{array}\right)(k=0,1,2,...,M\×N/({64}^2-1),s=64)---\left(1\right).$

3.如权利要求1所述的基于颜色滤波阵列特性的拼接图像篡改检测方法，其特征在于，所述第2步中估算原始图像模式时将I_CFA的像素分为M₁和M₂两类，其中M₁表示通过插值得到的像素值，M₂表示通过传感器直接获得的像素值，I_CFA(m,n)表示插值点(m,n)处的像素值，具体步骤如下：

第2.1步，对每一个图像块中插值点(m,n)处的像素值建立线性插值模型：

$I_{CFA}^k(m,n)=\underset{i,j={-N}_0}{\overset{N_0}{\Σ}}v(i,j)I_{CFA}^k(m+i,n+j)+r(m,n)---\left(2\right)$

其中，参数 $\stackrel{\‾}{v}=\left\{v(i,j)\right|-N_0\≤i,j\≤N_0\}(N_0\∈N^{+},v(0,0)=0),$ 参数r(m,n)是服从均值为0、方差为σ²正态分布的残余误差；

第2.2步，对参数进行初始化，令N₀＝1，即与其相邻的8个像素值相关，方差σ＝2，属于M₂的条件概率为P₀＝1/256，对每一个图像块利用EM算法估算出其插值系数，记为计算所有的平均值，记为

$m_{\stackrel{\→}{v}}=\frac{{64}^2}{M\×N}\underset{k=1}{\overset{M\×N/{64}^2}{\Σ}}{\stackrel{\→}{v}}^k---\left(3\right);$

第2.3步，利用构造最终的插值系数矩阵，记为H：

$H=\left\{\begin{array}{ccc}{m}_{\stackrel{\‾}{v}}\left(1\right)& m_{\stackrel{\‾}{v}}\left(\text{2}\right)& m_{\stackrel{\‾}{v}}\left(\text{3}\right)\\ m_{\stackrel{\‾}{v}}\left(\text{4}\right)& \text{0}& m_{\stackrel{\‾}{v}}\left(\text{5}\right)\\ m_{\stackrel{\‾}{v}}\left(\text{6}\right)& m_{\stackrel{\‾}{v}}\left(\text{7}\right)& m_{\stackrel{\‾}{v}}\left(\text{8}\right)\end{array}\right.\text{---}\left(4\right);$

第2.4步，记绿色分量I_CFA插值点(m,n)的邻域矩阵为

$I_{CFA}^{m,n}=\left[\begin{array}{ccc}{I}_{CFA}(m-1,n-1)& I_{CFA}(m-1,n)& I_{CFA}(m-1,n+1)\\ I_{CFA}(m,n-1)& I_{CFA}(m,n)& I_{CFA}(m,n-1)\\ I_{CFA}(m+1,n-1)& I_{CFA}(m+1,n)& I_{CFA}(m+1,n-1)\end{array}\right]---\left(5\right);$

第2.5步，利用最终的插值系数矩阵H和差值点(m,n)邻域矩阵得到原始图像模式I'_CFA内的像素值I'_CFA(m,n)：

${I^{\′}}_{CFA}(m,n)=\mathrm{\Σ}\mathrm{\Σ}H\·I_{CFA}^{m,n}---\left(6\right).$

4.如权利要求1所述的基于颜色滤波阵列特性的拼接图像篡改检测方法，其特征在于，所述第3步利用边缘检测算子进行篡改定位检测具体步骤如下：

第3.1步，定义新矩阵I_C，其元素为I_CFA与I'_CFA的对应元素差的平方：

$I_C=\left[\begin{array}{cccc}{(I_{CFA}^{11}-I_{CFA}^{\′11})}^2& {(I_{CFA}^{1\text{2}}-I_{CFA}^{\′1\text{2}})}^2& ...& {(I_{CFA}^{1\text{N}}-I_{CFA}^{\′1\text{N}})}^2\\ {(I_{CFA}^{21}-I_{CFA}^{\′21})}^2& {(I_{CFA}^{22}-I_{CFA}^{\′22})}^2& ...& {(I_{CFA}^{2N}-I_{CFA}^{\′2N})}^2\\ ...& ...& ...& ...\\ {(I_{CFA}^{M1}-I_{CFA}^{\′M1})}^2& {(I_{CFA}^{M2}-I_{CFA}^{\′M2})}^2& ...& {(I_{CFA}^{MN}-I_{CFA}^{\′MN})}^2\end{array}\right]\text{---}\left(7\right);$

第3.2步，对I_C进行二值化处理得到I'_C，然后利用Canny边缘检测算子对I'_C进行边缘检测，得到初步篡改定位结果I_L：

I_L＝E(I'_C,'canny')    (8)。

5.如权利要求4所述的基于颜色滤波阵列特性的拼接图像篡改检测方法，其特征在于，所述第3步还包括：

第3.3步，将初步篡改定位结果I_L使用形态学闭运算进行处理，得到最终的篡改定位结果I_Lend：

I_Lend＝imclose(I_L,SE)    (9)

其中，SE是结构元素。