潘瑞冬 廖琪男

摘要：為了有效實(shí)現(xiàn)對(duì)數(shù)字圖像的版權(quán)保護(hù)和提高算法的抗攻擊性能，提出了一種基于DCT直流系數(shù)和模運(yùn)算的數(shù)字圖像零水印算法。首先根據(jù)灰度水印的大小對(duì)載體圖像分塊，直接在空域中提取每一子塊的DCT直流系數(shù)作為載體圖像的特征圖像；然后用灰度水印與特征圖像相加的模運(yùn)算結(jié)果作為零水印。理論分析和實(shí)驗(yàn)結(jié)果表明，該文算法與基于DWT和XOR操作的同類(lèi)算法相比，不僅具備空域算法計(jì)算復(fù)雜度低的優(yōu)點(diǎn)，還具有更強(qiáng)得多的抗JPEG壓縮、抗圖像濾波處理和抗噪聲攻擊等能力。

關(guān)鍵詞：零水?。浑x散余弦變換；直流系數(shù)；模函數(shù)；空域

中圖分類(lèi)號(hào)：TP751 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1009-3044（2015）11-0210-03

1 引言

隨著互聯(lián)網(wǎng)和多媒體技術(shù)的飛速發(fā)展和廣泛應(yīng)用，數(shù)字水印技術(shù)也因此成為信息網(wǎng)絡(luò)安全技術(shù)研究的熱點(diǎn)。不可感知性和魯棒性是數(shù)字水印方案的最重要性質(zhì)。傳統(tǒng)的數(shù)字水印技術(shù)都要對(duì)載體圖像空域或變換域的信息進(jìn)行一定的修改， 從而造成載體圖像質(zhì)量下降。為了解決水印的不可感知性和魯棒性之間的矛盾，零水印方案被提了出來(lái)。所謂零水印是指利用載體圖像的重要特征來(lái)構(gòu)建水印，而不對(duì)載體圖像作任何修改的水印技術(shù)。

零水印算法可以分為兩大類(lèi)：空域方法[1 ，3]和變換域方法[2，4，5]。本文在深入研究已有零水印算法的基礎(chǔ)上，深入研究模運(yùn)算特性，提出簡(jiǎn)單而魯棒性很強(qiáng)的基于DCT直流系數(shù)和模操作的灰度零水印算法。

2 相關(guān)技術(shù)

2.1離散余弦變換

離散余弦變換（discrete cosine transform， DCT），是數(shù)字圖像處理以及信號(hào)處理常用的一種正交變換。設(shè)大小為M×N的圖像或圖像塊為f（x，y）， 其DCT定義為

[F（u，v）=x=0M-1y=0N-1c（u）c（v）f（x，y）cosπ（2x+1）u2Mcosπ（2Y+1）v2N] （1）

式中，u=0，1，2， …，M-1， v=0，1，2， …，N-1。

[c（u）= 1M u=0 2M u=1，2，…，M-1] （2）

[c（v）= 1N v=0 2N v=1，2，…，N-1] （3）

由式（1）可知，當(dāng)u=0， v=0， 即DCT直流（direct current， DC）系數(shù)可表示為

[F（0，0）=1MNx=0M-1y=0N-1f（x，y）] （4）

由式（4）可以看出，DC系數(shù)可以代表圖像或圖像塊平均亮度，可以直接在空域中通過(guò)簡(jiǎn)單的數(shù)學(xué)運(yùn)算求得，而無(wú)須經(jīng)過(guò)復(fù)雜的DCT得到。

2.2 Logistic混沌系統(tǒng)

Logisti映射是一個(gè)非常簡(jiǎn)單而在數(shù)字水印技術(shù)中最常用的一維混沌系統(tǒng)，其定義為

[xn+1=μxn（1-xn）] （5）

其中，0≤μ≤4稱(chēng)為分枝參數(shù)。當(dāng)xn∈（0，1），且3.5699456<μ≤4 時(shí)，Logistic 映射工作處于混沌狀態(tài)，此時(shí)生成的序列{xn ， n=0，1，2，…}具有良好的隨機(jī)統(tǒng)計(jì)特性，是較理想的混沌偽隨機(jī)序列。

3 本文零水印算法

設(shè)載體圖像和水印圖像均為256級(jí)灰度圖像，載體圖像I大小為M×N，水印圖像W大小為m×n，本文零水印算法使用在密鑰key=（μ，x0）下由式（5）產(chǎn)生的Logistic混沌序列為

[Xi|i=1，2，3，…，M×N] （6）

序列X的前m×n項(xiàng)序列用于置亂水印，后m×n項(xiàng)序列用于加密零水印信息。

3.1 水印置亂

為提高水印算法的魯棒性和安全性，對(duì)原水印W進(jìn)行混沌置亂加密。取式（6）的混沌序列X的前m×n項(xiàng)序列組成序列{aXi|i=1，2，…，m×n}，aX排序得到有序序列{tXi|i=1，2，…，m×n}，由tX在aX中的位置索引為序列值生成服從[1，m×n]的混沌隨機(jī)序列{sXi|i=1，2，…，m×n}，按光柵掃描順序，該序列值對(duì)應(yīng)W的像素位置（row，col）為

[row（i）=sX（i）ncol（i）=sX（i）-（row（i）-1）n] （7）

W從（row，col）變換到（i，j）得到置亂水印Wc為

[Wc（i，j）=W（row（（i-1）×n+j），col（（i-1）×n+j））] （8）

式中，i=1，2，…， m，j=1，2，…，n。

3.2 零水印構(gòu)建算法

3.2.1圖像特征信息的提取

根據(jù)水印圖像W大小把載體圖像劃分為m×n不重疊子塊，子塊大小為h×w， 則

[ h=Mm w=Nn] （9）

取子塊集B的每一子塊的DCT直流系數(shù)構(gòu)成載體圖像I的特征C，特征矩陣C根據(jù)式（4）直接從空域中計(jì)算，計(jì)算公式為

[C（u，v）=1hwi=1hj=1wB（u，v）（i，j）] （10）

式中，u=1，2，…，m，v=1，2，…，n。

3.2.2 零水印生成

由載體圖像特征C（i，j）與置亂后的水印像素Wc（i，j）通過(guò)模運(yùn)算構(gòu)建零水印信息F為

[F（i，j）=mod（C（i，j）+Wc（i，j），256）] （11）

式中，i=1，2，…，m，j=1，2，…，n。

3.2.3 零水印加密

為了進(jìn)一步提高零水印在注冊(cè)中心的安全性，用XOR操作對(duì)其加密。取混沌序列{Xi|i=1，2，…， M×N}的后m×n項(xiàng)組成序列{bXi|i=1，2，…， m×n}作為加密序列，將其變換為8位無(wú)符號(hào)整數(shù)，計(jì)算公式為

[bX（i）=round（bX（i）×10k）mod256] （12）

式中，round（）為取最近整數(shù)函數(shù)，i=1，2，…， m×n，k≥4。將bX（i）序列組成m×n矩陣sY為

[sY（i，j）=bX（（i-1）n+j）] （13）

式中，i=1，2，…，m，j=1，2，…，n。

利用XOR操作把零水印信息加密為

[F（i，j）=F（i，j）⊕sY（i，j）] （14）

經(jīng)過(guò)式（14）加密后的零水印的像素統(tǒng)計(jì)特性將被徹底改變。

3.3 水印檢測(cè)算法

水印的檢測(cè)過(guò)程是水印構(gòu)建過(guò)程的逆過(guò)程。在待檢測(cè)的載體圖像I提取圖像特征信息，并與在IPR水印數(shù)據(jù)庫(kù)中注冊(cè)的零水印信息進(jìn)行逆運(yùn)算得到檢測(cè)水印W′。

對(duì)于有特定含義的水印圖像而言，一般所提取的水印效果可由原水印W和檢出的水印W′進(jìn)行歸一化相關(guān)性計(jì)算來(lái)輔助度量，或作為魯棒性比較。歸一化相關(guān)計(jì)算公式為

[NC=i=1mj=1nW（i，j）W′（i，j）i=1mj=1nW2（i，j）i=1mj=1nW′2（i，j）] （15）

實(shí)驗(yàn)證明，對(duì)于灰度水印，歸一化系數(shù)NC不能較好反映檢測(cè)水印和原水印的相似程度，故本文用峰值信噪比（PSNR）來(lái)衡量檢測(cè)水印的視覺(jué)質(zhì)量及評(píng)價(jià)算法性能比較，其定義為

[PSNR=10×log10（2552/MSE）] （16）

式中，MSE是原水印圖像W與檢測(cè)水印圖像W′之間的均方差，計(jì)算公式為

[] [MSE=1m×ni=1mj=1n（W（i，j）-W′（i，j））2] （17）

式中，m和n分別為水印圖像的高和寬。

4 實(shí)驗(yàn)結(jié)果與分析

已通過(guò)大量實(shí)驗(yàn)證明本文算法的有效性和可行性。限于篇幅，在此僅給出代表性的實(shí)驗(yàn)結(jié)果。

實(shí)驗(yàn)平臺(tái)為Matlab2012b，載體圖像為8位 Lena 512×512標(biāo)準(zhǔn)灰度測(cè)試圖像，水印為8位logo 128×128灰度圖像，如圖1所示。logo為高亮圖像，式（18）和（19）中，a=0，b=220。Logistic系統(tǒng)密鑰key=（μ，x0）=（3.9，0.32579），式（12）中k=7。

（a）Lena圖像 （b）水印圖像 （a）零水印信息 （b）直方圖

圖1 實(shí)驗(yàn)圖像 圖2 零水印信息加密效果

4.1 零水印信息與安全性

零水印信息圖像及其灰度直方圖如圖2所示。

水印信息的安全性取決于密鑰和水印信息像素特征值統(tǒng)計(jì)。密鑰的安全性由Logistic系統(tǒng)保證，Logistic系統(tǒng)的密鑰key不僅密鑰空間足夠大，而且很敏感。從圖2實(shí)驗(yàn)效果表明不僅水印信息圖像加密效果好，而且其灰度直方圖在[0，255]范圍內(nèi)分布均勻，可以抵御來(lái)自像素特征值統(tǒng)計(jì)的攻擊。

4.2 零水印算法魯棒性實(shí)驗(yàn)

本實(shí)驗(yàn)為本文算法與文獻(xiàn)[9-10]為代表的基于DWT和XOR操作的同類(lèi)算法的魯棒性對(duì)比實(shí)驗(yàn)。文獻(xiàn)[9-10]算法在DWT二層低頻子帶構(gòu)建零水印。實(shí)驗(yàn)項(xiàng)目與數(shù)據(jù)如表1所示。

表1 本文算法與基于DWT算法的魯棒性實(shí)驗(yàn)對(duì)比結(jié)果

[實(shí)驗(yàn)項(xiàng)目＼&文獻(xiàn)[9-10]算法＼&本文算法＼&NC＼&PSNR（dB）＼&NC＼&PSNR（dB）＼&JPEG壓縮

質(zhì)量因子75＼&0.9921＼&19.5599＼&0.9999＼&40.8317＼&中值濾波3×3＼&0.9875＼&17.5778＼&0.9994＼&30.7544＼&增加對(duì)比度

灰度[0.1，0.9]到[0，1]＼&0.8716＼&7.5646＼&0.8988＼&8.7105＼&降低對(duì)比度

灰度[0，1]到[0.1，0.9]＼&0.8964＼&8.4354＼&0.9696＼&13.7843＼&圖像增亮

灰度[0，1]到[0.1，1]＼&0.8321

檢不出＼&6.3693

檢不出＼&0.9834＼&11.8441＼&圖像變暗

灰度[0，1]到[0，0.9]＼&0.8454＼&6.7098＼&0.7341＼&4.8707＼&高斯噪聲

均值0，方差0.01＼&0.9233＼&9.7279＼&0.9718＼&14.0959＼&裁剪四周25＼&0.9631＼&12.9594＼&0.9654＼&13.2491＼&]

在載體圖像沒(méi)有受到任何攻擊時(shí)，可以檢出無(wú)損水印圖像。與表2相對(duì)應(yīng)，經(jīng)JPEG壓縮攻擊和經(jīng)高斯噪聲攻擊后的檢測(cè)水印視覺(jué)效果分別如圖3和圖4所示。限于篇幅，其他魯棒性實(shí)驗(yàn)效果可通過(guò)表內(nèi)數(shù)據(jù)來(lái)間接比較。

實(shí)驗(yàn)表明，NC不能準(zhǔn)確表示檢測(cè)灰度水印的相似性，如本實(shí)驗(yàn)的抗JPEG壓縮攻擊實(shí)驗(yàn)，本文算法的NC與文獻(xiàn)[9-10]的NC相近，但本文算法檢出水印的視覺(jué)質(zhì)量和PSNR遠(yuǎn)遠(yuǎn)超過(guò)文獻(xiàn)[9-10]算法。

（a）文獻(xiàn)[9-10]算法， NC=0.9921 （b）本文算法， NC=0.9999

圖3 JPEG壓縮（質(zhì)量系數(shù)75%）攻擊檢出水印

（a）文獻(xiàn)[9-10]算法， NC=0.9233 （b）本文算法， NC=0.9718

圖4 高斯噪聲攻擊檢出水印

本次實(shí)驗(yàn)的圖像變暗處理，表2數(shù)據(jù)顯示檢出NC和PSNR都很小，且本文算法的NC和PSNR比文獻(xiàn)[9-10]算法的小，但本文算法檢出水印質(zhì)量卻明顯好很多（見(jiàn)圖5）。

由表2實(shí)驗(yàn)數(shù)據(jù)和檢出水印視覺(jué)質(zhì)量（如與表2對(duì)應(yīng)的圖3-5實(shí)驗(yàn)結(jié)果）表明，本文算法的抗JPEG壓縮、抗圖像濾波處理，抗高斯噪聲攻擊、抗圖像對(duì)比度、抗亮度處理和抗裁剪能力等能力遠(yuǎn)遠(yuǎn)強(qiáng)于文獻(xiàn)[9-10]同類(lèi)算法。檢出水印的視覺(jué)質(zhì)量和PSNR遠(yuǎn)遠(yuǎn)優(yōu)于文獻(xiàn)[9-10]算法。

（a）文獻(xiàn)[9-10]算法， NC=0.8454 （b）本文算法， NC=0.7341

圖5 圖像變暗攻擊檢出水印

5 結(jié)論

理論分析和實(shí)驗(yàn)結(jié)果表明，本文零水印算法與基于DWT和XOR操作的同類(lèi)算法相比，具有強(qiáng)得多的抗JPEG壓縮、濾波等圖像處理和抗高斯噪聲攻擊等方面的能力，具有更強(qiáng)的抗圖像對(duì)比度和亮度處理、抗裁剪和抗幾何攻擊等方面的能力。

參考文獻(xiàn)：

[1] 溫泉， 孫錟鋒， 王樹(shù)勛， 等. 零水印的概念與應(yīng)用[J]. 電子學(xué)報(bào)， 2003，31（2）：214-216.

[2] 付劍晶，王 珂.基于第一主成分方向穩(wěn)定性的圖像零水印[J]. 中國(guó)圖象圖形學(xué)報(bào)，2012，17（7）：756-769.

[3] 宋偉，侯建軍，李趙紅，黃亮.一種基于logistic混沌系統(tǒng)和奇異值分解的零水印算法[J].物理學(xué)報(bào)，2009，58（7）：4449-4556.

[4] 馬建湖，何甲興.基于小波變換的零水印算法[J].中國(guó)圖象圖形學(xué)報(bào)，2007，12（4）：581-585.

[5] 何冰，王晅. 基于圖像置亂和小波變換的零水印算法[J].計(jì)算機(jī)工程與設(shè)計(jì)，2009，30（1）：145-147，154.