[发明专利]一种基于DFT交叉相关分析的数码照片来源识别方法

权利要求书

1.一种DFT交叉相关分析的数码照片来源识别方法，其特征在于，包括以下步骤：

(1)首先对参考和测试照片集进行R、G和B颜色通道分割：照片I可以分割为Iⁱ，i∈{R，G，B}；选择一个BLS-GSM小波低通滤波器对Iⁱ，i ∈{R，G，B}进行滤波得到无PRNU的参考照片Iⁱ_rj，j＝1，2...n和测试照片Iⁱ_tk，k＝1，2...m；

(2)对应照片原通道Iⁱ和无PRNUIⁱ_rj，Iⁱ_tk，相减得到残差图像nⁱ_r，nⁱ_t；

(3)对每一通道计算参考PRNU w_ci，i ∈{R，G，B}和融合参考PRNU w_c；计算式如下：

wci=Σk=1LnikL;]]>

w_c＝0.3*w_cr+0.6*w_cg+0.1*w_cb.

(4)计算融合参考PRNU w_c和测试图像残差噪声n_t的交叉相关系数c_wn；计算式如下：

cwn(j,k)=Σx,ywc(x,y)nt*(x-j,y-k)]]>

=Σu,vWc(u,v)Nt*(u,v)exp[i2π(ujM+vkN)]j=1KM,k=1KN]]>

式中，N和M是照片的横向和纵向像素数，(*)表示取共轭，大写标示为对应小写表示的DFT变换；DFT变化式为：

Wc(u,v)=Σx,yWc(x,y)MNexp[-i2π(uxM+uyN)]]]>

(5)最后定位交叉相关系数最大值位置，并计算其到原点的欧氏距离；确定判定阈值；算式如下：

m(x₀，y₀)＝max||c_m(j，k)||j＝1K M，k＝1KN

ρ(x0,y0)=x02+y02>Th]]>

式中，(x₀，y₀)为交叉相关系数最大值位置，Th为判定阈值，ρ(x₀，y₀)＞Th表示不是此相机拍摄的。

2.根据权利要求1所述的方法，其特征在于，选择判定阈值Th的方法为：由中心极限定理可知，不是同一数码相机照片交叉相关性数位置的欧氏距离值近似服从高斯分布，n个测试照片其交叉相关性数位置欧氏距离值的累加分布密度函数、高斯分布密度函数和高斯分布函数的关系为：

c(z)=1-nz∫-∞+∞f(xz)Fn-1(x)dx,Z>1.]]>

取Th＝z，使：P_FA＝1-c(Th)最小。