【導(dǎo)讀】師的指導(dǎo)下進行的研究工作及取得的成果。盡我所知，除文中特別加。而使用過的材料。對本研究提供過幫助和做出過貢獻的個人或集體，均已在文中作了明確的說明并表示了謝意。環(huán)境下，保護版權(quán)和認證來源及完整性的新型技術(shù)。該算法選擇了檢測結(jié)果直觀、有特殊意義的二值圖像作為原始水印，見性和魯棒性，利用多分辨率分析思想進行水印的嵌入與提取?？s、噪聲、濾波、剪切等，均有較好的魯棒性。

　　

【正文】 止信息的丟失，我們要求采樣間隔 ? 滿足 Nyquist 采樣定理，采樣率大于等于該尺度下頻率通常的二倍。所以每當(dāng) m增加 1 時，尺度 a增加一倍，對應(yīng)的頻率減小一半，可見采樣率可以降低 一半而不致引起信息的丟失 (帶通信號的采樣率決定于其帶寬，而不是決定于其頻率上限 )。所以在尺度 j下，由于 ????????? tf0?的帶寬時 ??t? 的 aj0 倍，因此采樣間隔可以擴大 aj0 ，同時也不會引起信息的丟失。這樣， )(, ta?? 就改成 :a ? ? ? ?002020200 )( ???? ktaakata jjjjj ??? ???? () 記為 )(00, tkaj ??離散小波變換定義為 : WT ),( 00 ?kajf = )()()(00 , tdttf ka j ??? j=0,1,2...,k Z? () 在以上的尺度以及位移均離散化的小波序列，如果取離散柵格 a0 = 2 ,0? =0， 即相當(dāng)于連續(xù)小波只 在尺度 a 上進行量化，平移參數(shù) ? 仍然連續(xù)不被離散，我們稱這類小波為二進小波，表示為 : )(,2 tK??=2 )2(2kk t ?? ?? () 22 二進小波介于連續(xù)小波和離散小波之間，由于它只是對尺度參量進行離散化，在時間域上的平移量仍保持著連續(xù)的變化，所以二進小波具有連續(xù)小波變換的時移共變性，這個特點也是離散小波不具有的。也正因為如 此，它在奇異性檢測、圖像處理方面都十分有用。 令小波函數(shù)為 ? (t)，其傅立葉變換為 )(?? ，若存在常數(shù) A,B， 當(dāng) 0A? B? ，使得 BA zk kw ?? ?? 2)2(? () 此時， ? (t)才是一個二進小波，我們稱上式為二進小波的魯棒性條件。 定義函數(shù) f(t) )(2 RL? 的二進小波變換系數(shù)為 : WTK2(? )=f(t) )(,2 tk??=tkk dttf )2()(2 2 ? ??? () 其中 )(,2 tk??=22k? )2(kt ?? ? () 由前面的知識可得它的小波逆變換公式是存在的。 二進小波變換的重建公式為 : f(t)=???? dtWT KKzk R )()( ,22 ?? ?? () 其中， )(,2 tk??? 為??,2k(t)的對偶框架，其上、下界分別為 B1? ，和 A1? 多分辨率分析 多分辨率概念是由 Mallat 和 Meyer 于 1986 年提出來的，它可將此前所有正交 小波基的構(gòu)造統(tǒng)一起來，使小波理論產(chǎn)生突破性的進展。同時，在多分辨率理論分析的基礎(chǔ)上， Mallat 引入了一種計算柵格上小波變換的快速算法，即 Mallat算 法。它可以避免 a 值越大，對 ? (t)的采樣就得越密的缺點，這一算法在小波分析 中的地位很重要，相當(dāng)于快速傅立葉算法 (FFT)在經(jīng) 典傅立葉分析中的地位。 我們把平方可積的函數(shù) f(t) )(2 RL? 看成使某一逐級逼近的極限情況，每級逼 近都是用某一低通平滑函數(shù) ? (t)對 f(t)做平滑的結(jié)果，在逐級逼近時平滑函 23 數(shù) ) )(t? 也做逐級伸縮，這就是“多分辨率” ???，即用不同分辨率來逐級逼近待分析函數(shù) f(t)。 空間 L2 (R)中的多分辨率分析是指 L2 (R)中滿足下列條件的一個空間序列 ??jezjV ： ① 單調(diào) 性 :對任意 j? Z，有 V 1??jj V 。 ② 逼近性 : ??0??jez jV, )(2 RLVj j ??????。 ③ 伸縮性 :f(t) 1)2( ???? jj VtfV ，伸縮性體現(xiàn)了尺度的變換、逼近正交 小波函數(shù)的變化和空間的變化具有一致性。 ④ 平移不變性 :對任意 k?Z，有 jjj Vkt ??? )2( 2? 。 ⑤ Riesz基存在性 :存在 0)( Vt ?? ,使得jj Vkt ??????? ?? )2( 2? 構(gòu)成 jV 的 Riesz。 為了構(gòu)造正交小波基，引用尺度函數(shù)概念 : 定義函數(shù) )()( 2 RLt ?? 為尺度函數(shù)，若其經(jīng)過整數(shù)平移 k和尺度 j上的伸縮，得到一個尺度和位移均可以變化的函數(shù)集合 : )2(2 2 ktjjjk ?? ?? ?? () 稱每一個尺度 j上的平移系列 )(tjk? 所組成的空間 Vj 為尺度 j的尺度空間 ： Vj =span? ?)(tjk? k z? () 對于任意函數(shù) f(t) jV? ， 有 )2(2 2 ktjjjk ?? ?? ?? () 所以尺度函數(shù) )(t? 在不同尺度下其平移系列組成了一系列的尺度空間??jezjV 隨著尺度的 j增大，函數(shù) )(tjk? 的定義域變大，平移的間隔 ( 02?j )也 變大，所以它的線性組合不能表示函數(shù)小于該尺度的細微變化，所以其張成的尺度空間 24 只能包括大尺度的緩變信號。反之，如果 j減小， 函數(shù) )(tjk? 的定義域變小，平移的間隔 ( 02?j )也變小，則它的線性組合就能表示函數(shù)的更細微的變化，所以其張成 的尺度空間所包含的函數(shù)增多，包括小尺度信號和大尺度的緩變信號。 實驗環(huán)境 ： 可實現(xiàn)數(shù)字水印技術(shù)的高效實用工具 —— Matlab Matlab 簡介 Matlab 是當(dāng)前在國內(nèi)外十分流行的工程設(shè)計和系統(tǒng)仿真軟件包。它是MathWorks 公司于 1982 年推出的一套高性能的數(shù)值計算和可視化軟件，它集數(shù)值分析、矩陣運算、信 號處理和圖形顯示于一體，構(gòu)成了一人方便的、界面友好的用戶環(huán)境。 Matlab 的推出得到了各個領(lǐng)域?qū)＜?、學(xué)者的廣泛關(guān)注，其強大的擴展功能為各個領(lǐng)域的應(yīng)用提供了基礎(chǔ)。由各個專家學(xué)者相繼推出了 MATLAB 工具箱，其中的信號處理 (signal processing)、控制系統(tǒng) (control system)、神經(jīng)網(wǎng)絡(luò)(neural work)、圖像處理 (image processing)、魯棒控制 (robust control)、非線性系統(tǒng)控制設(shè)計 (nonlinear system control design)、系統(tǒng)辨識 (system identification)、最優(yōu)化 (optimization)、 模糊邏輯 (fuzzy logic)、小波(wavelet)、通信 (munication)、統(tǒng)計 (statistics)等工具箱，這些工具箱給各個領(lǐng)域的研究和工程應(yīng)用提供了有力的工具，借助于這些 “ 巨人肩上的工具 ” ，各個層次的研究人員可直觀、方便地進行分析、計算及設(shè)計工作，從而大大地節(jié)省了時間。 用 Matlab 研究數(shù)字水印的優(yōu)點 集成了 DCT、 DWT 等函數(shù)有豐富的小波函數(shù)和處理函數(shù)，這不僅方便了研究人員，而且使源程序簡潔明了、易實現(xiàn)。 強大的數(shù)學(xué)運算功能。能夠方便、高效地實現(xiàn)音頻、視頻中的大量矩陣運算。 提供了圖像處理工具箱、小波分析工具箱、數(shù)字信號處理工具箱。用來編制跨數(shù)字圖像處理技術(shù)、數(shù)字信號處理等多學(xué)科的數(shù)字水印技術(shù)是非常好的選擇。 MATLAB 與目前最強大的編程工具 —— Visual C++具有良好的接口。 25 第 4章 基于小波變換的數(shù)字水印算法 Amold變換是 Amold在遍歷理論研究中提出的一種變換，俗稱貓臉變換原意為cat mapping。設(shè) 想在平面單位正方形內(nèi)繪制一個貓臉圖像，通過如下變換 ????????39。39。yx = ??????1211 ??????yxmod1 () 這個貓臉圖像將由清晰變?yōu)槟：@就是 Arnold變換。注意到式 ()定義的Amold變換實際上是一種點的位置移動，并且這種變換是一一對應(yīng)的。此外，這個變換可以迭代地做下去。類似的變換還有面包師變換。 對于數(shù)字圖像來說，可 以將其看成是一個函數(shù)在離散網(wǎng)格點處的采樣值，這樣我們就得到了一個表示圖像的矩陣，矩陣中元素的值是對應(yīng)點處的灰度值或RGB顏色分量值。對于正方形數(shù)字圖像，我們有離散化的 Amold變換 : ????????39。39。yx = ??????1211 ??????yxmodN,x,y ? ?1...1,0 ?? N () 其中 N為圖像的高度和寬度。 對于數(shù)字化圖像而言，我們所說的位置移 動實際上是對應(yīng)點的灰度值或者 RGB顏色值的移動，即將原來點 (x,y)處象素對應(yīng)的灰度值或 RGB顏色值移動至變換后的點 (x39。 ,y39。 )處。如果我們對一個數(shù)字圖像迭代地使用離散化的 Amold變換，即將左端輸出的 (x39。 ,y39。 )T 作為下一次 Amold變換的輸入可以重復(fù)這個過程一直作下去當(dāng)?shù)侥骋?步時，如果出現(xiàn)的圖像符合我們對圖像的“雜亂無章”標準的要求，這即是一幅置亂了的圖像。 需要注意的是， Amold變換具有周期性，即當(dāng)?shù)侥骋徊綍r，將重新得到 原始圖像。 Dyson和 Falk分析了離散 Amold變換的周期性，給出了對于任意 N2, Amold 變換的周期性幾 T 2NN? /2，這也許是迄今為止最好的結(jié)果了。 依據(jù) Cox的觀點，水印信息應(yīng)該具有不可預(yù)測的隨機性，具有與噪聲相同的特性，這不僅可以提高水印的透明性，而且可以加強水印抗千擾的能力 ，提高水 26 印的魯棒性。但是無意義的偽隨機序列通常應(yīng)用價值不大，有意義水印可附帶許多證明信息，如原作者的個人標志，產(chǎn)品序列號等，在版權(quán)證明上顯然較無意義水印更具有直觀性和可驗證性。所以在本算法中，選取有特殊意義的二值圖像做為原始水印。 另外 ,為了提高水印的不可見性和魯棒性，在水印嵌入之前，我們先對原始水印進行 Amold變換預(yù)處理，處理之后的水印圖像各像素點變?yōu)殡s亂均勻分布，這樣不僅提高了水印的透明性，而且加強水印抗干擾的能力，提高了水印的魯棒性。當(dāng)嵌有水印的圖像在網(wǎng)絡(luò)中傳播時，難免會遇到有意或無意的干擾破壞， 印也會受到相應(yīng)損壞，當(dāng)經(jīng)過預(yù)處理的水印被提取出來并利用 Amold變換恢復(fù)出原始水印圖像時，被損壞的部分就被分散到了全圖，對人類視覺的影響也就不明顯了，這樣就相當(dāng)于加強了水印的魯棒性。同時，由于提取出的水印是一幅被置亂的圖像，只有原始嵌入者知道采用的何種變換及變換次數(shù)，從而利用 Arnold變換的周期性將之恢復(fù)成原始水印圖像，而非法攻擊者面對雜亂圖像，不僅不能夠得到有用信息，并且很難進行偽造，所以，預(yù)處理也增強了算法的安全性。 目前的小波域水印算法 ,對于水印嵌入位置的選擇有不同的意見。一種意見認為低頻子圖是圖像 的平滑部分，人眼對這部分的失真比較敏感，基于水印的不可感知性考慮，應(yīng)將水印數(shù)據(jù)隱藏在圖像的高頻部分亦即小波分解后的高頻系數(shù)中，而不應(yīng)在低頻系數(shù)嵌入水印。另一種意見則認為中高頻子圖的小波系數(shù)幅度 一般較小，常接近于 0，而低頻部分集中了圖像的大部分能量，系數(shù)的振幅比細節(jié)子圖的系數(shù)大得多，由人類視覺特性知，背景亮度越大，嵌入信號的 JND就越高，即低頻逼近子圖具有較大的感覺容量，相當(dāng)于一個強背景，可以容納更強或者更多的水印信息，只要迭加的水印信號低于 JND值，視覺系統(tǒng)就無法感覺到信號的存在。這樣在圖像有一定失真的情 況下，仍能保留主要成分，可保持原始載體圖像的主觀視覺質(zhì)量基本不變，于是提出水印嵌入低頻系數(shù)中。 (根據(jù)小波變換的特性和小波分解系數(shù)分析可知，各級小波子圖對視覺系統(tǒng)的影響是不同的，隨著分級的增加，其重要性也隨之增加，在同一尺度下，水平、垂直子圖的重要性稍高于對角子圖，人眼對水平、垂直分量上的變化比對角線分量上的變化要敏感。 ) 以前的很多算法不在低頻系數(shù)中加入水印，原因是避免出現(xiàn)方塊效應(yīng)，但經(jīng) 27 過實驗證明，不在低頻部分嵌入所有水印，只嵌入一部分水印，再在中頻部分嵌 入一部分水印，既能保證不可見性又有很好的魯棒性。 在小波域 ,為了使數(shù)字水印具有較好的魯棒性，用于嵌入水印的小波系數(shù)就應(yīng)該滿足以下兩個條件 :小波系數(shù)不應(yīng)該過多的被信號處理和噪聲干擾所改變 ； 具有較大的感覺容量，以便嵌入一定強度的水印后不會引起原始圖像視覺質(zhì)量的明顯改變 ???。 綜合考慮上述嵌入位置的探討以及小波分解系數(shù)的特點，本文將水印的嵌入 位置選擇為原始圖像經(jīng)過小波三級分解后的中頻和低頻分