程 莎，唐向宏，林軍海，陳宏炳

(杭州電子科技大學(xué)通信工程，浙江杭州310018)

0 引言

目前大多數(shù)水印算法存在水印不可見性與魯棒性難以均衡的問題，由于零水印算法［1］沒有改變圖像質(zhì)量，即保證了水印的不可見性，因此該算法研究的重點是如何提高其魯棒性。利用高階累積量、最高有效位以及圖像的二值化［1－3］來構(gòu)造零水印的算法，對常規(guī)攻擊(如加噪、濾波、JPEG壓縮)具有較好的魯棒性，但其對幾何攻擊(如旋轉(zhuǎn)、縮放、剪貼)的魯棒性能較差。為了提高對幾何攻擊的魯棒性，文獻4提出了基于小波提升與Sobel邊緣及膨脹的零水印算法，文獻5則利用對數(shù)極坐標映射(Log-Polar Mapping，LPM)將笛卡爾坐標系中的旋轉(zhuǎn)、縮放轉(zhuǎn)為循環(huán)平移來構(gòu)造零水印。這些算法均提高了水印抗幾何攻擊的能力，尤其是抗旋轉(zhuǎn)攻擊能力得到了較大改善，但是仍存在諸如對鏡像攻擊和噪聲攻擊的抵抗能力差、運算量太大等問題。本文利用對數(shù)極坐標映射，基于圖像的重要邊緣信息構(gòu)造零水印，以提高水印算法對鏡像、旋轉(zhuǎn)、縮放、加噪以及JPEG壓縮等攻擊的魯棒性，進一步加強了算法的綜合魯棒性能。

1 SUSAN邊緣檢測與對數(shù)極坐標映射特點

1.1 SUSAN邊緣檢測法

最小核值相似區(qū)(Smallest Univalue Segment Assimilating Nucleus，SUSAN)［6］算法具有抗擊噪聲能力強、計算量小等優(yōu)點［7］，其具體實現(xiàn)過程如下。

用一個圓形模板對圖像進行掃描比較，記I(r0)是模板中象素(核)的灰度值;I(r)為板內(nèi)其他任意象素的灰度值;D是以r0為中心的圓形模板;t是灰度差門限;g是區(qū)域門限。則:

比較g與n(r0)，若n(r0)＜g，則檢測到的象素位置r0可認為是個邊緣點。

1.2 對數(shù)極坐標映射

設(shè)圖像在笛卡爾坐標中的任意一點為(x，y)，相應(yīng)的極坐標變換為［8］:

式中，r為極徑，θ為極角，即 x=rcos θ，y=rsin θ。在復(fù)平面中 z=x+iy，則［8］:

令 D=lnz=lnr+i θ，對(x，y)進行尺度 k 倍同時旋轉(zhuǎn) Δ θ得到 z1，則［8］:

比較D和D1發(fā)現(xiàn)，在笛卡爾坐標系中的尺度和旋轉(zhuǎn)變換在對數(shù)極坐標系中分別變成了沿極徑r極角θ的平移運動［8］。

2 水印算法

整個水印系統(tǒng)的原理框圖如圖1所示:

圖1 水印系統(tǒng)原理框圖

2.1 水印嵌入算法

本文算法是基于SUSAN邊緣檢測法和對數(shù)極坐標映射實現(xiàn)的，具體嵌入步驟如下:

(1)為減小計算量，將大小為N×N的原始圖像進行8×8分塊，計算每一塊的灰度均值Mj，其中j=1，2，…，N×N/64，組成大小為(N/8)×(N/8)的近似于原圖像的低尺度圖像;

(2)選取適當?shù)膮?shù)對低尺度圖進行SUSAN邊緣檢測得到邊緣，再對圖像邊緣進行對數(shù)極坐標映射，并將其二值化處理得到二值LPM邊緣圖;

(3)將大小為(N/8)×(N/8)的二值水印圖像進行Arnold置亂得到加密圖像;

(4)將二值LPM邊緣圖和二值加密水印圖進行異或處理，得到最終的注冊圖。

2.2 水印提取算法

水印的提取與水印的構(gòu)造類似，是嵌入過程的逆過程。具體實現(xiàn)步驟如下:

(1)將待測圖像的低尺度圖進行邊緣提取、LPM和二值化得到二值LPM邊緣圖;

(2)將注冊圖像與二值LPM邊緣圖進行異或運算，再利用置亂密鑰得到重構(gòu)水印;

(3)由于本算法采用了對數(shù)極坐標映射，因此需通過循環(huán)遍歷搜索法［5］來檢測水印是否存在，主要步驟如下:首先設(shè)定閾值T，若原始水印和重構(gòu)水印的相似度值大于T，則判定水印存在;若相似度小于T，則將LPM邊緣圖沿著極角的方向平移k個單位，k=1，2，…，θ(θ為極角的采樣個數(shù))，重復(fù)(1)－(2)步。若經(jīng)過θ次遍歷，相似度值仍小于T，則水印不存在。

3 實驗仿真與分析

本文算法采用的原始圖像是512×512的256級灰度的Lena圖，水印圖像為64×64的二值圖像。在進行SUSAN邊緣檢測時，采用37象素模板，g為27，閾值t設(shè)為38。對數(shù)極坐標映射中的極軸極角的采樣個數(shù)均設(shè)為64。無干擾時的實驗結(jié)果如圖2所示。

圖2 無干擾時的實驗結(jié)果

對于水印魯棒性能的評價，采取相似度Sim作為衡量標準［4］。

S指的是原水印和提取水印象素值相同的個數(shù)，N為水印圖像的總象素。受攻擊圖像的質(zhì)量用峰值信噪比(PSNR)來衡量［2］:

式中，MSE指的是原圖像與處理圖像之間的均方誤差。

為驗證本算法對常見攻擊的魯棒性能，實驗中對Lena圖進行各種攻擊來進行測試，并與文獻4進行了性能對比實驗，其實驗結(jié)果如表1和圖3、4所示。

表1 各類攻擊下本文算法與文獻4魯棒性比較

圖3 JPEG壓縮實驗比較

圖4 高斯噪聲實驗比較

由實驗結(jié)果得知本文算法的綜合魯棒性更高，尤其是對于JPEG壓縮和高斯噪聲攻擊具有很好的魯棒性。由于零水印不需要嵌入到載體圖像當中，因此產(chǎn)生了對于內(nèi)容相近的圖像的版權(quán)歸屬問題。為驗證本文算法的是否存在誤判問題，采用4幅與Lena圖片相近的圖像進行測試，如圖5所示，仿真結(jié)果如表2所示。

圖5 Lena近似圖

表2 誤警實驗分析結(jié)果

從表2中可以看出對具有類似內(nèi)容的圖像本文提取水印的平均相似度值小于文獻2，即本算法誤警率相對較低，不容易產(chǎn)生誤判。

4 結(jié)束語

本文利用SUSAN邊緣檢測法和對數(shù)極坐標映射實現(xiàn)了零水印的構(gòu)造和嵌入。本算法的主要優(yōu)勢在于提高了對各類攻擊的抵抗力，特別是對于旋轉(zhuǎn)和鏡像攻擊具有良好的魯棒性。由于本算法是對低尺度近似圖進行對數(shù)極坐標映射，所以大大降低了計算量?？傊?，本算法對各種攻擊都具有良好的魯棒性能，并且對于相似圖像的誤警率較低，具有實際應(yīng)用意義。

［1］溫泉，孫錟鋒，王樹勛.零水印的概念與應(yīng)用［J］.電子學(xué)報，2003，31(2):214－216.

［2］牛萬紅，顏惠琴.一種基于MSB構(gòu)造密鑰的零水印算法［J］.計算機系統(tǒng)應(yīng)用，2009，18(12):66－69.

［3］葉登攀.基于二值化圖像構(gòu)造的零水印版權(quán)保護方案［J］.計算機應(yīng)用研究，2007，24(8):239－241.

［4］范禮，高鐵杠，楊群亭.基于小波提升與Sobel膨脹特征提取的零水印版權(quán)認證方案［J］.計算機應(yīng)用研究，2010，27(2):726－729.

［5］何冰.基于分塊DCT變換的抗旋轉(zhuǎn)，縮放攻擊零水印算法［J］.微計算機應(yīng)用，2010，31(7):1－10.

［6］Smith S M，Brady J M.SUSAN—A new approach to low level image processing［J］.International Journal of Computer Vision，1997，23(1):45 －78.

［7］馬桂珍，房宗良，姚宗中.SUSAN邊緣檢測算法性能分析與比較［J］.現(xiàn)代電子科技，2007，30(8):189－191.

［8］余博，郭雷，趙天云.基于對數(shù)極坐標變換的灰度投影穩(wěn)像算法［J］.計算機應(yīng)用，2008，28(12):3 127－3 128.