CN101452568A - 基于反卷积的篡改图像盲检测方法 - Google Patents
基于反卷积的篡改图像盲检测方法 Download PDFInfo
- Publication number
- CN101452568A CN101452568A CNA200810203998XA CN200810203998A CN101452568A CN 101452568 A CN101452568 A CN 101452568A CN A200810203998X A CNA200810203998X A CN A200810203998XA CN 200810203998 A CN200810203998 A CN 200810203998A CN 101452568 A CN101452568 A CN 101452568A
- Authority
- CN
- China
- Prior art keywords
- original
- tampered
- image
- circletimes
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
一种基于反卷积的篡改图像盲检测方法,属于图像信息安全领域。本发明利用反卷积和维纳滤波技术实现降晰篡改图像的盲检测,首先人为初步确定待检测图像的篡改区域位置,提取剩余未篡改区域的多个规则子图像,利用多个子图像构建一个近似满足完全卷积关系的数据块,采用交替迭代盲反卷积方法从该数据块中估计出未篡改区域的降晰函数,最后,利用估计的降晰函数对整幅图像做维纳滤波,并利用滤波结果辨别待检测图像是否经过篡改和具体的篡改区域。本方法简单、有效,在图像信息安全领域具有一定的应用前景。
Description
技术领域
本发明涉及一种基于反卷积的篡改图像盲检测方法。该方法是考虑到某些篡改图像的原始区域(未篡改区域)经过某种降晰,而篡改区域则为人为降晰的结果,原始区域和篡改区域有着不同的降晰函数。由此,可利用盲反卷积技术从原始区域的图像数据估计出原始降晰函数,基于此降晰函数,对整幅图像做维纳滤波,凸显出篡改区域。此方法在图像信息安全领域具有一定的应用前景。
背景技术
图像作为信息的一种极为重要的表达方式,广泛地应用在社会各个领域。但随着低成本、高性能及有着友好操作界面的图像处理软件的不断推出,数字化的图像信息面临着极大的安全隐患(即图像信息易被恶意篡改),由此,图像内容的真实性认证成为现代社会亟待解决的信息安全性问题。数字水印和数字签名曾作为一种有效手段被成功用于图像的安全认证,但由于其需要在未篡改的源图像中嵌入认证信息,使得其在某些情况下并不能得到有效应用。篡改图像盲检测是近年发展起来的一种新图像认证方法,由于其需要先验信息少,不需要人为在源图像中嵌入认证信息,因此适合更多场合的应用。
但是,篡改者篡改图像的方式多种多样,研究出适合所有篡改手段的盲检测方法几乎不可能。目前,篡改图像盲检测技术大体上可以分为两类:基于图像内容的盲检测技术和基于图像获取设备特性的盲检测技术。前者方法主要理论依据是考虑到篡改者在对图像进行篡改的同时,篡改手段会对图像的统计信息造成影响,如利用双谱分检测篡改过程中的双线性影响,利用期望最大化(EM)算法检测图像经历的重采样,利用篡改图像后的重压缩对图像统计数据造成的影响等等。基于图像内容盲检测的另一方面是检测篡改后图像的某种特性不一致性,如通过估计图像中不同部分的光照方向实现篡改图像的盲检测。基于图像获取设备的盲检测技术是利用图像获取设备会对图像数据的统计特性造成的变化或直接估计得到的设备特性参数实现篡改图像的盲检测,如通过估计数码相机插值部分引入的数据相关性,提取图像中遗留的数码相机模式噪声,计算数码相机的响应函数一致性,以及识别插值模式和插值算法及估计篡改滤波器来实现篡改盲检测等等。
发明内容
本发明的目的在于针对经过某种降晰的篡改图像,提出一种基于反卷积的篡改图像检测方法。利用盲反卷积技术,从未篡改区域估计出降晰函数,基于此估计的降晰函数,对整幅图像作维纳滤波,由于在篡改区域是人为降晰的结果,存在与未篡改区域不同的降晰函数,故可凸显出篡改部分。
为达到上述目的,本发明采用以下技术方案:
假设篡改图像为y,其中未篡改区域为yoriginal,篡改区域为ydoctored,且未篡改区域的降晰函数为horiginal,篡改区域的降晰函数为hdoctored,有下式关系成立:
式(1)(2)中xoriginal和xdoctored为未篡改区域和篡改区域的源图像数据。N、M为未篡改区域降晰函数的阶数,P、Q为篡改区域降晰函数的阶数。基于反卷积的篡改图像检测方法首先人为的估计图像中的可疑区域(即篡改区域)和原始区域(未篡改区域),利用未篡改区域的数据估计出降晰函数用此降晰函数对整幅图像y作维纳滤波,维纳滤波的结果凸显出的模糊区域即为篡改区域。
具体的篡改检测步骤如下:
1)人为辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal。
7)对第6步中维纳滤波结果逆傅立叶变换,得到维纳滤波空间域图像数据因为用估计的未篡改区域的降晰函数对整幅图像做维纳滤波,篡改区域具有不同的降晰函数,所以,若发现第1步中人为确定的可疑区域(篡改区域)存在模糊或者边界存在振铃效应,则此区域为篡改区域,否则,篡改的可能性不大。
8)若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。
现对第4步中用相邻子图像求均构造满足循环卷积关系方法作如下说明,首先,为方便解释,假设有二维矩阵c、b,阶数为C1×C2和B1×B2,二维卷积核a,阶数为A1×A2。三者满足线性卷积关系 且有C1=B1+A1-1、C2=B2+A2-1。把c划分为多个域,其中
Q01=c(0:A1-2,0:A2-2) (3)
Q02=c(C1-A1+2:C1,0:A2-2) (4)
Q03=c(C1-A1+2:C1,C2-A2+2:C2) (5)
Q04=c(0:A1-2:C1,C2-A2+2:C2) (6)
Q1=c(A1-1:C1-A1+1,0:A2-2) (7)
Q2=c(C1-A1+2:C1,A2-1:C2-A2+1) (8)
Q3=c(A1-1:C1-A1+1,C2-A2+2:C2) (9)
Q4=c(0:A1-2,A2-1:C2-A2+1) (10)
Q5=c(A1-1:C1-A1+1,A2-1:C2-A2+1) (11)
利用线性卷积构造循环卷积 区域划分如图1(a-2)所示,根据线性卷积和循环卷积的关系不难得到: 设从未篡改区域截取的子图像如图1(b)所示。对于每个子图像,根据卷积特性划分为四种不同区域 和显然,它们不仅不满足线性卷积关系,也不满足循环卷积关系,但它们和线性卷积有着很特殊的联系。以第一行子图像(1≤i≤K2)的为例分析,有下列关系存在:
直到
由上面分析 故
对式(16)求均
同理,对i列所有行相加求均有
式中误差error项相对项较小,而且随着K=K1×K2变大而减小。
上述第5步交替迭代估计未篡改区域的降晰函数作说明如下:构造迭代估计的代价函
式(22)中为循环卷积算子,右边二三两项为对原始图像和降晰函数的平滑约束项
具体迭代步骤如下:
②求解
③更新
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
②求解
③更新
④赋值J2′=J2,将步骤③中更新代入式(24)更新J2,求解
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
4)回到步骤2)步,循环30次结束;
本发明与现有技术相比较,具有如下显而易见的突出实质性特点和显著优点:本发明利用反卷积技术,从未篡改区域估计出降晰函数,基于此降晰函数,对整幅图像作维纳滤波,可凸显出篡改部分。本方法简单、有效,在图像信息安全领域具有一定的应用前景。
附图说明
图1为利用子图像构建满足循环卷积关系子块的示意图;(a)为线性卷积和循环卷积的关系示意图,(b)为子图像提取的示意图;
图2整个篡改图像检测程序示意图;
图3交替迭代估计降晰函数的程序示意图;
具体实施方式
本发明的一个优选实施例结合图详述如下:本发明旨在提供一种基于反卷积的图像盲检测方法,具体的检测流程如图2所示。该方法首先利用待检测图像,人为估计出篡改区域和未篡改区域,并从未篡改区域估计出未篡改区域经历的降晰函数,最终用该降晰函数对待检测图像作维纳滤波,实现篡改图像检测。具体的步骤为:
1.人为辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal。
②求解
③更新
④赋值J1′=J1,将步骤③中更新代入式(23)更新J1,求解
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对施加约束 即小于零值将其置零,大于255置为255。
①将 代入式(24)获得v(n)=-J2(n);
②求解
③更新
④赋值J2′=J2,将步骤③中更新代入式(24)更新J2,求解
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
4)回到步骤2)步,循环30次结束;
8.若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。
Claims (2)
1.一种基于反卷积的篡改图像盲检测方法,其特征在于首先人为的对图像中篡改区域进行估计,利用未篡改区域部分数据估计未篡改区域的降晰函数,然后利用估计的降晰函数对待检测图像进行维纳滤波,滤波图像中的模糊或周围有明显振铃效应的部分被认为是篡改区域;具体的操作步骤如下:
1)人为的辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal;
2)根据第一步得到的未篡改区域yoriginal大小,确定待提取子图像阶数为p、q,要求且待提取子图像个数为K个,N、M为未篡改区域降晰函数的阶数;3)在未篡改区域yoriginal提取规则子图像(1≤i≤K1,1≤j≤K2),要求规则子图像依次相邻,排列成K1行K2列的方格,且K=K1×K2;
4)将K个规则子图像数据相加求均,构造新的未篡改区域数据且 构造新数据和horiginal近似满足循环卷积关系,即 其中ω为频域表示,Xoriginal(ω)为未篡改区域的原始未降晰数据,Horiginal(ω)为未篡改区域降晰函数频域变换,IDFT(·)为逆傅立叶变换操作;
8)若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。
1)初时参数设置:设置参数代价函数中加权因子λ、γ,要求大致满足 为高通算子; 其中为各像素点周围5×5邻域的局部方差,而 w2初值设为全1,初始值设为 初始值为随机取值,同时补零扩展到大小,其中n为二维坐标(i,j)。
①将 代入下式
获得u(n)=-J1(n);
③更新
④赋值J1′=J1,将步骤③中更新代所列入J1(n)的求解式,更新J1,求解
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
①将 代入式
获得v(n)=-J2(n);
②求解
③更新
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
4)回到步骤2)步,循环30次结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA200810203998XA CN101452568A (zh) | 2008-12-04 | 2008-12-04 | 基于反卷积的篡改图像盲检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA200810203998XA CN101452568A (zh) | 2008-12-04 | 2008-12-04 | 基于反卷积的篡改图像盲检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101452568A true CN101452568A (zh) | 2009-06-10 |
Family
ID=40734783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA200810203998XA Pending CN101452568A (zh) | 2008-12-04 | 2008-12-04 | 基于反卷积的篡改图像盲检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101452568A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373115A (zh) * | 2016-08-22 | 2017-02-01 | 天津大学 | 一种平滑滤波操作类型及参数识别方法 |
CN106709915A (zh) * | 2017-01-06 | 2017-05-24 | 天津大学 | 一种图像重采样操作检测方法 |
CN106845540A (zh) * | 2017-01-06 | 2017-06-13 | 天津大学 | 一种图像重采样操作插值类型识别方法 |
CN108920700A (zh) * | 2018-07-17 | 2018-11-30 | 中国联合网络通信集团有限公司 | 一种虚假图片识别方法及装置 |
CN112861786A (zh) * | 2021-03-09 | 2021-05-28 | 南开大学 | 人脸篡改视频检测方法和人脸篡改视频检测装置 |
CN113674195A (zh) * | 2020-05-13 | 2021-11-19 | 中国移动通信集团有限公司 | 图像检测方法、装置、设备及存储介质 |
-
2008
- 2008-12-04 CN CNA200810203998XA patent/CN101452568A/zh active Pending
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373115A (zh) * | 2016-08-22 | 2017-02-01 | 天津大学 | 一种平滑滤波操作类型及参数识别方法 |
CN106373115B (zh) * | 2016-08-22 | 2019-02-05 | 天津大学 | 一种平滑滤波操作类型及参数识别方法 |
CN106709915A (zh) * | 2017-01-06 | 2017-05-24 | 天津大学 | 一种图像重采样操作检测方法 |
CN106845540A (zh) * | 2017-01-06 | 2017-06-13 | 天津大学 | 一种图像重采样操作插值类型识别方法 |
CN106845540B (zh) * | 2017-01-06 | 2020-02-07 | 天津大学 | 一种图像重采样操作插值类型识别方法 |
CN106709915B (zh) * | 2017-01-06 | 2020-02-07 | 天津大学 | 一种图像重采样操作检测方法 |
CN108920700A (zh) * | 2018-07-17 | 2018-11-30 | 中国联合网络通信集团有限公司 | 一种虚假图片识别方法及装置 |
CN108920700B (zh) * | 2018-07-17 | 2021-04-27 | 中国联合网络通信集团有限公司 | 一种虚假图片识别方法及装置 |
CN113674195A (zh) * | 2020-05-13 | 2021-11-19 | 中国移动通信集团有限公司 | 图像检测方法、装置、设备及存储介质 |
CN112861786A (zh) * | 2021-03-09 | 2021-05-28 | 南开大学 | 人脸篡改视频检测方法和人脸篡改视频检测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zeng et al. | Image splicing localization using PCA-based noise level estimation | |
Tang et al. | Median filtering detection of small-size image based on CNN | |
Huang et al. | Improved DCT-based detection of copy-move forgery in images | |
Singh et al. | Video frame and region duplication forgery detection based on correlation coefficient and coefficient of variation | |
US6266452B1 (en) | Image registration method | |
CN104933721B (zh) | 基于颜色滤波阵列特性的拼接图像篡改检测方法 | |
Kumar et al. | A fast DCT based method for copy move forgery detection | |
CN101452568A (zh) | 基于反卷积的篡改图像盲检测方法 | |
CN106097379B (zh) | 一种使用自适应阈值的图像篡改检测与定位方法 | |
Israni et al. | Edge detection of license plate using Sobel operator | |
Benrhouma et al. | Tamper detection and self-recovery scheme by DWT watermarking | |
CN106530200A (zh) | 一种基于深度学习模型的隐写图像检测方法及系统 | |
CN103747271B (zh) | 一种基于混合感知哈希的视频篡改检测方法和装置 | |
CN102542570B (zh) | 一种微波图像中人体隐藏危险物体自动检测方法 | |
CN108682007A (zh) | 基于深度随机森林的jpeg图像重采样自动检测方法 | |
CN106897694A (zh) | 用于国土资源监控的违建场景识别方法 | |
CN104636764B (zh) | 一种图像隐写分析方法以及其装置 | |
CN104036485A (zh) | 一种关于图像重采样篡改检测的方法 | |
CN109064376A (zh) | 基于Harris-Laplace与SIFT描述子的DFT域图像抗屏摄水印算法 | |
CN106056523A (zh) | 数字图像拼接篡改盲检测方法 | |
Yang et al. | Image tamper detection based on noise estimation and lacunarity texture | |
Zhu et al. | Remote sensing image watermarking based on motion blur degeneration and restoration model | |
CN102024244B (zh) | 基于图像特征区域的水印嵌入、检测方法及其装置 | |
CN103079066B (zh) | 脊波变换与dct变换相结合的图像水印嵌入和提取方法 | |
Zhang et al. | Image splicing localization using noise distribution characteristic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20090610 |