【正文】 f x y??? ? ? ? ? ? ? (2 , , ) ( ( , ) , ( , ) ) ( , )y j y yW f x y f u v u x v y f x y? ? ? ? ? ? ? （ 54） 任意的 22()f L R? 的二進制小波變換定義為如下函數(shù)族： ? ?( 2 , , ) ( 2 , , ) , ( 2 , , )j x j y jjZW f x y W f x y W f x y ?? （ 55） 為了保證小波變換的完備性和穩(wěn)定性，還必須滿足下面的條件：存在兩個正常數(shù)A和 B，對 ? ?2( , ) (0 , 0 )xyw w R? ? ?，使得 22( 2 , 2 ) ( 2 , 2 )xyj j j jx y x yjA w w w w B???? ? ?? ? ?? （ 56） 上式中， x? 和 y? 分別表示 x? 和 y? 的傅里葉變換。此外，因物體大小有差別，它們的邊緣也就有不同的尺度。小波變換能夠?qū)D像分解成多種尺度成分，并對大小不同的尺度采用相應(yīng)的時域或空域取樣步長，因此能夠捕捉到任何微小的細節(jié)。由于求導(dǎo)和求卷積是可以結(jié)合的，所以可以先用高斯函數(shù)濾波，然后再求導(dǎo)，即 ? ?( , ) ( ( , ) ) ( , ) ( , ) ( , )H x y G x y I x y G x y I x y? ? ? ? ? ? （ 416） 顯然，上式可以看成是先用高斯濾波器進行平滑，再求梯度。39。39。200( ) ( ) ( ) ( )( ) ( ) ( ) ()()WWwwwwG x f x d x G x f x d xJ f S N R f Lo c f n f x d xn f x d x??????? ? ? ?????????? （ 412） （ 3）在理想情況下，用濾波器對噪聲響應(yīng)的兩個最大值之間的距離來近似濾波器對一個邊緣點的長度。 39。39。39。0 0 0 0( ) ( 0 ) ( 0 ) ( ) ( 0 )G G G GH x H H x o x H x? ? ? ? 從而有： ? ?? ?239。 39。39。三準(zhǔn)則的數(shù)學(xué)表達式如下： （ 1）好的檢測結(jié)果。即標(biāo)示出的邊緣位置要和圖像上真正的邊緣位置充分接近。由于上面的原因， Laplace 算子很少直接用于檢測邊緣，而主要應(yīng)用于已知邊緣的像素后，確定 該像素是在圖像的暗區(qū)或明區(qū)一邊。因為它只用一個模板，并且不必綜合 各模板的值。依次用邊緣樣板去檢測圖像，與被檢測區(qū)域最為相似的樣板給出最大值，然后用這個最大值作為算子的輸出值 (, )Pi j ， 這樣可以將邊緣檢測出來。 它 對噪聲具有平 滑作用，能夠得到較為準(zhǔn)確的邊緣信息。因此， Sobel 算子的定義如下： ( , ) xys i j f f? ? ? ? 19 = ( 1 , 1 ) 2 ( 1 , ) ( 1 , 1 ) ( 1 , ) 2 ( 1 , ) ( 1 , 1 )f i j f i j f i j f i j f i j f i j? ? ? ? ? ? ? ? ? ? ? ? ? ?+ ( 1 , 1 ) 2 ( , 1 ) ( 1 , 1 ) ( 1 , 1 ) 2 ( , 1 ) ( 1 , 1 )f i j f i j f i j f i j f i j f i j? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 它的卷積算子為： 1 0 1: 2 0 21 0 1x f?????????? 1 2 1: 0 0 01 0 1y f? ? ???????? 選擇適當(dāng)?shù)拈T限 TH，如果 (, )Ri j TH，則 (, )ij 為階躍狀邊緣點， ? ?(, )Ri j 即為邊緣圖像。 4． 2． 2 Roberts 邊緣檢測算子 Roberts 邊緣檢測算子根據(jù)任意一對相互垂直方向上的差分可用來計算梯度的原理，采用對角線方向相鄰兩像素之差，即： ( , ) ( 1 , 1 )x f f i j f i j? ? ? ? ? ( , 1 ) ( 1, )y f f i j f i j? ? ? ? ? （ 43） 有了 yf? ， xf? 后，可以根據(jù)下式計算出 Roberts 的梯度幅度值 (, )Ri j ： 22( , ) xyR i j f f? ? ? ? 或 ( , ) xyR i j f f? ? ? ? （ 44） 它們的卷積算子為： 10:01x f ??? ???? 01:10y f ??? ????? （ 45） 選擇合適的門限 TH，如果 (, )Ri j TH，則 (, )ij 為階躍狀邊緣點， ? ?(, )Ri j 即為邊緣圖像。但是，用這種方法就必須使差分的方向與邊緣方向垂直，這就需要對圖像的不同方向都進行差分運算，實 際上增加了運算的復(fù)雜性。為了簡便，也可以將 ( , )exy 定義為偏導(dǎo)數(shù) ,xyff的絕對值之和，如下： ( , )exy ? ( , ) ( , )xyf x y f x y? （ 42） 人們根據(jù)這些理論提出了許多算法，現(xiàn)在比較常用的邊緣檢測方法有差分檢測法、Roberts 邊緣檢測算子、 Sobel 算子、 Prewitt【 14】 算子、 Robinson 算子、 Laplace 算子、Canny 算子 【 15】 和 LOG 算子 【 16】 等等。我們將邊緣定義為圖像中灰度發(fā)生急劇變化的區(qū)域邊界，它的變化情況可以用圖像灰度分布的梯度來反映，因此我們可以用局部圖像微分的方法來獲得邊緣檢測算子 。所以， 應(yīng)將局部邊緣聚合到一副圖像中，使其中只出現(xiàn)對應(yīng)于存在的物體或圖像部分的邊緣鏈。本章首先討論了傳統(tǒng)的基于分割的邊緣檢測方法；然后對幾種經(jīng)典的邊緣算子進行理論分析，并對各自的特點做出了比較和評價；最后介紹 Canny 連續(xù)準(zhǔn)則及其算法。它可以進一步細化為多個研究方向：圖片處理、模式識別、景物分析、圖像理解、光學(xué)處理等等。 總的來說，小波分析的應(yīng)用是與小波分析的理論研究緊密聯(lián)系在一起的。 對得到的新小波系數(shù) W? 進行逆小波變換，得到去噪后的圖像。 軟閾值化與硬閾值化是對超過閾值的小波系數(shù)進行縮減的兩種主要方法。如果一個信號能量集中于變換域少數(shù)小波系數(shù)上，那么它們的取值必然大于在變換域內(nèi)能量分散的大量信號和噪聲的小波系數(shù)，這種情況下就可以用閾值方法。 b：基于相關(guān)性：首先對含噪信號做小波變換，然后計算相鄰尺度間小波系數(shù)的相關(guān)性，根據(jù)相關(guān)性的大小區(qū)別小波系數(shù)的類型，做出取舍，最后進行重構(gòu)。 c：基波選擇的靈活性。 小波去噪 【 12】 的成功主要是因為小波變換有如下的特點： a：低熵性。而對于一個圖像，表現(xiàn)它的最主要的 部分是低頻部分，所以最簡單的壓縮方法是利用小波分解，去掉圖像中的高頻部分而保留低頻。它具有壓縮速度快、壓縮比高、壓縮后能保持信號與圖像的相關(guān)特征不變等特點，且在傳輸中可以抗干擾。壓縮就是想法去掉各種冗余，保留真正有用的信息。小波變換由于具有良好的時域局部化性能，有效的克服了 Fourier 變換在處理非平穩(wěn) 的復(fù)雜圖像信號時的局限性，因而在圖像壓縮領(lǐng)域 受 到了廣泛的重視?，F(xiàn)在，對于性質(zhì)穩(wěn)定不變的信號，處理的理想工具仍然是 Fourier 變換。 3． 2． 3 小波級數(shù) 首先，從正交小波的定義開始，記： 2, ( ) 2 ( 2 ) , ,kkkn t t n k n Z??? ? ? （ 317） 一個小波 2( ) ( )t L R? ? 稱為一個正交小波，如果 ? ?,kn? 是 2()LR的一個規(guī)范正交基，即： , , , , , , ( ) ( ) , , , ,k n l m k n l m k l n mt t dt k m n l Z? ? ? ? ? ????? ? ?? （ 318） 而且，每個 2()f L R? 都能寫成： ,( ) ( )k n k nknf t d t?????? ? （ 319） 13 上面的定理中，對 2, ( )f g L R? ， , ( ) ( )f g f t g t dt???? ?表示空間 2()LR的內(nèi)積，而()gt 表示函數(shù)的復(fù)共軛。因此，選取 0maa? ， 00mb nba? ，其中 ,mn取遍 Z，而 001, 0ab??是固定的。 0ab tbt a a b R aa??? ?? ? ? （ 314） 這時，用函數(shù) f(t)的連續(xù)小波變換就可以重構(gòu)函數(shù)，我們知道，這種重構(gòu) 是 冗余的。以后再講到小波 的時候，若無特別聲明，指的就是基小波。為了研究的需要，上式也常記為 ? ?( , )W f a b? 。由于函數(shù)族 ? ?, ()abt? 是由 ()t? 生成的，所以有時將其稱為母小波。 下面給出一些小波的例子。 下面 簡述一下小波在直觀上的特點。因為這些特點，小波分析可以探測正常信號中的瞬態(tài)，并展示其頻率成分，被稱為“數(shù)學(xué)顯微鏡”，廣泛應(yīng)用于各個時域分析領(lǐng)域?？偠灾?， 傳統(tǒng)傅里葉分析將一段較長的時間內(nèi)的信號看作從以前到未來信號一直按照此規(guī)律周期性變化，并分析信號的頻率成分 ；而Gabor 變換在要分析的信號上提取出信號中的每一個小段 (長度自己定 )，將此小段進行兩端周期性延拓，并對這樣的信號進行傳統(tǒng)傅里葉分析，得到此小段內(nèi)信號的頻率特性，平移原有分析信號中小段的位置，得到整個要分析信號在每個小范圍內(nèi)頻率成分。 2． 2． 2 離散 Gabor 變換 在基于計算的實際應(yīng)用中，我們使用的是離散變換而不是連續(xù)變換。窗函數(shù) ()gt 中心 *t 與半徑 g? 分別定義為： 2*221: ( )t t g t dtg???? ? （ 23） ? ? 1 22*221: ( ) ( )g t t g t d tg ???? ? ?? (24) 其中，窗函數(shù)的寬度為半徑的 2倍。再利用快速卷積的方法計算它們的卷積，提取卷積運算的峰值，得到在圖像 中字母“ a”的定位 （白色小點） ，程序執(zhí)行的結(jié)果如下： 圖 21 圖像 圖 22 字母“ a”的圖像 7 圖 23 快速卷積結(jié)果 圖 24 對字母“ a”的定位結(jié)果 2． 2 Gabor變換 本節(jié)介紹的 Gabor 變換是信號時 頻分析的一種重要工具。 MATLAB 中提供了函數(shù) conv2 來實現(xiàn)二維矩陣的卷積。 6 （ 2） 快速卷積 能夠?qū)崿F(xiàn)快速卷積是 Fourier 變換 的另一個重要應(yīng)用。 理由如下： （ 1） DFT的輸入 /輸出 均為離散值，非常適合于計算機的運算操作。對大多數(shù)工程信號來 說，這個條件是容易滿足的。本節(jié)主要講述 Fourier 分析的一些基礎(chǔ)知識 和它在圖像處理中的一些應(yīng)用 。比如，F(xiàn)ourier 變換可以使處理工作在頻率中進行，簡化了運算； Gabor 變換用開窗的方法作為 Fourier 變換的一種簡單局部化，使計算更為方便。 ( , )f xy 表示在特定點 (, )xy 處的函數(shù)值，用來表征圖像在該點相應(yīng)的顏色強度或者灰度。本章用改進的 B樣條函數(shù)作為小波函數(shù)，對圖像 進行多級的邊緣檢測。 本章首先研究了基于傳統(tǒng)算子 Sobel、 Laplace 等的邊緣檢測方法，對它們進行 了 理論分析，然后對各自的特點做出了比較和評價。 第三章是小波變換理論。主要闡述了課題背景和國內(nèi)外的研究現(xiàn)狀。在研究和分析了現(xiàn)有的圖 像邊緣檢測方法后，針對可能漏檢微弱邊緣和邊緣定位不夠準(zhǔn)確的不足，采用小波變換對圖像進行邊緣檢測。生物識別技術(shù)以其特有的穩(wěn)定性、唯一性和方便性，得到越來越廣泛的應(yīng)用。于是，將多種算法綜合起來的多特征融合的方法便成為研究的方向。基于子空間的特征提取是指將掌紋圖像經(jīng)過映射變換 或 矩陣運算，實現(xiàn)從樣本空間 到特征子空間的轉(zhuǎn)換。 掌紋可以被認為是無規(guī)則但在個體間獨一無二的一種紋理。現(xiàn)階段，國內(nèi)外主要有以下幾種掌紋識別算法： （ 1）基于點特 征和線特征的識別方法。除此之外，掌紋識別還是一種非侵犯性的識別方法，用戶比較容易接受，同時，對設(shè)備的要求也不是太高。 （ 3） 紋理特征，是指比紋線更短、更細的一些紋線， 并且是毫無規(guī)律的分布在手掌上。 掌紋 體現(xiàn)在圖像上的特征主要