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

摘要

本发明公开了一种DFT交叉相关分析的数码照片来源识别方法，包括以下步骤：(1)首先对参考和测试照片集进行R、G和B颜色通道分割；(2)对应照片原通道Ii和无PRNU Iirj，Iitk，相减得到残差图像nir，nit(3)对每一通道计算参考PRNU wci，i ∈{R，G，B}和融合参考PRNU wc；(4)计算融合参考PRNU wc和测试图像残差噪声nt的交叉相关系数cwn；(5)最后定位交叉相关系数最大值位置，并计算其到原点的欧氏距离；确定判定阈值。本发明使用的特征信息对7中型号不同品牌的数码相机拍摄的照片进行来源识别，平均识别率达到了96.3％。因此，本方法具有较高的应用价值。

主权项

一种DFT交叉相关分析的数码照片来源识别方法，其特征在于，包括以下步骤：(1)首先对参考和测试照片集进行R、G和B颜色通道分割：照片I可以分割为Ii，i∈{R，G，B}；选择一个BLS‑GSM小波低通滤波器对Ii，i ∈{R，G，B}进行滤波得到无PRNU的参考照片Iirj，j＝1，2...n和测试照片Iitk，k＝1，2...m；(2)对应照片原通道Ii和无PRNUIirj，Iitk，相减得到残差图像nir，nit；(3)对每一通道计算参考PRNU wci，i ∈{R，G，B}和融合参考PRNU wc；计算式如下： w ci = Σ k = 1 L n i k L ; wc＝0.3*wcr+0.6*wcg+0.1*wcb.(4)计算融合参考PRNU wc和测试图像残差噪声nt的交叉相关系数cwn；计算式如下： c wn ( j , k ) = Σ x , y w c ( x , y ) n t * ( x - j , y - k ) = Σ u , v W c ( u , v ) N t * ( u , v ) exp [ i 2 π ( uj M + vk N ) ] j = 1 KM , k = 1 KN 式中，N和M是照片的横向和纵向像素数，(*)表示取共轭，大写标示为对应小写表示的DFT变换；DFT变化式为： W c ( u , v ) = Σ x , y W c ( x , y ) MN exp [ - i 2 π ( ux M + uy N ) ] (5)最后定位交叉相关系数最大值位置，并计算其到原点的欧氏距离；确定判定阈值；算式如下：m(x0，y0)＝max||cm(j，k)||j＝1K M，k＝1KN ρ ( x 0 , y 0 ) = x 0 2 + y 0 2 > Th 式中，(x0，y0)为交叉相关系数最大值位置，Th为判定阈值，ρ(x0，y0)＞Th表示不是此相机拍摄的。