技術(shù)特征：

1.一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于包括以下步驟：

①獲取N幅原始數(shù)字照片圖像，將獲取的第n幅原始數(shù)字照片圖像記為I_org,n；并任意選取一幅數(shù)字照片圖像作為待識別照片圖像，記為I_test；其中，N≥2，1≤n≤N，獲取每幅原始數(shù)字照片圖像所采用的相機(jī)與I_test所聲明的來源相機(jī)為同一部相機(jī)，每幅原始數(shù)字照片圖像和I_test的尺寸大小相同，且寬度為col、高度為row；

②獲取每幅原始數(shù)字照片圖像的降噪圖像，將I_org,n的降噪圖像記為I'_org,n；然后根據(jù)每幅原始數(shù)字照片圖像及其降噪圖像，獲取每幅原始數(shù)字照片圖像的噪聲殘差圖像，將I_org,n的噪聲殘差圖像記為W_org,n，將W_org,n中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值記為W_org,n(i,j)，W_org,n(i,j)＝I_org,n(i,j)-I'_org,n(i,j)；接著根據(jù)所有原始數(shù)字照片圖像的噪聲殘差圖像，獲取相機(jī)的傳感器模式噪聲，記為R；再消除R中的CFA插值噪聲，得到相機(jī)的無CFA插值噪聲的傳感器模式噪聲，記為R'；其中，1≤i≤col,1≤j≤row，I_org,n(i,j)表示I_org,n中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，I'_org,n(i,j)表示I'_org,n中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值；

同樣，獲取I_test的降噪圖像，記為I'_test；然后根據(jù)I_test及I'_test，獲取I_test的噪聲殘差圖像，記為W_test，將W_test中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值記為W_test(i,j)，W_test(i,j)＝I_test(i,j)-I'_test(i,j)；接著將W_test作為I_test的傳感器模式噪聲；再消除W_test中的CFA插值噪聲，得到I_test的無CFA插值噪聲的傳感器模式噪聲，記為W'_test；其中，I_test(i,j)表示I_test中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，I'_test(i,j)表示I'_test中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值；

③隨機(jī)生成一幅寬度為col且高度為row的標(biāo)準(zhǔn)高斯白噪聲圖像，記為G；然后根據(jù)G在空域?qū)'進(jìn)行空域平滑濾波處理，獲得相機(jī)的最終的傳感器模式噪聲，記為R^*；同樣，根據(jù)G在空域?qū)'_test進(jìn)行空域平滑濾除處理，獲得I_test的最終的傳感器模式噪聲，記為W^*_test；

④采用SPCE方法計(jì)算W^*_test與R^*的相關(guān)性，記為SPCE(W^*_test,R^*)；

⑤用事先設(shè)置的I_test所聲明的來源相機(jī)下的相關(guān)性判別閾值對SPCE(W^*_test,R^*)進(jìn)行判決，若SPCE(W^*_test,R^*)大于或等于相關(guān)性判別閾值，則判定I_test來自其所聲明的來源相機(jī)；若SPCE(W^*_test,R^*)小于相關(guān)性判別閾值，則判定I_test來自非其所聲明的來源相機(jī)，即I_test來自其它相機(jī)。

2.根據(jù)權(quán)利要求1所述的一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于所述的步驟②中獲取每幅原始數(shù)字照片圖像的降噪圖像和獲取I_test的降噪圖像所采用的方法為小波降噪處理方法；所述的步驟②中采用最大似然估計(jì)法對所有原始數(shù)字照片圖像的噪聲殘差圖像進(jìn)行處理，獲得R；所述的步驟②中采用零均值化法對R進(jìn)行處理以消除R中的CFA插值噪聲，獲得R'；所述的步驟②中采用零均值化法對W_test進(jìn)行處理以消除W_test中的CFA插值噪聲，獲得W'_test。

3.根據(jù)權(quán)利要求1或2所述的一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于所述的步驟③中根據(jù)G在空域?qū)'進(jìn)行空域平滑濾波處理獲得R^*及根據(jù)G在空域?qū)'_test進(jìn)行空域平滑濾除處理獲得W^*_test的具體過程為：

③_1、對G的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；

同樣，對R'的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；

同樣，對W'_test的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；

③_2、計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的方差，得到對應(yīng)的方差圖像，記為varG，varG的寬度為col且高度為row；

同樣，計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的方差，得到對應(yīng)的方差圖像，記為varR，varR的寬度為col且高度為row；

同樣，計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的方差，得到對應(yīng)的方差圖像，記為varW，varW的寬度為col且高度為row；

然后計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值與中對應(yīng)位置的區(qū)域塊中的所有像素點(diǎn)的像素值之間的協(xié)方差，得到對應(yīng)的協(xié)方差圖像，記為covGR，covGR的寬度為col且高度為row；同樣，計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值與中對應(yīng)位置的區(qū)域塊中的所有像素點(diǎn)的像素值之間的協(xié)方差，得到對應(yīng)的協(xié)方差圖像，記為covGW，covGW的寬度為col且高度為row；

③_3、計(jì)算中的每個區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第一系數(shù)和第二系數(shù)，將中左上角像素點(diǎn)的坐標(biāo)位置為(i,j)的區(qū)域塊與中左上角像素點(diǎn)的坐標(biāo)位置為(i,j)的區(qū)域塊之間的第一系數(shù)和第二系數(shù)記為a^R(i,j)和b^R(i,j)，然后將中的所有區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第一系數(shù)構(gòu)成第一系數(shù)圖像，記為a^R，a^R的寬度為col且高度為row；并將中的所有區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第二系數(shù)構(gòu)成第二系數(shù)圖像，記為b^R，b^R的寬度為col且高度為row；其中，1≤i≤col,1≤j≤row，covGR(i,j)表示covGR中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，varG(i,j)表示varG中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，ξ為懲罰因子；

同樣，計(jì)算中的每個區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第一系數(shù)和第二系數(shù)，將中左上角像素點(diǎn)的坐標(biāo)位置為(i,j)的區(qū)域塊與中左上角像素點(diǎn)的坐標(biāo)位置為(i,j)的區(qū)域塊之間的第一系數(shù)和第二系數(shù)記為a^W(i,j)和b^W(i,j)，然后將中的所有區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第一系數(shù)構(gòu)成第一系數(shù)圖像，記為a^W，a^W的寬度為col且高度為row；并將中的所有區(qū)域塊與中對應(yīng)位置的區(qū)域塊之間的第二系數(shù)構(gòu)成第二系數(shù)圖像，記為b^W，b^W的寬度為col且高度為row；其中，covGW(i,j)表示covGW中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值；

③_4、對a^R的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；接著計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；

同樣，對b^R的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；接著計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；

同樣，對a^W的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；接著計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；

同樣，對b^W的四周邊界進(jìn)行擴(kuò)充，且將擴(kuò)充的像素點(diǎn)的像素值設(shè)定為0，將擴(kuò)充得到的寬度為col+2且高度為row+2的圖像記為然后采用尺寸大小為3×3的滑動窗口，在中將橫坐標(biāo)范圍為[1,col]且縱坐標(biāo)范圍為[1,row]中的每個像素點(diǎn)對應(yīng)滑動窗口的左上角進(jìn)行逐像素點(diǎn)滑動，將劃分為col×row個互相重疊的尺寸大小為3×3的區(qū)域塊；接著計(jì)算中的每個區(qū)域塊中的所有像素點(diǎn)的像素值的均值，得到對應(yīng)的均值圖像，記為的寬度為col且高度為row；

③_5、獲取R^*，將R^*中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值記為R^*(i,j)，其中，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，G(i,j)表示G中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值；

同樣，獲取W^*_test，將W^*_test中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值記為W^*_test(i,j)，其中，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值，表示中坐標(biāo)位置為(i,j)的像素點(diǎn)的像素值。

4.根據(jù)權(quán)利要求3所述的一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于所述的步驟③_3中取ξ＝0.75。

5.根據(jù)權(quán)利要求3所述的一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于所述的步驟④中其中，sign()為求符號函數(shù)，P_max表示W(wǎng)^*_test與R^*的相關(guān)函數(shù)圖像P中最大的像素值，Ne×Ne表示W(wǎng)^*_test與R^*的相關(guān)函數(shù)圖像P中以像素值為P_max的像素點(diǎn)為中心的窗口區(qū)域的尺寸大小，8≤Ne≤15，Ω表示W(wǎng)^*_test與R^*的相關(guān)函數(shù)圖像P中的所有像素點(diǎn)的坐標(biāo)位置的集合，Ω_Ne×Ne表示W(wǎng)^*_test與R^*的相關(guān)函數(shù)圖像P中以像素值為P_max的像素點(diǎn)為中心的窗口區(qū)域中的所有像素點(diǎn)的坐標(biāo)位置的集合，P(s)表示W(wǎng)^*_test與R^*的相關(guān)函數(shù)圖像P中坐標(biāo)位置為s的像素點(diǎn)的像素值。

6.根據(jù)權(quán)利要求5所述的一種使用空域平滑濾波的數(shù)字照片圖像來源識別方法，其特征在于所述的步驟⑤中相關(guān)性判別閾值的獲取過程為：

⑤_1、使用I_test所聲明的來源相機(jī)任意拍攝U幅同源圖像；然后按照步驟②和步驟③的過程，以相同的方式獲得每幅同源圖像的最終的傳感器模式噪聲；其中，U≥100，每幅同源圖像的寬度為col且高度為row；

⑤_2、選取與I_test所聲明的來源相機(jī)不同的V部相機(jī)，要求選取的V部相機(jī)的分辨率大于或等于I_test所聲明的來源相機(jī)的分辨率；然后使用選取的V部相機(jī)任意拍攝U幅異源圖像，共得到U幅異源圖像，再將每幅異源圖像按照步驟⑤_1獲得的同源圖像的尺寸大小截取中心區(qū)域作為新的異源圖像；接著按照步驟②和步驟③的過程，以相同的方式獲得每幅新的異源圖像的最終的傳感器模式噪聲；其中，1≤V≤9；

⑤_3、采用SPCE方法計(jì)算R^*與每幅同源圖像的最終的傳感器模式噪聲的相關(guān)性，共得到U個同源樣本相關(guān)性值；

同樣，采用SPCE方法計(jì)算R^*與每幅異源圖像的最終的傳感器模式噪聲的相關(guān)性，共得到U個異源樣本相關(guān)性值；

⑤_4、在區(qū)間[0,30]內(nèi)選取不同的實(shí)數(shù)值對步驟⑤_3中獲得的U個同源相關(guān)性值和U個異源相關(guān)性值進(jìn)行判決，獲得每個實(shí)數(shù)值下的假正率(False Positive Rate，F(xiàn)PR)，記為FPR，其中，N_FP表示將異源圖像判決為同源的數(shù)量，N_TN表示將異源圖像判決為異源的數(shù)量；

⑤_5、將FPR＝0.01時對應(yīng)的實(shí)數(shù)值作為相關(guān)性判決閾值。