CN101452568A - 基于反卷积的篡改图像盲检测方法 - Google Patents

基于反卷积的篡改图像盲检测方法 Download PDF

Info

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
Application number
CNA200810203998XA
Other languages
English (en)
Inventor
王睿
方勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CNA200810203998XA priority Critical patent/CN101452568A/zh
Publication of CN101452568A publication Critical patent/CN101452568A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

一种基于反卷积的篡改图像盲检测方法,属于图像信息安全领域。本发明利用反卷积和维纳滤波技术实现降晰篡改图像的盲检测,首先人为初步确定待检测图像的篡改区域位置,提取剩余未篡改区域的多个规则子图像,利用多个子图像构建一个近似满足完全卷积关系的数据块,采用交替迭代盲反卷积方法从该数据块中估计出未篡改区域的降晰函数,最后,利用估计的降晰函数对整幅图像做维纳滤波,并利用滤波结果辨别待检测图像是否经过篡改和具体的篡改区域。本方法简单、有效,在图像信息安全领域具有一定的应用前景。

Description

基于反卷积的篡改图像盲检测方法
技术领域
本发明涉及一种基于反卷积的篡改图像盲检测方法。该方法是考虑到某些篡改图像的原始区域(未篡改区域)经过某种降晰,而篡改区域则为人为降晰的结果,原始区域和篡改区域有着不同的降晰函数。由此,可利用盲反卷积技术从原始区域的图像数据估计出原始降晰函数,基于此降晰函数,对整幅图像做维纳滤波,凸显出篡改区域。此方法在图像信息安全领域具有一定的应用前景。
背景技术
图像作为信息的一种极为重要的表达方式,广泛地应用在社会各个领域。但随着低成本、高性能及有着友好操作界面的图像处理软件的不断推出,数字化的图像信息面临着极大的安全隐患(即图像信息易被恶意篡改),由此,图像内容的真实性认证成为现代社会亟待解决的信息安全性问题。数字水印和数字签名曾作为一种有效手段被成功用于图像的安全认证,但由于其需要在未篡改的源图像中嵌入认证信息,使得其在某些情况下并不能得到有效应用。篡改图像盲检测是近年发展起来的一种新图像认证方法,由于其需要先验信息少,不需要人为在源图像中嵌入认证信息,因此适合更多场合的应用。
但是,篡改者篡改图像的方式多种多样,研究出适合所有篡改手段的盲检测方法几乎不可能。目前,篡改图像盲检测技术大体上可以分为两类:基于图像内容的盲检测技术和基于图像获取设备特性的盲检测技术。前者方法主要理论依据是考虑到篡改者在对图像进行篡改的同时,篡改手段会对图像的统计信息造成影响,如利用双谱分检测篡改过程中的双线性影响,利用期望最大化(EM)算法检测图像经历的重采样,利用篡改图像后的重压缩对图像统计数据造成的影响等等。基于图像内容盲检测的另一方面是检测篡改后图像的某种特性不一致性,如通过估计图像中不同部分的光照方向实现篡改图像的盲检测。基于图像获取设备的盲检测技术是利用图像获取设备会对图像数据的统计特性造成的变化或直接估计得到的设备特性参数实现篡改图像的盲检测,如通过估计数码相机插值部分引入的数据相关性,提取图像中遗留的数码相机模式噪声,计算数码相机的响应函数一致性,以及识别插值模式和插值算法及估计篡改滤波器来实现篡改盲检测等等。
发明内容
本发明的目的在于针对经过某种降晰的篡改图像,提出一种基于反卷积的篡改图像检测方法。利用盲反卷积技术,从未篡改区域估计出降晰函数,基于此估计的降晰函数,对整幅图像作维纳滤波,由于在篡改区域是人为降晰的结果,存在与未篡改区域不同的降晰函数,故可凸显出篡改部分。
为达到上述目的,本发明采用以下技术方案:
假设篡改图像为y,其中未篡改区域为yoriginal,篡改区域为ydoctored,且未篡改区域的降晰函数为horiginal,篡改区域的降晰函数为hdoctored,有下式关系成立:
y original ( i , j ) = Σ k = 0 N - 1 Σ l = 0 M - 1 x original ( k , l ) · h original ( i - k , j - l ) - - - ( 1 )
y doctored ( i , j ) = Σ k = 0 P - 1 Σ l = 0 Q - 1 x doctored ( k , l ) · h doctored ( i - k , j - l ) - - - ( 2 )
式(1)(2)中xoriginal和xdoctored为未篡改区域和篡改区域的源图像数据。N、M为未篡改区域降晰函数的阶数,P、Q为篡改区域降晰函数的阶数。基于反卷积的篡改图像检测方法首先人为的估计图像中的可疑区域(即篡改区域)和原始区域(未篡改区域),利用未篡改区域的数据估计出降晰函数
Figure A200810203998D00063
用此降晰函数对整幅图像y作维纳滤波,维纳滤波的结果凸显出的模糊区域即为篡改区域。
具体的篡改检测步骤如下:
1)人为辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal
2)根据第一步得到的未篡改区域yoriginal大小,确定待提取子图像阶数为p、q,要求
Figure A200810203998D00065
待提取子图像个数为K个。
3)在未篡改区域yoriginal提取规则子图像
Figure A200810203998D00066
(1≤i≤K1,1≤j≤K2),要求规则子图像依次相邻,排列成K1行K2列的方格,且K=K1×K2
4)将K个规则子图像数据相加求均,构造新的未篡改区域数据
Figure A200810203998D00067
y ~ original ( k , l ) = Σ i = 1 K 1 Σ j = 1 K 2 y original ( i , j ) ( k , l ) , 构造新数据
Figure A200810203998D00069
和horiginal近似满足循环卷积关系,即 y ~ original = IDFT ( X original ( ω ) H original ( ω ) ) .
5)利用第4步中数据
Figure A200810203998D00071
采用迭代交替盲反卷积技术估计未篡改区域的降晰函数
6)利用第5步中估计的未篡改区域的降晰函数
Figure A200810203998D00073
对整幅图像做维纳滤波
Figure A200810203998D00074
γ为噪信比,可根据信噪比大小预设成固定值。
7)对第6步中维纳滤波结果
Figure A200810203998D00075
逆傅立叶变换,得到维纳滤波空间域图像数据
Figure A200810203998D00076
因为用估计的未篡改区域的降晰函数对整幅图像做维纳滤波,篡改区域具有不同的降晰函数,所以,若发现第1步中人为确定的可疑区域(篡改区域)存在模糊或者边界存在振铃效应,则此区域为篡改区域,否则,篡改的可能性不大。
8)若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。
现对第4步中用相邻子图像求均构造满足循环卷积关系方法作如下说明,首先,为方便解释,假设有二维矩阵c、b,阶数为C1×C2和B1×B2,二维卷积核a,阶数为A1×A2。三者满足线性卷积关系 c ( i , j ) = Σ k = 0 A 2 Σ l = 0 A 1 a ( k , l ) · b ( i - k , j - l ) , 且有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)
利用线性卷积构造循环卷积 c ~ = b ⊗ a ,
Figure A200810203998D00079
区域划分如图1(a-2)所示,根据线性卷积和循环卷积的关系不难得到: Q ~ 00 = Q 01 + Q 02 + Q 03 + Q 04 , Q ~ 1 = Q 1 + Q 3 , Q ~ 4 = Q 2 + Q 4 , Q ~ 5 = Q 5 . 设从未篡改区域截取的子图像如图1(b)所示。对于每个子图像,根据卷积特性划分为四种不同区域
Figure A200810203998D00085
Figure A200810203998D00086
Figure A200810203998D00087
Figure A200810203998D00088
显然,它们不仅不满足线性卷积关系,也不满足循环卷积关系,但它们和线性卷积有着很特殊的联系。以第一行子图像(1≤i≤K2)的
Figure A200810203998D000810
为例分析,有下列关系存在:
Figure A200810203998D000812
直到
式(12)假设
Figure A200810203998D000814
为第一行未截取子图像产生。将
Figure A200810203998D000815
(1≤i≤K2)相加得到
Figure A200810203998D000816
+ · · · · · · + ( Q 1 ( 1 , K 2 - 1 ) + Q 3 ( 1 , K 2 - 1 ) ) + Q 1 ( 1 , K 2 )
由上面分析 Q 1 ( 1 , i ) + Q 3 ( 1 , i ) = Q ~ 1 ( 1 , i ) ,
Figure A200810203998D000819
对式(16)求均
Figure A200810203998D000820
平均后不满足循环卷积
Figure A200810203998D000821
项较小,可看作误差项
Figure A200810203998D000822
同理,对i列所有行相加求均有
Figure A200810203998D000823
不难分析,对所有子图像
Figure A200810203998D000824
(1≤i≤K1,1≤j≤K2)相加求均有
对于
Figure A200810203998D00092
(1≤i≤K1,1≤j≤K2)来说,与循环卷积
Figure A200810203998D00093
是相等。结合式(18-20),得出结论
y ~ original ( k , l ) ≈ Σ i = 1 K 1 Σ j = 1 K 2 y original ( i , j ) ( k , l ) + error - - - ( 21 )
式中误差error项相对项较小,而且随着K=K1×K2变大而减小。
上述第5步交替迭代估计未篡改区域的降晰函数作说明如下:构造迭代估计的代价函
J ( x ^ original , h ^ original ) = 1 2 | | y ~ original - x ^ original ⊗ h ^ original | | 2
+ 1 2 λ | | c ⊗ x ^ original | | w 1 2 + 1 2 γ | | d ⊗ h ^ original | | w 2 2          (22)
式(22)中为循环卷积算子,右边二三两项为对原始图像和降晰函数的平滑约束项
| | c ⊗ x ^ original | | w 1 2 = Σ i , j = 0 p , q w 1 ( i , j ) · ( c ( i , j ) ⊗ x ^ ( i , j ) ) | | d ⊗ h ^ original | | w 1 2 = Σ i , j = 0 p , q w 2 ( i , j ) · ( d ( i , j ) ⊗ h ^ ( i , j ) )
其中 c = 0 0.25 0 0.25 - 1 0.25 0 0.25 0 , d = 1 0.5 0.5 0 为高通算子。 w 1 = 1 1 + μLoc _ var ( y ~ original ) , 其中
Figure A200810203998D000913
Figure A200810203998D000914
各像素点周围5×5邻域的局部方差,而 μ = 1000 max i , j ( Loc _ var ( y ~ original ) ) . w2初值设为全1,在迭代中更新。设n=(i,j)为二维坐标,求式(22)的梯度为
J 1 ( n ) = ∂ J ∂ x ^ original = h ^ original ( n ) ⊗ x ^ original ( n ) ⊗ h ^ original ( - n ) - y ~ original ( n ) ⊗ h ^ original ( - n ) - - - ( 23 )
+ λ · w 1 ( n ) · [ c ( n ) ⊗ x ^ original ( n ) ⊗ c ( - n ) ]
J 2 ( n ) = ∂ J ∂ h ^ original = h ^ original ( n ) ⊗ x ^ original ( n ) ⊗ x ^ original ( - n ) - y ~ original ( n ) ⊗ x ^ original ( - n ) - - - ( 24 )
+ γ · w 2 ( n ) · [ d ( n ) ⊗ h ^ original ( n ) ⊗ d ( - n ) ]
具体迭代步骤如下:
1)初时参数设置:代价函数的加权因子λ、γ的设置对算法影响不大,可大致满足 γ λ = Σ i , j p , q x ^ original ( i , j ) · max ( x ^ original ( i , j ) ) 即可。其它参数c、d、w1和w2如上所述,
Figure A200810203998D00102
初始值设为 x ^ original 0 ( n ) = y ~ original ( n ) ,
Figure A200810203998D00104
初始值为随机取值,同时补零扩展到
Figure A200810203998D00105
大小。
2)固定迭代更新
Figure A200810203998D00107
①将
Figure A200810203998D00108
Figure A200810203998D00109
代入式(23)获得u(n)=-J1(n);
②求解 α = | | J 1 | | 2 | | h ^ original ⊗ u | | 2 + λ | | c ⊗ u | | w 1 2 ;
③更新 x ^ original = x ^ original + α · u ;
④赋值J1′=J1,将步骤③中更新
Figure A200810203998D001012
代入式(23)更新J1,求解 β = | | J 1 | | 2 | | J 1 ′ | | 2 ;
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对
Figure A200810203998D001014
施加约束 0 ≤ x ^ original ( n ) ≤ 255 , 即小于零值将其置零,大于255置为255。
3)固定
Figure A200810203998D001016
迭代更新
Figure A200810203998D001017
①将
Figure A200810203998D001018
代入式(24)获得v(n)=-J2(n);
②求解 α ′ = | | J 2 | | 2 | | x ^ original ⊗ v | | 2 + γ | | d ⊗ v | | w 2 2 ;
③更新 h ^ original = h ^ original + α ′ · v ;
④赋值J2′=J2,将步骤③中更新代入式(24)更新J2,求解 β ′ = | | J 2 | | 2 | | J 2 ′ | | 2 ;
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对
Figure A200810203998D001024
施加约束 0 ≤ h ^ original ( n ) ≤ 1 , 即小于零值将其置零,同时进行归一化。更新 w 2 = 1 1 + ηLoc _ var ( h ^ original ) , 其中
Figure A200810203998D001027
各像素点周围3×3邻域的局部方差,而 η = 1000 max i , j ( Loc _ var ( h ^ original ) ) .
4)回到步骤2)步,循环30次结束;
本发明与现有技术相比较,具有如下显而易见的突出实质性特点和显著优点:本发明利用反卷积技术,从未篡改区域估计出降晰函数,基于此降晰函数,对整幅图像作维纳滤波,可凸显出篡改部分。本方法简单、有效,在图像信息安全领域具有一定的应用前景。
附图说明
图1为利用子图像构建满足循环卷积关系子块的示意图;(a)为线性卷积和循环卷积的关系示意图,(b)为子图像提取的示意图;
图2整个篡改图像检测程序示意图;
图3交替迭代估计降晰函数的程序示意图;
具体实施方式
本发明的一个优选实施例结合图详述如下:本发明旨在提供一种基于反卷积的图像盲检测方法,具体的检测流程如图2所示。该方法首先利用待检测图像,人为估计出篡改区域和未篡改区域,并从未篡改区域估计出未篡改区域经历的降晰函数,最终用该降晰函数对待检测图像作维纳滤波,实现篡改图像检测。具体的步骤为:
1.人为辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal
2.根据第一步得到的未篡改区域yoriginal大小,确定待提取子图像阶数为p、q,要求
Figure A200810203998D00112
Figure A200810203998D00113
待提取子图像个数为K个。
3.在未篡改区域yoriginal提取规则子图像
Figure A200810203998D00114
(1≤i≤K1,1≤j≤K2),要求规则子图像依次相邻,排列成K1行K2列的方格,且K=K1×K2
4.将K个规则子图像数据相加求均,构造新的未篡改区域数据
Figure A200810203998D00115
y ~ original ( k , l ) = Σ i = 1 K 1 Σ j = 1 K 2 y original ( i , j ) ( k , l ) , 构造新数据
Figure A200810203998D00117
Figure A200810203998D00118
近似满足循环卷积关系,即 y ~ original = IDFT ( X original ( ω ) H original ( ω ) ) .
5.利用第4步中数据
Figure A200810203998D00121
采用交替迭代盲反卷积技术估计未篡改区域的降晰函数
Figure A200810203998D00122
参见图3,估计
Figure A200810203998D00123
的具体步骤如下:
1)初时参数设置:设置参数λ、γ,要求大致满足 γ λ = Σ i , j p , q x ^ original ( i , j ) · max ( x ^ original ( i , j ) ) . c = 0 0.25 0 0.25 - 1 0.25 0 0.25 0 , d = 1 0.5 0.5 0 为高通算子。 w 1 = 1 1 + μLoc _ var ( y ~ original ) , 其中
Figure A200810203998D00128
各像素点周围5×5邻域的局部方差,而 μ = 1000 max i , j ( Loc _ var ( y ~ original ) ) . w2初值设为全1,
Figure A200810203998D001211
初始值设为 x ^ original 0 ( n ) = y ~ original ( n ) ,
Figure A200810203998D001213
初始值为随机取值,同时补零扩展到
Figure A200810203998D001214
大小。
2)固定迭代更新
Figure A200810203998D001216
①将
Figure A200810203998D001217
Figure A200810203998D001218
代入式(23)获得u(n)=-J1(n);
②求解 α = | | J 1 | | 2 | | h ^ original ⊗ u | | 2 + λ | | c ⊗ u | | w 1 2 ;
③更新 x ^ original = x ^ original + α · u ;
④赋值J1′=J1,将步骤③中更新代入式(23)更新J1,求解 β = | | J 1 | | 2 | | J 1 ′ | | 2 ;
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对施加约束 0 ≤ x ^ original ( n ) ≤ 255 , 即小于零值将其置零,大于255置为255。
3)固定
Figure A200810203998D001225
迭代更新
Figure A200810203998D001226
①将 代入式(24)获得v(n)=-J2(n);
②求解 α ′ = | | J 2 | | 2 | | x ^ original ⊗ v | | 2 + γ | | d ⊗ v | | w 2 2 ;
③更新 h ^ original = h ^ original + α ′ · v ;
④赋值J2′=J2,将步骤③中更新代入式(24)更新J2,求解 β ′ = | | J 2 | | 2 | | J 2 ′ | | 2 ;
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对
Figure A200810203998D00133
施加约束 0 ≤ h ^ original ( n ) ≤ 1 , 即小于零值将其置零,同时进行归一化。 w 2 = 1 1 + ηLoc _ var ( h ^ original ) , 其中
Figure A200810203998D00136
Figure A200810203998D00137
各像素点周围3×3邻域的局部方差,而 η = 1000 max i , j ( Loc _ var ( h ^ original ) ) .
4)回到步骤2)步,循环30次结束;
6.利用第5步中估计的未篡改区域的降晰函数
Figure A200810203998D00139
对整幅图像做维纳滤波γ为噪信比,可根据信噪比大小预设成固定值。
7.对第6步中维纳滤波结果
Figure A200810203998D001311
逆傅立叶变换,得到维纳滤波空间域图像数据
Figure A200810203998D001312
若发现第1步中人为确定的可疑区域(篡改区域)存在模糊或者边界存在振铃效应,则此区域为篡改区域,否则,篡改的可能性不大。
8.若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。

Claims (2)

1.一种基于反卷积的篡改图像盲检测方法,其特征在于首先人为的对图像中篡改区域进行估计,利用未篡改区域部分数据估计未篡改区域的降晰函数,然后利用估计的降晰函数对待检测图像进行维纳滤波,滤波图像中的模糊或周围有明显振铃效应的部分被认为是篡改区域;具体的操作步骤如下:
1)人为的辨别图像存在的可疑区域,即篡改区域ydoctored,一般篡改区域是为掩盖原始图像的部分数据,而达到篡改某种事实的目的,故篡改区域一般具有完整的内容特性,相应余下区域为未篡改区域yoriginal
2)根据第一步得到的未篡改区域yoriginal大小,确定待提取子图像阶数为p、q,要求
Figure A200810203998C00021
且待提取子图像个数为K个,N、M为未篡改区域降晰函数的阶数;3)在未篡改区域yoriginal提取规则子图像
Figure A200810203998C00022
(1≤i≤K1,1≤j≤K2),要求规则子图像依次相邻,排列成K1行K2列的方格,且K=K1×K2
4)将K个规则子图像数据相加求均,构造新的未篡改区域数据
Figure A200810203998C00023
y ~ original ( k , l ) = Σ i = 1 K 1 Σ j = 1 K 2 y original ( i , j ) ( k , l ) , 构造新数据
Figure A200810203998C00025
和horiginal近似满足循环卷积关系,即 y ~ original = IDFT ( X original ( ω ) H original ( ω ) ) , 其中ω为频域表示,Xoriginal(ω)为未篡改区域的原始未降晰数据,Horiginal(ω)为未篡改区域降晰函数频域变换,IDFT(·)为逆傅立叶变换操作;
5)利用第4步中数据
Figure A200810203998C00027
采用交替迭代盲反卷积技术估计未篡改区域的降晰函数
Figure A200810203998C00028
6)利用第5步中估计的未篡改区域的降晰函数
Figure A200810203998C00029
对整幅图像做维纳滤波
Figure A200810203998C000210
其中,σ为噪信比,可根据信噪比大小预设成固定值,
Figure A200810203998C000211
为估计的为篡改区域降晰函数的傅立叶变换,
Figure A200810203998C000212
为其共轭,Y(ω)为整幅待检测图像;
7)对第6步中维纳滤波结果
Figure A200810203998C00031
逆傅立叶变换,得到维纳滤波空间域图像数据
Figure A200810203998C00032
若发现第1步中人为确定的可疑区域存在模糊或者边界存在振铃效应,则此区域为篡改区域,否则,篡改的可能性不大;
8)若未发现篡改区域,可再次确定其它可疑区域,继续第1至第7步过程,进行篡改检测。
2.根据权利要求1所述的基于反卷积的篡改图像盲检测方法,其特征在于所述第5)步骤中的利用未篡改区域数据交替迭代估计的降晰函数
Figure A200810203998C00033
步骤如下:
1)初时参数设置:设置参数代价函数中加权因子λ、γ,要求大致满足 γ λ = Σ i , j p , q x ^ original ( i , j ) · max ( x ^ original ( i , j ) ) ; c = 0 0.25 0 0.25 - 1 0.25 0 0.25 0 , d = 1 0.5 0.5 0 为高通算子; w 1 = 1 1 + μLoc _ var ( y ~ original ) , 其中
Figure A200810203998C00038
各像素点周围5×5邻域的局部方差,而 μ = 1000 max i , j ( Loc _ var ( y ~ original ) ) ; w2初值设为全1,
Figure A200810203998C000311
初始值设为 x ^ original 0 ( n ) = y ~ original ( n ) ,
Figure A200810203998C000313
初始值为随机取值,同时补零扩展到
Figure A200810203998C000314
大小,其中n为二维坐标(i,j)。
2)固定
Figure A200810203998C000315
迭代更新
Figure A200810203998C000316
①将 h ^ original ( n ) , x ^ original ( n ) 代入下式
J 1 ( n ) = ∂ J ∂ x ^ original = h ^ original ( n ) ⊗ x ^ original ( n ) ⊗ h ^ original ( - n )
- y ~ original ( n ) ⊗ h ^ original ( - n ) + λ · w 1 ( n ) · [ c ( n ) ⊗ x ^ original ( n ) ⊗ c ( - n ) ]
获得u(n)=-J1(n);
②求解 α = | | J 1 | | 2 | | h ^ original ⊗ u | | 2 + λ | | c ⊗ u | | w 1 2 , 式中
Figure A200810203998C000322
为卷积算子;
③更新 x ^ original = x ^ original + α · u ;
④赋值J1′=J1,将步骤③中更新代所列入J1(n)的求解式,更新J1,求解
β = | | J 1 | | 2 | | J 1 ′ | | 2 ;
⑤更新u(n)=-J1(n)+β·u(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对
Figure A200810203998C00042
施加约束 0 ≤ x ^ original ( n ) ≤ 255 , 即小于零值将其置零,大于255置为255。
3)固定
Figure A200810203998C00044
迭代更新
Figure A200810203998C00045
①将 h ^ original ( n ) , x ^ original ( n ) 代入式
J 2 ( n ) = ∂ J ∂ h ^ original = h ^ original ( n ) ⊗ x ^ original ( n ) ⊗ x ^ original ( - n )
- y ~ original ( n ) ⊗ x ^ original ( - n ) + λ · w 2 ( n ) · [ d ( n ) ⊗ h ^ original ( n ) ⊗ d ( - n ) ]
获得v(n)=-J2(n);
②求解 α ′ = | | J 2 | | 2 | | x ^ original ⊗ v | | 2 + γ | | d ⊗ v | | w 2 2 ;
③更新 h ^ original = h ^ original + α ′ · v ;
④赋值J2′=J2,将步骤③中更新
Figure A200810203998C000412
代入所列J2(n)的求解式,更新J2,求解
β ′ = | | J 2 | | 2 | | J 2 ′ | | 2 ;
⑤更新v(n)=-J2(n)+β′·v(n);
⑥回到步骤②步,继续上述过程,循环10次;
更新结束对
Figure A200810203998C000414
施加约束 0 ≤ h ^ original ( n ) ≤ 1 , 即小于零值将其置零,同时进行归一化;更新 w 2 = 1 1 + ηLoc _ var ( h ~ original ) , 其中
Figure A200810203998C000417
Figure A200810203998C000418
各像素点周围3×3邻域的局部方差,而 η = 1000 max i , j ( Loc _ var ( h ~ original ) ) ;
4)回到步骤2)步,循环30次结束。
CNA200810203998XA 2008-12-04 2008-12-04 基于反卷积的篡改图像盲检测方法 Pending CN101452568A (zh)

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)

* Cited by examiner, † Cited by third party
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 中国移动通信集团有限公司 图像检测方法、装置、设备及存储介质

Cited By (10)

* Cited by examiner, † Cited by third party
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