張好好+吳游+于睿+馮金金

摘 要： 根據(jù)壓縮感知具有壓縮短、解碼復(fù)雜和魯棒性強(qiáng)的特征，提出一種新的基于壓縮感知的數(shù)字圖像水印算法。該算法有效克服了信號重建時(shí)稀疏度未知的問題，增強(qiáng)了數(shù)字水印系統(tǒng)的抗提取防破壞能力。在嵌入水印時(shí)，將水印的高頻系數(shù)用構(gòu)造的伯努利矩陣進(jìn)行觀測后，嵌入于圖像的低頻系數(shù)中；在提取水印時(shí)，利用正則化自適應(yīng)匹配追蹤（RAMP）算法對高頻系數(shù)進(jìn)行恢復(fù)，最后與低頻系數(shù)結(jié)合實(shí)現(xiàn)水印的重構(gòu)。實(shí)驗(yàn)結(jié)果表明：該算法在重構(gòu)精度和重構(gòu)時(shí)間上要高于同類算法，具有很好的魯棒性、透明性和安全性。

關(guān)鍵詞： 數(shù)字水?。?壓縮感知； RAMP算法； 水印重構(gòu)

中圖分類號： TN919?34； TP301.6 文獻(xiàn)標(biāo)識碼： A 文章編號： 1004?373X（2014）22?0010?04

A digital image watermark algorithm based on compressive perception

ZHANG Hao?hao， WU You， YU Rui， FENG Jin?jin

（School of Electronics and Information， Jiangsu University of Science and Technology， Zhenjiang 212003， China）

Abstract： Since compressive perception has the characteristics of short compression， decoding complexity and strong robustness， a new digital image watermark algorithm based on compressive perception is proposed. The algorithm overcame the shortcoming that the sparsity is unknown in signal reconstruction， and enhanced the function to resist illegal extraction and prevent damage of digital watermark system. When watermark is embedded， the high?frequency coefficients of the watermark are measured by the constructed Bernoulli matrix， and then the watermark is embedded in the low?frequency coefficients of the image. As for the watermark extraction， the regularized adaptive matching pursuit （RAMP） algorithm is used to recover the high?frequency coefficients， and make them combined with the preserved low?frequency coefficients to obtain the reconstructed watermark. Compared with the original compressed sensing algorithm， The experimental results demonstrate that the proposed algorithm is better than the original compressive perception algorithm in reconstruction accuracy and time consumption， and has good robustness， high transparency and security.

Keywords： digital watermark； compressive perception； RAMP algorithm； watermark reconstruction

0 引 言

隨著計(jì)算機(jī)網(wǎng)絡(luò)和信息技術(shù)的飛速發(fā)展，數(shù)字多媒體信息的傳輸和獲取變得越來越容易，由此引發(fā)的數(shù)字多媒體的版權(quán)、經(jīng)濟(jì)利益等問題受到人們的廣泛關(guān)注。數(shù)字水印技術(shù)[1]是一種有效解決版權(quán)和安全性問題的技術(shù)，被廣泛應(yīng)用在圖像、音頻、視頻等多媒體信息的處理中，并且以數(shù)字圖像水印為主，經(jīng)歷了從時(shí)域水印到變換域水印再到零水印，脆弱水印到魯棒性水印再到二重水印的發(fā)展過程。早期的水印算法主要應(yīng)用在空間域中，具有速度快的優(yōu)點(diǎn)，雖算法簡單，但魯棒性不高，且易影響載體圖像；而變換域水印具有容量大、抗干擾能力強(qiáng)、透明度高等特點(diǎn)，并且可以與現(xiàn)有的數(shù)據(jù)壓縮標(biāo)準(zhǔn)JPEG兼容，信號能擴(kuò)展到所有像素上，成為數(shù)字水印技術(shù)研究的熱點(diǎn)。雖然已經(jīng)進(jìn)行了大量的研究，但到目前為止水印的抗攻擊和防提取能力有限，這大大限制了數(shù)字水印技術(shù)在其他領(lǐng)域的應(yīng)用。因此，尋求一種具有更好的抗攻擊能力和防提取能力的數(shù)字水印新技術(shù)具有重要的研究意義。

近幾年新興的壓縮感知技術(shù)為其提供了這種可能，它突破奈奎斯特采樣定理的限制，使壓縮、采樣同步完成，通過合適的重構(gòu)算法，使用少量的觀測值便可恢復(fù)原始的稀疏信號，可有效減小系統(tǒng)存儲和降低硬件的設(shè)計(jì)，保護(hù)數(shù)字水印信息，比傳統(tǒng)的數(shù)字水印技術(shù)易操作和實(shí)現(xiàn)。本文把壓縮感知理論應(yīng)用于數(shù)字水印技術(shù)[2?3]中，結(jié)合人類視覺系統(tǒng)（HVS）的特性選擇嵌入位置，增強(qiáng)了水印的魯棒性和安全性，同時(shí)降低了存儲空間。

1 壓縮感知理論

壓縮感知[4]的核心內(nèi)容是已知信號進(jìn)行某一正交變換后得到的系數(shù)大部分是近似為零的，就可以以遠(yuǎn)少于奈奎斯特采樣定理所要求的采樣數(shù)來重構(gòu)出原始信號。它包含兩個(gè)要點(diǎn)：信號的稀疏性和觀測的不相關(guān)性，稀疏性是壓縮感知的先決條件，不相關(guān)性是壓縮感知觀測可行性的保證。稀疏性和不相關(guān)性分別體現(xiàn)在信號的稀疏表示和觀測矩陣的設(shè)計(jì)中，以下分別從三個(gè)方面來介紹壓縮感知的基本原理。

1.1 信號的稀疏表示

將信號[X∈RN]通過正交變換基計(jì)算信號的變換系數(shù)[θ=ψTX]，稱為信號的稀疏化?，F(xiàn)實(shí)生活中，大部分信號都不是稀疏的，但是它們在變換域是稀疏的。對于圖像信號，幾乎所有的像素都不為零，但經(jīng)過某個(gè)正交基（如DCT，DWT等）的變換，使得大多數(shù)系數(shù)的絕對值非常小，幾乎可忽略不計(jì)視為零。取絕對值較大的系數(shù)幾乎可以包含圖像的全部信息，這樣就可以將其看成是可壓縮的。只要是可壓縮的信號，就可以使用壓縮感知采樣理論。

1.2 觀測矩陣的設(shè)計(jì)

觀測矩陣的構(gòu)造直接影響到壓縮信號的重構(gòu)，不同的信號會選擇不同的稀疏表達(dá)基[5]。而針對每種不同的稀疏表達(dá)基[ψ]能否找到一個(gè)與其不相關(guān)的觀測矩陣就顯得尤為重要。用大小為M×N的采樣矩陣[φ]對變換系數(shù)[θ]進(jìn)行采樣得到[Y=φθ=φψTX]。其中觀測矩陣[φ]要求是一個(gè)與變換基[ψ]不相關(guān)、平穩(wěn)的矩陣。常用的高斯隨機(jī)矩陣滿足上述條件，但需要的存儲空間較大且計(jì)算復(fù)雜，因此本文選擇重構(gòu)信號概率較大且速度快的伯努利矩陣作為觀測矩陣。

1.3 信號的重構(gòu)

信號的重構(gòu)是通過求解一個(gè)欠定的線性方程組，從低維度的觀測值中恢復(fù)出高維度的原始信號的過程。將求解欠定方程組的問題轉(zhuǎn)化為求解0?范數(shù)意義的最優(yōu)化問題。通過公式計(jì)算可以得到X的精確或者近似X：

[minψTX0 s.t. Y=φψTX]

目前常用的重建算法有凸優(yōu)化算法和貪婪算法。由于貪婪算法具有算法簡單、重建效果好、計(jì)算量小等特點(diǎn)得到廣泛應(yīng)用，本文采用提出的RAMP算法來實(shí)現(xiàn)。

2 RAMP算法

RAMP算法把ROMP算法[6]正則化過程和SAMP算法的自適應(yīng)思想結(jié)合起來，有效克服了OMP及其延伸類算法只有首先對信號的稀疏度估計(jì)才能精確重建的局限，使得算法本身可以在迭代過程中通過可變步長自動進(jìn)行原子數(shù)目篩選來重建稀疏度未知的信號。算法中采用轉(zhuǎn)換階段的方式篩選候選集中的原子，把一次迭代過程分為多個(gè)階段，通過設(shè)置一個(gè)可變步長的方式取代所選原子數(shù)目，隨著步長的增加，實(shí)現(xiàn)在稀疏度未知的條件下逐步逼近稀疏度K取得原始信號的精確重構(gòu)。

下列步驟具體說明了該算法的過程：

（1） 余量[r0]=y，初始步長s[≠]0，階段j=1，迭代次數(shù)k=1，索引值集[Λ=?]，候選集[J=?]；

（2） 若[rk≤ε1]，則停止迭代，利用所得的原子進(jìn)行信號重構(gòu)；否則進(jìn)入下一步運(yùn)算；

（3） 用公式計(jì)算相關(guān)系數(shù)[ut]：

[ut=ujuj=rj，φj，j=1，2，…，N]

并從[ut]中找到s個(gè)最大值，將對應(yīng)的索引存入候選集J中，完成原子的第一次篩選；

（4） 第二次篩選原子：通過正則化將候選集J中索引值對應(yīng)的原子的相關(guān)系數(shù)結(jié)果存入集合[J0?J]中，其中原子的相關(guān)系數(shù)必須滿足式[u（i）≤2u（j），][其中i，j∈J0]；

（5） 更新支撐集[ΦΛ]， [Λ=Λ?J0]；

（6） 用公式[argminy-ΦΛX2]求出信號逼近[X]，同時(shí)用公式[rt=y-ΦΛX]更新余量；

（7） 如果[rt-rt-1≤ε2]時(shí)，j=j+1，s=s·j，返回步驟（3）； 否則返回步驟（2）。

其中[ε1]與[ε2]分別為迭代終止和階段轉(zhuǎn)換的閾值，在信號重構(gòu)過程中通過比較殘差和閾值的關(guān)系來判斷是否自動調(diào)整當(dāng)前步長，進(jìn)入下一個(gè)階段或者下一次迭代還是完成信號的重建，從而完成自適應(yīng)重建過程，避免了一般的貪婪算法由于選擇的稀疏度與信號的實(shí)際稀疏度不匹配現(xiàn)象的出現(xiàn)。

步驟（4）采用正則化過程進(jìn)行原子的二次篩選，可以保證未被選入支撐集的原子的能量一定遠(yuǎn)小于被選入原子的能量，同時(shí)也繼承正則化思想快速的特點(diǎn)，實(shí)現(xiàn)信號在稀疏度未知下高速、穩(wěn)定、精確地信號重建。

3 水印的嵌入與提取

本文采用256×256的Lena灰度圖作為原始載體圖像，標(biāo)有“ABCDEF”的64×64的二值圖像作為水印圖像，如圖4（a）、圖4（c）所示，選擇“db1”作為小波稀疏基，觀測矩陣為構(gòu)造的隨機(jī)伯努利矩陣，采用RAMP貪婪算法將稀疏的水印信號重構(gòu)出來，完成水印的嵌入與提取。

3.1 水印的嵌入

水印信息的具體嵌入步驟如下：

（1） 對原圖像X進(jìn)行小波多尺度變換[7]，得到含有{LH，HL，HH，LL} 多個(gè)高頻成分和低頻成分子帶系數(shù)X1；

（2）對水印圖像w進(jìn)行單層小波變換，得到 [{LH1，HL1，HH1，LL1}]四個(gè)小波子帶系數(shù)w1；

（3） 用服從伯努利分布的隨機(jī)觀測矩陣[?]分別對稀疏化的水印信息w1的高頻系數(shù)[{LH1，HL1，HH1}]進(jìn)行測量，而低頻[LL1]子帶系數(shù)則保持不變；

（4） 對觀測后的水印信號w2乘上嵌入強(qiáng)度系數(shù)q加到X1的低頻系數(shù)上，得到帶水印信息的信號X2，即[X2=X1+qw2]；

（5） 對X2采用小波逆變換重構(gòu)形成加入水印后的圖像X3。

整個(gè)嵌入過程如圖1所示。

圖1 數(shù)字水印的嵌入

3.2 水印的提取

水印信息的詳細(xì)提取過程如下：

（1） 對含水印的圖像X3用小波基實(shí)現(xiàn)稀疏化，得到P1；

（2） 分別對P1和稀疏化的原始載體圖像X1的對應(yīng)低頻系數(shù)利用加性逆規(guī)則求差值得到嵌入水印的低頻系數(shù)[LL1]；

（3） 對水印圖像w1的高頻[{LH1，HL1，HH1}]利用RAMP算法進(jìn)行重構(gòu)得到[{LH1，HL1，HH1}]；

（4） 最后[{LH1，HL1，HH1}]與[LL1]子帶結(jié)合一起進(jìn)行小波反變換得到水印圖像。

整個(gè)提取過程如圖2所示。

圖2 數(shù)字水印的提取

4 實(shí)驗(yàn)仿真分析

為了達(dá)到本文算法的預(yù)期效果，分別通過了3個(gè)實(shí)驗(yàn)進(jìn)行驗(yàn)證。實(shí)驗(yàn)一驗(yàn)證提出的RAMP算法的重構(gòu)效果，實(shí)驗(yàn)二是驗(yàn)證此算法下水印的嵌入與提取效果，實(shí)驗(yàn)三通過攻擊測試驗(yàn)證本文算法的魯棒性。

由于在Matlab 2010b下構(gòu)造的M行N列的伯努利矩陣[7]在仿真實(shí)驗(yàn)時(shí)，每次選擇的M行值進(jìn)行觀測時(shí)均不相同，這會影響到圖像重建的效果。因此，為了確保誤差比較小，在同一M值下進(jìn)行大量實(shí)驗(yàn)，得出重構(gòu)圖像的峰值信噪比（PSNR）相差會在±0.2 dB，重構(gòu)時(shí)間和互相關(guān)系數(shù)（NC）為±0.01。下面圖形的結(jié)果均為將程序運(yùn)行3次取其平均值繪制而成的。

實(shí)驗(yàn)一：對二值水印的壓縮感知編碼和重建過程進(jìn)行分析，圖3是單層小波變換后壓縮感知的處理結(jié)果，其中采樣率分別為0.3，0.5，0.7和0.9。可以看出，采樣率>0.5時(shí)，重構(gòu)效果比較好，若降低采樣率，重構(gòu)圖像將會變得很模糊。隨著采樣系數(shù)的增加，重建圖像的質(zhì)量越高，但運(yùn)行時(shí)間也越長。改進(jìn)的RAMP算法提高了運(yùn)算速度，重構(gòu)時(shí)間最長不超過0.2 s， PSNR最高可達(dá)34.957 5 dB。

圖3 水印圖像的重建效果比較

實(shí)驗(yàn)二：本文選擇采樣率為0.7在原始載體中嵌入二值水印圖像，對其不做任何形式的攻擊處理，并從中提取出水印信息[8]，實(shí)驗(yàn)結(jié)果如圖4所示。載有水印的圖像很好地被人類視覺系統(tǒng)接受，并沒有引起圖像的失真。分別對嵌入水印的宿主信號、載體圖像和原始水印、提取的水印的性能進(jìn)行評價(jià)，得到PSNR為36.139 6 dB，NC為0.948 1，滿足數(shù)字水印系統(tǒng)的評價(jià)標(biāo)準(zhǔn)，有較好的不可感知性和魯棒性。這說明本文提出的RAMP算法具有較好的重構(gòu)性能，達(dá)到本次實(shí)驗(yàn)的效果。

實(shí)驗(yàn)三：抗攻擊測試。對含有水印的載體圖像分別進(jìn)行JPEG壓縮、高斯低通濾波、疊加高斯噪聲和椒鹽噪聲四種攻擊方式來驗(yàn)證該算法的性能，攻擊后水印的提取結(jié)果如圖5所示。分別通過公式計(jì)算載體圖像的PSNR和水印的NC，與文獻(xiàn)[9]進(jìn)行比較，比較結(jié)果見表1。

通過圖5和表1的結(jié)果可以看出：該算法對嵌入圖像的水印有一定的魯棒性，特別是對JPEG壓縮和高斯低通濾波有較好的魯棒性。

圖4 水印的嵌入與提取

圖5 攻擊下提取的水印

表1 不同算法在攻擊下的性能對比

5 結(jié) 語

考慮到小波變換4個(gè)子帶系數(shù)的不同特性，只對高頻子帶進(jìn)行壓縮感知，結(jié)合RAMP算法完成水印信號的嵌入與提取。通過實(shí)驗(yàn)仿真得到本文方法較之傳統(tǒng)方法圖像重構(gòu)的PSNR平均提高2～4 dB。再加上本文提出的算法不需要預(yù)知稀疏度，算法的通用性得到增強(qiáng)，提高了算法效率[10]。同時(shí)，本文在抵御JPEG壓縮、低通濾波等常見攻擊方面擁有很強(qiáng)的魯棒性，保證了水印的安全性，具有較好的抗提取能力。為提高算法的總體性能，后續(xù)還將尋求更適合數(shù)字水印系統(tǒng)的采樣矩陣和重構(gòu)算法，努力實(shí)現(xiàn)水印的盲提取。

參考文獻(xiàn)

[1] 陳國法，郭樹旭，李楊，等.基于壓縮感知的數(shù)字圖像水印算法[J].現(xiàn)代電子技術(shù)，2012，35（13）：98?104.

[2] 魏豐，梁棟，張成，等.基于壓縮感知觀測值的數(shù)字圖像水印算法[J].安徽大學(xué)學(xué)報(bào)，2013，37（3）：61?68.

[3] LI Zhe?tao， XIE Jing?xiong， YU Zu?guo. Compressed sensing based on best wavelet packet basis for image processing [J]. Journal on Communications， 2013， 8（8）：1947?1950.

[4] 王澤宇.壓縮感知在數(shù)字水印中的應(yīng)用研究[D].北京：北方工業(yè)大學(xué)，2013.

[5] 張敏.基于壓縮感知的數(shù)字水印方法研究[D].武漢：武漢科技大學(xué)，2012.

[6] CHEN Lan?lan. Research on image digital watermarking based on wavelet transform [J]. Energy Procedia， 2011，13： 8527?8533.

[7] 朱延萬，趙擁軍，孫兵.一種改進(jìn)的稀疏度自適應(yīng)匹配追蹤算法[J].信號處理，2012，28（1）：80?86.

[8] 岑翼剛，陳曉方，岑麗輝，等.基于單層小波變換的壓縮感知圖像處理[J].通信學(xué)報(bào)，2010，31（8）：52?55.

[9] 李蘊(yùn)華.壓縮感知框架下基于ROMP算法的圖像精確重構(gòu)[J].計(jì)算機(jī)應(yīng)用，2011，31（10）：2714?2716.

[10] 吳曉歡，劉智威，朱衛(wèi)平.一種壓縮感知中的改進(jìn)的正則化牛頓算法[J].通信學(xué)報(bào)，2013，34（3）：11?13.