張 陽， 卿粼波， 何小海

（四川大學(xué) 電子信息學(xué)院， 成都 610065）

0 引 言

互聯(lián)網(wǎng)技術(shù)的迅速發(fā)展，極大的促進(jìn)了媒體信息的傳播，但也讓篡改者有機(jī)可乘。 針對創(chuàng)作者的版權(quán)保護(hù)，一種有效的方法是數(shù)字水印技術(shù)。 數(shù)字水印將版權(quán)信息嵌入到傳輸?shù)膱D像數(shù)據(jù)中，這些圖像數(shù)據(jù)與嵌入版權(quán)信息前肉眼觀察無差別，即便篡改者攻擊嵌入圖像也能大體提取出水印圖，能有效保護(hù)創(chuàng)作者版權(quán)。

現(xiàn)代數(shù)字水印技術(shù)發(fā)展迅速，出現(xiàn)了各種水印算法，傳統(tǒng)最低有效位（LSB）算法實(shí)現(xiàn)原理簡單，透明性好、但魯棒性差，嵌入水印的圖像稍受頻率域攻擊就無法提取準(zhǔn)確的水印。 變換域水印算法的出現(xiàn)大大改善了這一問題。 李文娜等應(yīng)用離散余弦變換（DCT）算法較為成熟，具有一定的抗攻擊性，但DCT 效應(yīng)導(dǎo)致圖像塊效應(yīng)嚴(yán)重，對于原圖質(zhì)量要求高的場景并不適用；以燕鯤鵬研究的分塊數(shù)字圖像水印算法為代表的離散小波變換（DWT）主流算法，對JPEG2000 等有損壓縮的攻擊具有更強(qiáng)的魯棒性，但提取水印需要原宿主圖或者原水印圖，無法實(shí)現(xiàn)水印盲提取，其透明性和魯棒性有進(jìn)一步提升的空間。 輪廓波（Contourlet）變換在小波變換的基礎(chǔ)進(jìn)一步發(fā)展，吳靜研究的基于Contourlet變換的圖像數(shù)字水印算法，相比小波變換有更加優(yōu)良的抗攻擊性，但仍然不能實(shí)現(xiàn)盲提取水印，并且對于原圖有更高要求。

本文在DWT－SVD 算法研究基礎(chǔ)上，通過在不同頻段嵌入水印，提取時(shí)利用各個頻段對原圖的估計(jì)影響因子進(jìn)行加權(quán)提取，融合成新的水印圖。 其中，給各個頻段進(jìn)行奇異值分解，利用奇偶嵌入，將水印的奇偶性嵌入奇異值中，再進(jìn)行逆奇異值分解（SVD），得到嵌入水印圖。 在對原圖影響可以接受的情況下，實(shí)現(xiàn)了更好的抗椒鹽攻擊效果，而且實(shí)現(xiàn)了水印圖像的盲提取。

1 理論基礎(chǔ)

1.1 算法框圖

傳統(tǒng)DWT 加水印算法，如圖1 所示。 由于圖像能量主要集中在低頻部分，傳統(tǒng)DWT 算法只將水印低頻二階小波變換低頻（LL2）分量信息加入到原圖低頻LL2 分量，并未考慮水印中存在的高頻分量的影響，提取的水印在魯棒性上存在一定缺陷。 因此，本文將水印高頻信息分別加入到原圖對應(yīng)高頻分量信息中，算法框架如圖2 所示。

圖1 傳統(tǒng)DWT 算法Fig.1 Traditional DWT algorithm

圖2 算法框架Fig.2 Algorithm block diagram

1.2 奇異值分解

奇異值分解應(yīng)用于正交矩陣中，是一種線性代數(shù)工具。 矩陣的SVD 分解定義：

、、都是二維矩陣，的大小；的大?。坏拇笮?。 奇異值公式（1）。

其中，

對于二維圖像（） 的SVD 分解，公式（2）：

其中，、、都是的二維矩陣，且是奇異值的對角矩陣。

圖像的奇異值矩陣具有較好的穩(wěn)定性，當(dāng)圖像被施加小的擾動，其奇異值變化不大，且在奇異值上嵌入水印不會損害圖像的幾何特性。 本文取原圖二維小波變換后的低中高頻8×8 塊作奇異值變換，得到矩陣，?。?，0） 作為計(jì)算后的奇異值，再嵌入水印值。

1.3 奇偶嵌入

奇偶嵌入的目的是將水印的奇偶性質(zhì)嵌入原圖中，從而可以對已經(jīng)嵌入圖像分析奇偶性，得到水印的值，實(shí)現(xiàn)水印盲提取。 本文采用二值圖像作為水印圖像，奇值為1，偶值為0。 把計(jì)算得到的奇異值設(shè)為，把二值水印圖像的值設(shè)為， 嵌入后的值設(shè)為，本文通過取整除法Math.floor（20） 計(jì)算出原圖奇異值的奇偶性，奇偶嵌入對照比見表1。

表1 奇偶嵌入對照表Tab.1 Comparison table of parity embedding

提取水印時(shí)根據(jù)各頻段計(jì)算奇異值，判斷Math.floor（20） 奇偶性得到原水印圖，實(shí)現(xiàn)水印圖像盲提取。

2 水印的嵌入與提取

2.1 水印的嵌入

實(shí)驗(yàn)選用512×512 的灰度圖線作為宿主圖像，64×64 的二值圖像作為水印圖像，水印嵌入過程圖，如圖3 所示。 水印嵌入算法具體描述如下：

圖3 水印嵌入過程圖Fig.3 Process diagram of watermark embedding

（1）已求得64×64 矩陣的Arnold 置換周期為48，對原始水印圖像置換24 次得到加密的水印圖像，并將其轉(zhuǎn)換為一個1×256 的行向量，長度為256；

（2）對宿主圖像進(jìn)行二階小波變換，分別對LL低頻段、LH中頻、HL中頻、HH高頻進(jìn)行處理，得到各頻段的256 三維塊，如原圖512×512 大小，那么LL可以處理為8×8×256 的三維塊，LH、HL、HH 都可以處理為16×16×256 的三維塊；

（3）遍歷水印行向量，對LL2 低頻段取8×8 塊（中高頻段取16×16 塊）求最大奇異值，作為參數(shù)，取對應(yīng)水印位置值，根據(jù)表1 計(jì)算得到該位置嵌入水印后的奇異值。 對HL1、LH1、HH1 3 個頻段的三維塊分別采用上述方法進(jìn)行嵌入處理。 如對HL1頻段，假如遍歷到第次，對一個16×16 的面進(jìn)行SVD 奇異值計(jì)算，取對角線第一個奇異值，也即最大奇異值作為， 再根據(jù)表1 計(jì)算得到該位置嵌入水印后的奇異值；

（4）對步驟（3）生成的奇異值進(jìn)行SVD 逆分解，得到8×8 的面或16×16 的面組成三維塊，再分別對應(yīng)還原為各個頻段，最后進(jìn)行小波逆變換，得到嵌入加密水印的圖像。

2.2 水印提取過程

以512×512 的含水印圖像進(jìn)行說明，默認(rèn)水印是64×64 大小的二值圖像，提取水印過程不需要水印原圖，水印提取過程如圖4 所示。 提取步驟如下：

圖4 水印提取過程Fig.4 Watermark extraction process

（1）對已嵌入圖像進(jìn)行二階小波變換，提取出LL、LH、HL、HH4 個頻段塊，并分別對這些頻段塊處理為三維塊，此例LL可以經(jīng)處理為8×8×256的三維塊，LH、HL、HH都可以處理為16×16×256的三維塊；

（2）從1 遍歷到256，加入遍歷到第次，對每個頻段三維塊第個面進(jìn)行奇異值計(jì)算，得到最大奇異值，通過取整除法Math.floor（20） 得出水印該位置的值，因此可以得到相應(yīng)LL、LH、HL、HH4 幅水印圖，，，；

（3）按照各個頻段加權(quán)決定最后水印圖，加權(quán)水印圖計(jì)算公式（3）：

（4）已求得64×64 大小矩陣置亂周期為48，將步驟（3）得到的加密圖繼續(xù)置亂24 次，得到可辨認(rèn)的水印圖。

3 實(shí)驗(yàn)結(jié)果及分析

實(shí)驗(yàn)選取64×64 像素的自制二值水印圖像water，512×512 像素 的宿主 圖像lenna、barbara、peppers，如圖5 所示。 分別對每張圖利用3 種不同算法嵌入水印，記錄嵌入水印前后的峰值信躁比（） 值；提取水印時(shí)，記錄提取水印后的歸一化相關(guān)系數(shù)（） 值。

圖5 水印原圖以及實(shí)驗(yàn)原圖Fig.5 Original watermarking and experiment

3.1 透明性測試

透明性指水印嵌入后的不可見性，也稱為透明性，一般用峰值信噪比（） 來衡量，表示宿主圖像的失真情況，其單位為dB（分貝）。 通過均方差（） 進(jìn)行定義，兩個的圖像和圖像，其均方差定義為公式（4），定義為公式（5）：

其中，max表示圖的圖像點(diǎn)顏色的最大值，本文中值為255。

結(jié)果表明，利用本文算法嵌入水印后的3 幅圖和DCT 變換、DWT 變換、Contourlet 變換＋SVD（簡稱Contourlet 算法）比較并無太大差異，都滿足透明性要求（30 dB），如圖6 所示。 分別對嵌入水印后的圖進(jìn)行參數(shù)為0.4 的低通濾波攻擊、參數(shù)為0.03 的椒鹽噪聲攻擊、參數(shù)為30 的壓縮、參數(shù)為1／4 面積的圖像裁剪，如圖7 所示。

圖6 含水印圖Fig.6 Image with watermark

圖7 被攻擊后圖像Fig.7 Image after attack

3.2 魯棒性測試

魯棒性一般用歸一化相關(guān)系數(shù)（） 衡量，的值越接近1，提取的水印與原始水印越相似，的計(jì)算公式（6）：

其中，（，） 表示原水印圖像，（，） 表示提取出的水印圖像。

圖7 中被攻擊圖像利用本文算法提取出的水印如圖8 所示。

圖8 提取水印圖Fig.8 Watermark extraction

實(shí)驗(yàn)針對每張測試圖，將本文算法與DCT 變換，DWT 變換以及Contourlet 變換作比較，并分別在3×3 模板的低通濾波器下（系數(shù)分別為0.4、0.25、0.1）攻擊、壓縮（質(zhì)量因子分別為90、60、30）攻擊、椒鹽噪聲（系數(shù)分別為0.03、0.02、0.01）攻擊，以及圖像中間1／4 剪切攻擊下對比各種算法的魯棒性。 測試得到的數(shù)值見表2～表4。

表2 lenna 圖提取水印NC 值Tab.2 NC value extraction of watermark from Lenna diagram

表3 barbara 圖提取水印NC 值Tab.3 NC value extraction of watermark from Barbara diagram

表4 peppers 圖提取水印NC 值Tab.4 NC value extraction of watermark from Peppers diagram

3.3 分析

對比各種算法的魯棒性，通過表2～表4 得知：

（1）不同圖片對于算法的測試結(jié)果基本無影響，提取水印的值隨攻擊強(qiáng)度的變化而變化；

（2）低通濾波系數(shù)越大，質(zhì)量因數(shù)越小，椒鹽噪聲越大，都會導(dǎo)致提取水印魯棒性變差；

（3）DCT 變換算法具有很好的透明性，但水印魯棒性差，稍加攻擊就無法提取水印，另外3 種算法犧牲了適量透明性（但仍然保持30），但水印魯棒性強(qiáng)；

（4）DWT 變換具有良好的透明性，抗剪切攻擊效果優(yōu)，但面對椒鹽噪聲攻擊，和質(zhì)量因數(shù)較低的壓縮攻擊，效果不能讓人滿意；

（5）Contourlet 變換有比較令人滿意的透明性和抗攻擊效果，但本文算法在透明性相差不大情況下，抗攻擊效果略勝一籌；

（6）本文算法在低通濾波，壓縮攻擊后提取效果和DWT 變換，Contourlet 變換相差不大的情況下，做到了在椒鹽噪聲攻擊下提取效果的提升。

4 結(jié)束語

本文提出了一種基于DWT－SVD 魯棒水印算法，在傳統(tǒng)離散小波變換基礎(chǔ)上，進(jìn)行奇異值分解，利用圖像奇異值良好的穩(wěn)定性，提高圖像的抗攻擊性；在圖像奇異值基礎(chǔ)上，進(jìn)行奇偶量化嵌入和提取水印，可實(shí)現(xiàn)盲提??；為了進(jìn)一步提高抗攻擊性，在圖像的各個頻段通過加權(quán)融合提取水印。 仿真和攻擊實(shí)驗(yàn)結(jié)果表明，與其它算法進(jìn)行比較分析，該算法的魯棒性有一定提高。