基于QR分解的彩色圖像自嵌入全盲水印算法

韓紹程， 張紅穎

（1. 中國(guó)民航大學(xué)基礎(chǔ)實(shí)驗(yàn)中心，天津 300300；2. 中國(guó)民航大學(xué)航空自動(dòng)化學(xué)院，天津 300300）

基于QR分解的彩色圖像自嵌入全盲水印算法

韓紹程1， 張紅穎2

（1. 中國(guó)民航大學(xué)基礎(chǔ)實(shí)驗(yàn)中心，天津 300300；2. 中國(guó)民航大學(xué)航空自動(dòng)化學(xué)院，天津 300300）

針對(duì)彩色圖像的版權(quán)保護(hù)問(wèn)題，基于 QR矩陣分解提出了一種自嵌入全盲水印算法。先將原始圖像的G通道分量進(jìn)行非下采樣剪切波變換，再對(duì)得到的低頻分量分塊QR分解，通過(guò)判斷各子塊R矩陣中第一行元素向量的l1范數(shù)與所有子塊R矩陣第一行元素l1范數(shù)均值之間的大小關(guān)系生成特征水印。然后對(duì)B通道分量DWT變換后的低頻分量進(jìn)行分塊QR分解，并通過(guò)修改該子塊QR分解后R矩陣中第一行最后一列元素來(lái)嵌入特征水印。特征水印的生成和嵌入在兩個(gè)通道內(nèi)獨(dú)立完成，水印檢測(cè)無(wú)需原始載體圖像，算法無(wú)需借助外加水印信息即可完成對(duì)圖像版權(quán)的鑒別。實(shí)驗(yàn)結(jié)果表明，該算法在經(jīng)歷添加噪聲、JPEG壓縮、縮放、剪切和行偏移等常見(jiàn)攻擊時(shí)，具有很強(qiáng)的魯棒性。

數(shù)字水??；自嵌入；全盲檢測(cè)；QR分解；剪切波變換

隨著互聯(lián)網(wǎng)和數(shù)字多媒體技術(shù)的發(fā)展與普及，圖像、音頻和視頻等數(shù)字化作品的非法復(fù)制和篡改也變得相對(duì)容易。數(shù)字水印技術(shù)作為一種可以對(duì)數(shù)字作品版權(quán)進(jìn)行保護(hù)的有效手段，受到廣泛關(guān)注[1]。根據(jù)數(shù)字水印可抵抗攻擊性能的強(qiáng)弱，一般將數(shù)字水印算法分成魯棒數(shù)字水印算法和脆弱數(shù)字水印算法兩大類(lèi)。魯棒圖像水印算法用于對(duì)原始圖像進(jìn)行版權(quán)保護(hù)；脆弱圖像水印算法常用于對(duì)原始圖像進(jìn)行內(nèi)容認(rèn)證[2]。自嵌入脆弱水印算法[3]作為脆弱水印技術(shù)中的一種特殊情況，是在嵌入端提取原始載體圖像的特征產(chǎn)生水印，并將其嵌入到原始載體圖像中以達(dá)到內(nèi)容完整性認(rèn)證的目的。

文獻(xiàn)[4-5]首次提出將脆弱水印技術(shù)中的“自嵌入”思想引入到魯棒水印技術(shù)中來(lái)，實(shí)現(xiàn)了“自嵌入全盲”水印算法，為數(shù)字水印技術(shù)的研究提出了一條新的思路。然而該類(lèi)算法大多以灰度圖像作為測(cè)試載體，原始特征水印的生成雖然借助“零水印[6]”的思想，但沒(méi)有對(duì)載體圖像做任何修改，只是將原始特征水印進(jìn)行“嵌入”時(shí)采用了系數(shù)比較、量化等方式，免不了對(duì)水印檢測(cè)端提取出的特征水印造成影響，這在一定程度上降低了算法的可靠性。

剪切波變換(shearlet transform)[7-8]作為一種新穎的多尺度幾何分析方法，相對(duì)于其他的多尺度分解工具，如曲線波變換(curvelet transform)和輪廓波變換(contourlet transform)等，對(duì)于圖像表征具有更優(yōu)越的性能和更高的計(jì)算效率[9]。目前，已有文獻(xiàn)將剪切波變換應(yīng)用到圖像去噪[10]、圖像融合[11]等領(lǐng)域，文獻(xiàn)[12]首次提出將剪切波變換應(yīng)用到數(shù)字水印技術(shù)中，并取得了較好的實(shí)驗(yàn)效果。此外，由于QR分解相對(duì)于SVD分解運(yùn)算復(fù)雜度低[13]，近幾年也在數(shù)字水印技術(shù)中得到了很好地應(yīng)用[13-15]。

本文將“自嵌入全盲”水印思想，應(yīng)用于彩色圖像的版權(quán)保護(hù)，結(jié)合彩色圖像具有的多通道特性，提出了一種基于QR分解的自嵌入全盲水印算法。利用QR矩陣分解理論，該算法先在載體圖像G分量子圖的非下采樣剪切波變換(non-subsampled shearlet transform, NSST)域下生成特征水印，然后在B分量子圖的小波變換域中完成以上特征水印的嵌入。實(shí)驗(yàn)表明，特征水印的生成和自嵌入算法在經(jīng)歷添加噪聲、JPEG壓縮、放縮、剪切和行偏移等常見(jiàn)攻擊操作后均表現(xiàn)出較強(qiáng)的魯棒性。

1 相關(guān)理論

1.1 剪切波變換

定義1. 對(duì)任意給定的ψ∈L2(R2)，二維具有合成膨脹的仿射系統(tǒng)定義為[8]：

對(duì)于大小為N×N的圖像，可以將其視為有限值序列構(gòu)成的函數(shù)，則可以表示為：其離散剪切波變換主要分成兩步進(jìn)行：①多尺度分解；②方向局部化[16]。NSST是在空域?qū)D像進(jìn)行剪切波變換時(shí)，通過(guò)非下采樣Laplacian金字塔變換與剪切波濾波器的不同組合來(lái)實(shí)現(xiàn)，其變換過(guò)程沒(méi)有對(duì)圖像進(jìn)行下采樣操作，因此對(duì)圖像進(jìn)行NSST后的結(jié)果仍然是一系列與原始圖像大小相同的子帶圖像，使得該變換具有平移不變性[17]。

1.2 圖像塊QR分解

定義 2. 如果實(shí)（復(fù)）非奇異矩陣A能夠化成正交（酉）矩陣Q與實(shí)（復(fù)）非奇異上三角矩陣R的乘積，即：

則稱(chēng)式(2)為A的QR分解。

定理. 設(shè)A是n階實(shí)（復(fù)）非奇異矩陣，則存在正交（酉）矩陣Q和實(shí)（復(fù)）非奇異上三角矩陣R使A有QR分解，且除去相差一個(gè)對(duì)角元素的絕對(duì)值（模）全等于 1的對(duì)角矩陣因子外，分解是唯一的[18]。

由以上定義及定理可知，對(duì)于雙精度圖像矩陣，無(wú)論行數(shù)和列數(shù)是否相等，均可以進(jìn)行QR分解，而且分解具有唯一性。通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，圖像QR分解之后的R矩陣第一行元素?cái)?shù)值較大，集中了圖像的絕大部分能量。文獻(xiàn)[13]指出圖像分塊QR分解后的R矩陣第一行元素中最后一列元素的改動(dòng)對(duì)圖像塊視覺(jué)質(zhì)量的改變影響是最小的，因此本文也在此位置嵌入原始特征水印來(lái)實(shí)現(xiàn)圖像版權(quán)的認(rèn)證。

2 自嵌入全盲彩色圖像水印算法描述

2.1 原始特征水印生成算法

在水印嵌入端，從彩色載體圖像中生成原始特征水印的具體步驟如下[12]：

(1) 將大小為N×N的RGB彩色載體圖像記為I1，從中提取出G通道分量記為I1_G；

(2) 對(duì)I1_G進(jìn)行尺度為d1的NSST，然后從其低頻分量中隨機(jī)抽取一幅子圖記為A，子圖大小為(N×N)/4；

(3) 將A劃分成互不重疊的子塊A(i,j)，子塊大小為n1×n1，子塊個(gè)數(shù)為K×K(K=N/2n1)，對(duì)每子塊進(jìn)行QR分解：

其中，qr表示QR分解，A(i,j)為A中與位置(i, j)相對(duì)應(yīng)的系數(shù)子塊，Q(i,j)和R(i,j)分別為A(i,j)經(jīng)QR分解后的Q矩陣和R矩陣；

(4) 計(jì)算每一個(gè)R(i,j)矩陣第一行元素的 l1向量范數(shù)F(i,j)，并將其掃描成一維向量，記為B(u)，這里u=1,2,…,K2：

其中，r(1,k)為R(i,j)矩陣中第一行元素值，k表示列數(shù)；

(5) 計(jì)算與所有系數(shù)子塊相對(duì)應(yīng)的R(i,j)矩陣中第一行元素l1向量范數(shù)的均值，記為Fave：

(6) 產(chǎn)生長(zhǎng)度為K2的正整數(shù)隨機(jī)序列P，并對(duì)其進(jìn)行升序排序，即[E X]=sort(P)，這里E為排序后的新序列，X用于存放序列P中的元素在E中的位置；

(7) 按照 X提供的映射地址提取與各子塊對(duì)應(yīng)的 B(u)，通過(guò)比較B(u)與 Fave的大小關(guān)系生成原始特征水?。?/p>

2.2 原始特征水印自嵌入算法

完成原始特征水印自嵌入的過(guò)程具體分解如下[13]：

(1) 從載體RGB圖像中提取出B通道分量，記為I1_B，并對(duì)I1_B進(jìn)行d2級(jí)DWT分解，獲得低頻分量I1_B_L；

(2) 將I1_B_L劃分成互不重疊的子塊A′(i,j)，子塊大小為 n2×n2，子塊個(gè)數(shù)為K×K，這里K=N/(n2·2d2)；

(3) 按照式(1)每一個(gè)圖像塊進(jìn)行QR分解：

(4) 通過(guò)修改分解后R′(i,j)矩陣中第一行最后一列元素r′(1,n2)，實(shí)現(xiàn)原始特征水印的嵌入：

① 根據(jù)要嵌入的水印信息，選取不同的修改幅度1T和2T：

其中，D為量化參數(shù)，文中通過(guò)實(shí)驗(yàn)獲得；

② 根據(jù)1T和2T得到C1和C2：

③ 按照以下條件來(lái)確定嵌入水印后r′(1,n2)的值，記為r′′(1,n2)：

④ 用r′′(1,n2)替換r′(1,n2)，并依次進(jìn)行逆QR分解運(yùn)算和IDWT得到含水印的B分量圖像，記為B′；

⑤ 最后將B′分量圖像與之前的R、G分量重新組合，構(gòu)成含水印的彩色圖像I2。

2.3 特征水印提取算法

將I2經(jīng)歷攻擊后的圖像記為I3，水印檢測(cè)時(shí)，從I3中提取特征水印的過(guò)程和從I1中生成原始特征水印的過(guò)程類(lèi)似[12]，步驟如下：

(1) 從I3中提取出I3_G，再將I3_G依次經(jīng)歷尺度為d1的NSST和低頻分量隨機(jī)子圖抽取，進(jìn)而得到不重疊子塊A*(i,j)；

(2) 對(duì)每一個(gè)系數(shù)子塊A*(i,j)進(jìn)行QR分解，并計(jì)算與每一個(gè) R*(i,j)相對(duì)應(yīng)的 F*(i,j)，同時(shí)獲得B*(u)；

(3) 按照式(5)計(jì)算與所有系數(shù)子塊相對(duì)應(yīng)F*(i,j)的均值，記為F*ave；

(4) 通過(guò)比較B*(u)和F*ave的大小關(guān)系提取特征水印：

2.4 認(rèn)證水印提取與版權(quán)鑒別

水印檢測(cè)時(shí)，從經(jīng)歷攻擊后的圖像I3中提取認(rèn)證水印的過(guò)程是特征水印自嵌入過(guò)程的逆過(guò)程，步驟如下[13]：

(1) 先從I3中提取出I3_B，對(duì)I3_B進(jìn)行d2級(jí)DWT分解獲得低頻分量I3_B_L；

(2) 將I3_B_L劃分成互不重疊的子塊A**(i,j)，對(duì)每一個(gè)A**(i,j)進(jìn)行QR分解，得到矩陣R**(i,j)；

(3) 利用式(12)，從矩陣R**(i,j)的第一行第n2列元素中提取認(rèn)證水?。?/p>

(4) 通過(guò)計(jì)算W**和 W*之間的歸一化相關(guān)系數(shù)(normalized correlation, NC)來(lái)鑒別篡改圖像的版權(quán)，NC定義如下：

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

3.1 實(shí)驗(yàn)參數(shù)說(shuō)明

采用Matlab7.6作為實(shí)驗(yàn)仿真環(huán)境，選擇大小為512×512的24位真彩色圖像Lena、Baboon和Tiffany為測(cè)試圖像，如圖1所示。文獻(xiàn)[19]提出采用結(jié)構(gòu)相似度(structural similarity index metric, SSIM)來(lái)評(píng)價(jià)含水印灰度圖像中水印的不可見(jiàn)性，針對(duì)彩色圖像，文獻(xiàn)[13]對(duì)其進(jìn)行了改進(jìn)，將 3種顏色通道的平均值作為最終水印不可見(jiàn)性的度量，文中選擇此方法來(lái)評(píng)價(jià)彩色圖像中水印的不可見(jiàn)性。SSIM越接近于1，表示含水印圖像與原始載體圖像越相似，水印的不可見(jiàn)性越好。實(shí)驗(yàn)中，NSST的分解尺度數(shù)d1=5，剪切矩陣向量為[70 70 100 100 100]，方向向量為[1 1 1 1 1]，進(jìn)行隨機(jī)子圖抽取和特征水印置亂時(shí)使用的密鑰與文獻(xiàn)[12]相同，DWT分解尺度數(shù)d2=1，圖像子塊的大小n1=n2=8，這樣保證K=512/2/8=32，由原始載體圖像B分量生成的原始特征水印的大小為32×32，共1 024位。自嵌入算法中，圖像子塊大小n2、量化步長(zhǎng)D和SSIM三者之間的關(guān)系相互制約又密不可分。這里分別通過(guò)實(shí)驗(yàn)獲得圖 1(a)～(c) 3幅圖像各自的量化步長(zhǎng)，D=100, 160, 100。此時(shí)，與圖1(a)～(c)相對(duì)應(yīng)的SSIM分別為0.972 5、0.980 3和0.978 1，可見(jiàn)原始特征水印自嵌入后，三幅含水印圖像均具有很好的不可見(jiàn)性。

3.2 算法性能評(píng)價(jià)

用峰值信噪比(peak signal noise ratio, PSNR)[20]來(lái)評(píng)價(jià)圖像經(jīng)歷攻擊后視覺(jué)上受損程度，用NC來(lái)衡量算法抵抗各種攻擊的魯棒性。表1給出了從經(jīng)歷不同攻擊后的3幅含水印圖像中提取出的特征水印W*與認(rèn)證水印W**和原始特征水印W之間的NC值關(guān)系。針對(duì)表1中給出的典型攻擊，采用3幅紋理不同的RGB彩色圖像測(cè)試本文算法，當(dāng)測(cè)試圖像嚴(yán)重受損時(shí)，W*與W**之間的NC值和W*與W之間的NC值均接近于1，表明文中特征水印的生成方法和特征水印的自嵌入認(rèn)證方法對(duì)常見(jiàn)攻擊均具有很強(qiáng)的魯棒性。同時(shí)，特征水印生成算法的抗攻擊魯棒性稍好于本文自嵌入認(rèn)證水印算法。

3.3 算法比較

將文獻(xiàn)[4-5]提出的自嵌入全盲水印算法和本文提出的算法進(jìn)行了比較。文獻(xiàn)[4-5]是以灰度圖像為測(cè)試載體，為了比較方便，將其算法移植到彩色圖像中，同樣分別是在G通道分量中生成原始特征水印，在B通道分量中完成特征水印的自嵌入。針對(duì)表1給出的典型攻擊方式，以Lena圖像進(jìn)行了驗(yàn)證，特征水印的大小均為32×32，與本文算法一致，同時(shí)在特征水印生成時(shí)采用本文算法中提到的置亂方法替換原文獻(xiàn)中的水印加密方法。這里，文獻(xiàn)[4]算法自適應(yīng)量化因子選取 0.046，對(duì)應(yīng)含水印圖像的SSIM為0.972 1，與本文算法獲得的含水印圖像的SSIM幾乎相同，視覺(jué)上與原始圖像基本沒(méi)有差異；文獻(xiàn)[5]算法量化步長(zhǎng)同樣使用原文獻(xiàn)中的數(shù)值55，對(duì)應(yīng)含水印圖像的SSIM為0.938 6，在視覺(jué)上與原始載體存在很大變化，這主要是由于對(duì)圖像分塊 SVD分解后的最大奇異值進(jìn)行了修改而引起的。使用 3種算法獲得的含水印圖像如圖 2所示，以 Lena為載體，3種算法的性能對(duì)比結(jié)果如表2所示。

圖1 載體測(cè)試圖像

表1 3幅含水印圖像抗攻擊性能評(píng)價(jià)

圖2 含水印Lena圖像

對(duì)比表2中數(shù)據(jù)可知，本文中采用的特征水印生成方法抵抗攻擊的魯棒性要稍好于文獻(xiàn)[4]和文獻(xiàn)[5]中提到的方法。本文中采用的基于 QR分解的原始特征水印自嵌入方法，在整體性能上稍好于文獻(xiàn)[4]中使用的方法，特別在抵抗 JPEG壓縮方面表現(xiàn)出明顯優(yōu)勢(shì)；與文獻(xiàn)[5]相比，抵抗添加噪聲、較大區(qū)域剪切和行偏移等攻擊的魯棒性明顯提高，抵抗 JPEG壓縮和縮放攻擊的魯棒性稍有遜色，但是，不能忽視的是，采用文獻(xiàn)[5]算法獲得的SSIM值僅為0.938 6，如果想獲得和本文算法具有相同視覺(jué)質(zhì)量的含水印Lena圖像，必須減小量化步長(zhǎng)D，隨著D的不斷減小，該算法在彩色圖像空間中的應(yīng)用將失去意義。綜合以上，本文算法在總體性能上稍強(qiáng)于文獻(xiàn)[4]算法，而明顯優(yōu)于文獻(xiàn)[5]算法。

本文提出的特征水印生成和自嵌入算法對(duì)各種攻擊具有頑強(qiáng)的魯棒性，其原因在于：①NSST可以對(duì)圖像實(shí)現(xiàn)最佳逼近描述，相比于其他變換域方法，圖像NSST分解后的低頻分量具有一定的穩(wěn)定性，而且其大小較原始圖像沒(méi)有減小，這為隨機(jī)子圖抽取和獲得較大尺寸的圖像分塊提供了可能性；②圖像QR分解后R矩陣中第一行元素?cái)?shù)值較大，代表了圖像的大部分能量，其向量范數(shù)在經(jīng)歷各種攻擊后具有較好的不變性；③隨機(jī)子圖抽取和特征水印生成時(shí)使用的置亂方法增強(qiáng)了算法的可靠性和安全性；④水印嵌入在圖像的彩色空間相對(duì)于灰度空間具有更強(qiáng)的隱蔽性；⑤特征水印自嵌入時(shí)，只修改了圖像QR分量后R矩陣中第一行向量的最后一個(gè)元素，通過(guò)實(shí)驗(yàn)可知，該元素的改變對(duì)原始圖像的視覺(jué)質(zhì)量影響最小。

表2 Lena圖像不同算法抗攻擊性能比較

4 結(jié) 束 語(yǔ)

本文利用NSST在圖像表征方面的顯著優(yōu)勢(shì)，結(jié)合彩色圖像的多通道特性，提出了一種基于 QR矩陣分解理論的自嵌入全盲檢測(cè)水印算法。先將從G通道分量NSST后的低頻子圖中提取出原始特征水印，再通過(guò)修改B通道分量小波變換低頻子圖系數(shù)嵌入水印，特征水印的提取和嵌入均沒(méi)有離開(kāi)QR矩陣分解理論，同時(shí)由于兩個(gè)通道相互獨(dú)立，原始特征水印的嵌入不會(huì)在水印檢測(cè)端對(duì)特征水印的提取帶來(lái)影響，有利于更有效地完成圖像的鑒別與保護(hù)。該算法幾乎免疫于JPEG壓縮和尺寸縮放兩種攻擊，并可以有效抵抗添加噪聲、濾波、剪切等常見(jiàn)的攻擊操作。本文首次以彩色圖像為載體，嘗試探索“自嵌入全盲”水印算法，下一步研究將集中在自嵌入水印算法的創(chuàng)新與優(yōu)化方面，例如針對(duì)圖像塊的紋理特征如何自適應(yīng)地選擇量化步長(zhǎng)。

[1] 李旭東. 抗亮度和對(duì)比度攻擊的 DCT域圖像水印算法[J]. 光電子·激光, 2013, 24(6): 1957-1963.

[2] 張憲海, 楊永田. 基于脆弱水印的圖像認(rèn)證算法研究[J].電子學(xué)報(bào), 2007, 35(1): 34-39.

[3] Qin Chuan, Chang Chinchen, Chen Peiyu. Self-embedding fragile watermarking with restoration capability based on adaptive bit allocation mechanism [J]. Signal Processing, 2012, (92): 1137-1150.

[4] 葉天語(yǔ). 自嵌入完全盲檢測(cè)頑健數(shù)字水印算法[J]. 通信學(xué)報(bào), 2012, 33(10): 7-15.

[5] 葉天語(yǔ). DWT-SVD域全盲自嵌入魯棒量化水印算法[J].中國(guó)圖象圖形學(xué)報(bào), 2012, 17(6): 644-650.

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

[7] Guo Kanghui, Labate D. Optimally sparse multidimensional representation using shearlets [J]. SIAM Journal on Mathematical Analysis, 2007, 39(1): 298-318.

[8] Easley G, Labate D, Lim W Q. Sparse directional image representations using the discrete shearlet transform [J]. Applied and Computational Harmonic Analysis, 2008, 25(1): 25-46.

[9] 王 雷, 李 彬, 田聯(lián)房. 基于平移不變剪切波變換的醫(yī)學(xué)圖像融合[J]. 華南理工大學(xué)學(xué)報(bào), 2011, 39(12): 13-19.

[10] Easley G R, Labate D, Colonna F. Shearlet-based total variation diffusion for denoising [J]. IEEE Transactions on Image Processing, 2009, 18(2): 260-268.

[11] Miao Qiguang, Shi Cheng, Xu Pengfei, et al. A novel algorithm of image fusion using shearlets [J]. Optics Communications, 2011, 284(6): 1540-1547.

[12] 韓紹程, 張兆寧, 張玉金. 基于非下采樣剪切波變換和 QR分解的魯棒零水印算法[J]. 光電子·激光, 2012, 23(10): 1957-1963.

[13] Su Qingtang, Niu Yugang, Wang Gang, et al. Color image blind watermarking scheme based on QR decomposition [J]. Signal Processing, 2014, 94: 219-235.

[14] Song Wei, Hou Jianjun, Li Zhaohong, et al. Chaotic system and QR factorization based robust digital image watermarking algorithm [J]. Journal of Central South University of Technology, 2011, 18(1): 116-124.

[15] Naderahmadian Y, Hosseini-Khayat S. Fast watermarking based on QR decomposition in wavelet domain [C]//Intelligent Information Hiding and Multimedia Signal Processing, Darmstadt, Germany, 2010: 127-130.

[16] 劉俊良, 雷 琳, 周石琳. 非下采樣形態(tài)學(xué) shearlet變換: 提高結(jié)構(gòu)細(xì)節(jié)捕捉的圖像表示新方法[J]. 信號(hào)處理, 2014, 30(2): 163-171.

[17] 高國(guó)榮, 許錄平, 馮冬竹. 利用非抽樣 shearlet域GSM 模型進(jìn)行圖像去噪[J]. 武漢大學(xué)學(xué)報(bào): 信息科學(xué)版, 2013, 38(7): 778-782, 812.

[18] 程云鵬. 矩陣論[M]. 西安: 西北工業(yè)大學(xué)出版社, 2003: 203-205.

[19] Wang Zhou, Bovik A C, Sheikh H R, et al. Image quality assessment: from error visibility to structural similarity [J]. IEEE Transactions on Image Processing, 2004, 13(4): 600-612.

[20] 韓紹程, 羅長(zhǎng)杰, 張兆寧, 等. 基于多小波變換和分塊SVD的彩色圖像水印算法[J]. 工程圖學(xué)學(xué)報(bào), 2010, 31(2): 128-133.

Self-Embedding Perfectly Blind Watermarking Algorithm Based on QR Decomposition for Color Images

Han Shaocheng1, Zhang Hongying2

(1. Basic Experimental Center, Civil Aviation University of China, Tianjin 300300, China; 2. Aeronautical Automation College, Civil Aviation University of China, Tianjin 300300, China)

A novel self-embedding perfectly blind watermarking algorithm is proposed based on QR matrix decomposition for copyright protection of digital color images. At first, G channel component of the original image is performed with non-subsampled shearlet transform, then the low-frequency component is performed with block-QR decomposition, the feature watermark is derived by judging the number relationship between vector l1norm from the first row elements of each sub-block′s R-matrix and the mean of vector l1norms from the first row elements of all sub-block′s R-matrix. Secondly, the low-frequency component is obtained after DWT from B channel of original image, then it is performed with block-QR decomposition. The feature watermark is embed by modifying the first row and the least column element of the R-matrix from every blocks. The generation and embedding of feature watermark are done in two channels independently, the watermark can be detected without the original image, and the algorithm can identify an image without the additional watermark information. Experimental results show that the proposed algorithm has strong robustness to resist various common attacks such as adding noise, JPEG compression, scaling, cropping and row shifting.

digital watermarking; self-embedding; perfectly blind detection; QR decomposition; shearlet transform

TN 911.73

A

2095-302X(2015)03-0345-07

2014-09-03；定稿日期：2014-11-24

國(guó)家自然科學(xué)基金資助項(xiàng)目(71171190，61379102，U1433105)；天津市青年科學(xué)基金資助項(xiàng)目(12JCQNJC00600)；中國(guó)民航大學(xué)科研基金資助項(xiàng)目(2012KYE02)；中央高校基本業(yè)務(wù)費(fèi)專(zhuān)項(xiàng)基金資助項(xiàng)目(3122014D030)

韓紹程(1981-)，男，天津人，實(shí)驗(yàn)師，碩士。主要研究方向?yàn)閳D像處理及數(shù)字水印技術(shù)。E-mail：schan@cauc.edu.cn