DWT域數字水印算法的FPGA實(shí)現

隨著(zhù)計算機網(wǎng)絡(luò )和通信技術(shù)的飛速發(fā)展，數字媒體（包括數字圖像、數字視頻、數字音頻）已得到了廣泛的應用，隨之而來(lái)的數字媒體的信息安全、知識產(chǎn)權保護和認證等問(wèn)題也變得日益突出。數字水印作為傳統加密方法的有效補充手段，是一種可以在開(kāi)放的網(wǎng)絡(luò )環(huán)境下保護版權和認證來(lái)源以及完整性的新技術(shù)[1]。

根據水印的實(shí)現過(guò)程，圖像水印算法可分為空域算法和變換域算法?？沼蛩惴ㄊ峭ㄟ^(guò)直接改變原始圖像的像素值來(lái)嵌入水印，通常具有較快的速度，但魯棒性差，且水印容量也會(huì )受到限制；變換域算法是通過(guò)改變某些變換系數來(lái)嵌入水印，通常具有很好的魯棒性和不可見(jiàn)性,其實(shí)現一般基于圖像變換(如DCT、DFT、DWT等)。大多數水印算法采用軟件實(shí)現,軟件實(shí)現具有易于應用、升級和適應性較好等特點(diǎn)，但存在速度受限、難以滿(mǎn)足實(shí)時(shí)處理的問(wèn)題。而采用硬件實(shí)現則可以克服這些問(wèn)題。本文基于上述考慮，利用FPGA設計了基于DWT(Discrete Wavelet Transform)域的數字水印算法，其中提升小波變換是該數字水印方案的關(guān)鍵單元之一，采用硬件描述語(yǔ)言Verilog HDL對算法進(jìn)行描述，并用ModelSim完成了功能和時(shí)序仿真。

1 離散小波變換原理

1.1 Mallat算法及提升格式

1987 年,Mallat將計算機視覺(jué)領(lǐng)域的多分辨思想引入了小波分析，提出了多分辨分析MRA(Multi-resolution Analysis)的理論，并給出了數學(xué)描述和一種子帶濾波器的離散小波變換與重構算法——Mallat算法。在圖像處理領(lǐng)域,需要處理的通常是二維圖形。因此,將Mallat算法擴展到二維空間, 適當地選取一組行和列變換正交的小波系數,對圖像(或分解后的低頻子圖)分別進(jìn)行行變換和列變換。然后,根據后續的具體應用對N次分解所得的圖像在不同的分辨率下進(jìn)行分析、處理或數據壓縮。二維Mallat算法的結構如圖1所示。

首先,用濾波器組g(n)和h(n)對圖像的每一行信號進(jìn)行濾波和2抽取,然后用相同的濾波器組對所得結果的每一列再次進(jìn)行濾波和2抽取。這樣,原始圖像就被分解成為4個(gè)子帶圖像,分別表示為L(cháng)L、LH、HL和HH。其中,LL為水平和垂直方向的低通信號;LH為水平方向低通,垂直方向高通;HL為水平方向高通,垂直方向低通;HH為水平和垂直方向的高通信號。

相對于Mallat算法而言，提升算法[2]是一種更快更有效的小波變換方法，標準的提升算法包含分裂(Split)、預測(Prediction)、更新(Update)三個(gè)步驟。分裂是將輸入數據分為偶數序列和奇數序列兩個(gè)部分;預測是用分裂的偶數序列預測奇數序列,得到的預測誤差為變換的高頻分量;更新是由預測誤差來(lái)更新偶數序列,得到變換的低頻分量。提升算法的結構框圖如圖2所示，其提升實(shí)現過(guò)程由式(1)和式(2)表示，其中濾波系數α=-1/2, β=1/4。


本文將采用5/3小波濾波器來(lái)實(shí)現小波變換。

1.2 邊界處理

由于圖像數據是有限長(cháng)的，因此離散小波變換必須對圖像數據進(jìn)行邊界延拓，在做小波提升算法時(shí)，同樣需要對其邊界數據進(jìn)行延拓，以保證邊緣數據的正確?；谫Y源消耗和功耗的考慮，這里采用內嵌對稱(chēng)延拓提升算法[3-4]，將延拓的過(guò)程融合到提升變換過(guò)程中。精簡(jiǎn)后的內嵌延拓方式如圖3所示，其中Xn表示圖像數據，dn表示高頻分量，cn表示低頻分量。

2 數字水印方案的設計

基于小波的圖像水印算法有很多，本文采用LSB方法[5]直接將水印信號嵌入到經(jīng)過(guò)小波變換的子圖中，這里所用到的水印信號是由線(xiàn)性移位寄存器產(chǎn)生的m序列。圖像經(jīng)過(guò)小波變換后，低頻子帶 LL 表示為 L(n)，數字水印序列設為a(n)，在低頻子帶中嵌入水印序列，算法如下：

其中α為常數，當水印序列為“1”時(shí)，對子帶數據加α；當水印序列為“0”時(shí)，對子帶數據減α。

整個(gè)設計過(guò)程的核心為提升小波變換的FPGA的實(shí)現。根據二維小波變換的特點(diǎn)，將總體設計分為水平一維DWT 單元、垂直一維DWT單元、外部存儲單元和控制單元?？傮w設計框圖如圖4所示。

2.1 二維DWT單元

二維提升式小波變換可以分解為兩個(gè)一維小波變換，先進(jìn)行行方向的離散小波變換，再進(jìn)行列方向離散小波變換，這樣就可以完成二維提升式小波變換。在本設計中采用的是并行的二維提升離散小波變換，根據提升式小波變換的特點(diǎn)，當列變換處理完三行數據后，即可以進(jìn)行行方向的變換。因此可以根據二維小波變換[6]這個(gè)特點(diǎn)，使用四行緩存來(lái)緩存列方向小波變換后的系數。在第三行緩存存滿(mǎn)后就開(kāi)始從左向右進(jìn)行“行方向”小波變換。因為在進(jìn)行行方向小波變換的同時(shí)，仍然在進(jìn)行列方向的小波變換，所以需要第四行緩存繼續緩存列方向變換得到的小波系數以避免數據的丟失。通過(guò)以上的處理，實(shí)現了并行的二維離散小波變換，極大地提高了小波變換的速度。在這整個(gè)架構包含了列方向小波變換模塊及其控制模塊、行方向小波變換模塊及其控制模塊，以及用于緩存一維小波系數的緩存組。其中列方向變換模塊在列方向變換控制模塊的控制下，從系統外部讀入影像數據進(jìn)行列方向小波變換處理，并根據列變換的控制信號將處理后得到的小波高、低頻系數存儲到系統內部緩存中。行方向變換模塊在行方向變換控制模塊的控制下，從緩存器中取出列方向小波變換處理后的小波系數，進(jìn)行行方向小波變換處理，得到二維小波變換后的LL、LH、HL、HH小波系數，最后輸出到系統外。

2.2 存儲單元

在變換過(guò)程中，首先對圖像數據進(jìn)行行方向的一維提升式小波變換，將變換后的數據存放在存儲器或臨時(shí)寄存器中，行方向變換結束后，再對該存儲器中的數據進(jìn)行列方向變換。這里使用的是FIFO存儲器。FIFO是一種先進(jìn)先出的存儲器，就好像一個(gè)單向管道，數據只能按照固定的方向從管道一頭進(jìn)來(lái)，再按照相同的方向從管道另一頭出去，最先進(jìn)來(lái)的數據必定最先出去。

2.3 控制單元

控制單元接收外部控制信號(芯片使能、數據時(shí)鐘等),由內部計數器進(jìn)行計時(shí)、產(chǎn)生水平一維DWT 單元、垂直一維 DWT 單元和存儲單元的控制信號(模塊使能、起始終止信號等),對外部數據或前一級的LL輸出數據進(jìn)行鎖存、處理、緩沖、輸出等操作。

3 實(shí)驗結果

在整個(gè)電路設計過(guò)程中,采用Verilog HDL語(yǔ)言對二維小波變換結構進(jìn)行了RTL設計，并且采用ModelSim仿真對設計進(jìn)行驗證，其結果與 Matlab仿真結果一致，說(shuō)明二維提升小波變換結構的正確。圖5為α＝2時(shí)嵌入水印后的圖片和算法仿真結果。由圖可知，嵌入水印后對圖的影響不大，基本不能被人的眼睛所感知。但當α>8時(shí)，圖像質(zhì)量還是有了一定的影響,這是因為L(cháng)L子帶包含了圖像的主要信息。

本文實(shí)現了一種基于DWT域的數字水印算法，其中DWT變換算法的實(shí)現是該設計中數字水印方案的關(guān)鍵算法之一。采用ModelSim仿真對設計進(jìn)行驗證，其結果與Matlab仿真結果一致，說(shuō)明了該設計的正確性。實(shí)驗結果表明,該方案利用FPGA實(shí)現水印算法比用軟件實(shí)現法速度更快，性能更穩定。