摘 要： 提出基于Logistic混沌映射加密與NSCT?DWT?SVD結(jié)合的關(guān)于YUV復(fù)合數(shù)字水印算法。首先將RGB載體圖像轉(zhuǎn)換成YUV彩色模型圖像，提取Y分量進(jìn)行非采樣Contourlet變換（NSCT）。將產(chǎn)生的低頻子帶進(jìn)行小波變換，并對(duì)小波變換的低頻系數(shù)進(jìn)行奇異值分解（SVD），然后對(duì)水印信息進(jìn)行Logistic混沌映射加密，確定嵌入強(qiáng)度，將加密的水印信息疊加到彩色圖像的奇異值上。仿真實(shí)驗(yàn)表明，該算法在保證水印嵌入信息量的前提下，滿足水印信息的隱蔽性，同時(shí)該算法對(duì)常規(guī)攻擊的穩(wěn)健性較強(qiáng)，尤其是抵抗旋轉(zhuǎn)攻擊和JPEG壓縮攻擊時(shí)，仍能提取出較清晰的水印圖像。

關(guān)鍵詞： Logistic混沌映射加密； YUV； DWT； NSCT； SVD； 多重?cái)?shù)字水印算法

中圖分類號(hào)： TN911.73?34； TP391.41 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2016）10?0037?05

NSCT?DWT?SVD color watermarking algorithm based on Logistic chaotic encryption

MA Ting1， CHEN Nongtian2

（1. School of Computer Science， Civil Aviation Flight University of China， Guanghan 618307， China；

2. Aviation Engineering Institute， Civil Aviation Flight University of China， Guanghan 618307， China）

Abstract： The compound digital watermarking algorithm related to YUV based on combination of Logistic chaotic mapping encryption and NSCT?DWT?SVD is proposed. The RGB carrier image is converted into YUV color model image， and then the Y component is extracted to perform NSCT. The wavelet transform for the produced low?frequency sub?band is conducted， and then SVD is performed for low?frequency coefficient of the wavelet transform. After that the Logistic chaotic mapping encryption for watermark information is carried out to determine the embedding strength， and superimpose the encrypted watermarking information on the singular value of a color image. The simulation experiment results show that the algorithm can meet the concealment demand of watermark information on the premise of ensuring the amount of watermarking embedding information. Meanwhile， the algorithm has strong robustness to resist common attack. Especially while the algorithm is resisting the rotation attack and JPEG compression attack， it still can extract the clear watermarking image.

Keywords： Logistic chaotic mapping encryption； YUV； DWT； NSCT； SVD； multi?digital watermarking algorithm

0 引 言

目前，計(jì)算機(jī)網(wǎng)絡(luò)技術(shù)發(fā)展迅速，人們可以通過網(wǎng)絡(luò)找到所需要的信息，也可以通過網(wǎng)絡(luò)傳播信息，在獲取信息如此便易的現(xiàn)狀下，信息的安全性自然就變得尤為重要，高魯棒性數(shù)字水印算法的研究也迫在眉睫。現(xiàn)在，大部分的算法都是針對(duì)灰色圖像，但在實(shí)際應(yīng)用的過程當(dāng)中，彩色圖像應(yīng)用更為廣泛；因此本文針對(duì)彩色圖像，提出一個(gè)具有隱蔽性又具有高魯棒性的彩色數(shù)字水印算法[1?7]。

離散小波變換（Discrete Wavelet Transform，DWT）是空間域算法中應(yīng)用最為廣泛的算法，它具有較好的時(shí)頻局部性和多尺度分析性，是MPEG和JPEG壓縮標(biāo)準(zhǔn)的核心算法。通過已有關(guān)于DWT水印算法[8?13]的研究發(fā)現(xiàn)，DWT在方向上的捕捉還不足，2005年，Do和Cunha提出的非采樣Contourlet變換NSCT[14?15]不僅具有Contourlet的優(yōu)點(diǎn)，還具備了平移不變性的特性。奇異值分解（Singular Value Decomposition，SVD）是一種數(shù)字矩陣變換方法，由于奇異值對(duì)常見的攻擊具有特有的穩(wěn)定性，所以在數(shù)字水印算法中引用SVD，可以加強(qiáng)算法對(duì)抗常見攻擊的能力，尤其是最難抵抗的幾何攻擊。文獻(xiàn)[2?3]是DWT與SVD結(jié)合的關(guān)于彩色圖像的水印算法。文獻(xiàn)[2]分析了單純基于DWT和單純基于SVD算法的不足，提出了一種將DWT與SVD結(jié)合的關(guān)于RGB彩色圖像的水印算法；實(shí)驗(yàn)證明，此算法對(duì)抗各種攻擊具有較好的魯棒性。文獻(xiàn)[3]針對(duì)全息水印抵抗濾波以及幾何攻擊不足的缺點(diǎn)，將SVD和全息算法的優(yōu)點(diǎn)結(jié)合起來，首先對(duì)彩色圖像的R分量進(jìn)行一級(jí)DWT，再對(duì)分解的低頻進(jìn)行SVD，在奇異值上嵌入全息水印，提出了一種DWT與SVD結(jié)合的關(guān)于RGB彩色圖像的全息水印算法。文獻(xiàn)[5?6]是基于Contourlet變換的關(guān)于彩色圖像的水印算法，文獻(xiàn)[5]將全息算法和Contourlet變換特點(diǎn)相結(jié)合，將全息水印嵌入到對(duì)Y分量進(jìn)行Contourlet變換后的子帶中，提出基于Contourlet變換關(guān)于YUV彩色圖像的全息水印算法。文獻(xiàn)[6]則是將Contourlet變換與SVD結(jié)合的關(guān)于YUV彩色圖像的水印算法。實(shí)驗(yàn)證明文獻(xiàn)[5?6]所提出算法的穩(wěn)健性比小波算法有所提高。文獻(xiàn)[7]則是基于NSCT關(guān)于RGB彩色圖像的水印算法，算法僅采用了NSCT，在分解后的低高頻部分采用不同的嵌入強(qiáng)度，在RGB各顏色分量中嵌入水印，NSCT算法同小波變換的算法相比較，魯棒性提高得很明顯。

綜上所述，本文提出了基于Logistic混沌映射加密的NSCT?DWT?SVD彩色水印算法，算法首先將RGB彩色圖像轉(zhuǎn)換成YUV彩色圖像，在分離的Y分量上嵌入水印，為了算法的防偽安全，算法先對(duì)水印圖像進(jìn)行Logistic混沌加密，通過算法的仿真實(shí)驗(yàn)結(jié)果表明，該算法不但確保了嵌入的水印信息量，而且提高了算法的穩(wěn)健性，尤其是抗幾何攻擊和JPEG壓縮攻擊的能力。

1 基本理論

1.1 YUV顏色模型

YUV顏色模型是與RGB不同的的一種顏色色度空間。顏色空間是一個(gè)三維坐標(biāo)系統(tǒng)，RGB和YUV都是色彩空間，用于表示顏色；YUV（亦稱YCrCb）是被歐洲電視系統(tǒng)所采用的一種顏色編碼方法（屬于PAL），是PAL和SECAM模擬彩色電視制式采用的顏色空間。YUV主要用于優(yōu)化彩色視頻信號(hào)的傳輸，使其向后兼容老式黑白電視。與RGB視頻信號(hào)傳輸相比，它最大的優(yōu)點(diǎn)在于只需占用極少的帶寬（RGB要求三個(gè)獨(dú)立的視頻信號(hào)同時(shí)傳輸）。YUV與RGB兩者可以相互轉(zhuǎn)化。

在YUV空間中，每一個(gè)顏色有一個(gè)亮度信號(hào)Y，即灰階值，是一個(gè)基帶信號(hào)。兩個(gè)色度信號(hào)U和V，作用是描述影像色彩及飽和度，用于指定像素的顏色。亮度Y是通過RGB輸入信號(hào)來創(chuàng)建的，方法是將RGB信號(hào)的特定部分疊加到一起。色度U和V則定義了顏色的兩個(gè)方面：色調(diào)與飽和度，分別用Cr和Cb來表示。其中，Cr反映了RGB輸入信號(hào)紅色部分與RGB信號(hào)亮度值之間的差異；而Cb反映的是RGB輸入信號(hào)藍(lán)色部分與RGB信號(hào)亮度值之同的差異。通過運(yùn)算，Y，U，V三分量可以還原出R（紅），G（綠），B（藍(lán)）。亮度信號(hào)是強(qiáng)度的感覺，它和色度信號(hào)斷開，采用YUV色彩空間的重要性是它的亮度信號(hào)Y和色度信號(hào)U，V是分離斷開的。這樣的話強(qiáng)度就可以在不影響顏色的情況下改變。如果只有Y信號(hào)分量而沒有U，V信號(hào)分量，那么這樣表示的圖像就是黑白灰度圖像。該模型將亮度和色調(diào)分開，適用于圖像處理中對(duì)光照比較敏感的情況。彩色電視采用YUV空間正是為了用亮度信號(hào)Y解決彩色電視機(jī)與黑白電視機(jī)的相容問題，使黑白電視機(jī)也能接收彩色電視信號(hào)。

YUV與RGB相互轉(zhuǎn)換的公式如下（RGB取值范圍均為0～255）︰

YUV轉(zhuǎn)換成RGB：

2 水印算法

2.1 水印嵌入過程

首先對(duì)水印信息進(jìn)行Logistic混沌加密，再將原始RGB載體圖像轉(zhuǎn)換成YUV色度模型圖像，對(duì)轉(zhuǎn)換后提取的Y分量進(jìn)行3層NSCT，將得到與原圖像大小相同的低頻部分進(jìn)行二級(jí)小波分解，加大了嵌入信息量，為提高算法的魯棒性，再次對(duì)分解的低頻部分進(jìn)行SVD，在確定嵌入強(qiáng)度的基礎(chǔ)上，在奇異值上加入加密的水印信息，重構(gòu)含水印信息的低頻分量，對(duì)各系數(shù)一同進(jìn)行小波反變換，以及NSCT逆變換，得到Y(jié)分量，與U和V分量重構(gòu)YUV色度圖像，再將其轉(zhuǎn)換成RGB水印圖像，見圖4。

水印嵌入步驟如下：

Step1：讀取原始RGB彩色載體圖像和水印圖像，將RGB彩色載體圖像轉(zhuǎn)換成YUV彩色模型圖像，提取Y亮度分量，對(duì)其進(jìn)行3層NSCT，對(duì)水印圖像進(jìn)行Logistic變換加密；

Step2：對(duì)NSCT變換后的與原載體圖像大小相同的低頻子帶進(jìn)行兩級(jí)小波變換分解；

Step3：對(duì)小波變換分解后的低頻子帶矩陣進(jìn)行SVD分解：[A=USVT]，得到奇異值區(qū)間[S=[λ1，λ2，…，λt]]；

Step4：確定嵌入強(qiáng)度，將置亂后的水印信息迭加到分解的奇異值矩陣上得到[S′]：[S′=S+αW]；

Step5：將含水印的奇異值矩陣[S′]再次進(jìn)行SVD：[S′=U1S1V1T]，得到奇異值矩陣[S′]，將其替代原來的奇異值矩陣[S=[λ1，λ2，…，λt]]，與對(duì)應(yīng)的酉矩陣重構(gòu)低頻子帶部分；

Step6：采用逆DWT變換對(duì)系數(shù)進(jìn)行重構(gòu)，再對(duì)重構(gòu)的系數(shù)連同其他子帶進(jìn)行NSCT逆變換，得到Y(jié)UV含水印圖像，再將其轉(zhuǎn)變成RGB含水印圖像。

2.2 水印提取過程

水印嵌入流程的逆過程則是水印提取的流程，水印提取流程如圖5所示。

水印提取過程如下：

Step1：將RGB載體圖像和含水印圖像進(jìn)行轉(zhuǎn)換，轉(zhuǎn)換成YUV模型圖像，提取Y分量，對(duì)它們分別進(jìn)行三層NSCT，對(duì)得到的低頻子帶矩陣進(jìn)行兩級(jí)小波變換分解；

Step2：對(duì)小波分解后的低頻子帶矩陣[A]和[A′]分別進(jìn)行SVD：[A=USVT]，[A′=U′S1′V′T]；

Step3：將得到的奇異值矩陣相減得到提取的水印奇異值；

Step4：最后得到提取的置亂水印信號(hào)：[W′=][（U1S1′VT1-S）α]；

Step5：對(duì)水印信號(hào)進(jìn)行解密，提取出水印信息。

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

本文實(shí)驗(yàn)采用相關(guān)系數(shù)NC如式（7）所示，峰值信噪比（Peak Signal to Noise Ratio，PSNR）如式（8）所示，作為衡量水印穩(wěn)健性和透明性的評(píng)判標(biāo)準(zhǔn)。

[NC=i=1mj=1na（i，j）*b（i，j）i=1mj=1n[a（i，j）]2] （7）

[PSNR=10lg（m?n?max（max（A?A′）））j=1ni=1m（A-B）2] （8）

實(shí)驗(yàn)選用512×512的RGB模型Peppers圖像作為原始載體圖像，64×64的“飛行學(xué)院”作為水印圖像，NSCT變換中的方向?yàn)V波器選用“pfilter”基函數(shù)，塔形濾波器選用“dfilter”基函數(shù)，圖6為原始RGB圖像、RGB圖像轉(zhuǎn)換成YUV模型圖像及嵌入水印后圖像，圖7為原始水印圖像及提取的水印圖像。

從圖6看出該算法具有較強(qiáng)的透明性，嵌入水印后的圖像幾乎沒有失真，圖像質(zhì)量高，即從人的感官上，嵌入水印后的圖像看不出有任何的改變。經(jīng)式（8）計(jì)算，原始載體圖像與嵌入水印后圖像的PSNR值為52.157 9 dB，也說明本文算法具有較好的隱蔽性。本算法由于采用NSCT，在同等規(guī)則下能夠在原始載體圖像中嵌入更多的水印信息量，圖像在感觀上仍未受到影響，同時(shí)算法對(duì)嵌入強(qiáng)度沒有要求，增加了嵌入強(qiáng)度嵌入水印后，仍能保證水印的不可見性。

從圖7可看出，提取的水印圖像與原始水印圖像相比完全沒有變化，NC值達(dá)到了1。為了更好地評(píng)估本算法的穩(wěn)健性，給出了圖像進(jìn)行噪聲攻擊、濾波攻擊、剪切攻擊以及壓縮攻擊、旋轉(zhuǎn)攻擊后的攻擊測(cè)試，如圖8所示。圖8（a）～（g）是在含水印圖像經(jīng)過常規(guī)攻擊后提取的水印信息圖像，通過圖中所顯示經(jīng)過攻擊后提取的信息，在此也直接可以識(shí)別出水印的內(nèi)容“飛行學(xué)院”， 圖8（d）～（g）提取的水印信息雖含有些許噪音，但水印的內(nèi)容仍然在圖中清晰可見，為了更直觀地說明本文算法的魯棒性。表1列出了各種攻擊下的NC值，可以看出，本文算法在抵抗常見的攻擊NC值均達(dá)到了0.98以上，尤其是對(duì)抗較難抵抗的旋轉(zhuǎn)攻擊，算法的NC值也達(dá)到0.989 8，從圖8（e）也清晰可見水印信息的內(nèi)容。

表2給出了本文算法與文獻(xiàn)[7]水印算法的抗常規(guī)攻擊的魯棒性測(cè)試比較分析表。從表1、表2中的實(shí)驗(yàn)數(shù)據(jù)結(jié)果顯示可知，文獻(xiàn)[7]在抵抗椒鹽噪聲攻擊時(shí)NC值在0.92以上，在其他各攻擊上也表現(xiàn)出了較好的穩(wěn)健性，通過比較可知本文算法抵抗常見攻擊的魯棒性均有明顯的提高，尤其是受到旋轉(zhuǎn)幾何攻擊與JPEG壓縮攻擊后，本文算法仍表現(xiàn)出了很好的穩(wěn)健性，NC值均達(dá)到0.989 8和0.999 1。由此可見，本文提出的算法對(duì)抗常見攻擊的抵抗能力上提高了很多，尤其是抵抗旋轉(zhuǎn)幾何攻擊的能力以及JPEG壓縮攻擊。

旋轉(zhuǎn)和高斯濾波攻擊后的圖像

4 結(jié) 語

根據(jù)目前應(yīng)用廣泛的彩色圖像，將RGB圖像首先轉(zhuǎn)換成YUV圖像，再將NSCT，DWT以及SVD的特點(diǎn)結(jié)合起來，提出一種魯棒性高，隱蔽性好的多重彩色數(shù)字水印算法。為了確保水印的安全性，首先對(duì)水印圖像進(jìn)行Logistic混沌加密，最后通過在Matlab 2010上進(jìn)行仿真實(shí)驗(yàn)。通過實(shí)驗(yàn)結(jié)果可知，本文所提出的水印算法，不但確保了水印的嵌入信息量，同時(shí)滿足了水印信息的隱蔽性；并且通過與文獻(xiàn)[7]的NC值相比較，算法抵抗常見攻擊的能力均有大幅提高，尤其是抵抗旋轉(zhuǎn)攻擊和JPEG壓縮攻擊的能力。

參考文獻(xiàn)

[1] 李毓彬，孫劉杰，包觀笑.一種基于QR碼的彩色圖像數(shù)字水印[J].包裝工程，2015（7）：91?95.

[2] 褚靜，徐安成，張美鳳.DWT和SVD相融合的彩色圖像水印算法[J].電視技術(shù)，2013，37（17）：29?32.

[3] 李晨璐，孫劉杰.一種新的彩色圖像變換域全息水印算法[J].包裝工程，2013（9）：84?88.

[4] 謝斌，劉珊，任克強(qiáng).基于DCT的自適應(yīng)多重彩色圖像盲水印算法[J].電視技術(shù)，2014，38（9）：21?24.

[5] 孫劉杰，徐卓.一種基于Contourlet變換的彩色圖像全息水印算法[J].包裝工程，2013（9）：10?13.

[6] 羅可，馮喬生，章秀君.基于CT?SVD彩色圖像盲水印方案[J].微型機(jī)與應(yīng)用，2010，29（3）：22?25.

[7] 梁棟，殷兵，于梅，等.基于非抽樣Contourlet變換的彩色圖像數(shù)字水印算法[J].光學(xué)學(xué)報(bào)，2008，28（8）：1469?1474.

[8] 周中原，孫劉杰，唐波，等.一種抗旋轉(zhuǎn)的全息水印算法[J].包裝工程，2013（19）：95?100.

[9] 陳偉琦，李倩.基于DWT?SVD的圖像雙零水印算法[J].計(jì)算機(jī)工程與科學(xué)，2014，36（10）：1991?1996.

[10] 張雷洪，唐波，李碧程，等.基于DWT?SVD的強(qiáng)魯棒全息水印算法研究[J].包裝工程，2013（21）：105?109.

[11] 王曉紅，許詩(shī)旸.基于W?SVD的強(qiáng)魯棒性復(fù)制圖像水印算法[J].包裝工程，2012（21）：113?119.

[12] 鄧小亞.小波變換和奇異值分解的魯棒水印算法[J].激光雜志，2015，36（4）：86?89.

[13] 朱光，張軍亮.基于SVD和小波包分解的自適應(yīng)魯棒水印算法[J].計(jì)算機(jī)應(yīng)用研究，2013，30（4）：1230?1233.

[14] ZHOU J P， CUNHA A L， DO M N. Nonsubsampled contourlet transform： construction and application in enhancement [C]// Proceedings of 2005 IEEE International Conference on Image Processing. Genoa： IEEE， 2005： 469?472.

[15] DA CUNHA A L， ZHOU J P， DO M N. The nonsubsampled contourlet transform： theory， design， and applications [J]. IEEE transactions on image processing， 2006， 15（10）： 3089?3101.