印刷品數字水印檢測器的實(shí)現

最近幾年，許多研究人員放棄了傳統密碼學(xué)的技術(shù)路線(xiàn)，嘗試用各種信號處理方法對聲像數據進(jìn)行隱藏加密，并將該技術(shù)用于制作多媒體的數字水印.數字水印(Digital Watermark)技術(shù)是指用信號處理的方法在數字化的多媒體數據中嵌入隱蔽的標記，這種標記通常是不可見(jiàn)的，只有通過(guò)專(zhuān)用的檢測器或閱讀器才能提取。數字水印是信息隱藏技術(shù)的一個(gè)重要研究方向。本文采用高端CMOS圖像傳感器采集圖像，基于DCT變換和相關(guān)算法，結合可編程邏輯器件FPGA和數字信號處理器DSP,實(shí)現了印刷品數字水印的快速檢測。

　　1 基本原理

　　1.1 色彩空間的選擇

　　本系統中含有數字水印的印刷品由計算機生成，通過(guò)激光打印機打印，再由CMOS傳感器采集。整個(gè)過(guò)程中，圖像信號經(jīng)歷了DAC變換和ADC變換，并且在不同設備上傳輸。為了保證圖像的失真最小，必須選擇一種獨立于設備的色彩空間來(lái)表示圖像。本系統采用CIE-XYZ色彩空間。XYZ顏色模型將彩色光表示為：

　　C=X(x)+Y(y)+Z(z)

　　其中，x、y、z是XYZ色彩空間的基色量，X、Y、Z為三色比例系數。

　　1.2 水印的添加和提取

　　所要添加的水印信息可以是一段長(cháng)度為n的由-1和1表示的二進(jìn)制數碼{ki},如：-1-1-1-1-1-1-1-1 1 1 1 1 1 1 1 1 -1-1-1-1-1-1-1-1.將ki與偽隨機序列M序列按位與，具體實(shí)現是ki序列長(cháng)度取M序列長(cháng)度的整數倍，從而得到一個(gè)具有良好自相關(guān)特性的序列，這就是水印hi.宿主圖像可以是任意一幅色彩豐富的彩色圖片。首先，將宿主圖像轉化為CIE-XYZ色彩空間方式，取出亮度分量，將其進(jìn)行16×16 DCT變換，選出每個(gè)塊的前10個(gè)較大的DCT系數Wi;然后將hi序列按照如下公式計算得到加有水印的DCT系數序列：

　　Wi′=Wi(1+аhi)

　　其中，a控制了水印的強度。最后對處理過(guò)的DCT系數進(jìn)行反DCT變換后與圖像的另外兩個(gè)基色量合成并打印，即可得到添加了水印的圖像。

　　水印的提取過(guò)程與水印的添加過(guò)程大體相反。為了保證圖像的大小與原圖相同，將一幅含有水印的印刷品在合適的燈光照明和適當的物距下，通過(guò)高端CMOS圖像傳感器及其外圍電路采集到存儲器，將圖像從RGB空間轉換為CIE-XYZ空間，提取其中的亮度分量;然后對其進(jìn)行16×16 DCT變換，選取前10個(gè)較大的DCT系數，用已知hi序列與之進(jìn)行相關(guān)運算。由于經(jīng)M序列調制的hi具有很好的自相關(guān)特性，當圖像中含有水印時(shí)，相關(guān)運算的值超過(guò)一定的閾值即表示此圖像含有水印，反之則表示這幅圖片中不含有水印信息。通過(guò)上述方法，可以在需要受到保護的圖像中添加水印信息，然后用本論文闡述的機讀設備進(jìn)行檢測，從而將不受保護的圖像檢測出來(lái)。對于含有水印信息圖像的復制印刷品，由于其水印能量被削弱，因此檢測得出的相關(guān)值偏低，也能將其淘汰。經(jīng)試驗得出，正品的相關(guān)值都在40%以上，而正品的復制品都在15%以下，不含水印信息的相同印刷品的相關(guān)值都在10%以下。

　　數字水印系統原理如圖1所示。

　　1.3 算法在FPGA及DSP上的實(shí)現

　　由于本文采用CMOS圖像傳感器所采集圖像的色彩空間屬于RGB空間，因此，必須先將其轉換為獨立于設備的CIE-XYZ空間，使后續處理得以繼續。

　　FPGA(Field-Programmable Gate Array)，即現場(chǎng)可編程門(mén)陣列，它是在PAL、GAL、CPLD等可編程器件的基礎上進(jìn)一步發(fā)展的產(chǎn)物。它是作為專(zhuān)用集成電路(ASIC)領(lǐng)域中的一種半定制電路而出現的，既解決了定制電路的不足，又克服了原有可編程器件門(mén)電路數有限的缺點(diǎn)。

　　本文所采用的數字水印算法包含了圖像的2-D DCT運算，計算量極大。為了保證計算的實(shí)時(shí)性，本文采用具有并行處理結構的FPGA實(shí)現。首先改進(jìn)2-D DCT算法，使之適合FPGA的并行處理，常用的方法是將二維運算轉化為兩次一維運算。

　　因此，二維DCT運算可由兩個(gè)一維運算得到，并且，公式中的余弦值可以事先計算好并做為運算參數保存。一維運算的實(shí)現采用了DA(Distributed Arithmetic)算法進(jìn)一步提高運算速度，具體原理如下：

　　Ai(i=0,1,2,……，15)表示某一行某一像素的亮度值，xk表示余弦變換的多項式系數。對于式(8)，宜采用16個(gè)相同的模塊進(jìn)行并行運算，所得結果為Z0,Z1,Z2,……，Z15,然后將Zy帶入式(4)，進(jìn)行與上一步相同的運算，得出DCT系數的值。以上運算用按位與運算和加法運算，可以避免對乘法器資源的占用，而且可以用FPGA中的查找表(LUT)實(shí)現。在電路設計上可采用乒乓操作分時(shí)復用一維運算模塊，以提高器件的利用率。

　　在得到圖片的DCT系數矩陣以后，取出每個(gè)16×16塊的前10個(gè)較大的DCT系數，形成一個(gè)一維序列，與前面提到的水印序列hi進(jìn)行相關(guān)運算。相關(guān)運算值超過(guò)一個(gè)由試驗得出的閾值則表明檢測到水印。