抗強(qiáng)剪切攻擊的四元數(shù)彩色圖像零水印算法

何　冰

(1. 渭南師范學(xué)院物理與電氣工程學(xué)院，陜西 渭南 714099；2. 陜西省X射線檢測與應(yīng)用研究開發(fā)中心，陜西 渭南 714099)

抗強(qiáng)剪切攻擊的四元數(shù)彩色圖像零水印算法

何冰1,2

(1. 渭南師范學(xué)院物理與電氣工程學(xué)院，陜西 渭南 714099；2. 陜西省X射線檢測與應(yīng)用研究開發(fā)中心，陜西 渭南 714099)

針對現(xiàn)有彩色圖像零水印算法對存在信息量丟失的各種攻擊(如剪切、行列移除、涂抹等)的不足，提出一種基于非負(fù)矩陣分解(NMF)和四元數(shù)實(shí)矩陣的彩色圖像零水印算法。該算法首先提取彩色圖像的3個(gè)分量(R、G、B)并將各分量進(jìn)行NMF分解；其次將分解后的各系數(shù)矩陣在版權(quán)中心進(jìn)行注冊，同時(shí)將原始彩色圖像使用四元數(shù)實(shí)矩陣來表征，并計(jì)算實(shí)數(shù)字矩陣的Krawtchouk不變矩；最后，在版權(quán)中心注冊Krawtchouk不變矩特征向量；由剩余未剪切部分的彩色圖像和注冊的系數(shù)矩陣共同恢復(fù)出完整的原始圖像。實(shí)驗(yàn)結(jié)果表明：在彩色圖像剪切或涂抹面積超過80%時(shí)，水印信息仍可以正確地被提取；使用四元數(shù)彩色圖像整體法表征圖像在圖像類聚以及水印誤檢率方面的性能都要優(yōu)于彩色圖像三通道分別處理或直接灰度化實(shí)現(xiàn)的效果。

四元數(shù)；零水印；彩色圖像；不變矩；非負(fù)矩陣分解

零水印技術(shù)因其很好解決了傳統(tǒng)水印技術(shù)中不可見性和魯棒性之間的矛盾，近幾年來在信息安全領(lǐng)域中成為研究的熱點(diǎn)。溫泉等[1]在2003年首次提出零水印的概念，同時(shí)利用高階累計(jì)量構(gòu)造零水印特征，通過仿真實(shí)驗(yàn)證明所提出算法的有效性和實(shí)用性。此后零水印方面相關(guān)的文獻(xiàn)和學(xué)術(shù)成果不斷增多[2-4]，馬建糊和何甲興[5]提出了一種基于小波變換的零水印算法，該算法利用一次小波變換后的低頻系數(shù)與二值水印圖像相異或運(yùn)算來構(gòu)造零水印信息，并將其在可信第三方(certificate authority，CA)中心進(jìn)行注冊；陳偉琦和李倩[6]采用雙零水印嵌入方式，同時(shí)利用(singular value dfecomposition，SVD)空間向量來表征圖像的算法，解決了目前零水印算法虛警率高的問題；張春凱等[7]設(shè)計(jì)了一種基于 Krawtchouk矩和NSCT(non-subsampled contourlet transform)變換的魯棒性零水印算法，對圖像進(jìn)行 NSCT分解后，計(jì)算其低頻成分的Krawtchouk低階矩不變量來構(gòu)造特征向量，并仿真實(shí)驗(yàn)來證明算法的有效性；馮銀波和陳善學(xué)[8]將彩色圖像在四元數(shù)域內(nèi)分別進(jìn)行分塊離散傅里葉變換、分塊離散余弦變換以及奇異值分解之后產(chǎn)生兩組二值序列，將二值序列與版權(quán)標(biāo)志信息相結(jié)合構(gòu)造零水印信息，實(shí)驗(yàn)結(jié)果表明，所設(shè)計(jì)的算法對常規(guī)攻擊以及部分組合攻擊具有較強(qiáng)的魯棒性。吳偉民等[9]利用細(xì)胞自動機(jī)變換將圖像分離成低頻子帶和高頻子帶，然后對低頻子帶圖像分塊后進(jìn)行奇異值分解，最后利用分解后的奇異值矩陣零水印信息的設(shè)計(jì)。

以上零水印算法對于常規(guī)的信號處理(如：濾波、JPEG有損壓縮、一般的幾何變換等)具有一定的有效性，而對于信息量受損或丟失的信號處理(如：剪切、涂抹、行列移除)的研究國內(nèi)相關(guān)文獻(xiàn)報(bào)道還比較少。文獻(xiàn)[10]提出一種根據(jù)人類視覺模型，使用混沌系統(tǒng)及提升小波變換的抵抗剪切攻擊的彩色圖像盲數(shù)字水印算法，實(shí)驗(yàn)結(jié)果證明在抗剪切方面的性能較為顯著；文獻(xiàn)[11]提出一種基于 Radon變換不變矩和小波提升的水印算法，實(shí)驗(yàn)結(jié)果證明對幾何變換具有很好的魯棒性；文獻(xiàn)[12]提出一種擴(kuò)頻水印算法，該算法在水印嵌入前，先對水印圖像進(jìn)行面包師變換，使得水印圖像的像素能夠均勻分布于每個(gè)區(qū)域，由此獲得了一組具有可以抵抗剪切攻擊能力的水印序列。文獻(xiàn)[13]利用Patchwork方法同時(shí)結(jié)合相應(yīng)的量化索引調(diào)制，提出一種抗剪切攻擊的數(shù)字水印方法，該方法對剪切及壓縮、噪聲等圖像處理具有良好的效果。文獻(xiàn)[10-13]所提出的水印算法，在抗剪切攻擊方面的性能均得到了一定的提高和改善；然而這三種算法都具有相同的不足之處：①隨著水印嵌入數(shù)據(jù)量的增加，嵌入水印后的圖像其不可感知性隨之下降，即水印的魯棒性和不可見性之間的矛盾是傳統(tǒng)水印算法均需要解決的問題；②這 3種算法對于圖像信息量損失或丟失百分比在 50%以下的剪切攻擊，水印檢測器可以正確地提取出完整的水印信息，但對于某些強(qiáng)剪切攻擊(圖像信息量丟失百分比超過75%)，均無法正確地提取出原始水印信息。近幾年來，相關(guān)學(xué)者將非負(fù)矩陣分解(nonnegative matrix factorization，NMF)應(yīng)用到圖像處理、圖像識別領(lǐng)域。文獻(xiàn)[14]提出一種基于 NMF的譜聚類集成(synthetic aperture radar，SAR)圖像分割算法，使用NMF方法合并相關(guān)圖像分割后的個(gè)體，實(shí)驗(yàn)結(jié)果表明，該方法優(yōu)于 K-means方法、基于 Nystrom逼近的方法和Meta-clustering方法，具有一定的實(shí)用性。雖然NMF[15-16]在圖像處理方面有一定的應(yīng)用，但將其應(yīng)用到零水印算法中的相關(guān)文獻(xiàn)鮮見報(bào)道[17-18]。

通過以上分析，本文提出一種基于NMF和四元數(shù)實(shí)矩陣的彩色圖像零水印算法。該算法根據(jù)NMF部分感知全局的特性，利用剪切后剩余部分圖像及之前保存過的相應(yīng)系數(shù)矩陣可以重構(gòu)出完整的原始圖像，解決了強(qiáng)剪切帶來的問題；同時(shí)引入四元數(shù)理論，將其應(yīng)用到彩色圖像中，由于彩色圖像3個(gè)通道(R、G、B)之間存在著非常緊密的光譜聯(lián)系，若將彩色圖像每一像素點(diǎn)處作為一個(gè)完整的體系來處理，那么每個(gè)像素點(diǎn)三個(gè)通道間的光譜聯(lián)系就會貫穿在整個(gè)圖像的運(yùn)算和處理之中，本文使用四元數(shù)實(shí)矩陣來表征彩色圖像，通過對零水印算法的設(shè)計(jì)，有效地提高了水印在提取階段的檢測精度。

1　非負(fù)矩陣分解

NMF為當(dāng)前的圖像分解提供了一條新的路徑。其基本原理可以描述為：任意給定一個(gè)非負(fù)矩陣 Vn×m，NMF可以尋找出一個(gè)非負(fù)矩陣Wn×r和另外一個(gè)非負(fù)矩陣Hr×m，使得這樣一個(gè)非負(fù)矩陣就分解為兩個(gè)非負(fù)矩陣的乘積。其中，W稱為基矩陣；H稱為系數(shù)矩陣；r稱為基矩陣維數(shù)，且滿足r<

1.1乘性迭代算法

Lee和Seung[15]首次提出了使用乘性迭代的方法來解決NMF問題，采用KL散度，其目標(biāo)函數(shù)形式為：

其中，F(xiàn)表示目標(biāo)函數(shù)；i、j表示二維矩陣的行和列數(shù)。

泊松噪聲模型下的迭代規(guī)則：

其中，α,u,i,k,v 分別表示二維矩陣的行或列數(shù)(即二維矩陣中的坐標(biāo)位置)。

該算法與梯度下降、交替最小平方算法相比，需要更多的迭代步驟，因此收斂速度較慢。因?yàn)镹MF本身不是一種全局最優(yōu)化的方法，所以就有可能得不到全局最優(yōu)的解。

1.2梯度下降算法

迭代算法如下：

其中，f(Wk,Hk)表示二維矩陣的基矩陣和系數(shù)矩陣；εW和εH是步長；?Wf(Wk,Hk)和?Hf(Wk,Hk)表示對目標(biāo)函數(shù)f(Wk,Hk)求偏導(dǎo)。梯度下降算法εW和εH的選擇與其收斂性有極大的關(guān)系。此算法對其初始值特別敏感，如果隨機(jī)取初始值，往往得不到令人滿意的結(jié)果，因此步長的選擇和初始值的選取決定著梯度下降算法的設(shè)計(jì)優(yōu)劣。

1.3交替最小平方算法

該算法的理論基礎(chǔ)在于：在解決最優(yōu)化KL散度目標(biāo)函數(shù)的問題中，如果只針對W矩陣或H矩陣來說其是凸函數(shù)，但同時(shí)對W矩陣和H矩陣而言并不是凸函數(shù)。所以可以在固定一個(gè)矩陣的條件下，對另外一個(gè)矩陣進(jìn)行最小平方計(jì)算。其迭代算法如下：

交替最小平方算法中，為了保證數(shù)據(jù)的非負(fù)性，將負(fù)值強(qiáng)制為0。這在一定程度上增加了數(shù)據(jù)的稀疏性，使該算法的計(jì)算量小了一些。

為了驗(yàn)證以上3種不同NMF算法的收斂性和分解后的重構(gòu)誤差，本文以64×64Lena灰度圖像為例進(jìn)行說明，實(shí)驗(yàn)結(jié)果見圖 1、2。實(shí)驗(yàn)環(huán)境：CPU Intel pentium4 3.2 GHz、內(nèi)存2 G、Windows 7操作系統(tǒng)。

圖1　3種NMF算法的收斂速度曲線(mm：乘性迭代；prob：梯度下降算法；als：交替最小平方算法)

圖2　3種NMF算法的分解后的重構(gòu)誤差曲線(mm：乘性迭代；prob：梯度下降算法；als：交替最小平方算法)

圖3　基矩陣維數(shù)r取不同值的誤差曲線

由圖1～3可知：NMF中，在相同迭代次數(shù)下采用交替最小平方算法其收斂速度效率明顯優(yōu)于乘性算法和梯度算法；然而對分解后的灰度 Lena圖像進(jìn)行重構(gòu)時(shí)，其交替最小平方算法產(chǎn)生的重構(gòu)誤差最大，而采用梯度算法產(chǎn)生的重構(gòu)誤差最小。基矩陣維數(shù)r的取值也非常重要，若r取值較小，計(jì)算運(yùn)行速度會加快，但經(jīng)過NMF分解后的矩陣重構(gòu)誤差會增大；若要減小分解后的重構(gòu)誤差，則需增大基矩陣維數(shù)r的取值，此時(shí)計(jì)算速度則會降低；因此NMF算法中收斂速度和矩陣的重構(gòu)誤差是一對矛盾，近年來關(guān)于NMF的應(yīng)用中，對于灰度圖像或彩色圖像分解后的圖像重構(gòu)及分解效率方面的文獻(xiàn)報(bào)道較少。本文采用乘性NMF算法對二值圖像、灰度圖像分別在基矩陣維數(shù)r取不同值時(shí)，對其圖像重構(gòu)誤差和計(jì)算時(shí)間效率進(jìn)行了仿真實(shí)驗(yàn)分析，這為本文將NMF算法應(yīng)用到彩色圖像零水印中提供了參考依據(jù)，實(shí)驗(yàn)結(jié)果見圖4～7。

圖4　26個(gè)大寫字母組成的二值圖像測試樣本(單個(gè)字母大小128×128)

圖5　26個(gè)大寫字母二值圖像的重構(gòu)結(jié)果

圖6　字母E的重構(gòu)結(jié)果

圖7　灰度圖像baboon的重構(gòu)結(jié)果

使用峰值信噪比PSNR、均方誤差MSE及絕對差ε對原始圖像和重構(gòu)后的圖像進(jìn)行度量，則：

式(5)～(8)中，I(x,y)和I′(x,y)分別表示原始圖像和重構(gòu)后的圖像在(x,y)處的灰度值。

以上實(shí)驗(yàn)結(jié)果表明對大小為 128×128的二值圖像，基矩陣維數(shù)r>10時(shí)；大小為256×256的灰度圖像，基矩陣維數(shù)r>20時(shí)，其圖像重構(gòu)誤差和計(jì)算時(shí)間效率均達(dá)到了理想效果(對于一幅256×256的圖像，重構(gòu)后的PSNR大于20 dB，及重構(gòu)時(shí)間小于 10 ms在圖像識別中是完全可以接受的)，那么將其應(yīng)用到彩色圖像零水印中其不可見性和水印的嵌入、檢測計(jì)算時(shí)間也是完全可以滿足的。

2　四元數(shù)彩色圖像表示

目前的彩色圖像數(shù)字水印算法，大部分處理過程都是先將彩色圖像轉(zhuǎn)換為單通道顏色模式即灰度圖像，再基于灰度的方法設(shè)計(jì)水印的嵌入和提取操作。這種方法損失了彩色圖像大量的色彩信息，同時(shí)又沒有利用到圖像各顏色通道的空間關(guān)系信息，影響到水印檢測時(shí)的精度。四元數(shù)是一般復(fù)數(shù)的推廣，其相關(guān)理論已經(jīng)成功應(yīng)用于彩色圖像去噪、圖像分割以及相關(guān)模式識別領(lǐng)域。采用四元數(shù)表示彩色圖像其優(yōu)勢在于：在多維的空間上將彩色圖像三個(gè)通道的像素作為一個(gè)矢量進(jìn)行整體的處理，而非將三個(gè)通道分別處理或直接進(jìn)行灰度化，考慮了各個(gè)彩色圖像每個(gè)通道之間的相關(guān)性，將四元數(shù)理論應(yīng)用于彩色圖像零水印算法中，有助于對圖像進(jìn)行更精確的表述。

一般一個(gè)四元數(shù)q可以表示為：

其中a,b,c,d∈R；i,j,k分別代表四元數(shù)的 3個(gè)虛部的單位向量，并且滿足以下的運(yùn)算法則：

其中，a為四元數(shù)q的實(shí)部；i×b+j×c+k×d 為q的虛部；令q的實(shí)部為0，3個(gè)虛部分量分別表示彩色圖像R、G、B三通道即：a=0,b=R, c=G,d=B ，則一幅彩色圖像IRGB(m,n)可以表示為：

任意的一個(gè)四元數(shù)q=a+ib+jc+kd=(a+ib)+(c+id)j=u+vj=u+j×conj(v)其矩陣形式為：

其中，conj表示四元數(shù)的共軛；u,v,conj(u),conj(v)對應(yīng)的矩陣形式為：conj(u)=[a?b;ba]；conj(v)=[c?d;dc]，將4個(gè)矩陣形式代入到Q矩陣中，得到四元數(shù)實(shí)矩陣W的表示形式：

根據(jù)式(12)RGB彩色圖像的實(shí)矩陣表示為：

圖 8～9給出了四元數(shù)實(shí)矩陣整體彩色圖像表示法及彩色圖像轉(zhuǎn)換為灰度圖像表示法的圖像聚類實(shí)驗(yàn)結(jié)果(以 Radon變換不變矩和 Krawtchouk不變矩為特征描述質(zhì)，Radon變換不變矩詳見文獻(xiàn)[11])，實(shí)驗(yàn)結(jié)果說明，實(shí)矩陣整體彩色圖像表示法在類間距的聚類效果明顯優(yōu)于彩色圖像轉(zhuǎn)換為灰度圖像表示法(類間距越大、類內(nèi)距越小，聚類效果越好)，為彩色圖像零水印采用四元數(shù)實(shí)矩陣整體表征圖像提供了依據(jù)。

圖8　16幅不同的彩色樣本圖像

圖9　以Radon不變矩均值和Krawtchouk不變矩均值為特征值下的16幅不同彩色圖像的聚類圖

3　本文算法的實(shí)現(xiàn)過程

3.1零水印生成及注冊過程框圖

零水印生成及注冊過程具體步驟如圖10所示：

步驟1. 分別提取原始彩色載體圖像的R、G、 B三分量并對其進(jìn)行非負(fù)矩陣分解；

步驟2. 將NMF分解后的三分量系數(shù)矩陣H分別在第三方認(rèn)證機(jī)構(gòu)(CA中心)進(jìn)行注冊；

步驟3. 使用四元數(shù)實(shí)矩陣整體法來表征原始彩色載體圖像，得到四元數(shù)數(shù)字矩陣；

步驟 4. 提取四元數(shù)數(shù)字矩陣的 Krawtchouk低階不變矩(Krawtchouk低階不變矩詳見文獻(xiàn)[7])作為特征向量并在CA中心進(jìn)行注冊；

步驟 5. CA中心對特征向量及用戶信息加蓋時(shí)間戳后，宣布原始彩色載體圖像已在版權(quán)保護(hù)之下。

3.2零水印提取與檢測過程

給定原始圖像矩陣V，其矩陣大小為m×n，經(jīng)過NMF分解后得到基矩陣W(大小為n×r)和系數(shù)矩陣H(大小為r×m)，r為基矩陣維數(shù)。具體模型如下：

圖10　水印生成及注冊過程

其中，vi為V的第i列；hi為H的第i列；由式(15)可知原始圖像矩陣的每一列與分解后的系數(shù)矩陣每一列一一對應(yīng)。

根據(jù)以上公式推論可知：原始圖像矩陣 V遭受剪切攻擊時(shí)，仍可以利用矩陣 V中剩余未被剪切部分恢復(fù)出完整基矩陣 W；然后根據(jù)保留的系數(shù)矩陣H，重構(gòu)原始圖像矩陣。為使基矩陣W存在唯一的解，系數(shù)矩陣V的最小秩應(yīng)等于NMF分解的維數(shù)，即T≥r。

零水印提取和檢測過程具體步驟如圖11所示：

步驟1. 在第三方版權(quán)認(rèn)證中心(CA)按照時(shí)間戳提取保存的系數(shù)矩陣HR,HG,HB。

步驟2. 對待檢測彩色圖像是否受到剪切或涂抹攻擊進(jìn)行預(yù)判，若遭受剪切或涂抹攻擊則從未受損彩色局部圖像提取VR,VG,VB矩陣(原始彩色圖像R、G、B三分量的局部矩陣)，否則執(zhí)行步驟4。

步驟3. 根據(jù)VR,VG,VB數(shù)字矩陣，從系數(shù)矩陣HR,HG,HB得到與之對應(yīng)的局部矩陣,,，然后按照式(17)、(14)完成對受損彩色圖像的重構(gòu)。

步驟4. 將待檢測彩色圖像或重構(gòu)圖像用四元數(shù)實(shí)矩陣表示，計(jì)算該實(shí)數(shù)字矩陣 10個(gè)低階Krawtchouk不變矩，并將其作為特征向量B。

步驟5. 求出特征向量B與零水印算法注冊階段在CA中心的特征向量A的絕對差之和，即：

若d≥ε并且時(shí)間戳與CA中心提供的信息不符(ε為經(jīng)驗(yàn)閾值系數(shù)，實(shí)驗(yàn)中取值為0.01)則驗(yàn)證結(jié)束，說明待檢測圖像中不含有水印信息；否則，待檢測彩色圖像中水印信息存在。

圖11　水印提取過程

4　實(shí)驗(yàn)結(jié)果和分析過程

為了驗(yàn)證所設(shè)計(jì)零水印算法的有效性，本文以Baboon彩色圖像(大小256×256)為宿主圖像(見圖12)，通過Matlab 8.0仿真軟件完成以下3組實(shí)驗(yàn)結(jié)果。實(shí)驗(yàn) 1是原始彩色圖像遭受常見典型攻擊后的實(shí)驗(yàn)結(jié)果，見表1；實(shí)驗(yàn)2是對Baboon彩色圖像進(jìn)行不同比例的剪切和涂抹攻擊后的實(shí)驗(yàn)，結(jié)果見表2；實(shí)驗(yàn)3是為了驗(yàn)證對于其他圖像是否會出現(xiàn)誤判的情況所作的誤檢率驗(yàn)證實(shí)驗(yàn)(采用四元數(shù)整體法和直接灰度化方法)，宿主圖像采用圖8的16幅圖像，實(shí)驗(yàn)結(jié)果見圖13，由實(shí)驗(yàn)結(jié)果可知使用四元數(shù)整體法的彩色圖像零水印算法相比直接灰度化或單通道(R、G、B取其中一個(gè)通道)的算法在零水印檢測以及誤檢率方面性能都有一定的提高(采用四元數(shù)整體法其絕對差d的數(shù)值明顯高于直接灰度化處理的方法，說明其區(qū)分度相對較大，水印檢測時(shí)不易產(chǎn)生誤判；其中圖 8中第3行第1列的彩色圖像采用直接灰度化的方法會導(dǎo)致誤檢，而四元數(shù)方法其區(qū)分度依然很大)。

圖12　Baboon原始彩色圖像

表1　常見典型攻擊測試結(jié)果

表2　抗剪切、涂抹攻擊測試結(jié)果

圖13　采用四元數(shù)整體法和直接灰度化方法的誤檢率曲線圖

5　結(jié)　論

為了提高零水印算法抵抗存在信息量丟失的各種攻擊(剪切、涂抹、行列移除等)的能力，本文將NMF和四元數(shù)理論相結(jié)合構(gòu)造出一種抗剪切攻擊的零水印算法。算法利用NMF的局部感知全局的特性，對受損后的圖像進(jìn)行重構(gòu)，并通過對二維圖像大量的仿真實(shí)驗(yàn)尋找非負(fù)矩陣分解后基矩陣維數(shù)r的選擇及迭代次數(shù)之間的折中；由實(shí)驗(yàn)結(jié)果可知使用四元數(shù)整體法的彩色圖像零水印算法相比直接灰度化或單通道(R、G、B取其中一個(gè)通道)的算法在零水印檢測以及誤檢率方面性能都有一定的提高。下一步的重點(diǎn)是研究NMF算法和四元數(shù)整體表示圖像在計(jì)算量上的相關(guān)快速算法，以滿足實(shí)時(shí)性的要求，從而進(jìn)一步對零水印系統(tǒng)進(jìn)行優(yōu)化。

[1] 溫泉, 孫錟鋒, 王樹勛. 零水印的概念與應(yīng)用[J].電子學(xué)報(bào), 2003, 31(2): 214-216.

[2] 吳偉民, 丁冉, 林志毅, 等. 基于混沌的醫(yī)學(xué)圖像篡改定位零水印算法[J]. 計(jì)算機(jī)應(yīng)用研究, 2014, 31(12): 3685-3688.

[3] 曲長波, 楊曉陶, 袁鐸寧. 小波域視覺密碼零水印算法[J]. 中國圖象圖形學(xué)報(bào), 2014, 19(3): 367-371.

[4] 謝勇, 張金龍, 張雯. 一種基于奇異值的抗打印/掃描的彩色圖像零水印方案[J]. 包裝學(xué)報(bào), 2014, 16 (3): 11-13.

[5] 馬建糊, 何甲興. 基于小波變換的零水印算法[J]. 中國圖象圖形學(xué)報(bào), 2007, 12(4): 582-585.

[6] 陳偉琦, 李倩. 基于DWT-SVD的圖像雙零水印算法[J]. 計(jì)算機(jī)工程與科學(xué), 2014, 36(10): 1992-1996.

[7] 張春凱, 楊德志, 叢佩麗. 基于 krawtchouk矩和NSCT變換的魯棒性零水印算法[J]. 計(jì)算機(jī)應(yīng)用與軟件, 2014, 31(3): 282-285.

[8] 馮銀波, 陳善學(xué). 基于四元數(shù)域的彩色圖像雙重零水印算法[J]. 計(jì)算機(jī)應(yīng)用與軟件, 2014, 31(9): 264-267.

[9] 吳偉民, 丁冉, 林志毅, 等. 基于細(xì)胞自動機(jī)與奇異值分解的零水印算法[J]. 計(jì)算機(jī)應(yīng)用, 2014, 34(6): 1690-1693.

[10] 趙玉霞. 基于混沌系統(tǒng)與提升小波的抗剪切攻擊的彩色圖像盲水印算法[J]. 工程圖學(xué)學(xué)報(bào), 2010, 24(4): 216-219.

[11] 何冰. 基于Radon變換不變矩和小波提升的水印算法[J]. 計(jì)算機(jī)與數(shù)字工程, 2011, 39(2): 125-128.

[12] 張鑫, 徐光憲, 付曉. 基于面包師變換的抗剪切擴(kuò)頻水印算法研究[J]. 計(jì)算機(jī)應(yīng)用研究, 2012, 29(6): 2246-2248.

[13] Lin C, Pan J S, Liao B Y. Parity modulationbased watermarking resisting to cropping [C]//Proc of ICICIC 2007. Kumamoto: IEEE, 2007: 294.

[14] 鄧曉政, 焦李成, 盧山. 基于非負(fù)矩陣分解的譜聚類集成 SAR圖像分割[J]. 電子學(xué)報(bào), 2011, 12(39): 2906-2909.

[15] Lee D D, Seung H S. Learning the parts of objects by non-negative matrix factorization [J]. Nature, 1999, 401: 788-791.

[16] Zhang T P, Fang B, Tang Y Y, et al. Topology preserving nonnegative matrix factorization for face recognition [J]. IEEE Trans on Image Processing, 2008, 17(4): 574-584.

[17] 韓紹程, 張兆寧. 一種基于子采樣的混合域穩(wěn)健零水印算法[J]. 圖學(xué)學(xué)報(bào), 2013, 34(3): 21-24.

[18] 韓紹程, 張紅穎. 基于 QR分解的彩色圖像自嵌入全盲水印算法[J]. 圖學(xué)學(xué)報(bào), 2013, 34(3): 346-351.

A Zero Color Image Watermarking Resisting Strong Cropping Attacks Using Quaternion Representation

He Bing1,2
(1. Department of Physics and Electronic Engineering, Wei Nan Normal University, Weinan Shaanxi 714099, China; 2. Center of X Ray Detection and Application of ShaanXi, Weinan Shaanxi 714099, China)

For some zero color image watermarking methods are insufficient for several attacks of information loss, a zero color image watermarking algorithm based on nonnegative matrix factorization and real matrix of quaternion is proposed. Firstly, three components (R, G, B) are extracted from color image and each component is decomposed. Secondly, the coefficient matrix decomposed is registered in the copyright center and the real matrix of quaternion is used to represent the original color image, also the krawtchouk invariant moments of the real matrix is calculated. Finally, some krawtchouk invariant moments are used to design and construct zero watermarking information. The experimental results show that original color image under strong shear can be recovered by residual shear part. When the area of shearing color image exceed 80%, watermark information can still be extracted correctly. This performance in image clustering and watermark error detection rate by using the quaternion to represent the whole color image is better than usingeach of three channels of color image processing or gray image directly.

quaternion; zero watermarking; color image; invariant moments; nonnegative matrix factorization

TP 391

A

2095-302X(2015)06-0909-11

2015-05-05；定稿日期：2015-09-08

陜西省教育廳科研計(jì)劃資助項(xiàng)目(14JK1248)；渭南師范學(xué)院第二批特色學(xué)科建設(shè)資助項(xiàng)目(14TSXK06)；渭南市科研發(fā)展計(jì)劃資助項(xiàng)目(2015KYJ-2-1)；渭南師范學(xué)院科研計(jì)劃資助項(xiàng)目(15YKS010)

何冰(1982–)，男，陜西合陽人，講師，碩士。主要研究方向?yàn)橹悄軋D像數(shù)據(jù)處理、信息安全研究。E-mail：hebing126@126.com