劉斌 熊靜雯

摘? 要： 針對張量積小波數(shù)字圖像水印算法在幾何攻擊下魯棒性較差的問題，提出一種將不可分小波、奇異值分解以及Arnold置亂變換相結合的數(shù)字圖像水印算法。該算法首先構造了二維不可分小波濾波器組，其次對水印圖像進行Arnold置亂變換，對原始圖像用不可分小波進行多尺度分析，然后分別對其低頻部分和置亂后的水印進行奇異值分解，最后在兩者的奇異值矩陣中嵌入水印。實驗結果表明，該算法能實現(xiàn)水印的嵌入和準確提取，并且水印的不可見性好，與張量積小波相比，該方法在對濾波、噪聲、剪切、旋轉、加噪的攻擊下，魯棒性較強。

關鍵詞： 圖像水印; 不可分小波; 奇異值分解; 多尺度分析; 水印嵌入; 版權保護

中圖分類號： TN918?34; TP391? ? ? ? ? ? ? ? ? ? 文獻標識碼： A? ? ? ? ? ? ? ? ? ? ?文章編號： 1004?373X（2020）11?0087?05

Non?separable wavelet construction and its application in image watermarking

LIU Bin， XIONG Jingwen

（School of Computer Science and Information Engineering， Hubei University， Wuhan 430062， China）

Abstract： In view of the poor robustness of tensor product wavelet digital image watermarking algorithm under the geometric attack， a digital image watermarking algorithm in combination with the non?separable wavelet， singular value decomposition （SVD） and Arnold scrambling transform is proposed. In the algorithm， a filter bank of two?dimensional non?separable wavelet is constructed first. And then， the watermarking image is subjected to the Arnold scrambling transform， and the non?separable wavelet is used to perform multi?scale analysis on the original image. The low?frequency part and the scrambled watermarking are subjected to the SVD. In the end， the watermarking is embedded in the singular value matrix of the two. The experimental results show that the algorithm can accurately embed and extract watermark， which is of good invisibility performance. In comparison with the tensor product wavelet， the proposed method is of better robustness under the attack of filtering， noise， shearing， rotation and noise adding.

Keywords： image watermarking; non?separable wavelet; SVD; multi?scale analysis; watermarking embedment; copyright protection

0? 引? 言

隨著信息技術的飛速發(fā)展，數(shù)字多媒體與人們的生活已經息息相關，人們能夠很容易地將信息進行復制、存儲、修改與轉發(fā)等，正是因為這個原因使得信息也很容易以非法途徑傳播。因此，針對數(shù)字信息的版權保護等問題，數(shù)字圖像水印技術等被提出。

數(shù)字水印技術是將一些標識信息（如文本、圖像及標識符等）通過某種方式嵌入數(shù)字載體（語音、圖像及視頻等）中，但不影響原載體的使用價值，同時，嵌入的信息不易被人眼所感知，且能夠通過其隱藏的水印信息實現(xiàn)內容認證、版權保護以及信息隱秘傳輸?shù)龋渲杏糜诎鏅啾Ｗo是目前一個非常重要的研究方向。目前，國內外對其進行了大量研究，常見的大多分為兩類：一種是空間域算法，典型的有LSB算法，這種算法雖然實現(xiàn)起來簡單，但是極易受到攻擊，例如，加噪和壓縮都能使水印失效;另一種是變換域算法，典型的有基于離散小波變換（DWT）[1?2]的水印算法和離散余弦變換（DCT）[3?4]的水印算法，前者擁有多尺度分析的特點，可以抵抗加噪、壓縮等攻擊，但因為二維張量積小波是由一維小波通過張量積生成的，存在著方向上信息提取不足的缺陷。后者的主要思想是通過分塊進行DCT將水印嵌入，能抵抗一部分攻擊，但在受到旋轉時水印信息損失較為嚴重。

文獻[1]是將宿主圖像進行張量積小波分解，得到低頻子圖并對其進行分塊，對每一塊做SVD變換得到[U]，[S]，[V]，將置亂后的水印以一定的強度多次重復地嵌入[S]中，并對其再次進行SVD變換得到[U1]，[S1]，[V1]，最后用[U]，[S]，[V]做逆奇異值變換等處理得到含水印圖像，該算法對于剪切、噪聲、濾波、壓縮都有較好的抵抗力，但是對于質量因子為4以下的壓縮攻擊魯棒性較差。文獻[2]采用將彩色RGB的宿主圖像轉換到Lab彩色空間，在其亮度分量[L]上進行張量積小波分解，得到1個低頻和3個高頻分量，分別對其和水印圖像做SVD變換得到相應的[S]，然后用加性原則以不同的強度將水印的奇異值矩陣[S]嵌入到其余的奇異值矩陣中，最后做相應的逆變換得到嵌入水印的圖像。文獻[3]首先把宿主圖像變換為整數(shù)DCT域，將系數(shù)劃分為非重疊的塊，利用非線性混沌映射完成塊的選擇，然后用所選擇的塊構造循環(huán)矩陣，最后用循環(huán)矩陣計算奇異值將水印嵌入。文獻[4]采取對水印圖像做置亂變換，將宿主圖像分為8×8的子塊，對每個子塊進行DCT變換，選擇每一個子塊的中頻系數(shù)，與在Zig?Zag排序中相鄰的6個系數(shù)的平均值做比較，以相應的公式嵌入中頻系數(shù)上，最后對每個子塊做DCT反變換得到水印圖像，該算法魯棒性較好，但是對于旋轉抗攻擊效果較差。文獻[5]選取宿主圖像DWT分解后的LH和HL子圖，然后將其8×8分塊做DCT變換，對AC分量根據(jù)Zig?Zag排序選取系數(shù)放到矩陣[A]中，再對[A]進行SVD奇異值變換得到[S]，在[S]上以加性原則嵌入水印，最后做相應的逆變換得到含水印圖像，該算法矯正了幾何攻擊使其對幾何攻擊有較好的抵抗力。文獻[6]在嵌入水印時結合了DWT，DCT和SVD，該算法在確保水印抗噪聲和濾波攻擊的同時也提高了抗旋轉、剪切攻擊的能力。文獻[7]采取SSVD和DWT結合的數(shù)字水印算法，先對水印圖像進行SSVD以提取其主分量，然后對宿主圖像進行一級DWT變換，取其低頻分塊進行SVD，將水印主成分以加性嵌入原則嵌入其奇異值矩陣的最大值位置，該算法對于剪切攻擊的魯棒性較差。

目前還沒有人將二維不可分小波用于數(shù)字圖像水印，因此結合上述各種技術的優(yōu)點，本文針對于大多數(shù)算法抗旋轉能力較差的特點，提出了將不可分小波、奇異值分解和Arnold置亂變換相結合的水印算法，實驗結果表明，該算法在濾波、壓縮、剪切、噪聲，特別是旋轉攻擊都具有更好的抵抗力。

1? 二維不可分小波變換

目前二維小波變換可以分為兩種：一種是張量積小波變換;另外一種是不可分小波變換[8]。前者由一維小波通過張量積生成，即使二維張量積小波延續(xù)了其部分良好特性，但其對圖像的處理依然存在著一些不足和缺陷。例如，大部分的二維張量積小波基的緊支撐正交的性質是以舍棄其對稱性作為代價換來的，正是因為其不具備對稱的特性，使得圖像分解得到的各個子圖會存在相位失真，并且這種相位失真將無法避免地產生圖像細節(jié)信息上的失真，與此同時，二維張量積小波變換只強調垂直、水平以及對角線方向上的信息，不具備各向同性，導致其不能捕獲各個方向上的奇異性。故本文選用二維小波變換中具備各向同性的不可分小波變換，以捕獲所有方向上的奇異性。

因為Daubechies[9]證實了小波分解的實現(xiàn)是用信號或圖像和離散濾波器進行卷積完成的，所以利用小波對圖像進行處理的核心問題就是小波低通濾波器和高通濾波器的構造。Chen等是將MRA作為基礎，采取的方法是矩陣擴充，提出了一種高維不可分小波濾波器組的一般構造方法[10]，其在理論上來說可以構造出無數(shù)組濾波器。本文在此基礎上給出二維不可分小波濾波器組的具體構造形式。

二維不可分小波變換可以由二維不可分濾波器組實現(xiàn)。設不可分小波變換過程中的抽樣矩陣[S=][2002]，不可分小波濾波器組為[mj（x，y）]，[j=0，1，2，3]，則具有緊支撐、正交性、對稱性的二維四通道不可分濾波器組的形式可構造[11]如下：

[（m0（x，y），m1（x，y），m2（x，y），m3（x，y））=? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 14（1，x，y，xy）j=1k（UjD（x2，y2）UTj）V] （1）

式中：[m0（x，y）]，[m1（x，y）]，[m2（x，y）]，[m3（x，y）]分別為濾波器[H0]，[H1]，[H2]，[H3]所對應的頻域形式。設[x=e-iω1]，[y=e-iω2]，[Ujj=1，2，…，k]為正交對稱陣，[D（x，y）=Diag（1，x，y，xy）]，[V2=（V0，V1，V2，V3）2]為正交陣，[V1，V2，V3]為[4×1]向量，[V0=（1，1，1，1）T]。為得到[6×6]具有對稱性不可分小波濾波器組，取[k=2]，構造：

[V=1-11-1111111-1-11-1-11]，取[P=10100101010-110-10]，[Aj=cos αj-sin αj00sin αjcos αj0000cos? βj-sin? βj00sin? βjcos? βj][j=1，2，…，k]。按[Uj=12×P×Aj×PT]構造[Ujj=1，2，…，k]，可以驗證[Ujαj，βj]是中心對稱正交陣，[V2]是酉矩陣。在此方法中，[αj，βj]的不同取值可以得到不同的具有正交、緊支撐且對稱的濾波器組[H0]，[H1]，[H2]，[H3]，其中，[H0]為低通濾波器，[H1]，[H2]，[H3]為高通濾波器。這里給出[α1=π7]，[β1=π]，[α2=π12]，[β2=π24]時的一組濾波器在時域下的形式：

很顯然，它們是不可分的，即不能被分解為兩個一維小波濾波器的張量積，且具有對稱的性質（關于濾波器中心對稱或者反對稱）、緊支撐的性質，也可以驗證得知它具有正交性，是一組能夠完全重構濾波器。

2? 基于不可分小波變換的水印算法

為了增強數(shù)字圖像水印對各種攻擊的抵抗效果，本文算法先將水印圖像進行Arnold變換[12]和SVD分解[13]。選擇宿主圖像進行多尺度不可分小波后的低頻部分進行SVD分解，在其奇異值矩陣上嵌入水印。通過其嵌入過程的逆過程可以提取出水印。

2.1? 水印的嵌入過程

設待嵌入水印的宿主圖像[I]的大小為[M×N]，原始水印圖像[W]的大小為[m×n]。

1） 對原始水印圖像進行Arnold置亂變換，得到置亂后的水印圖像[W]。

2） 對宿主圖像[I]用本文方法構造的低高通濾波器進行[n1]級不可分小波變換，得到[3n1+1]個子圖。

3） 將第[n1]級低頻子圖[LLn]作為嵌入水印的區(qū)域，對[W]以及[LLn]根據(jù)式（2）和式（3）做奇異值分解，得到奇異值矩陣[SW]和[SLLn]。

[SVD（LLn）=[ULLn? SLLn? VLLn] ] （2）

[SVD（W）=[UW? SW? VW] ] （3）

4） 使用加性原則在奇異值矩陣上嵌入水印得到[S′LLn]，嵌入的強度為[λ]，[λ]根據(jù)嵌入水印后圖片的魯棒性以及可視性而定。嵌入公式如下：

[S′LLn=SLLn+λSW] （4）

5） 用[S′LLn]根據(jù)式（5）進行逆奇異值分解得到[LLn′]，最后對其進行[n1]級不可分小波逆變換得到含水印的圖像[IW]。

[LLn′=ULLn×S′LLn×VTLLn] （5）

2.2? 水印的提取過程

1） 對宿主圖像[I]和含水印圖像[IW]用本文方法構造的低高通濾波器進行[n1]級不可分小波變換。

2） 分別對經過不可分小波變換后的低頻子圖[LLn]和[LLn′]進行SVD奇異值分解，得到[SLLn]和[S′LLn]。

3） 用式（6）進行加性嵌入過程的逆過程，得到水印的奇異值矩陣[S′W]。

[S′W=（S′LLn-SLLn）λ] （6）

4） 用[S′W]進行逆奇異值分解得到[W′r]，最后將[W′r]用相同的置亂次數(shù)進行反Arnold置亂變換得到提取出的水印圖像[Wr]。

3? 仿真實驗與測試

3.1? 實驗結果分析

宿主圖像的大小是256×256的灰度Lena圖像和512×512的灰度Mandrill圖像，如圖1a）和圖2a）所示，水印圖像Word1和Word2為64×64的二值圖像如圖1b）和圖2b）所示，嵌入水印的圖像如圖1c）和圖2c）所示。在無攻擊時，對圖1c）和圖2c）做水印提取，提取出的水印與原始水印的NC均為1。

圖3a）～圖3k）前兩幅圖像是本文算法的攻擊后圖像以及提取水印圖像，后兩幅圖像是文獻[4]的攻擊后圖像以及提取水印圖像。

圖4a）～圖4n）前兩幅圖像是本文算法的攻擊后圖像以及提取水印圖像，后兩幅圖像是文獻[14]的攻擊后圖像以及提取水印圖像。

3.2? 測試指標評價

為了驗證魯棒性和不可見性的優(yōu)劣，則要討論含水印圖像的質量即可見性和提取出水印與原始水印的相似度。其中，含水印圖像的質量一般采用峰值信噪比PSNR來衡量嵌入水印的圖像與原始圖像之間的相似度[15]，通常PSNR越大，含水印圖像與原始圖像的差異就越小，含水印圖像的質量越高，水印透明性就越好，且當PSNR超過30 dB時，人眼基本不能分辨出兩幅圖像之間的差異。其公式如下：

[PSNR=10×logM×N×2552i=0M-1j=0N-1I（i，j）-Iw（i，j）2] （7）

對于提取水印與原始水印的相似性，一般采用歸一化相關系數(shù)NC來評價[16]。其中，NC越趨近于1，提取的水印與原水印越接近，表明水印效果越好。其計算公式如下：

[NC=i=0m-1j=0n-1W（i，j）×Wr（i，j）i=0m-1j=0n-1W2（i，j）×i=0m-1j=0n-1W2r（i，j）] （8）

表1列出了Mandrill圖像在本文算法和文獻[4]算法在不同攻擊下提取出水印和原始水印之間的NC值對比，表2列出了Lena圖像在本文算法和文獻[14]算法在不同攻擊下提取出水印和原始水印之間的NC值對比。

為了進一步說明本文算法嵌入水印具有較強的魯棒性，將本文算法與文獻[4]算法進行對比可知：用512×512的灰度級Mandrill圖像，用構造的6×6不可分濾波器組進行3層分解，水印圖像Word2的Arnold置亂次數(shù)為10，嵌入強度[λ]為7，嵌入水印后的圖像的PSNR為32.544。將本文算法與文獻[14]算法進行對比可知：用256×256的灰度級Lena圖像，用構造的6×6不可分濾波器組進行2層分解，水印圖像Word1的Arnold置亂次數(shù)為10，嵌入強度[λ]為7，嵌入水印后的圖像的PSNR為37.349 4。為了對比的可行性和對等性，文獻[14]用db3作為小波基，嵌入強度為7。

根據(jù)表1和表2中數(shù)據(jù)以及圖3和圖4中受到攻擊后提出水印的效果，可以看到，文獻[4]對旋轉以及JPEG壓縮因子為3以下的攻擊抵抗力較差，文獻[14]在均值濾波、中值濾波、椒鹽噪聲、剪切以及旋轉180°的攻擊下魯棒性較差。而本文算法在大多數(shù)攻擊下的NC值都在0.9以上，特別是在旋轉攻擊的魯棒性上有了很大的提升。除此之外，本文算法在遭受其他攻擊時也明顯高于文獻[4，14]的NC值，因此可以說明，在抵抗以上攻擊時，本文算法的水印魯棒性強于文獻[4，14]中的算法。

4? 結? 論

本文通過分析現(xiàn)在已提出的數(shù)字水印的相關算法中存在的缺陷與問題，提出了不可分小波與SVD奇異值分解以及Arnold變換相結合的算法。實驗數(shù)據(jù)表明，該算法在遭受濾波、噪聲、旋轉、剪切、JPEG壓縮攻擊后所提取出的水印圖像具有較好的魯棒性，特別是在對于旋轉攻擊的抵抗力明顯強于文獻[4，14]算法，且其NC值大部分都在0.9以上。綜合實驗結果表明，本文提出的算法提高了水印在幾何攻擊下的魯棒性，且對其他攻擊也有較好的抵抗力，是一種行之有效的方法。

參考文獻

[1] 姚軍財.基于小波變換和奇異值分解的圖像水印算法研究[J].光學技術，2017，43（5）：439?444.

[2] 付文亭，鄧體俊，陳海生.小波變換系數(shù)奇異值分解數(shù)字水印防偽算法[J].包裝工程，2015，36（21）：128?132.

[3] SINGH S P， BHATNAGAR G. A new robust watermarking system in integer DCT domain [J]. Journal of visual communication and image representation， 2018， 53： 86?101.

[4] 張曉琪，唐天國，胡振，等.基于DCT的數(shù)字圖像水印技術研究[J].電子制作，2018（7）：50?53.

[5] FAZLI S， MOEINI M. A robust image watermarking method based on DWT， DCT， and SVD using a new technique for correction of main geometric attacks [J]. Optik?international journal for light and electron optics， 2016， 127（2）： 964?972.

[6] SINGH A K， DAVE M， MOHAN A. Hybrid technique for robust and imperceptible image watermarking in DWT?DCT?SVD domain [J]. National academy science letters， 2014， 37（4）： 351?358.

[7] GUO J M， PRASETYO H. False?positive?free SVD?based image watermarking [J]. Journal of visual communication and image representation， 2014， 25（5）： 1149?1163.

[8] 劉斌，高強.基于二維不可分小波變換的矩不變量[J].電子與信息學報，2016，38（8）：2085?2090.

[9] LIBAL U， HASIEWICZ Z. Risk upper bound for a NM?type multiresolution classification scheme of random signals by Daubechies wavelets [J]. Engineering applications of artificial intelligence， 2017， 62： 109?123.

[10] LIU Bin， LIU Weijie. The lifting factorization of 2D 4?channel nonseparable wavelet transforms [J]. Information sciences， 2018， 456： 113?130.

[11] 劉斌，鄭凱凱.基于四通道不可分小波的均值漂移目標跟蹤方法[J].量子電子學報，2018，35（1）：13?22.

[12] 李琪，廖鑫，屈國慶，等.基于Arnold變換的數(shù)字圖像自適應隱寫算法[J].通信學報，2016，37（6）：192?198.

[13] GUO Qiang， ZHANG Caiming， ZHANG Yunfei， et al. An efficient SVD?based method for image denoising [J]. IEEE tran?sactions on circuits & systems for video technology， 2016， 26（5）： 868?880.

[14] 劉琰，周理.基于小波變換域的數(shù)字圖像嵌入和提取方法[J].沈陽工業(yè)大學學報，2019，41（1）：68?72.

[15] 李國慶，趙洋，劉青萌，等.多層感知分解的全參考圖像質量評估[J].中國圖象圖形學報，2019，24（1）：149?158.

[16] 曹軍梅.一種用于MP3文件版權保護的魯棒水印算法[J].現(xiàn)代電子技術，2016，39（4）：6?8.