楊江濤，馬喜宏*，鄔　琦

(1.中北大學電子測試技術(shù)國家重點實驗室，太原030051; 2.中北大學儀器科學與動態(tài)測試教育部重點實驗室，太原030051; 3.中北大學計算機與控制工程學院，太原030051)

?

數(shù)字圖像復制粘貼篡改盲取證算法研究與實現(xiàn)

楊江濤1，2，馬喜宏1，2*，鄔琦3

(1.中北大學電子測試技術(shù)國家重點實驗室，太原030051; 2.中北大學儀器科學與動態(tài)測試教育部重點實驗室，太原030051; 3.中北大學計算機與控制工程學院，太原030051)

摘要:針對數(shù)字圖像檢測中一類常見的復制－粘貼圖像篡改，提出了一種基于小波變換和奇異值分解的檢測算法。該算法利用小波變換提取圖像的低頻分量，對低頻分量分塊提取奇異值特征，然后將特征矢量進行按行字典排序，并且配合圖像塊的偏移位置信息，進行圖像復制偽造區(qū)域的檢測和定位。實驗表明該算法大大減小了特征向量的維數(shù)，從而提高了相似塊的匹配檢測效率。為了更方便快捷的檢測圖像是否被惡意篡改，設(shè)計了簡單明了的系統(tǒng)檢測界面，只需載入待檢測的圖像并輸入相應(yīng)的參數(shù)就能進行檢測，最后將檢測結(jié)果返回給界面，而且系統(tǒng)完成了篡改檢測算法的DSP硬件實現(xiàn)，該算法將有利于推動數(shù)字圖像取證技術(shù)的理論研究與應(yīng)用推廣的發(fā)展。

關(guān)鍵詞:數(shù)字圖像;篡改取證;小波變換;奇異值分解; DSP

圖像數(shù)據(jù)作為信息傳遞的重要載體，被廣泛應(yīng)用在新聞報道、智能信息收集、犯罪調(diào)查、安全監(jiān)控、醫(yī)療和科學研究等領(lǐng)域［1－2］。但隨著成本、高性能、友好操作界面的圖像處理軟件的不斷推出，數(shù)字化的圖像信息面臨著極大的安全隱患(即圖像信息易被惡意篡改)，因此數(shù)字圖像真實性認證已經(jīng)成為目前的研究熱點［3－4］。

為了進一步減少運算量，本文從縮小圖像尺寸和降低特征空間維數(shù)兩方面考慮，提出了一種基于小波和奇異值分解的圖像復制區(qū)域檢測算法。該算法利用小波變換提取的圖像低頻分量作為分析對象，并使用奇異值分解提取小波低頻圖像特征進行降維。實驗表明該算法能夠定位出復制和粘貼的圖像區(qū)域，并顯著地減少了運算量，提高了檢測效率。為了更方便快捷的檢測圖像是否被惡意篡改，設(shè)計了簡單明了的系統(tǒng)檢測界面，只需載入待檢測的圖像并輸入相應(yīng)的參數(shù)就能進行檢測，最后將檢測結(jié)果返回給界面，而且系統(tǒng)完成了篡改檢測算法的DSP［5］硬件實現(xiàn)，該算法將有利于推動數(shù)字水印技術(shù)的理論研究與應(yīng)用推廣的發(fā)展。

1　數(shù)字圖像復制粘貼篡改盲取證算法

1.1一級小波窮舉搜索法

1.1.1算法思想

該算法［6］首先將圖像進行不重疊分塊，把每一個圖像塊當成一個模板，然后對圖像塊進行一級小波變換，將經(jīng)過一級小波變換的圖像塊遍歷剩余的部分(圖像塊同樣要經(jīng)過一級小波變換)，看是否存在跟模板完全一樣的圖像塊，若存在，則標記這兩個圖像塊的位置，繼續(xù)搜索，直到找到所有的相同塊。

首先將圖像分為ms×ns的不重疊塊B，然后對小塊做DWT，把DWT后圖像小塊當成一個模板，遍歷圖像剩余的部分(圖像塊同樣要經(jīng)過一級小波變換)，看是否存在跟模板完全一樣的圖像塊，若存在，則標記這兩個圖像塊的位置，繼續(xù)搜索，直到找到所有的相同塊。

1.1.2算法步驟

(1)不重疊快的劃分和特征向量的提取

將圖像分為ms×ns的不重疊塊，其中塊大小的選擇要求小于篡改塊。假設(shè)圖像大小為mc×nc，重疊塊Bi大小為ms×ns，對于圖像塊的劃分，分塊越大檢測的速度越快，在同一閾域T前提下漏檢的可能性也越大;與之相對應(yīng)分塊劃分過小，檢測的準確性雖然提高了，但降低了檢測的速度，同時可能會出現(xiàn)大量的虛警現(xiàn)象，所以關(guān)于塊的大小的選擇要根據(jù)人工懷疑的篡改區(qū)域大體估計，并通過反復試驗，適當選擇適當?shù)姆謮K大小。對每個小塊做Haar離散小波變換(DWT):［Ai，Hi，Vi，Di］= dwt2 (Bi，’harr’)，圖像DWT后所對應(yīng)的各子頻帶的值為相應(yīng)的小波系數(shù)。低頻逼近子帶表示小波分解在最大尺度、最小分辨率下對原始圖像的最佳逼近，它的統(tǒng)計特征和原圖像相似，含有最重要的信息，所以我們采用低頻子帶Ai作為小塊Bi的特征向量。

(2)重疊快的劃分和特征向量的提取

將圖像按照單像素位移被劃分為多個重疊塊Ci，圖像塊的大小和重疊塊的大小一樣，同為ms×ns大小。對每個小塊做Haar離散小波變換(DWT):［Ai，Hi，Vi，Di］=dwt2(Ci，’harr’)。

(3)塊匹配檢測

每一個DWT重疊塊和所有DWT非重疊塊做相似性檢測。常用的塊相似性檢測方法有基于直方圖統(tǒng)計的檢測方法，基于像素統(tǒng)計的檢測方法和基于像素統(tǒng)計的檢測方法。本文采用Pearson相關(guān)系數(shù)檢測方法。

1.2DCT水印算法

1.2.1算法思想

該算法［7］首先將水印圖像置亂，然后將水印嵌入到原始圖像中。檢測時，若提取出的水印和原始水印不一致，則說明圖像被篡改過。

1.2.2算法步驟

(1)水印圖像置亂，本文采用的置亂算法為Arnold變換。

(2)圖像分為8×8大小塊并進行DCT變換，把水印嵌入圖像塊中。

(3)檢測時，提取水印，反置亂得到提取水印，若得到的水印和原始水印相同則圖像未被篡改，否則為篡改圖像［8］。

1.3LSB水印算法

1.3.1算法思想

該算法［9］的主要思想是圖像2×2不重疊分塊，提取圖像分塊的灰度均值并量化為8 bit的精確值，經(jīng)混沌logistic映射加密和確定嵌入位平面后，嵌入到原始圖像所有像素的2個最低有效位中;在圖像需要認證時，提取圖像中水印信息和重新計算圖像分塊灰度均值，經(jīng)加密和確定嵌入位平面后生成參考水印，比較提取的水印信息與參考水印，若相同則通過認證，反之判定為篡改。在認證檢測過程中，不需要原始圖像和任何有關(guān)水印的附加信息，保證了水印的安全性和保密性，實現(xiàn)對圖像內(nèi)容完整性認證和篡改的精確定位。

1.3.2算法步驟

(1)原始圖像分2×2不重疊分塊，提取圖像塊的灰度均值并量化為8 bit的精度值。

(2)求分塊高灰度均值二值化。

(3)經(jīng)混沌logistic映射加密和確定嵌入為平面后，嵌入原始圖像所有像素的2個最低有效位中。

(4)在圖像需要認證時，提取圖像中水印信息和重新計算圖像分塊灰度均值，經(jīng)加密和確定嵌入位平面后生成參考水印。

(5)比較提取的水印信息與參考水印，若相同則通過認證，反之判定為篡改。定位篡改位置。

1.4基于小波和奇異值分解算法

1.4.1算法思想

本文算法首先對檢測圖像進行小波變換，對小波低頻分量進行固定大小的滑窗操作，對取出的圖像塊進行奇異值分解［10］，然后對圖像奇異值特征矩陣進行按行字典排序。根據(jù)相似的圖像塊則奇異值矢量相近，排序后兩個矢量會在排序矩陣中相鄰，遍歷排序后矩陣，并且配合圖像塊的偏移頻率信息，檢測出復制偽造區(qū)域。

1.4.2算法步驟

(1)設(shè)Imc×nc為待檢測的灰度圖像，其中mc×nc為圖像的大小。

(2)對圖像I進行j級小波變換，變換后的小波低頻圖像為Am×n，其中m=［mc/2j］，n =［nc/2j］，［］為下取整操作。

(3)設(shè)定窗口塊w的大小為ms×ns(這里假設(shè)設(shè)定分塊的大小小于要檢測的復制圖像區(qū)域的大小)，對小波低頻圖像進行一次移動一個像素的滑窗操作，則共有Nw=(m－ms+1)×(n－ns+1)個塊。

(4)對每個圖像塊進行奇異值分解，得到降維后的圖像特征描述，特征空間由ms×ns維下降到了r維，其中r=min(ms，ns)。

(5)對得到的Nw×r維的特征矩陣，按行進行字典排序，得到排序后的特征矩陣S。向量si表示S中的一行，i=1，…，Nw，(xi，yi)表示與向量si對應(yīng)圖像塊的左上角坐標值。

(6)遍歷S，對S中相鄰的兩行向量si，sy，計算它們對應(yīng)的圖像塊的坐標的偏移矢量(Δx，Δy)。

(7)對滿足Δx＞|ms|orΔy＞|ns|的向量si，sy，統(tǒng)計其偏移頻率C(Δx，Δy)，得到偏移頻率矩陣C(初始化為零):C(Δx，Δy)= C(Δx，Δy)+1

(8)對得到偏移頻率矩陣C進行大于閾域T判斷，判別出圖像中是否存在復制圖像區(qū)域。如果存在，進入下一步，如果不存在，則算法結(jié)束。

(9)對偏移頻率大于閾值T的那些偏移矢量對應(yīng)的塊位置進行標識:如果這些具有相同偏移矢量的塊是相鄰的，則組成了復制粘貼的區(qū)域;如果標識塊是孤立的，則可能是噪聲。

2　DSP硬件平臺設(shè)計

本系統(tǒng)主要由圖像采集和圖像處理兩大部分組成。由CMOS傳感器負責圖像的采集，CPLD進行邏輯控制，SRAM進行數(shù)據(jù)存儲，DSP進行數(shù)據(jù)處理，將數(shù)據(jù)采集和數(shù)據(jù)處理集成到一個系統(tǒng)中，從而發(fā)揮了不同類型可編程芯片的，提高了系統(tǒng)的性能。系統(tǒng)的硬件框圖如圖1所示。

圖1　系統(tǒng)硬件設(shè)計框圖

系統(tǒng)圖像數(shù)據(jù)的采集、處理過程如下:首先，由DSP啟動對圖像數(shù)據(jù)的采樣，通過擴展的I/O口，發(fā)送控制信號到CPLD;控制信號主要由開始采集、采集完成等幾個狀態(tài)組成。CPLD一旦接收到開始采集信號，就立刻向CMOS傳感器發(fā)送采集控制信號。為了實現(xiàn)SRAM中數(shù)據(jù)地址的同步，將CMOS傳感器輸出的數(shù)字信號送往CPLD進行緩沖，然后送到SRAM中供DSP讀取。DSP通過查詢和中斷方式，監(jiān)控CPLD發(fā)出的中斷信號，一旦接收到CPLD發(fā)出的采集完成中斷信號，DSP在通知CPLD停止采集的同時，延時一段時間后從SRAM讀取數(shù)據(jù)，然后在內(nèi)部進行數(shù)據(jù)處理，將處理結(jié)果通過顯示器件顯示或者根據(jù)用戶的要求，將放在存儲器中的結(jié)果通過McBSP、USB接口送往PC機。

2.1圖像數(shù)據(jù)格式轉(zhuǎn)換

為了彌補CCS環(huán)境圖像顯示能力的不足，系統(tǒng)對圖像進行的性能評估工作主要在MATLAB平臺下完成。DSP平臺主要用來實時采集圖像并提取水印信息，圖像在CCS與MATLAB之間通過.dat圖像數(shù)據(jù)文件進行傳輸。

MATLAB程序可以直接讀取圖像文件，圖像數(shù)據(jù)經(jīng)過處理后，通過格式轉(zhuǎn)換程序，將圖像數(shù)據(jù)封裝成CCS可識別的.dat數(shù)據(jù)文件。CCS可以通過File－＞Data－＞Load很方便地將文件中的數(shù)據(jù)導入到DSP的存儲空間，利用DSP對圖像數(shù)據(jù)進行快速處理。然后通過File－＞Data－＞Save將存儲空間中的圖像數(shù)據(jù)導出至.dat數(shù)據(jù)文件，在MATLAB中通過格式轉(zhuǎn)換程序?qū)⑽募械膱D像數(shù)據(jù)轉(zhuǎn)換成MATLAB可識別的unit8格式圖像數(shù)據(jù)，進而對圖像進行顯示輸出和質(zhì)量評估。

2.2視頻驅(qū)動程序設(shè)計

TI公司為C64X系列DSP提供了一種類/微型驅(qū)動模型，同時為開發(fā)DSP的外設(shè)驅(qū)動程序提供了驅(qū)動程序開發(fā)包DDK，而且提供專門針對視頻設(shè)備的FVID模型。FVID模型將微型驅(qū)動也分為兩層結(jié)構(gòu):上層為通用視頻端口層;下層為編解碼芯片微驅(qū)動層，它們通過外部設(shè)備控制接口EDC實現(xiàn)對外圍編解碼芯片的操作。

系統(tǒng)在DSPBIOS的用戶驅(qū)動配置中添加了3個自定義視頻輸入輸出設(shè)備驅(qū)動，并在其屬性中添加GIO的驅(qū)動實例。在任務(wù)管理中添加任務(wù)，并綁定函數(shù)作為視頻驅(qū)動的處理程序。主程序中視頻驅(qū)動處理任務(wù)函數(shù)的處理過程可描述為:在初始化部分建立視頻采集和顯示設(shè)備的實例句柄，利用FVID相關(guān)函數(shù)打開并初始化設(shè)備以取得句柄，同時分別申請一個輸入和輸出緩沖區(qū)。進入任務(wù)循環(huán)后，通過DMA操作將輸入緩沖區(qū)中的數(shù)據(jù)讀出來，對其進行處理，并將數(shù)據(jù)送輸出緩沖區(qū)進行顯示，然后向視頻口驅(qū)動程序切換輸入和輸出緩沖區(qū)完成更新。操作完成后釋放句柄，關(guān)閉驅(qū)動設(shè)備。

2.3無線串口通信程序設(shè)計

系統(tǒng)串口模塊采用TL16C752B通用異步收發(fā)器進行設(shè)計。對TL16C752B訪問之前需要配置一些控制寄存器對其進行初始化。主程序中首先定義串口結(jié)構(gòu)體和串口句柄，使用初始化函數(shù)打開一個串口通道，并返回該通道句柄，此后便可利用此句柄進行數(shù)據(jù)的接收和發(fā)送。數(shù)據(jù)收發(fā)可以通過直接讀寫模塊寄存器完成，也可以利用McBSP模塊的CSL函數(shù)來實現(xiàn)。

為了解決DSP硬件與上位機進行串口通信受限于串口線長度的問題。系統(tǒng)使用2塊DRF1605H模塊，實現(xiàn)計算機與DSP開發(fā)板之間的無線串口通信。DRF1605H模塊是利用TI的CC2530F256設(shè)計的，它運行ZigBee2007協(xié)議。系統(tǒng)配置上位機ZigBee模塊作為主節(jié)點，DSP硬件ZigBee模塊作為從節(jié)點，為它們配置相同的無線頻率及網(wǎng)絡(luò)ID，從而組成一個ZigBee網(wǎng)絡(luò)。節(jié)點上電后即可自動組網(wǎng)，主節(jié)點自動給從節(jié)點分配網(wǎng)絡(luò)地址，采用數(shù)據(jù)透明傳輸模式，主節(jié)點從計算機串口接收到的數(shù)據(jù)會自動發(fā)送給從節(jié)點，從節(jié)點由DSP串口接收到的數(shù)據(jù)會自動發(fā)送給主節(jié)點，主從節(jié)點之間的通信就像用一條串口線直接相連一樣方便。

3　仿真實驗與結(jié)果分析

3.1系統(tǒng)測試實驗

(1)測試一級小波窮舉搜索法如圖2所示。

用例描述:在基于一級小波窮舉搜索法界面，輸入正確參數(shù);

前置條件:進入一級小波窮舉搜索法界面;

輸入:要檢測的圖像，圖像塊長為10，寬為15，閾值為0.95;

測試結(jié)果:正確顯示執(zhí)行結(jié)果。

(2)測試DCT水印算法如圖3所示。

用例描述:在嵌入DCT水印算法界面，載入正確圖像并點擊嵌入水印按鈕;

前置條件:進入嵌入DCT水印算法界面;

輸入:要嵌入水印的圖像和水印圖像;

測試結(jié)果:正確顯示執(zhí)行結(jié)果。

圖2　一級小波窮舉搜索法測試圖

圖3　DCT水印法測試圖

(3)測試LSB水印算法如圖4所示。

用例描述:在LSB水印算法界面，載入圖像并點擊嵌入水印按鈕;

前置條件:進入LSB水印算法界面;

輸入:要檢測的圖像;

測試結(jié)果:正確顯示執(zhí)行結(jié)果。

(4)測試基于小波和奇異值分解算法如圖5所示。

用例描述:在基于小波和奇異值分解算法界面，輸入?yún)?shù)正確;

前置條件:進入基于小波和奇異值分解算法界面;

輸入:要檢測的圖像，小波級數(shù)為3，閥值為100;

測試結(jié)果:正確顯示執(zhí)行結(jié)果。

圖4　LSB水印法測試圖

圖5　基于小波和奇異值分解算法測試圖

3.2測試結(jié)果分析

基于以上算法，系統(tǒng)在MATLAB平臺下，對彩色圖像嵌入水印后，對含水印圖像進行各種攻擊處理后提取水印，對提取結(jié)果進行理論分析對比。同時利用DSP系統(tǒng)采集經(jīng)打印機輸出的含水印圖像，在DSP平臺下進行水印提取，然后將圖像數(shù)據(jù)傳輸給MATLAB程序，在MATLAB平臺上對圖像進行質(zhì)量評估，以驗證本水印算法有效抵抗打印掃描過程中各種圖像攻擊處理的魯棒性。

本文提出的算法是基于小波和奇異值分解的圖像復制區(qū)域檢測算法，而所有的基于圖像塊匹配的算法時間消耗都比較高，該算法利用小波變換提取的圖像低頻分量作為分析對象，并使用奇異值分解提取小波低頻圖像特征進行降維。通過表1可以比較得出，在相同實驗條件下，本文算法的圖像塊的特征維數(shù)為6，而窮舉搜索法、DCT水印算法和LSB水印算法中的圖像特征維數(shù)分別是無、64維和32維。因此，本文提出的算法減少了提取的特征維數(shù)和運算量，大大降低了算法的復雜度。

表1　本文算法與其他算法的比較結(jié)果

4　結(jié)論

圖像部分區(qū)域復制－粘貼篡改是常見的篡改方式，在日常生活新聞媒體醫(yī)學法律等方面都有應(yīng)用。本文提出了一種基于小波和奇異值分解的圖像復制區(qū)域檢測算法，該算法利用小波變換提取的圖像低頻分量作為分析對象，并使用奇異值分解提取小波低頻圖像特征進行降維。實驗表明該算法能夠定位出復制和粘貼的圖像區(qū)域，并顯著地減少了運算量，提高了檢測效率。為了更方便快捷的檢測圖像是否被惡意篡改，設(shè)計了簡單明了的系統(tǒng)檢測界面，只需載入待檢測的圖像并輸入相應(yīng)的參數(shù)就能進行檢測，最后將檢測結(jié)果返回給界面。本系統(tǒng)將圖像篡改檢測算法在DSP平臺上進行實現(xiàn)，促進了圖像數(shù)字圖像取證技術(shù)理論與實際應(yīng)用的有效結(jié)合。系統(tǒng)進一步完善后如能實現(xiàn)產(chǎn)品應(yīng)用推廣，將有力推動數(shù)字圖像取證技術(shù)產(chǎn)品化的發(fā)展。

參考文獻:

［1］李亞標，王寶光，李溫溫.基于小波變換的圖像紋理特征提取方法及其應(yīng)用［J］.傳感技術(shù)學報，2009，22(9):1308－1311.

［2］歐紅玉，陳曦，宋燕輝，等.基于LBP的圖像復制篡改檢測［J］.計算機應(yīng)用與軟件，2013，30(9):170－172.

［3］葉鴻敏，沈永增，張敏捷.基于二代小波提升算法的快速圖像邊緣檢測［J］.計算機測量與控制，2007，15(7):866－869.

［4］崔宇寅.數(shù)字圖像盲取證技術(shù)研究及其應(yīng)用［D］.上海:華東理工大學，2014.

［5］李志勇，張鵬.基于DSP的偏振圖像快速融合研究［J］.電子器件，2010，33(4):498－501.

［6］龔薇，斯科，葉秀清，等.一種強魯棒性的實時圖像增強算法［J］.傳感技術(shù)學報，2007，20(9):2024－2028.

［7］林克正，李紹華，李東勤.基于DWT和DCT的自適應(yīng)水印算法［J］.微計算機信息，2007，23(6):284－286.

［8］古麗江·庫爾班.基于DCT編碼圖像分割算法研究［J］.電子器件，2014，37(1):168－170.

［9］何冰.基于LSB的抗旋轉(zhuǎn)攻擊魯棒性數(shù)字水印算法［J］.系統(tǒng)仿真技術(shù)，2010，6(4):304－307.

［10］蔣華，張敏.基于提升小波和奇異值分解的灰度水印算法［J］.計算機應(yīng)用研究，2009，26(8):3028－3030.

楊江濤(1988－)，男，山西運城人，碩士研究生，研究方向為慣性測試技術(shù)及控制;

鄔　琦(1989－)，女，山西大同人，碩士研究生，研究方向為計算機仿真與圖形圖像處理，wuqi6567@163.com。

馬喜宏(1964－)，女，漢族，山西大同人，優(yōu)異級高級工程師，碩士生導師，研究方向為微系統(tǒng)集成技術(shù);

收稿日期:2014－05－23修改日期:2014－06－26

中圖分類號:TP391

文獻標識碼:A

文章編號:1005－9490(2015)04－0961－06