陳 曉 鋒

(浙江大學(xué)寧波理工學(xué)院信息科學(xué)與工程學(xué)院 浙江 寧波 315000)

0 引 言

在無線移動(dòng)信道(開放信道)上的多媒體數(shù)據(jù)傳輸面臨諸如版權(quán)保護(hù)和有限傳輸帶寬等安全問題。數(shù)字水印一般是隱含著著作權(quán)者的相關(guān)信息，用于版權(quán)歸屬設(shè)定，用戶對(duì)數(shù)據(jù)信息的處理需要首先對(duì)水印代碼進(jìn)行處理。數(shù)字水印具有兩種屬性：永久嵌入屬性和易于讀取屬性，同時(shí)不能被人體視覺輕易感知到，并且嵌入的水印與原始信息之間不存在較為明顯的干擾。為了安全起見，嵌入水印不可被用戶輕易地刪除。雖然水印需要進(jìn)行永久性標(biāo)識(shí)，但相對(duì)而言不可對(duì)原始數(shù)據(jù)信息的訪問產(chǎn)生影響[1-2]。

當(dāng)前的數(shù)字水印研究中已有很多成熟的算法。例如文獻(xiàn)[3]將圖像中常用的數(shù)字水印算法拓展到音視頻等多種格式數(shù)據(jù)文件中，采用的基本方式是利用擴(kuò)頻方式進(jìn)行數(shù)字水印的顯頻分量隱藏，問題是算法簡(jiǎn)單的進(jìn)行移植操作，其算法適應(yīng)性相對(duì)較差。文獻(xiàn)[4]采用傅里葉算法對(duì)版權(quán)圖像的數(shù)字水印問題進(jìn)行研究，將數(shù)字水印按照非遞增序列的形式進(jìn)行嵌入，缺點(diǎn)是算法存在較大的失真問題，水印與原始圖像之間存在較大的干擾問題。文獻(xiàn)[5]采用圖像融合方式進(jìn)行數(shù)字水印的嵌入，其可有效防止用戶的隨意復(fù)制，缺點(diǎn)是水印的嵌入分辨率較差。文獻(xiàn)[6]采用小波算法對(duì)數(shù)字水印問題進(jìn)行研究，提出了具有可變換分辨率算法的數(shù)字水印檢測(cè)方法，缺點(diǎn)是算法的計(jì)算復(fù)雜度過高。文獻(xiàn)[7]設(shè)計(jì)了一種具有互補(bǔ)特性的數(shù)字水印檢測(cè)方法，算法的優(yōu)點(diǎn)是不需要對(duì)圖像的原始數(shù)據(jù)信息進(jìn)行參照，缺點(diǎn)是算法的水印檢測(cè)精度偏低。文獻(xiàn)[8]基于離散余弦變變換對(duì)數(shù)字水印算法的頻域檢測(cè)方案進(jìn)行設(shè)計(jì)，可實(shí)現(xiàn)對(duì)水印圖像的信息進(jìn)行調(diào)制，對(duì)于水印的魯棒攻擊性具有較強(qiáng)的抵抗能力。

上述算法在進(jìn)行數(shù)字水印算法進(jìn)行設(shè)計(jì)過程中雖然針對(duì)某一方面進(jìn)行了針對(duì)性的設(shè)計(jì)，但是其優(yōu)化解具有局限性，是一種局部極優(yōu)的數(shù)字水印檢測(cè)方法。本文針對(duì)數(shù)字水印檢測(cè)算法設(shè)計(jì)了一種最優(yōu)化計(jì)算模型，是一種多準(zhǔn)則極值優(yōu)化問題。近些年，在多準(zhǔn)則極值優(yōu)化問題研究中，已有很多成熟的計(jì)算方法，例如權(quán)值法、灰度值法、熵值法、均衡法等。這些方法對(duì)于多準(zhǔn)則極值問題，具有一定的優(yōu)化效果，但是問題是算法的主觀性過強(qiáng)，特別是算法的收斂性在理論上缺乏足夠的收斂依據(jù)。對(duì)此，本文利用分段線性坐標(biāo)優(yōu)化算法對(duì)多極值目標(biāo)優(yōu)化問題進(jìn)行研究，實(shí)現(xiàn)對(duì)數(shù)字水印檢測(cè)的最優(yōu)化實(shí)現(xiàn)。

從上述研究成果可知，以往的數(shù)字水印研究成果側(cè)重的是算法的實(shí)現(xiàn)，對(duì)于主機(jī)信號(hào)干擾HSI和衰減退化影響檢測(cè)器性能，考慮不是非常細(xì)致。為解決該問題，本文采用擴(kuò)頻SS水印和壓縮感知CS進(jìn)行算法融合，提出了一種優(yōu)化問題模型，以最小化失真(包括CS重建和水印檢測(cè))作為優(yōu)化指標(biāo)，采用檢測(cè)概率為特征的檢測(cè)可靠性約束，進(jìn)行最優(yōu)測(cè)量次數(shù)(提供帶寬)和嵌入強(qiáng)度的優(yōu)化。

1 基于退化壓縮感知的擴(kuò)頻水印檢測(cè)

在原始變換空間Ψ∈RN×N中，將長(zhǎng)度為N的具有零均值和單位方差的獨(dú)立同分布(i.i.d)高斯模式P(∈RN)加到主機(jī)數(shù)據(jù)X∈RN×N中，如傅里葉、小波變換(這里使用)、離散余弦變換(DCT)等。SS水印可表示為：

(1)

圖1 CS-SS水印圖像傳輸和水印檢測(cè)的系統(tǒng)模型

在檢測(cè)器側(cè)接收的CS測(cè)量y可表示為：

(2)

式(2)中的水印檢測(cè)可以被建模為兩個(gè)假設(shè)過程，H0(水印缺失)和H1(水印存在)，如式(3)所示，均服從高斯分布。

(3)

式中:總干擾協(xié)方差矩陣是C可表示為：

C=E[yyH/H0]=E[(hφx+η)(hφx+η)H]=

(4)

式中:E(<φx,η>)=E(<η,φx>)=0；操作符〈p,q〉表示向量p和q的內(nèi)積?；讦咋誋=IM,總干擾協(xié)方差矩陣C可改寫為：

(5)

基于式(5)，在H0和H1情況下，y的概率密度函數(shù)(PDF)可表示為：

(6)

式中：det(·)表示矩陣的行列式。測(cè)試統(tǒng)計(jì)量t(y)可推導(dǎo)為：

(7)

(8)

(9)

(10)

(11)

式中:對(duì)于已知的信道條件(退化函數(shù))和給定的pfa,pd值取決于CS測(cè)量m的數(shù)目和水印強(qiáng)度α，即：

(12)

2 水印檢測(cè)可靠性優(yōu)化

由于CS和水印檢測(cè)過程所導(dǎo)致的失真現(xiàn)象，可表示為如下模型：

(13)

計(jì)算失真(e)的平均“L2”范數(shù)，具體計(jì)算過程如下：

(14)

(15)

(16)

理論上，具有檢測(cè)可靠性約束的優(yōu)化框架可以表述為：

(17)

(18)

(19)

式中:d2=z+α2+M2。將約束優(yōu)化問題轉(zhuǎn)化為無約束優(yōu)化問題，形式為：

(20)

當(dāng)▽d2|M=0時(shí)，可得到d2的最小值，即d2的梯度。因此，可得：

(21)

求解三次方程并放棄虛根，得到最優(yōu)值M*的取值形式為：

M*= (b/2+((-a3/27+b/2)2-a6/729)1/2-

a3/27)1/3-a/3+a2/(9×(b/2+((b/2-

a3/27)2-a6/729)1/2-a3/27)1/3)

(22)

(23)

由此可得到目標(biāo)函數(shù)(d2)的二階導(dǎo)數(shù)，計(jì)算形式如下：

(24)

式中:▽2d2|M是正的，這清楚地表明目標(biāo)函數(shù)的值是最優(yōu)M*值的最小值。因此，在給定的pfa=ξ中，pd的最佳值可以通過在式(13)中代入α*和M*來獲得。

3 分段線性坐標(biāo)優(yōu)化求解

3.1 算法描述

式(19)退化壓縮感知的擴(kuò)頻水印檢測(cè)模型的求解，可采用坐標(biāo)法進(jìn)行優(yōu)化。一般來說，對(duì)于給定的水印檢測(cè)模型網(wǎng)格頂點(diǎn)，主頂點(diǎn)位置上的擴(kuò)頻水印檢測(cè)模型的值是最大的，而位于網(wǎng)匯聚點(diǎn)位置上的點(diǎn)最小，算法計(jì)算過程為：

Step1算法參數(shù)進(jìn)行初始操作，設(shè)定t=1、k=1、i=1。

Step2設(shè)定分段線性坐標(biāo)優(yōu)化頂點(diǎn)νt=(i1,i2,i3),式中i1∪i2∪i3={1,2,3},參數(shù)i是分段線性坐標(biāo)的圖標(biāo)號(hào)。

Step3查找xs=i,xs-1=max{Ns/xs},xs-2=max{Ns/(xs,xs-1)},…,x4=max{Ns/(xs,…,x5)}。采取升序方式對(duì){Ns/(xs,xs-1…,x4)}進(jìn)行序列的排列，且有j1

Step4查找xs=i,xs-1=min{Ns/xs},xs-2=min{Ns/(xs,xs-1)},…,x4=min{Ns/(xs,…,x5)}。采取升序方式對(duì){Ns/(xs,xs-1,…,x4)}進(jìn)行序列的排列，且有j1

Step5對(duì)于構(gòu)建的坐標(biāo)優(yōu)化網(wǎng)絡(luò)，對(duì)f(p1)和f(pst1)進(jìn)行對(duì)比，可對(duì)算法下一步的搜索方向進(jìn)行確定，對(duì)于設(shè)定參數(shù)y*，若滿足f(p1)-y*≤y*-f(p1)，按照從Sink節(jié)點(diǎn)到主頂點(diǎn)的逆向順序進(jìn)行數(shù)字水印檢測(cè)模型的優(yōu)化，并執(zhí)行Step 6。若不滿足，則采取相反的方向進(jìn)行算法目標(biāo)函數(shù)的搜索優(yōu)化，并執(zhí)行Step 7。該算法過程的優(yōu)點(diǎn)是算法的執(zhí)行效率較高。

Step6采取降序方式對(duì)目標(biāo)輸入xk排序，式中k∈Nk,jk>jk-1>…>j1。針對(duì)設(shè)定的序列jk?jk-1?…?j1進(jìn)行坐標(biāo)xk轉(zhuǎn)換，將產(chǎn)生k-1個(gè)坐標(biāo)p2,p3,…,pk，得到算法的上層代碼，然后跳轉(zhuǎn)到Step 8。

Step7采取升序方式對(duì)xk排序，k∈Nk，jk>jk-1>…>j1。針對(duì)設(shè)定的序列jk?jk-1?…?j1進(jìn)行坐標(biāo)xk轉(zhuǎn)換，將產(chǎn)生k-1個(gè)坐標(biāo)pst2,pst3,…,pstk，進(jìn)而可得到算法的下層代碼。

Step8基于求解的坐標(biāo)，形式如下：對(duì)于優(yōu)化目標(biāo)f(pn)=f(pn-1)-Δn-1，根據(jù)函數(shù)取值進(jìn)行排列，式中μ(λ)為碼pn-1中jλ位置，Δn-1=(jn-jn-1)(cn-cμ(n-1))。

Step9按照順序?qū)ο铝袟l件進(jìn)行檢測(cè)：

(a) 若滿足水印檢測(cè)模型條件f(pm)≥y*，則跳轉(zhuǎn)Step 10，繼續(xù)目標(biāo)算法優(yōu)化。

(b) 若滿足水印檢測(cè)模型條件f(pm)

(c) 若滿足水印檢測(cè)模型條件f(pm)

(d) 若滿足t+1≤6，則轉(zhuǎn)Step 2；否則，滿足該約束下對(duì)算法進(jìn)行終結(jié)。

Step10假定k=k+1，若滿足k

Step11采取降序方式對(duì)xk排序，k∈Nk,jk>jk-1>…>j1。針對(duì)設(shè)定的序列jk?jk-1?…?j1進(jìn)行坐標(biāo)xk轉(zhuǎn)換，可產(chǎn)生k-1個(gè)坐標(biāo)q2,q3,…,qk。

Step12基于坐標(biāo)f(qn)=f(qn-1)-δn-1,δn-1=(jn-jn-1)(cn-cμ(n-1))，對(duì)函數(shù)排列值進(jìn)行查找，μ(λ)是函數(shù)目標(biāo)排列碼qn-1中jλ所處的位置區(qū)域。

Step13若滿足水印檢測(cè)模型條件f(qn)≤y*，則保留坐標(biāo)頂點(diǎn)qn，跳轉(zhuǎn)Step 10。

3.2 改進(jìn)坐標(biāo)優(yōu)化算法

對(duì)于多準(zhǔn)則極值優(yōu)化問題，對(duì)分段線性坐標(biāo)優(yōu)化求解算法進(jìn)行設(shè)計(jì)：

Step1對(duì)于目標(biāo)模型輸入，設(shè)定其參數(shù)約束以及目標(biāo)準(zhǔn)則。

Step2生成水印檢測(cè)模型的組合元素。

Step3對(duì)于設(shè)定的約束k，基于改進(jìn)坐標(biāo)法對(duì)組合問題排列進(jìn)行求解，可得到可行解集D*子集Di?X,i∈Nk。

Step4查找可行解D*=D1∩D2∩…∩Dk。

Step7對(duì)于點(diǎn)x∈D*位置進(jìn)行函數(shù)f*值優(yōu)化求解。

Step8對(duì)選取優(yōu)化方向與Step 7所得數(shù)值進(jìn)行對(duì)比，可得目標(biāo)極值。

Step9對(duì)數(shù)字水印檢測(cè)模型的線性分式函數(shù)近求解。終結(jié)算法。

數(shù)字水印檢測(cè)模型的線性分式函數(shù)求解算法如圖2所示。

圖2 算法流程圖

4 實(shí)驗(yàn)結(jié)果

實(shí)驗(yàn)平臺(tái)：CPU i7-6500k 2.9 GHz，RAM 6 GB DDR4 2 400 MHz，系統(tǒng)為Win10旗艦，仿真軟件是MATLAB 2013b。選擇lena(萊娜)測(cè)試圖作為載體，Punjabi Logo作為水印,大小是M/3×N/3，圖像大小M×N=300×300。

4.1 水印檢測(cè)：求最佳參數(shù)設(shè)置

這里設(shè)定的最佳參數(shù)主要是M*和α。利用Daubechies 9/7小波對(duì)256×256，8位/像素萊娜圖像(圖3(a))進(jìn)行小波變換，將其分解為3級(jí)。假定小波系數(shù)的分布遵循GSM模型。水印檢測(cè)在方差小于閾值(1 000)的小波系數(shù)集合中。圖3(b)和圖3(c)分別給出用于78%的CS測(cè)量的水印圖像(PSNR=44.86 dB)和重建水印圖像(PSNR=32.96 dB)。利用CS測(cè)量的最佳值(M*=3 204)，水印檢測(cè)強(qiáng)度(α*=1.5)，可獲得具有零均值和方差0.7(h～CN(0,0.7))的復(fù)正態(tài)分布的瑞利衰落的水印檢測(cè)強(qiáng)度，而具有零均值和方差0.25的加性高斯噪聲(η～N(0,0.25))。

(a) 原始圖像 (b) 水印圖像 (c) 重構(gòu)圖像圖3 實(shí)驗(yàn)圖像

最佳參數(shù)M*和α的選取結(jié)果如圖4所示。

(a) 最優(yōu)CS測(cè)量中的歸一化平均誤差與水印檢測(cè)強(qiáng)度(M*=3 204)

(b) 歸一化平均誤差與CS測(cè)量在最佳水印檢測(cè) 強(qiáng)度上的比較(α*=1.5)圖4 測(cè)量結(jié)果對(duì)比

圖5 檢測(cè)概率對(duì)比

根據(jù)圖5結(jié)果可知，文獻(xiàn)[6]中的工作提出了一種魯棒的和分層的CS水印方案，其中水印的檢測(cè)可以在加密域或解密域中完成。本文算法相對(duì)于文獻(xiàn)[6]算法可實(shí)現(xiàn)顯著性能改進(jìn)，約提升5.43%，因?yàn)楸疚牟捎玫氖亲顑?yōu)檢測(cè)器設(shè)計(jì)。

4.2 抵御攻擊測(cè)試

為了驗(yàn)證所提出的算法對(duì)于數(shù)字水印檢測(cè)過程中高頻分量影響的魯棒性，利用噪聲點(diǎn)對(duì)數(shù)字水印的嵌入過程進(jìn)行模擬攻擊測(cè)試。選取的模擬攻擊噪聲主要有4種類型，分別是：斑點(diǎn)、椒鹽、壓縮和濾波等。其中：斑點(diǎn)噪聲的參數(shù)設(shè)置是mean=0,var=0.4；椒鹽噪聲的參數(shù)設(shè)置是var=0.4；壓縮攻擊過程的參數(shù)設(shè)置是Q=20；濾波攻擊過程的參數(shù)設(shè)置是[50，50]。對(duì)比算法選取FrWT和FrDCT兩種算法，分別見參考文獻(xiàn)[14-15]。實(shí)驗(yàn)結(jié)果如圖6所示。

(a) 斑點(diǎn)攻擊

(b) 椒鹽攻擊

(c) 壓縮攻擊

(d) 濾波攻擊圖6 水印檢測(cè)算法性能對(duì)比

根據(jù)圖6的實(shí)驗(yàn)結(jié)果可知，當(dāng)采用lena圖像作為水印圖像時(shí)，圖6(a)、(b)以及(d)所示的椒鹽攻擊、斑點(diǎn)攻擊、濾波攻擊中，其類似于真實(shí)環(huán)境下的人造邊緣攻擊，是一種高頻的攻擊形式。實(shí)驗(yàn)結(jié)果顯示，本文算法相對(duì)于文獻(xiàn)[14-15]可得到更加清晰的數(shù)字水印圖像，而文獻(xiàn)[14-15]兩種算法檢測(cè)到水印圖像存在非常明顯的高頻干擾，水印圖像不夠清晰，甚至無法進(jìn)行準(zhǔn)確的識(shí)別。在壓縮攻擊算法中，本文算法可以提取到清晰的水印，而FrWT算法提取的水印雖可識(shí)別，但是存在一定的模糊性，效果相對(duì)較差，F(xiàn)rDCT提取到的水印難以識(shí)別。

圖6實(shí)驗(yàn)對(duì)比結(jié)果通過視覺方式對(duì)算法的水印檢測(cè)性能進(jìn)行了驗(yàn)證和對(duì)比，是一種定性分析結(jié)果。下面從定量角度對(duì)算法的水印檢測(cè)性能進(jìn)行實(shí)驗(yàn)驗(yàn)證和對(duì)比，結(jié)果如表1所示。

表1 水印檢測(cè)性能量化對(duì)比

根據(jù)表1實(shí)驗(yàn)結(jié)果可知，對(duì)于選取的集中攻擊類型，本文算法對(duì)于水印檢測(cè)的效果均要優(yōu)于選取的文獻(xiàn)[14-15]兩種對(duì)比算法。根據(jù)以上實(shí)驗(yàn)結(jié)果，本文所提算法在水印圖像檢測(cè)過程中相對(duì)于選取的文獻(xiàn)[14-15]兩種對(duì)比算法具有更高的抗攻擊能力和魯棒性。

5 結(jié) 語

本文提出一種基于搜索優(yōu)化改進(jìn)坐標(biāo)法的退化壓縮感知水印檢測(cè)算法，以最小化失真作為優(yōu)化指標(biāo)，采用檢測(cè)概率為特征的檢測(cè)可靠性約束，進(jìn)行最優(yōu)測(cè)量次數(shù)(提供帶寬)和嵌入強(qiáng)度的優(yōu)化。采用圖論處理方式進(jìn)行坐標(biāo)搜索算法的優(yōu)化過程改進(jìn)，獲得一組能夠符合約束限制的解集，實(shí)現(xiàn)對(duì)水印檢測(cè)算法性能的提升。存在的問題：1) 本文算法側(cè)重于過程的實(shí)現(xiàn)和實(shí)驗(yàn)的驗(yàn)證，在理論分析方面稍顯不足；2) 算法在實(shí)驗(yàn)室條件下得到了很好的驗(yàn)證，對(duì)于真實(shí)應(yīng)用場(chǎng)景的效果，有待于進(jìn)一步驗(yàn)證。這兩方面是今后研究的重點(diǎn)。