劉 鵬 ，魏 娟，楊?lèi)?ài)麗

(安陽(yáng)師范學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，河南 安陽(yáng) 455000)

一種隨機(jī)選取嵌入位置的DWT域盲水印

劉 鵬 ，魏 娟，楊?lèi)?ài)麗

(安陽(yáng)師范學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，河南 安陽(yáng) 455000)

提出了一種基于離散小波變換域的盲水印算法.水印嵌入前，利用隨機(jī)的狀態(tài)對(duì)水印圖像進(jìn)行置亂.對(duì)載體圖像進(jìn)行離散小波變換后，利用隨機(jī)狀態(tài)在分解后的低頻子帶中量化系數(shù)的思想隨機(jī)地嵌入水印信息.由于置亂和選取的位置具有隨機(jī)性，所以該算法有效地增加了攻擊者提取水印的難度，提高了算法的安全性.仿真結(jié)果表明，該算法具有較強(qiáng)的魯棒性，可以抵抗諸如噪聲、壓縮、幾何裁剪以及中值濾波等常見(jiàn)的圖像攻擊.

離散小波變換(DWT)；圖像置亂；盲水??；量化水印嵌入

數(shù)字水印作為保護(hù)數(shù)字圖像版權(quán)和完整性的有效手段，主要分為空間域算法和變換域算法.一般地，變換域算法比空間域算法的魯棒性能更高.通常使用的變換有離散余弦變換（DCT）、離散傅立葉變換（DFT）以及離散小波變換（DWT）[1]，由于小波變換域數(shù)字水印方法兼具空間域方法和DCT變換域方法的優(yōu)點(diǎn)，因此基于離散小波變換的數(shù)字水印算法已經(jīng)成為當(dāng)前研究的熱點(diǎn)[2-4].

魯棒性被認(rèn)為是水印算法的難點(diǎn)和關(guān)鍵，為增加水印算法的魯棒性，本文提出一種隨機(jī)選取嵌入位置的DWT域盲數(shù)字水印算法.該算法根據(jù)水印信息，在載體圖像的小波域中的低頻子帶隨機(jī)地找取嵌入的位置，然后利用量化思想進(jìn)行水印的嵌入.由于通過(guò)密鑰控制使水印的嵌入位置隨機(jī)變化，因此可增強(qiáng)水印的魯棒性.

1 隨機(jī)選取嵌入位置的DWT域盲水印算法

水印嵌入前，先對(duì)水印信息進(jìn)行置亂.置亂不改變圖像的大小，只改變像素的位置，可以加強(qiáng)水印的魯棒性，這里的置亂是由生成的隨機(jī)序列確定的.下面給出置亂的具體方法、水印的嵌入規(guī)則、量化嵌入的思想及水印嵌入的具體步驟.

1.1 圖像的置亂

用MATLAB命令rand('state',t)[5]設(shè)置一隨機(jī)狀態(tài)t，再用命令randperm(n)輸出從1到n的隨機(jī)排序的數(shù).隨機(jī)序列由狀態(tài)t唯一確定.

設(shè)要置亂的圖像矩陣f的大小是m×n，按照列增長(zhǎng)的方式將矩陣f變換成長(zhǎng)度為m×n的向量f1，用上述生成隨機(jī)序列的方法生成一個(gè)長(zhǎng)度為m×n的隨機(jī)序列f2，使f1中的元素變換后的位置是f2中相對(duì)應(yīng)的元素的值，再把變換位置后的列向量f1變回矩陣的形式，這個(gè)矩陣就是f置亂后的圖像矩陣.二值水印圖像w用這種方法置亂，有很好的隨機(jī)性.隨機(jī)序列輸出時(shí)的狀態(tài)t可用作密鑰.圖1是水印圖像，圖2是狀態(tài)t=4置亂后的圖像.

圖1 水印圖像

圖2 置亂后的圖像

1.2 嵌入位置的選取

本文在嵌入水印前對(duì)原始圖像進(jìn)行了離散小波變換(DWT).小波變換是一種信號(hào)空間(時(shí)間)尺度分析方法，具有多分辨率分析的特點(diǎn)，可以很好地表現(xiàn)信號(hào)的時(shí)(空)頻域局部特性.利用小波變換可以提高水印在視覺(jué)上的隱蔽性和穩(wěn)健性.二維離散小波變換將原圖像分解為4個(gè)子帶，分別為低頻子帶LL、中頻子帶HL、LH和高頻子帶HH.由于低頻部分集中了圖像的大部分能量，在這些地方嵌入的水印不易丟失，因此本文選取在低頻子帶嵌入水印.

設(shè)經(jīng)過(guò)DWT處理后圖像I的大小是m×n，水印大小是m1×n1，在某一狀態(tài)t下用上述的隨機(jī)序列生成法生成一個(gè)長(zhǎng)度為m×n的隨機(jī)序列I2，把序列中前m1×n1個(gè)最大的數(shù)的位置的元素值改為1，其它位置的元素值改為0，得到I3，最后再將I3變換成m×n的矩陣，矩陣中元素值為1的位置就是水印信息要嵌入的位置.例如，圖3為Barbara圖像，圖4中的白色斑點(diǎn)部分是狀態(tài)t=4下所選取的嵌入水印的位置.

圖3 Barbara圖像

圖4 水印嵌入的位置

1.3 量化分析

量化思想的主要目的是為了實(shí)現(xiàn)水印盲檢測(cè).其主要思想是根據(jù)水印信息的不同將原始載體數(shù)據(jù)量化到不同的量化區(qū)間，而檢測(cè)時(shí)根據(jù)數(shù)據(jù)所屬的量化區(qū)間來(lái)識(shí)別水印信息.

設(shè)量化步長(zhǎng)為a，二值水印圖像為w，待修改的水印信息為g，修改后的水印信息為g′，令m=[g/a]（中括號(hào)為取整運(yùn)算），r=g-ma，則量化公式為：

1.4 水印的嵌入

水印的嵌入模型如圖5，具體過(guò)程如下：

第1步：按照1.1中的方法對(duì)水印圖像w進(jìn)行置亂得到w′，水印的置亂狀態(tài)t1作為密鑰；

第2步：將原始圖像f進(jìn)行n級(jí)小波變換，選取低頻區(qū)域，根據(jù)水印圖像的大小按照2.2所述的方法在變換后的低頻區(qū)域X中選取嵌入水印的位置，其中選取嵌入位置的狀態(tài)t2作為密鑰；

第3步：利用2.3所述的量化思想，將w′嵌入到X選取的位置中得到X′；

第4步：將X′插回到原始圖像分解的系數(shù)列中的位置上去，然后將新的系數(shù)列進(jìn)行逆小波變換，得到含水印圖像.

2 關(guān)于隨機(jī)選取嵌入位置的DWT域盲水印的提取

水印的提取過(guò)程是嵌入的逆過(guò)程，由于利用了量化的思想，因此該算法屬于盲水印.提取時(shí)需要用到的密鑰為水印置亂狀態(tài)t1，嵌入位置選取時(shí)的狀態(tài)t2以及量化步長(zhǎng)a.提取模型如圖6.

設(shè)w′為提取出的水印信息，g′是量化參數(shù)g量化后的值， r′=g ′? [g ′/a ]ia ，則水印提取方法為：

3 仿真結(jié)果與分析

本文采用512×512的灰度圖像作為載體圖像如圖7-a，64×64的二值圖像作為水印圖像如圖7-b.利用嵌入算法進(jìn)行水印嵌入，嵌入時(shí)對(duì)載體圖像進(jìn)行n=2級(jí)的DWT變換，選取置亂狀態(tài)為t1=4，位置狀態(tài)為t2=4，量化步長(zhǎng)a=0.21.含水印圖像如圖7-c，提取出的水印圖像如圖7-d.

圖7 含水印圖像及提取出的水印

表1 含水印圖像抵抗各種攻擊后的PSNR值與NC

嵌入信息后的圖像相對(duì)于原圖像有較高的峰值信噪比PSNR= 42. 3748dB，并且通過(guò)圖7-c可知，該算法具有良好的不可見(jiàn)性.另外，對(duì)含水印圖像分別進(jìn)行噪聲、濾波、縮放、裁剪、JPEG有損壓縮等各種攻擊能力測(cè)試，得到攻擊后的PSNR值和相似度(NC)見(jiàn)表1.結(jié)果表明，該算法具有良好的不可見(jiàn)性和魯棒性.

4 結(jié)語(yǔ)

本文利用離散小波變換來(lái)嵌入水印，由于嵌入前對(duì)水印圖像進(jìn)行了置亂，并且嵌入過(guò)程中選取的水印嵌入位置是隨機(jī)的，因此增強(qiáng)了算法的隱蔽性以及水印的安全性.嵌入水印的圖像不可見(jiàn)性和魯棒性都很好.仿真結(jié)果表明，該算法可以有效地抵抗常見(jiàn)的攻擊如JPEG壓縮、高斯噪聲、椒鹽噪聲、裁剪、中值濾波等.

[1] 張永庫(kù). 數(shù)字水印技術(shù)及典型的攻擊分析[J]. 信息技術(shù), 2006(11): 156-158.

[2] 馮象初, 甘小冰, 宋國(guó)鄉(xiāng). 數(shù)值泛函與小波理論[M]. 西安: 電子科技大學(xué)出版社, 2003.

[3] 孫圣和, 陸哲明, 牛夏牧. 數(shù)字水印技術(shù)及應(yīng)用[M]. 北京: 科學(xué)出版社, 2004.

[4] MEMON N, WONG P W. Protecting digital content[J]. Communications of the ACM, 1998, 41(7): 35-43.

[5] 趙 翔, 郝 林. 數(shù)字水印綜述[J]. 計(jì)算機(jī)工程與設(shè)計(jì), 2006, 27(11): 1946-1950.

[6] 金 聰. 數(shù)字水印理論與技術(shù)[M]. 北京: 清華大學(xué)出版社, 2008.

[7] 王秋生. 變換域數(shù)字水印嵌入算法研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2001.

[8] ANDERSON R J, PETITCOLAS F. On the limits of steganography[J]. IEEE Journal of Selected Areas in Communications, 1998, 16(4): 474-481.

（責(zé)任編輯：陳 丹）

A DWT Based on Watermarking Algorithm by Randomly Selecting Embedding Positions

LIU Peng, WEI Juan, YANG Ai-li
(School of Mathematics and Statistics, Anyang Normal University, Anyang 455000, China)

A blind watermarking algorithm based on discrete wavelet transform (DWT) was proposed. Before embedding, it scrambles the watermark by random state, then chooses the embedding position randomly in the wavelet low-frequency sub-band of the cover image according to the watermark, then embeds the watermark image with the quantization thought. Since using secret keys that make the embedding positions change randomly, it greatly enhances the security of the watermark. The simulation result shows that this algorithm is robust to common watermarking attacks.

DWT; Scrambling; Watermarking; Robustness; Watermark image embedding with the quantization thought

TP391

A

1009-2854(2010)02-0025-03

2009-11-12;

2010-01-21

劉 鵬(1982— ), 男, 河南焦作人, 安陽(yáng)師范學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院助教.