崔漢國， 劉健鑫， 李 彬

(1. 海軍工程大學(xué)動力工程學(xué)院，湖北 武漢430033；2. 海軍蚌埠士官學(xué)校機電系，安徽 蚌埠 233012)

基于DD-DT CWT和SIFT的體數(shù)據(jù)數(shù)字水印算法

崔漢國1， 劉健鑫2， 李 彬1

(1. 海軍工程大學(xué)動力工程學(xué)院，湖北 武漢430033；2. 海軍蚌埠士官學(xué)校機電系，安徽 蚌埠 233012)

為了確保三維體數(shù)據(jù)數(shù)字模型信息在網(wǎng)絡(luò)傳輸過程中的安全性、完整性，對體數(shù)據(jù)數(shù)字水印算法進行研究。提出基于雙密度雙樹復(fù)小波變換的體數(shù)據(jù)數(shù)字水印嵌入算法，算法提高了水印的嵌入容量，在兼顧水印不可見性的同時提高了水印抵抗壓縮、噪聲等攻擊的能力；提出基于尺度不變特征變換的體數(shù)據(jù)數(shù)字水印提取算法，實現(xiàn)了水印嵌入和提取時的同步性，提高了水印抵抗旋轉(zhuǎn)、縮放等幾何攻擊的能力。實驗證明：算法不可見性好，實現(xiàn)速度快，在盲檢測下能夠抵抗壓縮、濾波、噪聲、平移、旋轉(zhuǎn)、縮放、剪切等攻擊，具有較好的魯棒性。關(guān) 鍵 詞：數(shù)字水??；體數(shù)據(jù)；雙密度雙樹復(fù)小波變換；尺度不變特征變換；奇異值分解中

三維體數(shù)據(jù)模型應(yīng)用廣泛，在網(wǎng)絡(luò)環(huán)境下實現(xiàn)三維體數(shù)據(jù)信息的安全性保護、完整性檢測已成為目前亟需解決的問題。文獻[1]基于三維離散余弦變換(three dimensional discrete cosine transform，3D-DCT)技術(shù)和DCT雙極性系數(shù)量化技術(shù)提出了一種三維體數(shù)據(jù)模型的魯棒盲提取數(shù)字水印算法，實現(xiàn)了盲提取，能夠抵抗剪切、加噪、濾波等常見攻擊，但經(jīng)實驗驗證依然難以抵抗較大強度的幾何攻擊。文獻[2]對離散小波變換(discrete wavelet transform，DWT)域內(nèi)的數(shù)字水印算法進行了研究，但由于DWT不具備平移不變性、缺乏方向選擇性，導(dǎo)致載體數(shù)據(jù)在嵌入水印后系數(shù)之間的能量分布會發(fā)生較大變化，降低了水印的性能。文獻[3-5]結(jié)合量化調(diào)制策略與雙樹復(fù)小波變換(dual tree complex wavelet transform，DT CWT)技術(shù)，將水印嵌入到關(guān)鍵熵圖像區(qū)域中，增強了數(shù)字水印抵抗壓縮、噪聲等攻擊的能力，但是存在水印嵌入容量有限與無法抵抗幾何攻擊的不足。為了提高水印抵抗平移、旋轉(zhuǎn)等幾何攻擊的能力，人們提出基于特征點的第二代數(shù)字水印方法，常用的特征點主要有墨西哥帽小波特征點[6]、SUSAN特征點[7]、Harris-Laplace特征點[8]等，但只能抵抗單一的幾何變換攻擊，不能抵抗聯(lián)合幾何變換攻擊。

本文針對三維體數(shù)據(jù)模型自身的特點，提出基于雙密度雙樹復(fù)小波變換(double-density dual-tree complex wavelet transform，DD-DT CWT)的體數(shù)據(jù)數(shù)字水印嵌入算法和基于尺度不變特征變換(scale invariant feature transform，SIFT)的體數(shù)據(jù)數(shù)字水印提取算法。

1 基于 DD-DT CWT的體數(shù)據(jù)數(shù)字水印嵌入算法

1.1 三維體數(shù)據(jù)模型的DD-DTCWT

DD-DT CWT[9]結(jié)合了雙密度小波變換(double density discrete wavelet transform，DD DWT)和DT CWT的特點，具有平移不變性、抗混疊特性、不隨尺度變化及良好的方向性等特點。本文將DD-DT CWT應(yīng)用到體數(shù)據(jù)數(shù)字水印的嵌入過程中，將傳統(tǒng)小波理論信息描述的主方向提高到 16個，使在同一尺度內(nèi)相鄰小波間頻帶間隔更小，并且每個主方向有兩個小波描述，對細節(jié)信息描述更加詳細，可以更精確地檢測到載體圖像的紋理和邊緣，提高了載體圖像的分解與重構(gòu)精度，從而能更好地提取原始載體圖像的信息。

本文在對載體圖像進行分塊的基礎(chǔ)上，利用DD-DT CWT處理后產(chǎn)生的4個低頻子帶和32個高頻子帶，在多個層次、多個方向上嵌入水印信息，以提高水印的嵌入容量，并選擇小波變換后±45°方向子帶系數(shù)嵌入水印信息，以實現(xiàn)在提高水印的嵌入容量、兼顧水印不可見性的同時提高水印抵抗壓縮、噪聲等一般信號處理的能力。

1.2 基于DD-DT CWT的體數(shù)據(jù)數(shù)字水印嵌入算法

本文算法步驟如下：

(1) 對三維體數(shù)據(jù)模型進行預(yù)處理，離散采樣三維體數(shù)據(jù)信息，得到二維斷層切片圖像，并對圖像進行灰度化處理，得到灰度圖像序列。

(2) 對載體灰度圖像進行分塊，設(shè)每塊子圖f0i的大小為N1×N1(N1≥8)，通過N1控制水印的嵌入容量和實現(xiàn)效率。由于在低頻子帶中加入水印信息容易引起載體圖像的較大失真，在高頻子帶中加入水印信息無法抵抗壓縮、濾波、噪聲等攻擊，而±45°子帶系數(shù)的變化對重構(gòu)圖像的視覺影響較小，且該方向高頻子帶的能量較其他方向高頻子帶高，信息在壓縮、濾波等后不易損失，因此，本文在對f0i進行4級DD-DT CWT時，將水印信息嵌入第4層±45°方向的小波高頻子帶Mk。

(3) 對Mk進行奇異值分解(singular value decomposition，SVD)，即Mk=USVT，得到兩個正交矩陣U、V及一個對角矩陣S，其中，奇異值矩陣S=diag(λ1,λ2,λ3,…,λn)(λ1≥λ2≥λ3≥…≥λn)是非負對角陣。為了避免水印嵌入后原始載體圖像失真，本文將水印嵌入奇異值矩陣S從λ2開始的R個奇異值中，產(chǎn)生新的矩陣Sw，即Sw=S+αw2(l)，其中，常數(shù)α ＞0是水印強度因子，用于調(diào)節(jié)水印的疊加強度，α 越大水印魯棒性越好，但是原始圖像的失真越明顯、水印的透明性越差。本文α 取經(jīng)驗值0.05。

(4) 進行奇異值逆變換得到新的含有水印信息的復(fù)小波中頻子帶M′k= USwVT。

綜上，DNV對DP-ER附加標志電力系統(tǒng)的要求為：DP-ER附加標志電力系統(tǒng)應(yīng)采用母聯(lián)閉合型母線結(jié)構(gòu)、n+1和n-1冗余設(shè)計原則以及冗余發(fā)電機熱備用保護模式原則。

(5) 對M′k進行雙密度雙樹復(fù)小波逆變換，得到新的含有水印信息的灰度斷層載體圖像序列。

(6) 對含有水印信息的二維斷層切片圖像進行三維重建，得到含水印信息的三維體數(shù)據(jù)模型。

2 基于SIFT的體數(shù)據(jù)數(shù)字水印提取算法

2.1 基于SIFT的體數(shù)據(jù)水印重定位

首先構(gòu)造灰度斷層載體圖像的高斯差分尺度空間(difference of gaussian，DOG)，檢測尺度空間中的極值點并組成一個具有尺度不變性的SIFT候選特征點集合，再從中選出對比度較高的非邊緣點作為SIFT特征點，以增強匹配的穩(wěn)定性和抵抗噪聲的能力，然后確定特征點的方向，實現(xiàn)特征點的旋轉(zhuǎn)不變性，并生成特征點描述符，使其具有亮度不變性，根據(jù)相似性判定度量實現(xiàn)基于SIFT的載體斷層載體圖像特征點的匹配，然后以匹配點對的位置參數(shù)為參照信息，計算出載體圖像遭受攻擊的變化參數(shù)，并進行反變換失真校正，實現(xiàn)對經(jīng)過旋轉(zhuǎn)、縮放和平移攻擊后載體圖像的校正，完成基于SIFT的水印重定位。

2.2 基于SIFT的體數(shù)據(jù)數(shù)字水印提取算法

具體算法步驟如下：

(2) 對含水印載體數(shù)據(jù)進行匹配與校正。

(3) 按照與水印嵌入時相同的規(guī)則對載體灰度圖像進行分塊，對每塊子圖f1i進行4級DD-DT CWT，從第4層±45°方向的小波高頻子帶中提取小波系數(shù)。

(4)對中頻子帶小波系數(shù)進行奇異值分解，得到一組呈單調(diào)遞減排列的奇異值序列S′，

(5) 與原始奇異值序列S進行比較，并提取水印。

在實際的應(yīng)用中，提取出的水印值可能會與嵌入的水印值不完全相同，而在大多數(shù)情況下，需關(guān)注的是水印信息是否存在及能否有效提取，這就需要給出一個判決標準，即閾值。本文采用歸一化互相關(guān)系數(shù) NC定量評價和分析提取水印與原始水印的相似度。

其中，M1和M2代表水印圖像的尺寸，w(i,j)和w′(i,j)分別表示原始水印和提取的水印點(i,j)處的灰度值。若NC大于閾值，則判定待檢測圖像中有相似水印存在且水印提取有效，否則則判定待檢測圖像中不含有水印或水印提取無效。本文閾值取經(jīng)驗值0.6。

3 實驗結(jié)果分析

為了測試和分析本文算法的性能，本文實驗以Stenford數(shù)據(jù)庫中的CT數(shù)據(jù)和二維圖像作為原始三維體數(shù)據(jù)模型和原始水印(圖1所示)，采用本文算法嵌入水印后的三維體數(shù)據(jù)模型如圖2所示。

圖2 嵌入水印后體數(shù)據(jù)模型

實驗結(jié)果分析：

(1) 不可見性測試。本文采用峰值信噪比(peakSignal to noise ratio，PSNR)衡量嵌入水印后的灰度斷層切片圖像的質(zhì)量。比較圖1和圖2可以看出，嵌入水印前后的三維體數(shù)據(jù)模型不存在較大視覺差異，對其中的50幅水印灰度斷層切片圖像進行PSNR值量化測試，得到的PSNR值在30 dB與50 dB之間浮動，表明增加水印嵌入容量后，圖像依然保持了較好的質(zhì)量，本文水印算法滿足水印的不可見性條件。

(2) 抗攻擊測試。表1給出了本文算法與對比算法在抗壓縮攻擊、抗高斯濾波攻擊、抗椒鹽噪聲攻擊、抗幾何變換攻擊(剪切、平移、縮放、旋轉(zhuǎn))、抗圖像處理聯(lián)合攻擊和抗幾何變換聯(lián)合攻擊等方面的實驗結(jié)果，表明本文算法具有較好的魯棒性且優(yōu)于對比算法。

表1 攻擊實驗的相關(guān)值表

4 結(jié) 論

本文算法針對三維體數(shù)據(jù)模型自身的特點，提高了水印的不可見性和水印的嵌入容量，具有更好的抗壓縮、高斯濾波、椒鹽噪聲、幾何變換(剪切、平移、縮放、旋轉(zhuǎn))、圖像處理聯(lián)合和幾何變換聯(lián)合攻擊的能力，實現(xiàn)了體數(shù)據(jù)水印信息嵌入和提取時的同步性，具有較好的魯棒性。

[1] 劉 旺, 趙春暉. 體數(shù)據(jù)模型盲提取數(shù)字水印嵌入算法[J]. 電子測量與儀器學(xué)報, 2009, 23(10)： 13-18.

[2] 車生兵, 黃 達, 李 光. 基于視覺特性的半脆弱水印算法[J]. 通信學(xué)報, 2007, 28(10)： 134-140.

[3] 蘭紅星, 陳松喬, 胡愛娜, 等. 基于小波域的第二代數(shù)字水印算法的研究[J]. 電子學(xué)報, 2007, 35(9)： 1799-1803.

[4] 劉金華, 佘 堃. 基于關(guān)鍵熵的雙樹復(fù)小波域盲圖像水印算法[J]. 光電子·激光, 2011, 22(5)： 757-762.

[5] Kingsbury N．Complex wavelets forShift invariant analysis and filtering ofSignals [J]．Applied and Computational Harmonic Analysis, 2001, 10(3)： 234-253.

[6] 尚吉揚, 張 宇, 陳 馳, 等. 小波改進最大平均相關(guān)高度法實現(xiàn)畸變目標識別[J]. 儀器儀表學(xué)報, 2011, 32(9)： 2057-2065.

[7] 劉 全, 張 樂, 張永剛, 等. 基于角點特征的幾何同步數(shù)字水印算法[J]. 通信學(xué)報, 2011, 32(4)： 25-31.

[8] 支力佳, 張少敏, 趙大哲, 等. 融合多種特征點信息的最小生成樹醫(yī)學(xué)圖像配準[J]. 計算機研究與發(fā)展, 2011,48(3)： 501-507.

[9]Selesnick I W. The double-density dual-tree DWT [J]. IEEE Transactions onSignal Processing, 2004, 52(5)：1304-1314.

Digital Watermarking Algorithm for Volume Data Based on DD-DT CWT andSIFT

Cui Hanguo1, Liu Jianxin2, Li Bin1
(1. College of Power Engineering, Naval University of Engineering, Wuhan Hubei430033, China; 2. Department of Electromechanical Engineering, Bengbu Naval Petty Officer Academy, Bengbu Anhui 233012, China)

To ensure theSecurity and integrity of volume data during transmission in the network, a watermarking embedding algorithm for volume data based on double-density dual-tree complex wavelet transform was proposed, the watermarking embedding capacity and robustness to imageManipulations attacks were improved, and watermarking extraction algorithm for volume data based onScale invariant feature transform was proposed, the original volume dataModel and the watermarkedModel were registered accurately, and the occurrence ofSynchronization error was avoided, the robustness to geometric attacks was improved. Experimental resultsShow that the presented algorithm is efficient and robust toMany attacksSuch as compression, wave filtering, noise, rotation, translation, uniformScaling and cropping under blind detection.

digital watermarking; volume data; double-density dual-tree complex wavelet transform;Scale invariant feature transform;Singular value decomposition

A

2095-302X(2015)02-0148-04

2014-06-26；定稿日期：2014-08-15

湖北省自然科學(xué)基金資助項目(2012FB06904)

崔漢國(1964–)，男，江蘇無錫人，教授，博士。主要研究方向為計算機圖形學(xué)、虛擬現(xiàn)實。E-mail：cuihanguo@163.com

TP 391