基于計(jì)算全息的強(qiáng)魯棒性自適應(yīng)水印算法

白韜韜，劉 真，盧 鵬

（上海理工大學(xué) 出版印刷與藝術(shù)設(shè)計(jì)學(xué)院，上海200093）

1 引 言

計(jì)算機(jī)技術(shù)和網(wǎng)絡(luò)通信的飛速發(fā)展，使得數(shù)字產(chǎn)品的獲取、使用、傳播變得異常便捷。隨著版權(quán)保護(hù)問題的日益突出，數(shù)字音頻、視頻、圖像的版權(quán)保護(hù)也引起了越來(lái)越多的重視，成為當(dāng)今社會(huì)和廣大學(xué)者的研究熱點(diǎn)。數(shù)字水印技術(shù)［1－5］為解決版權(quán)保護(hù)問題提供了一個(gè)有效的途徑，已經(jīng)廣泛應(yīng)用到了版權(quán)保護(hù)領(lǐng)域。目前，大多數(shù)水印算法都是針對(duì)數(shù)字圖像的傳播、保護(hù)進(jìn)行設(shè)計(jì)，且只能有效抵抗一些常見的圖像處理攻擊，對(duì)印刷稿圖像的版權(quán)無(wú)法進(jìn)行有效的保護(hù)。但是印刷圖像的版權(quán)保護(hù)同等重要，如何設(shè)計(jì)一種算法使得數(shù)字圖像即可以有效抵抗打印掃描攻擊，又可以有效抵抗一般圖像處理攻擊是該領(lǐng)域的難點(diǎn)問題。

近年來(lái)，結(jié)合現(xiàn)代信息光學(xué)技術(shù)的數(shù)字全息水印日益受到學(xué)者關(guān)注。2002年，Takai等［6］首次將數(shù)字全息技術(shù)引入數(shù)字水印領(lǐng)域，提出了數(shù)字全息水印技術(shù)。雖然該方法結(jié)合了全息圖像良好的抵抗裁切攻擊的性能，但是由于其選擇在空間域?qū)⒃紙D像和全息水印直接疊加，并且水印嵌入前要對(duì)原始圖像做低通濾波處理，細(xì)節(jié)信息大量丟失，導(dǎo)致含水印圖像視覺質(zhì)量差，且安全性不高，抵抗其他攻擊的能力也有待進(jìn)一步提高。隨后，許多學(xué)者［7－10］對(duì)全息水印及其嵌入方式進(jìn)行了 大 量、深 入 的 研 究。2007 年，Sun 等［11－12］提出了加密同軸全息水印技術(shù)，該種技術(shù)可以抵抗圖像的打印掃描攻擊，并取得了初步的成果。但其得到的含水印圖像視覺效果較差。并且由于同軸全息解密生成共軛像，使得從全息圖恢復(fù)出的圖像背景疊加了嚴(yán)重的噪聲，對(duì)水印圖像的辨別產(chǎn)生影響。2008年，Sun等［13］提出了彩色印刷圖像的計(jì)算全息數(shù)字水印。但是其并沒有說(shuō)明算法抵抗其他圖像處理攻擊的魯棒性。2009 年，Huang等［14］提出了一種新的生成計(jì)算全息的方式，即共軛對(duì)稱延拓傅立葉計(jì)算全息，并在2011年將生成的全息水印嵌入到圖像中［15］，實(shí)現(xiàn)了含水印圖像經(jīng)打印掃描后的水印提取。但是水印的安全性和抵抗其他圖像處理攻擊的能力需要進(jìn)一步提高。

本文提出了一種基于計(jì)算全息的強(qiáng)魯棒性自適應(yīng)水印算法，充分利用了Contourlet變換在圖像處理中的獨(dú)特優(yōu)點(diǎn)——多方向性和各向異性來(lái)自適應(yīng)隱藏兩個(gè)水印信息，設(shè)計(jì)水印算法。實(shí)驗(yàn)證明，本算法具有較高的安全性和不可見性，可以抵抗一般的圖像處理攻擊，并且經(jīng)過(guò)打印掃描后提取出的水印信息清晰。

2 基于計(jì)算全息的強(qiáng)魯棒性自適應(yīng)水印算法

2.1 共軛對(duì)稱延拓傅立葉計(jì)算全息

本文利用共軛對(duì)稱延拓傅立葉計(jì)算全息方法生成全息水印圖像。設(shè)原始實(shí)值水印圖像為已經(jīng)歸一化的f（x，y），圖像的大小為mpixel×npixel.待記錄的復(fù)數(shù)圖像φ（x，y）＝f（x，y）exp［j2πp（x，y）］，其中p（x，y）為［0，1］范圍內(nèi)均勻分布的隨機(jī)數(shù)。

將復(fù)數(shù)圖像φ（x，y）進(jìn)行共軛對(duì)稱延拓，得到M pixel×N pixel大小的復(fù)數(shù)圖像，其中M＝2×（m＋1），N＝n＋1。如式（1）：

式中：上標(biāo)“＊”表示對(duì)復(fù)數(shù)做共軛運(yùn)算，Φ（x，y）為共軛對(duì)稱延拓后的復(fù)數(shù)圖像。對(duì)Φ（x，y）做二維離散傅立葉變換，得到實(shí)值的全息圖像.設(shè)得到的全息圖像為H（x，y），則：

式中：F 表示傅立葉變換。由傅立葉變換的性質(zhì)，即共軛對(duì)稱的復(fù)數(shù)矩陣經(jīng)過(guò)快速傅立葉變換得到實(shí)值矩陣，Φ（x，y）經(jīng)過(guò)DFT 后得到實(shí)值圖像H（x，y）。此實(shí)數(shù)圖像記錄了復(fù)值圖像Φ（x，y）的振幅和相位信息。將此圖像經(jīng)過(guò)適當(dāng)編碼即可得到實(shí)值數(shù)字全息圖。

2.2 Contourlet變換

2002年，Do和Vetterli［16］提出了一種“真正”的圖像二維表示法：Contourlet變換。Contourlet變換，也稱金字塔型方向?yàn)V波器組PDFB（Pyramidal Directional Filter Bank）。PDFB 首先利用拉普拉斯金字塔（Laplacian Pyramid，LP）變換對(duì)原始圖像進(jìn)行多尺度分解，以“捕獲”二維圖像信號(hào)中存在的奇異點(diǎn)。然后由方向?yàn)V波器組（Directional Filter Bank，DFB）將分布在同方向上的奇異點(diǎn)連成線，形成輪廓段。PDFB 通過(guò)雙濾波器組將多尺度分析和方向分析分開進(jìn)行，這種互補(bǔ)的結(jié)構(gòu)使其具有較優(yōu)的非線性逼近性能。Contourlet變換不僅繼承了小波變換的多分辨率和時(shí)頻局部性，而且允許每個(gè)尺度上有不同數(shù)目的方向，提供了更高水平的方向性和各向異性。圖1和圖2分別給出了Contourlet對(duì)圖像進(jìn)行分解的過(guò)程示意圖和Contourlet變換頻域子帶分布示意圖。圖3 展示了應(yīng)用Contourlet變換對(duì)lena圖像進(jìn)行二級(jí)分解的結(jié)果，高頻和中頻分別分解為4方向和8方向。

圖1 Contourlet變換分解Fig.1 Contourlet transform decomposition

圖2 Contourlet變換尺度和方向分布Fig.2 Frequency subbands distribution of Contourlet transform

圖3 灰度圖像lena的Contourlet分解結(jié)果Fig.3 Contourlet decomposition results of the gray image“l(fā)ena”

2.3 算法分析

由于全息圖像本身信息量大，一般二值水印嵌入方法無(wú)法嵌入大信息量的全息水印，使得利用全息圖像作為水印的算法較少，并且多數(shù)全息水印算法無(wú)法有效抵抗JPEG 壓縮等常見的圖像處理攻擊。為了增強(qiáng)算法的魯棒性，本文嵌入了兩個(gè)全息水印圖像，進(jìn)一步增加了水印的嵌入量，所以本文選擇具有很高冗余性的Contourlet變換域嵌入水印信息。Contourlet變換冗余度為4／3，這種性質(zhì)非常適合隱藏大信息量的全息水印。目前，基于Contourlet變換的數(shù)字水印算法［17－20］也提出了很多，但現(xiàn)有大多數(shù)算法都是簡(jiǎn)單利用LP分解出的低頻信息更加穩(wěn)定這一特點(diǎn)來(lái)設(shè)計(jì)水印算法，并沒有充分利用Contourlet變換在圖像處理中特有的優(yōu)勢(shì)：更高水平的方向性和各向異性這一特點(diǎn)，并且應(yīng)用此種變換結(jié)合全息技術(shù)設(shè)計(jì)的數(shù)字水印算法并不多見。本文正是考慮到Contourlet變換可以將LP分解出的中高頻信息進(jìn)行多方向分解，提供更多水印嵌入位置這一特性來(lái)隱藏水印。在嵌入位置選擇上，本文進(jìn)行了巧妙的設(shè)計(jì)，由于分解出的各個(gè)中頻子帶信息與全息圖像的相似度不同，故當(dāng)全息水印替換不同中頻信息完成水印嵌入后，逆變換得到的含水印圖像與原始圖像的PSNR 值也會(huì)產(chǎn)生很大的變化。所以本文將得到的PSNR 最大值對(duì)應(yīng)的圖像作為最終的含水印圖像，可以在不損失相似度的情況下進(jìn)一步提高嵌入強(qiáng)度，實(shí)現(xiàn)了水印的自適應(yīng)嵌入，增強(qiáng)算法的魯棒性。本文算法充分利用了Contourlet變換多方向性和各向異性的特點(diǎn)，將全息水印以本文設(shè)計(jì)方法自適應(yīng)嵌入到原始圖像穩(wěn)定的中頻信息中，使其可以有效抵抗JPEG 壓縮、低通濾波攻擊等一般圖像處理攻擊和打印掃描攻擊。

3 水印的嵌入與提取

3.1 全息水印的預(yù)處理

利用2.1節(jié)共軛對(duì)稱延拓計(jì)算全息生成兩個(gè)初始全息水印圖像，大小均為M pixel×N pixel。針對(duì)生成的全息圖像安全性差的問題，本文采用一種簡(jiǎn)單、有效的雙水印預(yù)處理加密方法。首先，利用隨機(jī)分散矩陣將兩幅全息圖像隨機(jī)分散合成一個(gè)圖像。然后，將圖像分割成兩個(gè)大小為M pixel×N pixel的圖像，作為加密后的待嵌入全息水印。具體步驟如下：

（1）利用共軛對(duì)稱延拓計(jì)算全息生成兩個(gè)大小為M pixel×N pixel的全息水印圖像。

（2）隨機(jī)生成一個(gè)大小為M×2 N，只包含0，1的矩陣作為隨機(jī)分散矩陣，其中0，1 數(shù)量各占一半。

（3）當(dāng)隨機(jī)分散矩陣中某一位置的值為0時(shí)，在合成圖像對(duì)應(yīng)位置記錄第一個(gè)全息圖像的像素值，反之，記錄第二個(gè)全息圖像的像素值，依次類推。即兩個(gè)全息圖像被隨機(jī)分散到一個(gè)合成圖像中。

（4）將合成全息圖像分割成兩個(gè)大小為M pixel×N pixel的圖像，即為加密后的待嵌入全息水印。

（5）記錄隨機(jī)分散矩陣，作為解密密鑰。

3.2 全息水印的嵌入

首先，將原始圖像進(jìn)行二級(jí)Contourlet變換，將中頻信息進(jìn)行8方向分解，得到8個(gè)中頻子帶圖像。其次，將兩個(gè)加密后的待嵌入全息水印自適應(yīng)的替換原始圖像的中頻子帶，實(shí)現(xiàn)水印的嵌入。最后，做二級(jí)Contourlet逆變換，得到含水印圖像。具體步驟如下：

（1）將4 M pixel×8 N pixel大小原始圖像進(jìn)行二級(jí)Contourlet變換，兩級(jí)分解分別為4方向和8方向，得到8個(gè)大小為M pixel×N pixel（或N pixel×M pixel）的中頻子帶圖像。

（2）將雙全息水印自適應(yīng)的替換原始圖像中頻信息，以實(shí)現(xiàn)水印信息的嵌入。為縮小所替換中頻信息的范圍，分別計(jì)算待嵌入水印與8個(gè)中頻子帶的PSNR 值，并記錄由大到小前3 個(gè)PSNR 值對(duì)應(yīng)的中頻系數(shù)編號(hào)。

（3）將水印信息按所得編號(hào)隨機(jī)替換兩個(gè)中頻子帶圖像，實(shí)現(xiàn)水印的嵌入。具體公式如下：

式中：LLi為原始圖像的中頻子帶圖像，i為步驟2 中記錄的編號(hào)。watermarkj為全息水印，j＝1，2。當(dāng)j不同時(shí)，i的編號(hào)也有可能相同。k為嵌入系數(shù)。

（4）做二級(jí)Contourlet逆變換，得到含水印圖像。計(jì)算其與原始圖像的PSNR，并記錄全息水印替換的中頻子帶編號(hào)。

（5）重復(fù)步驟3～4，得到PSNR 值最大時(shí)對(duì)應(yīng)的中頻子帶編號(hào)。對(duì)應(yīng)的含水印圖像即為自適應(yīng)嵌入水印得到的含水印圖像。

（6）記錄所替換中頻子帶的編號(hào)，作為提取水印的密鑰。

3.3 全息水印的提取

水印提取過(guò)程為水印嵌入過(guò)程的逆過(guò)程，本文的水印檢測(cè)無(wú)需原始圖像，為盲水印。具體步驟如下：

（1）將待檢測(cè)圖像進(jìn)行二級(jí)Contourlet變換，兩級(jí)分解分別為4方向和8方向，得到8個(gè)大小為M pixel×N pixel（或N pixel×M pixel）的中頻子帶。

（2）根據(jù)存儲(chǔ)的中頻子帶的編號(hào)提取出水印信息并合并。

（3）利用存儲(chǔ)的隨機(jī)分散矩陣，恢復(fù)出兩個(gè)原始全息水印圖像。

（4）對(duì)恢復(fù)的全息圖像做逆傅立葉變換得到水印信息。

4 仿真實(shí)驗(yàn)

為了驗(yàn)證本算法的有效性，實(shí)驗(yàn)采用1 024 pixel×1 024pixel大小的8bit灰度lena圖像作為實(shí)驗(yàn)的原始圖像，如圖4（a）所示。原始水印圖像采用127pixel×127pixel大小的二值圖像，如圖4（b）和4（c）所示。利用共軛對(duì)稱延拓計(jì)算全息生成兩幅全息水印，用3.1節(jié)的預(yù)處理方式，將兩幅全息水印加密，得到如圖5（a）和5（b）所示的加密圖像。利用正確的分散矩陣對(duì)加密全息圖像解密，并做逆傅立葉變換得到原始水印，如圖5（c）和5（d）所示。

圖4 原始圖像和水印圖像Fig.4 Original image and watermark images

圖5 加密全息圖及其再現(xiàn)Fig.5 Encryption and reconstruction of holograms

本文用峰值信噪比（PSNR）來(lái)客觀評(píng)價(jià)兩幅圖像的相似程度。利用3.2節(jié)的方法將雙水印嵌入到原始圖像中。水印嵌入系數(shù)k 取為0.4，由3.2節(jié)步驟2計(jì)算兩個(gè)水印分別與8個(gè)方向子帶的PSNR 值，如表1所示。得到由大到小前3個(gè)PSNR 值對(duì)應(yīng)的中頻系數(shù)編號(hào)均為4、6、8。將雙水印對(duì)應(yīng)替換不同的中頻系數(shù)（所要替換中頻系數(shù)編號(hào)均為4、6、8），不同替換組合得到的含水印圖像與原始圖像PSNR 值如表2所示。

由表2知，當(dāng)兩個(gè)全息水印的嵌入強(qiáng)度確定時(shí)，替換不同中頻系數(shù)，得到的PSNR 值從37.402 5變化到38.385 7，變化范圍較大。需要說(shuō)明的是，如果替換與水印圖像相似度更小的中頻子帶，PSNR 值的變化范圍更大。本文正是利用Contourlet變換提供的更高水平的多方向性設(shè)計(jì)隱藏水印，使得在嵌入強(qiáng)度確定的情況下，找到最優(yōu)的嵌入位置，使得PSNR 提高。達(dá)到了在不犧牲相似度的前提下，提高水印嵌入強(qiáng)度的目的，嵌入強(qiáng)度的提高使其抵抗攻擊能力也進(jìn)一步增強(qiáng)。

本文實(shí)驗(yàn)中，兩個(gè)水印所替換的中頻系數(shù)其編號(hào)分別為6 和8，含水印圖像和原始圖像的PSNR 值為38.385 7，客觀說(shuō)明了兩幅圖像的相似度很高。圖6所示為含水印圖像和用正確的密鑰提取出的水印圖像。從圖中可以看出提取的水印信息清晰。本文仿真實(shí)驗(yàn)中，提取出的水印圖像均未做任何處理。為了說(shuō)明本算法除了可以抵抗打印掃描攻擊外，還具有更好的抵抗一般圖像處理攻擊的魯棒性，本算法將與文獻(xiàn)［15］算法進(jìn)行對(duì)比實(shí)驗(yàn)。需要說(shuō)明的是，文獻(xiàn)［15］算法是將一個(gè)水印信息嵌入在原始圖像中，且得到的含水印圖像與原始圖像的PSNR 值為37.5dB，小于本算法的PSNR。為保證嵌入的水印圖像相同，對(duì)比實(shí)驗(yàn)中采用的水印圖像為雙水印的拼接圖像。為了對(duì)比實(shí)驗(yàn)直觀，從文獻(xiàn)［15］提取的水印再裁切成兩個(gè)水印進(jìn)行效果對(duì)比。

表1 兩個(gè)全息圖與8個(gè)方向子帶的PSNR 值Tab.1 PSNR between two holograms and eight directional subbands

表2 原始圖像與替換不同方向子帶所得水印圖像的PSNR 值Tab.2 PSNR between original image and watermarked images

圖6 含水印圖像和提取的水印圖像Fig.6 Watermarked image and extracted watermark images

4.1 雙水印預(yù)處理安全性測(cè)試

由于選用的全息圖生成方式幾乎沒有安全性，故采取一種簡(jiǎn)單的雙水印預(yù)處理方法提高其安全性。加密后的全息圖像如圖5（a）和（b）所示，用錯(cuò)誤的隨機(jī)分散矩陣恢復(fù)的圖像，經(jīng)過(guò)逆傅立葉變換得到的圖像如圖7所示，均為無(wú)意義的圖像。

由實(shí)驗(yàn)可知，本文使用的雙水印預(yù)處理方法簡(jiǎn)單、安全性好。因?yàn)閮煞既D的信息被隨機(jī)分散，只占每個(gè)加密全息圖的一部分；加密圖像所含的信息是兩幅原始全息圖的部分信息隨機(jī)組合，這樣會(huì)對(duì)全息圖的逆傅立葉變換解密產(chǎn)生相互干擾。所以，即使正確提取出嵌入的水印信息，如果沒有密鑰矩陣，也無(wú)法恢復(fù)出原始的雙水印信息。

圖7 錯(cuò)誤的隨機(jī)分散矩陣恢復(fù)的圖像Fig.7 Reconstruction images using wrong random scatter matrix

4.2 JPEG 壓縮攻擊實(shí)驗(yàn)

JPEG 圖像壓縮是最常見的圖像處理攻擊之一，所以能否有效抵抗JPEG 壓縮攻擊也是評(píng)價(jià)一個(gè)算法魯棒性的重要指標(biāo)。如圖8所示，標(biāo)號(hào)a～f和m～r的圖像是本算法提取出來(lái)的水印圖像。對(duì)應(yīng)的壓縮質(zhì)量因子分別為20、30、40、50、60、80。標(biāo)號(hào)g～l和s～x圖像為文獻(xiàn)［15］算法提取的水印圖像。

圖8 JPEG 壓縮攻擊實(shí)驗(yàn)Fig.8 Experiment for attack of JPEG compression

由圖8 可以看出，文獻(xiàn)［15］算法在抵抗JPEG 壓縮攻擊方面存在很大缺陷，只有當(dāng)壓縮因子為50及以上時(shí)，提取的水印信息才有一定的主觀辨識(shí)性，但是也沒有本算法提取的水印清晰。由于Contourlet域的中低頻具有很強(qiáng)的穩(wěn)定性，故本算法可以有效抵抗JPEG 壓縮攻擊，當(dāng)壓縮因子在30 及以上時(shí)，提取出的水印信息均很清晰，視覺辨識(shí)度很高，當(dāng)壓縮因子為20時(shí)，提取的水印信息也有一定的視覺辨識(shí)度。說(shuō)明本算法抵抗JPEG 壓縮攻擊的魯棒性強(qiáng)。

4.3 噪聲、低通濾波攻擊

水印算法對(duì)噪聲攻擊和對(duì)濾波攻擊的魯棒性也是評(píng)價(jià)算法優(yōu)劣的重要指標(biāo)。為了檢驗(yàn)本文提出算法對(duì)噪聲攻擊和濾波攻擊的魯棒性強(qiáng)弱，本文選用的噪聲攻擊有加性噪聲：高斯噪聲和椒鹽噪聲，也有乘性噪聲：speckle噪聲；低通濾波攻擊有高斯濾波，中值濾波和均值濾波。選用的噪聲攻擊其均值都為0，高斯噪聲的方差為0.002，椒鹽噪聲的方差為0.01，speckle噪聲的方差為0.03；高斯濾波選用7×7大小模版，中值濾波和均值濾波均選用3×3大小模版。如圖9所示，標(biāo)號(hào)a～f和m～r的圖像是本算法從受到上述攻擊的對(duì)應(yīng)含水印圖像中提取出來(lái)的水印圖像。標(biāo)號(hào)g～l和s～x圖像為文獻(xiàn)［15］算法對(duì)應(yīng)提取的水印圖像。

圖9 噪聲和濾波攻擊實(shí)驗(yàn)Fig.9 Experiment for attack of noise and low－pass filtering

由圖9可知，本文算法和文獻(xiàn)［15］算法對(duì)加性噪聲攻擊都有很好的抵抗能力，但本算法抵抗spckle乘性噪聲的能力更強(qiáng)，提取的水印更加清晰；本文算法抵抗濾波攻擊的能力強(qiáng)于文獻(xiàn)［15］算法，提取的水印背景噪聲較少。這是因?yàn)槲墨I(xiàn)［15］算法，將水印隱藏在圖像小波變換的對(duì)角線中頻系數(shù)中，當(dāng)受到濾波攻擊時(shí)，小波中頻信息的穩(wěn)定較差而受到破壞，導(dǎo)致提取的水印背景有很明顯的、不均勻的對(duì)角線噪聲。實(shí)驗(yàn)說(shuō)明本算法可以有效抵抗噪聲和低通濾波攻擊，提取的水印信息均清晰、辨識(shí)度高。

4.4 裁切攻擊

由于全息圖的不可撕毀性，該算法對(duì)裁切攻擊也有很強(qiáng)的抵抗能力。本文將含水印圖像分別進(jìn)行裁切10%和30%處理。如圖10所示，編號(hào)（a）和（b）為從裁切10%的含水印圖像中提取出的水印信息，編號(hào)（c）和（d）為從裁切30%的含水印圖像中提取出的水印信息。由圖可以看出，提取出的水印信息清晰，本算法可以有效抵抗裁切攻擊。

圖10 剪切攻擊實(shí)驗(yàn)Fig.10 Experiment for attack of cropping

4.5 打印掃描攻擊

打印掃描實(shí)驗(yàn)中采用的打印機(jī)型號(hào)為HPZ3200彩色打印機(jī)，打印分辨率為150dpi；采用的掃描儀型號(hào)為EPSON V500，掃描分辨率為600dpi。本算法生成的圖像經(jīng)打印掃描后得到電子圖像，從中提取出的水印信息如圖11（a）和圖11（b）所示。從圖中可以看出，本算法得到的水印信息清晰，視覺辨識(shí)度高，說(shuō)明本算法除了有很強(qiáng)的抵抗一般圖像處理攻擊的魯棒性外，還可以有效抵抗打印掃描攻擊。

圖11 打印掃描攻擊實(shí)驗(yàn)Fig.11 Experiment for print－scan attack

圖12 打印掃描和JPEG 壓縮混合攻擊實(shí)驗(yàn)Fig.12 Experiment for hybrid attack of print－scan and JPEG compression

JPEG壓縮是最常見的圖像處理攻擊，圖像很可能經(jīng)過(guò)壓縮后用于打印輸出，所以算法對(duì)JPEG壓縮和打印掃描混合攻擊的魯棒性也是一個(gè)重要的指標(biāo)。圖12（a）和圖12（b）所示為本算法得到的含水印圖像，經(jīng)過(guò)壓縮因子為30的JPEG 壓縮攻擊和打印掃描攻擊后提取出的水印信息。文獻(xiàn)［15］生成圖像經(jīng)過(guò)相同操作得到的水印如圖12（c）和圖12（d）所示。從圖中可以看出，文獻(xiàn)［15］算法提取出的水印信息已經(jīng)無(wú)法辨認(rèn)，說(shuō)明該算法無(wú)法有效抵抗JPEG 壓縮和打印掃描混合攻擊，而本文提出的算法提取出的水印信息依然清晰可辨，說(shuō)明本算法可以有效抵抗這種常見的組合攻擊，魯棒性很強(qiáng)。

5 結(jié) 論

基于CGH 和Contourlet變換提出了一種基于計(jì)算全息的強(qiáng)魯棒性自適應(yīng)水印算法。利用Contourlet變換冗余信息量大有利于水印隱藏、在圖像處理中的獨(dú)特優(yōu)勢(shì)——多方向性和各向異性，以及其中、低頻信息更加穩(wěn)定這些特點(diǎn)，將雙全息水印自適應(yīng)的隱藏到圖像Contourlet變換的中頻部分，使得在PSNR 值保持相對(duì)較高的情況下，水印嵌入的強(qiáng)度增強(qiáng)，提高了算法抵抗攻擊的能力。實(shí)驗(yàn)結(jié)果表明：利用本算法得到的含兩個(gè)全息水印的圖像與原始圖像的PSNR 值達(dá)到了38.385 7，說(shuō)明本算法有很強(qiáng)的不可見性。本算法不僅可以從打印掃描攻擊后的圖像中提取出清新可辨認(rèn)的水印信息，也對(duì)JPEG 壓縮攻擊的魯棒性強(qiáng)，壓縮因子達(dá)到30時(shí)，提取的水印依然可以清晰辨認(rèn)。同時(shí)對(duì)剪切、噪聲和低通濾波攻擊均有很強(qiáng)的抵抗能力。并且本算法可以抵抗常見的JPEG 壓縮和打印掃描的混合攻擊，提取出的水印信息清晰可辨認(rèn)。在水印提取時(shí)，不需要借助原始圖像，屬于盲水印算法。本文所提出的數(shù)字水印算法不僅可以保護(hù)數(shù)字圖像，也可以對(duì)印刷圖像進(jìn)行有效的版權(quán)保護(hù)。

［1］ 李新娥，班皓，任建岳，等.一種多傳感器圖像融合與數(shù)字水印技術(shù)［J］.液晶與顯示，2014，29（2）：286－292.Li X E，Ban H，Ren J Y，et al.Method of multisensory image fusion and digital watermarking［J］.Chinese Journal of Liquid Crystals and Displays，2014，29（2）：286－292.（in Chinese）

［2］ 申靜.一種利用人眼視覺特性的數(shù)字圖像水印技術(shù)［J］.液晶與顯示，2011，26（2）：234－240.Shen J.Digital image watermarking technology based on human vision characteristics［J］.Chinese Journal of Liquid Crystals and Displays，2011，26（2）：234－240.（in Chinese）

［3］ Papakostas G A，Tsougenis E D，Koulouriotis D E.Moment－based local image watermarking via genetic optimization［J］.Applied Mathematics and Computation，2014，227：222－236.

［4］ Ishtiaq M，Jaffar M A，Choi T S.Optimal sub－band and strength selection for blind watermarking in wavelet domain［J］.Imaging Science Journal，2014，62（3）：171－177.

［5］ Chen C H，Tang Y L，Wang C P，et al.A robust watermarking algorithm based on salient image features［J］.Optik，2014，125（3）：1134－1140.

［6］ Takai N，Mifune Y.Digital watermarking by a holographic technique［J］.Applied Optics，2002，41（5）：865－873.

［7］ 馬建設(shè)，任振波，蘇萍，等.多重計(jì)算全息水印技術(shù)［J］.光學(xué)精密工程，2013，21（8）：2111－2120.Ma J S，Ren Z B，Su P，et al.Multiple watermarking technique based on computer－generated hologram［J］.Optics and Precision Engineering，2013，21（8）：2111－2120.（in Chinese）

［8］ Martinez C，Laulagnet F，Lemonnier O.Gray tone image watermarking with complementary computer generated holography［J］.Optics Express，2013，21（13）：15438－15451.

［9］ Li J Z.Robust image watermarking scheme against geometric attacks using a computer－generated hologram ［J］.Applied Optics，2010，49（32）：6302－6312.

［10］ Giuseppe S S，Michele D S.Holographic watermarking for authentication of cut images［J］.Optics and Lasers in Engineering，2011，49（12）：1447－1455.

［11］ 孫劉杰，莊松林.雙隨機(jī)相位加密同軸傅里葉全息水印防偽技術(shù)［J］.光學(xué)學(xué)報(bào)，2007，27（4）：621－624.Sun L J，Zhuang S L.Forgery prevention based on inline fourier holographic watermark with double random phase encrypt ion［J］.Acta Optica Sinica，2007，27（4）：621－624.（in Chinese）

［12］ 孫劉杰，莊松林.加密同軸全息數(shù)字水?。跩］.光學(xué)精密工程，2007，15（1）：131－137.Sun L J，Zhuang S L.Digital watermarking of encrypted in－line holography［J］.Optics and Precision Engineering，2007，15（1）：131－137.（in Chinese）

［13］ 孫云峰，翟宏琛，楊曉蘋，等.傅立葉計(jì)算全息數(shù)字水印在彩色圖像印刷防偽中的應(yīng)用［J］.光電子·激光，2008，19（7）：952－955.SunY F，Zhai H C，Yang X P，et al.Application of Fourier CGH digital watermarking technique in color image Forgery－Prevention printing［J］.Journal of Optoelectronics·Laser，2008，19（7）：952－955.（in Chinese）

［14］ 黃素娟，王朔中，于瀛潔.共軛對(duì)稱延拓傅里葉計(jì)算全息［J］.物理學(xué)報(bào)，2009，58（2）：952－958.Huang S J，Wang S Z，Yu Y J.Computer generated holography based on Fourier transform using conjugate symmetric extension［J］.Acta Physica Sinica，2009，58（2）：952－958.（in Chinese）

［15］ 黃素娟，王杜瑤，劉曉靜，等.基于小波變換的硬拷貝全息水?。跩］.光電子·激光，2011，22（9）：1415－1420.Huang S J，Wang D Y，Liu X J，et al.Hardcopy hologram watermarking based on discrete wavelet transform［J］.Journal of Optoelectronics·Laser，2011，22（9）：1415－1420.（in Chinese）

［16］ Do M N，Vetterli M.Contourlets：A new directional multiresolution image representation［C］.Signals，Systems and Computers，2002，1：497－501.

［17］ Ranjbar S，Zargari F，Ghanbari M.A highly robust two－stage Contourlet－based digital image watermarking method［J］.Signal Processing：Image Communication，2013，28（10）：1526－1536.

［18］ Qin，C，Wen，X B.A novel digital watermarking algorithm in Contourlet domain［J］.Journal of Information and Computational Science，2014，11（2）：519－526.

［19］ Mo J，Ma Z F，Huang Q L，et al.An adaptive watermarking scheme using SVD in Contourlet domain［J］.Advances in Information Sciences and Service Sciences，2012，4（15）：221－232.

［20］ Jiang A P，Yao D.A Contourlet watermarking algorithm based on fractal dimension and one－dimensional logistic chaotic encryption［J］.International Journal of Advancements in Computing Technology，2012，4（19）：317－323.