任克強(qiáng)，王法印

(江西理工大學(xué) 信息工程學(xué)院，江西 贛州 341000)

?

基于BCH和JND的小波域數(shù)字圖像水印算法

任克強(qiáng)，王法印

(江西理工大學(xué) 信息工程學(xué)院，江西 贛州 341000)

為了提高數(shù)字圖像水印的魯棒性，提出一種結(jié)合BCH碼和JND模型的小波域數(shù)字圖像水印算法。該算法采用Arnold變換和BCH碼對(duì)二值水印圖像進(jìn)行置亂加密以及編碼，對(duì)載體圖像分塊后再進(jìn)行二維離散小波變換，將加密并且編碼后的水印信息嵌入到載體圖像各子塊的二維離散小波變換的低頻子帶中，嵌入強(qiáng)度根據(jù)載體圖像的JND值自適應(yīng)的選取。實(shí)驗(yàn)結(jié)果表明，水印的透明性很好，在遭受不同類型和強(qiáng)度的水印攻擊時(shí)，具有較強(qiáng)的魯棒性。

數(shù)字圖像水?。蛔钚】捎X察誤差；BCH碼；離散小波變換；Arnold變換

數(shù)字水印技術(shù)是在不影響原始數(shù)字媒體質(zhì)量的條件下，通過嵌入算法在載體中添加一定數(shù)量的標(biāo)志性信息，并且保證這些信息具有不易察覺和破壞，但易于被原始載體所有者提取和識(shí)別。作為版權(quán)保護(hù)和信息安全領(lǐng)域的一個(gè)重要分支，數(shù)字水印技術(shù)被普遍用于文本、音頻、視頻和圖像等數(shù)字媒體的版權(quán)保護(hù)[1-4]。

當(dāng)前，數(shù)字水印算法主要集中在兩個(gè)方面：一種是空間域，基于該區(qū)域算法是直接對(duì)載體中的像素值進(jìn)行修改完成水印信息的嵌入；另一種是變換域，先將圖像做某種變換，然后再進(jìn)行水印嵌入[5]。文獻(xiàn)[6]提出了一種抗旋轉(zhuǎn)、剪切和縮放等幾何型攻擊的LSB數(shù)字水印算法，通過修改載體中每個(gè)像素值二進(jìn)制表示的最低位進(jìn)行水印嵌入，雖然水印嵌入的容量較大，但受到攻擊時(shí)魯棒性不理想。文獻(xiàn)[7]基于抗亮度攻擊的基本思想，利用圖像亮度的變換公式，對(duì)離散余弦變換之后的直流系數(shù)進(jìn)行修正來嵌入水印，在抵御高斯、中值以及平均濾波、JPEG壓縮和重采樣等攻擊時(shí)，表現(xiàn)出良好的魯棒性。文獻(xiàn)[8]提出一種基于改進(jìn)的脊波變換數(shù)字水印算法，選擇恰當(dāng)?shù)南禂?shù)，利用奇偶量化水印信息，達(dá)到在載體中嵌入水印的目的。文獻(xiàn)[9]利用人類視覺的固有屬性以及離散小波變換(Discrete Wavelet Transform，DWT)的優(yōu)良特性將水印嵌入載體圖像中，具有較好的透明性和魯棒性。

本文結(jié)合最小可覺察誤差(Just Noticeable Distortion，JND)模型、DWT和BCH碼的特點(diǎn)，提出一種數(shù)字圖像水印算法。算法對(duì)Arnold置亂加密后的水印圖像進(jìn)行BCH編碼，將載體圖像進(jìn)行分塊二維DWT變換，計(jì)算各子塊DWT低頻子帶的JND值，并將置亂加密和編碼的水印信息嵌入到載體圖像中，嵌入強(qiáng)度根據(jù)載體圖像的JND值自適應(yīng)選取。

1　相關(guān)理論

1.1JND模型

JND模型是根據(jù)人類視覺系統(tǒng)(Human Visual System，HVS)特性創(chuàng)建的一種模型。根據(jù)HVS的掩蓋特性，在含水印的載體圖像中，水印的透明性會(huì)受到頻率、背景亮度和紋理特性的影響。

人眼對(duì)不同頻帶信號(hào)中的頻率敏感性不同，對(duì)低頻信號(hào)的頻率較為敏感，對(duì)高頻信號(hào)的不同尺度、方向的頻率不敏感，特別是對(duì)高頻中45°方向的頻率更加不敏感。對(duì)DWT中4個(gè)子帶的頻率敏感性如下

(1)

式中：l為分解層數(shù)(l=0,1,2,3)；s為DWT分解方向(s=LL,HL,LH,HH)。

背景亮度敏感性為

j+1+y·2l-3)

(2)

式中：(x,y)為l層中的坐標(biāo)。

紋理掩蔽特性為

x·2l-3,1+j+y·2l-3]}

(3)

根據(jù)式(1)、式(2)和式(3)可以得到在第l層坐標(biāo)(x,y)處的JND值

(4)

JND值的求解過程結(jié)合了HVS對(duì)頻率和背景亮度的敏感性以及紋理的掩蓋特性，可為水印嵌入的強(qiáng)度提供依據(jù)。

1.2DWT

目前大多數(shù)數(shù)字水印算法是基于變換域的，在變換域中，小波變換具有良好的空域分解，可以模擬HVS特性，并且具有多尺度的處理能力，因而被廣泛應(yīng)用于圖像處理中。

DWT是將時(shí)間域上連續(xù)的小波變換進(jìn)行方向尺度和平移的離散化過程，將目標(biāo)圖像做高通、低通和降頻處理后分解成低頻子帶LL、水平子帶HL、垂直子帶LH和45°對(duì)角子帶HH。LL為圖像的邊緣和亮度部分，與原圖像最為接近，包含了圖像的大部分能量；其他3個(gè)子帶為載體的紋理、細(xì)節(jié)信息。載體圖像經(jīng)過兩次二維DWT，得到一、二級(jí)子帶的分布如圖1所示。

圖1　兩級(jí)小波分解后的子帶分布

2　本文算法

2.1水印預(yù)處理

1)水印圖像置亂加密

Arnold變換首先由弗拉基米爾·阿諾爾德應(yīng)用于一幅貓臉的圖像中，所以又稱貓臉變換。設(shè)載體圖像尺寸為M×M，則二維Arnold變換為

(5)

式中：n為變換次數(shù)。

對(duì)應(yīng)的Arnold逆變換為

(6)

對(duì)載體圖像進(jìn)行n次Arnold變換后可將圖像轉(zhuǎn)換為雜亂無章，相當(dāng)于對(duì)圖像進(jìn)行了加密。本文的水印圖像為64×64的二值圖像，如圖2a所示，經(jīng)過30次Arnold變換后得到圖2b所示的加密水印。

圖2　原始水印和加密水印

對(duì)于置亂加密后的水印圖像，想要從中恢復(fù)出原始水印圖像必須知道M和n的取值，否則難以得到原始水印圖像，從而可以有效提升水印的安全性。

2)水印信息糾錯(cuò)編碼

數(shù)字水印系統(tǒng)中的水印信息嵌入和提取相當(dāng)于數(shù)字通信系統(tǒng)完成一次通信過程，通信系統(tǒng)通過對(duì)待發(fā)送的消息進(jìn)行糾錯(cuò)編碼可以有效降低碼字在接收端發(fā)生錯(cuò)誤的概率。在數(shù)字圖像水印算法中，可以借鑒數(shù)字通信系統(tǒng)中糾錯(cuò)編碼的思想，利用糾錯(cuò)碼對(duì)水印信息進(jìn)行編碼，以提高水印的魯棒性。本文采用本原BCH(63，16，11)對(duì)置亂并且加密之后的水印信息進(jìn)行BCH編碼。

2.2水印嵌入

本文把置亂加密并且編碼之后的水印嵌入到載體圖像DWT的低頻子帶中，嵌入強(qiáng)度根據(jù)載體圖像的JND值自適應(yīng)選取，水印的嵌入步驟如下。

步驟1：將載體圖像裁剪成大小為8×8的子塊，并且兩兩之間均不重疊，共64×64塊。然后分別對(duì)每個(gè)子塊進(jìn)行兩級(jí)DWT變換，得到低頻子帶矩陣LL2。

步驟2：對(duì)水印圖像進(jìn)行Arnold置亂加密，將置亂加密后的水印降為一維序列并對(duì)其進(jìn)行BCH編碼得到序列M_BCH。

步驟3：計(jì)算經(jīng)兩級(jí)DWT變換得到的低頻子帶LL2的JND值，并計(jì)算每個(gè)子塊的平均值得到矩陣AVE_LL。

步驟4：將矩陣AVE_LL轉(zhuǎn)換成一維向量并進(jìn)行降序排列，并按從大到小選取低頻子帶進(jìn)行水印嵌入

(7)

式中：I(i)為原始載體圖像低頻子帶LL2的像素值；JND(i)為I(i)對(duì)應(yīng)的JND值；M_BCH(i)為序列M_BCH中第i個(gè)值。

步驟5：記錄含水印的位置信息，存入矩陣Loa_M中，將含水印的低頻子帶與其他未作修改的子帶經(jīng)過兩級(jí)IDWT變換，生成已嵌入水印之后的載體圖像。

2.3水印提取

從含水印載體圖像之中提取水印的步驟如下。

步驟1：將含水印的載體圖像裁剪成大小為8×8的子塊，并且兩兩之間均不重疊，并通過兩級(jí)DWT得到低頻子帶LL2′。

步驟2：由矩陣Loa_M找到含水印的像素值，提取水印如下

(8)

步驟3：對(duì)提取出的水印序列M_BCH′進(jìn)行BCH解碼以及Arnold逆變換，獲得提取的水印圖像。

3　實(shí)驗(yàn)結(jié)果及分析

實(shí)驗(yàn)軟硬件環(huán)境為：Intel(R) Core(TM) i5-5200U CPU @ 2.20 GHz，內(nèi)存8 Gbyte；Windows 10專業(yè)版；MATLAB R2015a。載體圖像是尺寸為512×512的Lena.jpg灰度值圖像，水印為圖2a所示的二值圖像。

衡量一種數(shù)字水印算法的性能優(yōu)劣一般包括主觀評(píng)價(jià)和客觀評(píng)價(jià)。主觀評(píng)價(jià)受測(cè)試者的主觀因素影響較大，不同的測(cè)試者可能對(duì)同一圖像產(chǎn)生完全不同的判斷，因此一般使用客觀評(píng)價(jià)來鑒別一種水印算法的優(yōu)劣。客觀評(píng)價(jià)一般采用峰值信噪比(Peak Signal to Noise Ratio，PSNR)和歸一化相關(guān)系數(shù)(Normalized Correlation Coefficient，NC)進(jìn)行評(píng)判。

(9)

式中：M和N為載體圖像的尺寸；I和I′分別表示水印嵌入前后的載體圖像；Imax為I中所有像素的最大值；I(x,y)和I′(x,y)分別為兩幅載體圖像坐標(biāo)(x,y)處的像素點(diǎn)的大小值。PSNR值越大，表明載體圖像嵌入水印或遭受攻擊后的失真越小，水印的透明性越好。

(10)

式中：m和n為水印圖像的尺寸；W(i,j)為原始水印的像素值；W′(i,j)為提取水印的像素值。NC值越大，表明提取出的水印圖像與原始的水印圖像越相似，水印的魯棒性越好。

圖3所示為原始載體圖像和已嵌入水印的載體圖像。當(dāng)PSNR>30 dB時(shí)，圖像通常表現(xiàn)出較高的視覺質(zhì)量[10]；當(dāng)PSNR>36 dB時(shí)，HVS則難以辨別其差異。含水印圖像的PSNR為44.957 8 dB，遠(yuǎn)遠(yuǎn)大于HVS對(duì)圖像分辨的閾值30 dB，說明本文水印的透明性很好。

圖3　原始圖像和含水印圖像

對(duì)含水印圖像分別進(jìn)行中值濾波(3×3)、椒鹽噪聲(0.01)、JPEG壓縮(20%)和剪切(25%)攻擊后，提取的水印如圖4所示。

圖4　攻擊后提取的水印圖像

從圖4可以看出，本文水印遭受上述攻擊后，可以較好地提取水印圖像，提取出來的水印識(shí)別度較高。

為進(jìn)一步地驗(yàn)證比較水印的透明性和魯棒性，對(duì)含水印圖像分別進(jìn)行不同攻擊類型、強(qiáng)度下的攻擊，提取水印、計(jì)算PSNR和NC值并與文獻(xiàn)[9]進(jìn)行比較，實(shí)驗(yàn)結(jié)果如表1所示。由表1可知，遭受不同類型和強(qiáng)度的攻擊后，本文和文獻(xiàn)[9]均表現(xiàn)出較好的透明性和魯棒性。除了遭受剪切(右上1/4)、剪切(右下1/4)以及高斯噪聲(0.01)攻擊后，本文的PSNR值小于文獻(xiàn)[9]外，本文的其他PSNR值均大于文獻(xiàn)[9]，說明本文水印的透明性總體上比文獻(xiàn)[9]更好。尤其是本文的NC值均大于文獻(xiàn)[9]，這是由于本文采用了糾錯(cuò)編碼技術(shù)，使得本文水印的魯棒性得到有效提升。

表1不同攻擊類型和強(qiáng)度的實(shí)驗(yàn)結(jié)果

攻擊類型和強(qiáng)度本文文獻(xiàn)[9]PSNR/dBNCPSNR/dBNC椒鹽噪聲0.0126.36550.998125.34800.93200.0222.36580.948422.27140.86170.0321.63270.902820.61030.83330.0419.38060.850818.40760.7656剪切左上1/411.97410.897311.55170.8499左下1/413.56250.91429.84160.9045右上1/410.85110.915611.94950.8186右下1/410.87810.911012.10940.9090中間1/411.65530.920711.06360.8420高斯噪聲0.0119.63820.913419.95870.80680.00820.89010.907119.95270.82060.00620.92470.850319.99670.81030.00421.03480.821019.98940.8177中值濾波3×335.71970.9572—0.9457JPEG壓縮Q=10021.03480.821019.98940.8177Q=8043.32391.000025.31030.9907Q=6039.21970.990625.17320.9416Q=4037.89110.923124.93710.8062縮放3/438.09360.9501—0.9457

4　結(jié)束語

提出一種基于BCH碼和JND模型的小波域數(shù)字圖像水印算法，算法通過Arnold變換來增強(qiáng)水印的安全性，利用HVS、JND模型以及DWT的特性來保證水印的透明性和魯棒性，并引入BCH碼來進(jìn)一步強(qiáng)化水印的魯棒性。通過仿真比較實(shí)驗(yàn)，該算法在抵御不同水印攻擊類型和攻擊強(qiáng)度時(shí)，表現(xiàn)出良好的透明性和魯棒性。

[1]WANG T，LI H. A novel scrambling digital image watermark algorithm based on double transform domains[J]. Mathematical problems in engineering，2015(2)：1-13.

[2]朱新山，陳硯鳴，董宏輝，等.基于雙域信息融合的魯棒二值文本圖像水印[J].計(jì)算機(jī)學(xué)報(bào)，2014，37(6)：1352-1364.

[3]王友衛(wèi)，劉元寧，朱曉東.用于圖像內(nèi)容認(rèn)證的半脆弱水印新算法[J].浙江大學(xué)學(xué)報(bào)(工學(xué)版)，2013，47(6)：969-976.

[4]方旺盛，李玉南，張蓉，等.基于整數(shù)小波變換的大容量彩色水印嵌入技術(shù)[J].江西理工大學(xué)學(xué)報(bào)，2011，32(3)：6-10.

[5]KASHYAP N，SINHA G R. Image watermarking using 3-level discrete wavelet transform(DWT)[J]. International journal of modern education and computer science，2012，4(3)：50-56.

[6]白冬慧，劉彥隆.抗幾何攻擊的最低有效位數(shù)字水印算法[J].電視技術(shù)，2013，37(5)：29-32.

[7]李旭東.抗亮度和對(duì)比度攻擊的DCT域圖像數(shù)字水印算法[J].光電子·激光，2013，24(6)：1184-1190.

[8]高虎明，賈麗媛，劉美玲.基于改進(jìn)脊波變換的抗攻擊數(shù)字水印算法[J].計(jì)算機(jī)應(yīng)用研究，2014，31(9)：2750-2753.

[9]葉闖，沈益青，李豪，等.基于人類視覺特性(HVS)的離散小波變換(DWT)數(shù)字水印算法[J].浙江大學(xué)學(xué)報(bào)(理學(xué)版)，2013，40(2)：152-156.

[10]姜傳賢，楊鐵軍，董明剛，等.基于線性空間隱藏模型的可逆圖像水印算法[J].自動(dòng)化學(xué)報(bào)，2014，40(10)：2324-2333.

任克強(qiáng)(1959— )，教授，碩士生導(dǎo)師，主要研究方向?yàn)樾畔㈦[藏、無線傳感器網(wǎng)絡(luò)等；

王法印(1989— )，碩士研究生，研究方向?yàn)樾畔㈦[藏。

責(zé)任編輯：時(shí)雯

Digital image watermarking algorithm in wavelet domain based on BCH and JND

REN Keqiang， WANG Fayin

(SchoolofInformationEngineering，JiangxiUniversityofScience&Technology，JiangxiGanzhou341000，China)

In order to improve the robustness of digital image watermarking， a digital image watermarking algorithm in wavelet domain based on BCH code and JND model is proposed. Arnold transform and BCH code are used to scrambling encryption and encoding of the binary watermarking image. The carrier image is divided to blocks，then each sub-block is carried on two-dimensional discrete wavelet transform respectively. The encrypted and encoded watermarking are embedded into the low frequency subband of two-dimensional discrete wavelet transform of each sub-block in the carrier image， the embedding strength is adaptively selected according to JND value of carrier image. The experimental results show that the watermarking transparency is very good， and has stronger robustness when it is attacked by different types and intensity.

digital image watermarking； just noticeable distortion； BCH code； discrete wavelet transform； Arnold transform

TN911.73； TP391

A

10.16280/j.videoe.2016.09.004

江西省教育廳青年科學(xué)基金項(xiàng)目(GJJ11132)；江西省研究生創(chuàng)新基金項(xiàng)目(YC2013-S199)

2016-01-21

文獻(xiàn)引用格式：任克強(qiáng)，王法印. 基于BCH和JND的小波域數(shù)字圖像水印算法[J]. 電視技術(shù)，2016,40(9)：22-25.

REN K Q， WANG F Y.Digital image watermarking algorithm in wavelet domain based on BCH and JND[J]. Video engineering，2016,40(9)：22-25.