小波變換和數(shù)字圖像處理中應(yīng)用

1、關(guān)于小波變換和數(shù)字圖像處理中的應(yīng)用第一張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的快速實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第二張，PPT共一百零八頁，創(chuàng)作于2022年6月 Fourier變換一直是信號處理領(lǐng)域中應(yīng)用最廣泛、效果最好的一種分析手段，是時域到頻域互相轉(zhuǎn)化的工具，從物理意義上講，傅里葉變換的實質(zhì)是把對原函數(shù)的研究轉(zhuǎn)化為對其傅里葉變換的研究。但是

2、傅里葉變換只能提供信號在整個時間域上的頻率，不能提供信號在某個局部時間段上的頻率信息。8.1 從傅里葉變換到小波變換的時頻分析法8.1.1 傅里葉變換第三張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.1 傅里葉變換 傅里葉變換：對于時域的常量函數(shù)，在頻域?qū)⒈憩F(xiàn)為沖擊函數(shù)，表明具有很好的頻域局部化性質(zhì)。 傅里葉變換反傅里葉變換第四張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.1 傅里葉變換x=sin(2*pi*50*t)+sin(2*pi*300*t);%產(chǎn)生50HZ和300HZ的信號f=x+3.5*randn(1,length(t);%在信號中加入白噪聲時間第五張，PPT共一百零八

3、頁，創(chuàng)作于2022年6月 由于傅立葉變換無法作局部分析，為此，人們提出了短時傅里葉變換（STFT）的概念，即窗口傅里葉變換。 短時傅里葉變換將整個時間域分割成一些小的等時間間隔，然后在每個時間段上用傅里葉分析，它在一定程度上包含了時間頻率信息，但由于時間間隔不能調(diào)整，因而難以檢測持續(xù)時間很短、頻率很高的脈沖信號的發(fā)生時刻。8.1.2 短時傅里葉變換第六張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.2 短時傅里葉變換基本思想是：把信號劃分成許多小的時間間隔，用傅立葉變換分析每一個時間間隔，以便確定該時間間隔存在的頻率。 STFT的處理方法是對信號施加一個滑動窗(反映滑動窗的位置)后，再作傅

4、立葉變換。即： 時限頻限第七張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.2 短時傅里葉變換第八張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.2 短時傅里葉變換 短時傅里葉變換的分析特點(a)頻率變化的影響 (b) 基本分析單元的特點第九張，PPT共一百零八頁，創(chuàng)作于2022年6月小波起源： 1984年Morlet提出;1985年Meyer構(gòu)造出小波;1988年，Daubechies證明了離散小波的存在;1989年，Mallat提出多分辨分析和二進(jìn)小波變換的快速算法;1989年Coifman、 Meyer引入小波包;1990年崔錦泰等構(gòu)造出樣條單正交小波基;1994年Swelden

5、s提出二代小波提升格式小波(Lifting Scheme)。小波定義：“小”是指在時域具有緊支集或近似緊支集，“波”是指具有正負(fù)交替的波動性，直流分量為0。小波概念：是定義在有限間隔而且其平均值為零的一種函數(shù)。8.1.3 小波變換第十張，PPT共一百零八頁，創(chuàng)作于2022年6月持續(xù)寬度相同振蕩波波與小波的差異：第十一張，PPT共一百零八頁，創(chuàng)作于2022年6月用鏡頭觀察目標(biāo) (待分析信號)。 代表鏡頭所起的作用(如濾波或卷積)。 相當(dāng)于使鏡頭相對于目標(biāo)平行移動。 的作用相當(dāng)于鏡頭向目標(biāo)推進(jìn)或遠(yuǎn)離。 小波變換的粗略解釋 8.1.4 小波變換的時頻分析第十二張，PPT共一百零八頁，創(chuàng)作于2022年

6、6月尺度a較大距離遠(yuǎn)視野寬概貌觀察尺度a較小距離近視野窄細(xì)節(jié)觀察分析頻率低分析頻率高由粗到精多分辨分析品質(zhì)因數(shù)保持不變第十三張，PPT共一百零八頁，創(chuàng)作于2022年6月小波變換的時頻分析特點： 小波變換的分析特點(a) 尺度a不同時時域的變化 (b)尺度a不同時頻域的變化第十四張，PPT共一百零八頁，創(chuàng)作于2022年6月小波變換的多分辨分析特性：不同a值下小波分析區(qū)間的變化不同a值下分析小波頻率范圍的變化第十五張，PPT共一百零八頁，創(chuàng)作于2022年6月頻窗時窗小波變換的時頻局部特性： 第十六張，PPT共一百零八頁，創(chuàng)作于2022年6月8.1.5 連續(xù)小波變換尺度因子 的作用是將基本小波 做伸

7、縮， 越大 越寬。 小波的位移與伸縮 第十七張，PPT共一百零八頁，創(chuàng)作于2022年6月設(shè) ，當(dāng) 滿足允許條件時：8.1.5 連續(xù)小波變換稱 為一個“基小波”或“母小波”。小波變換的含義是：把基本小波(母小波)的函數(shù) 作位移后，再在不同尺度下與待分析信號作內(nèi)積，就可以得到一個小波序列。第十八張，PPT共一百零八頁，創(chuàng)作于2022年6月連續(xù)情況時，小波序列為： (基本小波的位移與尺度伸縮)其中 為尺度參量， 為平移參量。離散的情況，小波序列為 ：第十九張，PPT共一百零八頁，創(chuàng)作于2022年6月根據(jù)容許條件要求，當(dāng)=0時，為使被積函數(shù)是有效值，必須有 ，所以可得到上式的等價條件為：此式表明 中不

8、含直流，只含有交流，即具有震蕩性，故稱為“波”，為了使 具有局部性，即在有限的區(qū)間之外很快衰減為零，還必須加上一個衰減條件： 第二十張，PPT共一百零八頁，創(chuàng)作于2022年6月衰減條件要求小波具有局部性，這種局部性稱為“小”，所以稱為小波。對于任意的函數(shù) 的連續(xù)小波變換定義為：逆變換為： 是尺度因子， 反映位移。 第二十一張，PPT共一百零八頁，創(chuàng)作于2022年6月線性 設(shè)： 平移不變性 若 ，則伸縮共變性 如果 的CWT是 則 的CWT是冗余性(自相似性) 由連續(xù)小波變換恢復(fù)原信號的重構(gòu)公式不是唯一的8.1.6 連續(xù)小波的性質(zhì)第二十二張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1

9、從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的快速實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第二十三張，PPT共一百零八頁，創(chuàng)作于2022年6月8.2 小波變換分類 小波函數(shù)中 三個變量均為連續(xù)變量，稱為連續(xù)小波?？梢詫?三個變量施加不同的離散化條件，并相應(yīng)地對小波及小波變換進(jìn)行分類。 其中，最重要的兩種分類： 離散小波及離散小波變換 二進(jìn)小波及二進(jìn)小波變換第二十四張，PPT共一百零八頁，創(chuàng)作于2022年6月8.2.1

10、 離散小波變換如果設(shè)定 ，則 對于任意函數(shù) ，定義相應(yīng)的離散小波變換為：如果這時 構(gòu)成空間 的一組規(guī)范正交基，對于任一函數(shù) 的反演式為一展開式： 第二十五張，PPT共一百零八頁，創(chuàng)作于2022年6月8.2.2 二進(jìn)小波及二進(jìn)小波變換在連續(xù)小波變換中，令參數(shù) 而參數(shù) 仍取連續(xù)值，則有二進(jìn)小波：這時， 的二進(jìn)小波變換定義為第二十六張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的快速實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波

11、變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第二十七張，PPT共一百零八頁，創(chuàng)作于2022年6月 多分辨分析是小波分析中最重要的概念之一，它將一個函數(shù)表示為一個低頻成分與不同分辨率下的高頻成分，并且多分辨分析能提供一種構(gòu)造小波的統(tǒng)一框架，提供函數(shù)分解與重構(gòu)的快速算法。由理想濾波器引入多分辨率分析的概念：8.3 小波變換的多分辨分析特性第二十八張，PPT共一百零八頁，創(chuàng)作于2022年6月多分辨分析定義： 空間 中的一系列閉子空間 ，稱為 的多分辨率分析或逼近，若下列條件滿足：單調(diào)性： ，對任意 逼近性： 伸縮性： 平移不變性： Riesz基：存在 ，使 構(gòu)成 的Riesz基，即 是線

12、性無關(guān)的，且存在常數(shù) 與 ，滿足 使得對任意的 ，總存在序列 使得 且 ，稱 為尺度函數(shù)，并稱 生成 的一個多分辨分析 。第二十九張，PPT共一百零八頁，創(chuàng)作于2022年6月 是一個無限維向量空間，稱為平方可積空間，將 用它的子空間 ， 表示，其中 稱為尺度空間， 稱為小波空間。 尺度空間的遞歸嵌套關(guān)系： 小波空間 是 和 之間的差，即 ，它捕捉由 逼近 時丟失的信息。推出： 多分辨率的空間關(guān)系圖第三十張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的實現(xiàn)8.6

13、 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第三十一張，PPT共一百零八頁，創(chuàng)作于2022年6月兩尺度方程 若 是尺度函數(shù)，它生成 的多分辨分析 ，則必然存在系數(shù)序列 ，使得以下尺度關(guān)系成立:這就是兩尺度方程，必須滿足下列條件：定義函數(shù) 為尺度函數(shù)，若其經(jīng)過整數(shù)平移 和尺度 上的伸縮，得到一個尺度和位移均可變化的函數(shù)集合： 第三十二張，PPT共一百零八頁，創(chuàng)作于2022年6月 和 的基本性質(zhì)是兩尺度差分方程：兩尺度方程的頻域表示為： 第三十三張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變

14、換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的快速實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第三十四張，PPT共一百零八頁，創(chuàng)作于2022年6月8.5.1 Mallat算法與塔式分解 系數(shù)分解的快速算法： Mallat小波快速分解算法的流程圖第三十五張，PPT共一百零八頁，創(chuàng)作于2022年6月 系數(shù)重構(gòu)的快速算法： Mallat小波快速重構(gòu)算法的流程圖第三十六張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換

15、到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的快速實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第三十七張，PPT共一百零八頁，創(chuàng)作于2022年6月8.6.1 二維小波變換的實現(xiàn)假定二維尺度函數(shù)可分離，則有 其中 、 是兩個一維尺度函數(shù)。若 是相應(yīng)的小波，那么下列三個二維基本小波： 與 一起就建立了二維小波變換的基礎(chǔ)。第三十八張，PPT共一百零八頁，創(chuàng)作于2022年6月8.6.2 圖像小波變換的正變換正變換 圖像小波分解的正變換可以依

16、據(jù)二維小波變換按如下方式擴展，在變換的每一層次，圖像都被分解為4個四分之一大小的圖像。第三十九張，PPT共一百零八頁，創(chuàng)作于2022年6月8.6.3 圖像小波變換的逆變換逆變換在每一層(如最后一層)都通過在每一列的左邊插入一列零來增頻采樣前一層的4個陣列(即4個分解圖像)；接著用重構(gòu)低通濾波器h和重構(gòu)高通濾波器g來卷積各行，再成對地把這幾個的陣列加起來；然后通過在每行上面再插入一行零來將剛才所得兩個陣列(圖像)的大小增頻采樣為NN；再用h和g與這兩個陣列的每列進(jìn)行卷積。這兩個陣列的和就是這一層次重建的結(jié)果。 第四十張，PPT共一百零八頁，創(chuàng)作于2022年6月 對于二維圖像信號，在每一層分解中，

17、由原始圖像信號與一個小波基函數(shù)的內(nèi)積后再經(jīng)過在x和y方向的二倍間隔抽樣而生成四個分解圖像信號。對于第一個層次(j=1)可寫成：8.6.4 二維小波變換的Mallat算法第四十一張，PPT共一百零八頁，創(chuàng)作于2022年6月 將上式內(nèi)積改寫成卷積形式，則得到離散小波變換的Mallat算法的通用公式： 二維小波變換Mallat算法的通用公式：第四十二張，PPT共一百零八頁，創(chuàng)作于2022年6月8.6.5 二維Mallat多分辨率分解與重構(gòu)第四十三張，PPT共一百零八頁，創(chuàng)作于2022年6月圖像的Mallat快速塔式分解實驗第四十四張，PPT共一百零八頁，創(chuàng)作于2022年6月8.6.6 多孔算法第四十

18、五張，PPT共一百零八頁，創(chuàng)作于2022年6月多孔算法的分解實驗第四十六張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第四十七張，PPT共一百零八頁，創(chuàng)作于2022年6月 邊緣像素實質(zhì)上是局部圖像范圍內(nèi)灰度的急劇變化點(奇異點)，圖像邊緣就是二維圖像中奇異點的集合。邊緣點在頻域表現(xiàn)為高頻信號，而圖像噪聲也多為高頻信號，這

19、使得兩者難以區(qū)分。邊緣檢測的目的就是既要將高頻信號從圖像中分離出來，又要區(qū)分邊緣與噪聲，準(zhǔn)確地標(biāo)定邊緣的位置。 8.7 小波變換在圖像邊緣檢測中的應(yīng)用第四十八張，PPT共一百零八頁，創(chuàng)作于2022年6月傳統(tǒng)的邊緣檢測方法為什么要用小波來進(jìn)行邊緣檢測？邊緣檢測的分類局部模極大值點Canny準(zhǔn)則第四十九張，PPT共一百零八頁，創(chuàng)作于2022年6月小波多尺度局部模極大值邊緣檢測的原理小波多分辨率邊緣檢測的具體實現(xiàn)小波函數(shù)的選取自適應(yīng)閾值的選取利用邊緣信息進(jìn)行目標(biāo)定位仿真實驗第五十張，PPT共一百零八頁，創(chuàng)作于2022年6月8.7.1小波多尺度局部模極大值邊緣檢測的原理 假設(shè) 是二維平滑函數(shù)，且滿足

20、可把它沿 兩個方向上的導(dǎo)數(shù)作為基本小波：對于一幅圖像 ，其小波變換為： 第五十一張，PPT共一百零八頁，創(chuàng)作于2022年6月整個圖像的二進(jìn)小波變換即矢量：模值為 ：相角為 ：第五十二張，PPT共一百零八頁，創(chuàng)作于2022年6月8.7.2 小波多分辨率邊緣檢測的具體實現(xiàn)搜尋模極大值： 第五十三張，PPT共一百零八頁，創(chuàng)作于2022年6月噪聲的濾除：(1)閾值法 硬閾值，軟閾值，自適應(yīng)閾值；(2)多分辨率分割 利用模極大值在各個尺度的傳播特性去除噪聲； 利用小波變換尺度間相關(guān)性去除噪聲； 平移不變量的小波去噪方法；第五十四張，PPT共一百零八頁，創(chuàng)作于2022年6月邊緣跟蹤算法的四個約束條件：a)

21、“方向不變性”原則； b)角度平滑條件； c)幅值最鄰近條件； d)“互認(rèn)”原則。 邊緣點的鏈接：第五十五張，PPT共一百零八頁，創(chuàng)作于2022年6月三種小波多尺度局部模極大值 邊緣檢測方法的比較 方法一：小波變換模極大值用于邊緣檢測的原始方法具體步驟： 利用多孔算法對原圖像 進(jìn)行保持圖像大小不變的小波變換，生成水平分量 和垂直分量 ；計算出梯度向量的模值計算出梯度向量的相角按照相角量化方法劃分為8個方向，求出不同方向的局部模極大值點；對各尺度的邊緣圖像進(jìn)行閾值處理；鏈化模極大值點，去除長度小于一定閾值的邊緣鏈，就得到各 個尺度上的邊緣圖像。第五十六張，PPT共一百零八頁，創(chuàng)作于2022年6月

22、方法二：小波變換模極大值用于邊緣檢測的簡化方法 具體步驟：對所給圖像的每一行進(jìn)行小波變換，求出 ；對所給圖像的每一列進(jìn)行小波變換，求出 ； 計算出梯度向量的模值 計算出梯度向量的相角 將相角 按8方向量化，求出 在不同方向的局部模極大值點；對各尺度的邊緣圖像進(jìn)行閾值處理；鏈化模極大值點，去除長度小于一定閾值的邊緣鏈，就得到各 個尺度上的邊緣圖像。第五十七張，PPT共一百零八頁，創(chuàng)作于2022年6月方法三：小波變換模極大值的多尺度綜合方法具體步驟：求出各尺度的模圖像簇 和相角圖像簇 ；對各尺度的邊緣圖像進(jìn)行閾值處理；將相角 按8方向量化，求出 在不同方向的局部模極大值點；由粗到細(xì)的邊緣鏈接：對經(jīng)

23、過閾值處理的最粗分辨率級上的模極大值開始，鏈接模極大值圖像中模值相近，相角相似的非零像素點，刪除長度小于鏈長閾值的邊緣鏈，得到大尺度下單像素寬的圖像邊緣 ；第五十八張，PPT共一百零八頁，創(chuàng)作于2022年6月針對尺度 的每一個邊緣像素，搜索 尺度下以這點為中心的3X3鄰域，將該鄰域內(nèi)出現(xiàn)的所有可能邊緣點均標(biāo)記為候選邊緣點，得到 尺度下的候選邊緣點圖像 非候選邊緣點標(biāo)記為零；將 尺度下的候選邊緣點圖像 中模值相近，相角相似的非零像素點鏈接，刪除長度小于閾值的邊緣鏈，得到 尺度下單像素寬的圖像邊緣 ；重復(fù)步驟，直到 為止，邊緣圖像即為綜合后形成的邊緣圖像，也就是我們最終所要得到的邊緣圖像。第五十九

24、張，PPT共一百零八頁，創(chuàng)作于2022年6月8.7.3 小波函數(shù)的選取 在實際邊緣檢測中，小波函數(shù)的選取直接關(guān)系到邊緣檢測的結(jié)果，小波變換相當(dāng)于對圖像進(jìn)行帶通濾波，在一定程度上減少了噪聲對圖像的影響，但同時也去掉了一些模糊邊緣。這就要求尋找一種具有好的去噪特性同時又能精確地提取邊緣的小波函數(shù)，而且這種小波在滿足Canny準(zhǔn)則的同時應(yīng)具有限緊支撐特性。在關(guān)于邊緣提取的文獻(xiàn)中，主要采用的是正交小波，雙正交小波和B樣條小波作為邊緣提取的小波函數(shù)。通常選取高斯函數(shù)作為平滑函數(shù)。三次B樣條函數(shù)已很接近高斯函數(shù)，能滿足大多數(shù)的應(yīng)用，而它的緊支性優(yōu)于高斯函數(shù)。在有噪聲環(huán)境下進(jìn)行邊緣檢測時，需要平衡噪聲抑制能

25、力與邊緣定位精度，此時三次B樣條函數(shù)是最優(yōu)的。第六十張，PPT共一百零八頁，創(chuàng)作于2022年6月其中， 為閾值初值， 為比例系數(shù)，N為采樣點數(shù)。根據(jù) 和 的大小來決定窗口內(nèi)均值對閾值的作用。8.7.4 自適應(yīng)閾值選取選擇一個窗口在小波系數(shù)上滑動，窗口大小可以取3232或1616，將尺度下小于閾值 的梯度值置為0，自適應(yīng)閾值為：第六十一張，PPT共一百零八頁，創(chuàng)作于2022年6月8.7.5 利用邊緣信息進(jìn)行目標(biāo)定位 對經(jīng)過綜合后得出的圖像邊緣，根據(jù)形心公式就可以計算出圖像的形心坐標(biāo)，判斷出目標(biāo)在視場中的位置，實現(xiàn)目標(biāo)定位。 第六十二張，PPT共一百零八頁，創(chuàng)作于2022年6月8.7.6 實驗結(jié)果

26、（一）五個尺度的模極大值提取第一個尺度的噪聲很多第六十三張，PPT共一百零八頁，創(chuàng)作于2022年6月多尺度鏈接后的模極大值第六十四張，PPT共一百零八頁，創(chuàng)作于2022年6月實驗二:閾值的選取對邊緣檢測結(jié)果影響的實驗 用固定閾值在一階尺度提取的邊緣用自適應(yīng)閾值在一階尺度提取的邊緣用固定閾值在二階尺度提取的邊緣用自適應(yīng)閾值在二階尺度提取的邊緣 第六十五張，PPT共一百零八頁，創(chuàng)作于2022年6月實驗三：三種模極大值邊緣檢測方法的性能比較 邊緣檢測效果抑制噪聲能力(加入均值為0,方差為0.001的高斯噪聲)計算量模極大值邊緣檢測的簡化方法相對最差 抑制噪聲能力最弱最小模極大值邊緣檢測的原始算法相對

27、較好 抑制噪聲能力較強稍大模極大值邊緣檢測多尺度綜合法效果最好 抑制噪聲能力最強最大第六十六張，PPT共一百零八頁，創(chuàng)作于2022年6月實驗四：基于三次B樣條小波的邊緣檢測實驗 抑制噪聲方面(加入均值為0，方差為0.001的高斯噪聲)提取邊緣能力方面（對比度為2%的低對比度邊緣）計算量方面(基于雙DSP平臺)Roberts最差，無法分離邊緣與噪聲最差，邊緣定位不準(zhǔn)確 1.34ms Sobel性能優(yōu)于Roberts優(yōu)于Roberts 1.62ms拉普拉斯高斯算子優(yōu)于Roberts和Sobel，能部分檢測出邊緣效果優(yōu)于Roberts和Sobel 2.69ms小波模極大值方法優(yōu)于上面三種方法，能檢測

28、出主要邊緣很好，少部分提取的不準(zhǔn)確 18.9ms小波多尺度綜合方法效果最好，檢測出的邊緣受噪聲影響很小最好，幾乎能提取出全部微弱邊緣信息 59.2ms第六十七張，PPT共一百零八頁，創(chuàng)作于2022年6月由于小波變換具有多分辨分析特性和時頻局部化能力，在邊緣檢測、去噪和圖像增強等方面都具有很強的優(yōu)勢；更適合用來檢測受噪聲污染嚴(yán)重的模糊圖像和低對比度圖像，尤其對微弱目標(biāo)，它首先能抑制噪聲、增強對比度，然后利用多尺度的模極大值方法有效檢測出目標(biāo)邊緣，從而實現(xiàn)3對比度下的目標(biāo)精確定位；這種算法要進(jìn)行多尺度運算，所以計算量很大。 小 結(jié)第六十八張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從

29、傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪與增強中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第六十九張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8 小波變換在圖像去噪與增強中的應(yīng)用高頻系數(shù)置零的線性去噪方法小波系數(shù)硬閾值去噪方法小波系數(shù)軟閾值去噪方法小波系數(shù)自適應(yīng)閾值去噪方法基于小波模極大值的去噪方法基于信號奇異性的去噪方法基于小波系數(shù)相關(guān)性的去噪方法基于小波變換的圖像增強方法第七十張，PPT共一百零八頁，創(chuàng)作于2

30、022年6月8.8.1 高頻系數(shù)置零去噪1.對噪聲圖像進(jìn)行二維離散小波分解，分解層數(shù)一般取2或3 層，分解過程如下圖。2.對每一層的高頻系數(shù)LH，HL，HH，置零。3.對小波系數(shù)進(jìn)行重構(gòu)。這是最簡單的利用小波變換性質(zhì)的去噪方法。第七十一張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8.2 基于小波硬閾值的去噪方法1.首先將圖像信號求小波變換。2.除了最粗尺度信號外，將各細(xì)節(jié)信號作閾值處理，當(dāng)某位置小波變換值大于閾值時，保留原值，否則置零。即3.利用小波變換重構(gòu)，求出信號的濾波值。第七十二張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8.3 基于小波軟閾值的去噪方法1.首先對圖像信號進(jìn)行小波

31、變換，得出帶有噪聲的小波系數(shù)。2.將各細(xì)節(jié)信號作閾值處理，當(dāng)某位置小波變換值大于閾值時，作下面運算，下式中sgn(x)代表符號函數(shù) 否則置零。3.利用小波變換重構(gòu)，求出信號的濾波值。第七十三張，PPT共一百零八頁，創(chuàng)作于2022年6月軟硬閾值濾波器第七十四張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8.4 非線性軟閾值去噪方法第七十五張，PPT共一百零八頁，創(chuàng)作于2022年6月其中，噪聲方差的估計為 ， MAD為圖像中位值，n為信號采樣點數(shù)。閾值選取1.2.3.4.第七十六張，PPT共一百零八頁，創(chuàng)作于2022年6月對含有高斯噪聲的Lena圖像利用硬閾值法、軟閾值法和自適應(yīng)閾值法去噪后實驗

32、結(jié)果 第七十七張，PPT共一百零八頁，創(chuàng)作于2022年6月基于小波變換閾值方法去噪的不足 只適合高斯白噪聲的去噪，對椒鹽噪聲效果不明顯。原因在 于閾值法保留的是大于閾值的小波系數(shù)，而椒鹽噪聲在圖像 上表現(xiàn)為或者是灰度值特別大的白像素，或者特別小的黑像 素，椒鹽噪聲的小波系數(shù)都很大，所以不能用閾值分離出來。 閾值方法的去噪效果依賴于信噪比的大小，它特別適合信噪 比高的圖像去噪。 在圖像信號不連續(xù)點處會有偽吉普斯現(xiàn)象。 閾值方法的關(guān)鍵在于閾值的選取，而選擇一種普適性很好的 閾值選取方法是很困難的，事實上，人們已經(jīng)證明在均方誤 差意義上閾值方法能得到原始圖像信號的最優(yōu)估計，然而在 實際應(yīng)用中還是需要

33、根據(jù)具體的情況和經(jīng)驗來對一些閾值進(jìn) 行改進(jìn)。第七十八張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8.5 基于多尺度模極大值的小波去噪方法 根據(jù)圖像邊緣點與噪聲點具有不同的奇異性，則它們小波變換后的模極大值在不同尺度下的傳播特性也不同這一特性進(jìn)行圖像去噪。當(dāng)信號 在 處Lipschitz指數(shù)為 時， 反映了函數(shù)在點 處的奇異性大小。是當(dāng)信號在 處的奇異性 時，表明小波系數(shù)模極大值將隨尺度 的增大而增大；當(dāng) 時，則隨 的增大而減小。通常的圖像具有時域上的相關(guān)性，因而 ；而噪聲由于在該點不可導(dǎo)， 。因此，邊緣點的模值和噪聲點的模值隨尺度的變化具有不同的規(guī)律，所以可以利用這一特性將噪聲分離出來。第

34、七十九張，PPT共一百零八頁，創(chuàng)作于2022年6月基于多尺度模極大值的小波去噪方法的具體實現(xiàn) 通過分析圖像和噪聲在小波域中對應(yīng)的系數(shù)模極大值在不同尺度上的分布情況，來研究圖像和噪聲的突變點在不同尺度上的傳遞特性。在從低到高的分解尺度中，圖像突變點對應(yīng)的小波系數(shù)極大值具有傳遞性，而噪聲突變點不具有這種傳遞性。根據(jù)這一性質(zhì)，在確定出各尺度小波系數(shù)極大值的基礎(chǔ)上，由粗到精的跟蹤不同尺度上的小波系數(shù)極大值，并依據(jù)其不同尺度間的傳遞性，識別信噪屬性，剔除噪聲部分對應(yīng)的小波系數(shù)極大值，從而抑制噪聲，提高圖像質(zhì)量。第八十張，PPT共一百零八頁，創(chuàng)作于2022年6月模極大值方法去噪的過程1.小波系數(shù)極大值的確

35、定 2.圖像極大值跟蹤3.噪聲極大值濾除 第八十一張，PPT共一百零八頁，創(chuàng)作于2022年6月 模極大值多尺度鏈接方法對高斯噪聲 圖像去噪的實驗結(jié)果 對含有高斯噪聲（ 0.01）的Lena圖像利用多尺度模極大值鏈接方法去噪后圖像第八十二張，PPT共一百零八頁，創(chuàng)作于2022年6月 模極大值多尺度鏈接方法對椒鹽噪聲 圖像去噪的實驗結(jié)果對含有椒鹽噪聲（ 0.005）的圖像利用多尺度模極大值鏈接方法去噪后圖像第八十三張，PPT共一百零八頁，創(chuàng)作于2022年6月含有椒鹽噪聲的圖像、均值濾波去噪圖像、中值濾波去噪圖像、硬閾值去噪圖像、軟閾值去噪圖像和模極大值奇異性去噪圖像的實驗結(jié)果對比基于多尺度模極大值

36、奇異性對椒鹽噪聲的去噪實驗結(jié)果第八十四張，PPT共一百零八頁，創(chuàng)作于2022年6月模極大值多尺度鏈接方法對混和噪聲圖像的去噪實驗含有高斯和椒鹽混合噪聲圖像、硬閾值去噪圖像、軟閾值去噪圖像、自適應(yīng)閾值去噪圖像、中值濾波去噪圖像和模極大值多尺度鏈接去噪圖像第八十五張，PPT共一百零八頁，創(chuàng)作于2022年6月8.8.6 基于小波系數(shù)相關(guān)性的去噪方法 信號與噪聲的小波變換在各尺度下的不同傳播特性表明，信號的小波變換在各尺度間有較強的相關(guān)性，而且在邊緣處具有很強的相關(guān)性；而噪聲的小波變換在各尺度間確沒有明顯的相關(guān)性，而且噪聲的小波變換主要集中在小尺度各層次中。根據(jù)信號與噪聲的小波變換在不同尺度間的上述不

37、同特點，可以通過將相鄰尺度的小波系數(shù)直接相乘來增強信號，抑制噪聲。由于噪聲主要分布在小尺度上，所以這種現(xiàn)象在小尺度上非常明顯。 第八十六張，PPT共一百零八頁，創(chuàng)作于2022年6月首先定義 為尺度 上 點處的相關(guān)系數(shù)。為使相關(guān)系數(shù)與小波系數(shù)具有可比性，定義規(guī)范化相關(guān)系數(shù)：其中 和 分別表示對應(yīng)于尺度的小波系數(shù)與相關(guān)系數(shù)的能量。顯然，在尺度 下，小波系數(shù)與規(guī)范化相關(guān)系數(shù)具有相同的能量，這為它們之間提供了可比性。第八十七張，PPT共一百零八頁，創(chuàng)作于2022年6月基于小波系數(shù)相關(guān)法去噪的實驗結(jié)果（一）（ 0.01）Lena圖像和基于小波系數(shù)相關(guān)法去噪后的圖像第八十八張，PPT共一百零八頁，創(chuàng)作于2

38、022年6月小 結(jié)對小波系數(shù)使用閾值的去噪方法速度快，實現(xiàn)簡單，所以一直被廣泛使用在高斯噪聲去除上，其中的硬閾值方法容易產(chǎn)生偽吉普斯效應(yīng)，軟閾值方法克服了這個缺點，圖像比較平滑，去噪效果也比硬閾值好，而自適應(yīng)閾值使用的是局部閾值，比軟、硬閾值都好。模極大值多尺度鏈接和小波系數(shù)相關(guān)去噪兩種方法雖然速度都比較慢，但是它們都比閾值方法效果好，對各種噪聲效果都比較穩(wěn)定，特別是對混雜噪聲的去除方面。但是在椒鹽噪聲圖像的去噪方面，只有模極大值鏈接的方法有效果，而且比傳統(tǒng)的中值濾波方法要保留更多的邊緣信息。在對同時帶有高斯噪聲和椒鹽噪聲的圖像去噪時，模極大值鏈接的方法去噪效果非常突出，而且穩(wěn)定。第八十九張，

39、PPT共一百零八頁，創(chuàng)作于2022年6月 8.8.7 基于小波的圖像增強方法 小波多分辨分析由于它能多尺度多角度提取信號特征，并在不同尺度上讓噪聲和信號明顯地區(qū)分開來，所以它在圖像去噪和增強方面有很大優(yōu)勢?；谛〔ǘ喾直娣治龅膱D像增強，就是突出圖像的邊緣細(xì)節(jié)，盡可能的消除負(fù)面因素，從而達(dá)到增強圖像的目的。因此，基于小波分析的圖像增強可以轉(zhuǎn)化為兩步。 第一步是區(qū)分噪聲和圖像的邊緣細(xì)節(jié)； 第二步是根據(jù)需要對圖像的邊緣細(xì)節(jié)適度增強。 通過一個合適的增益函數(shù)完成上面所述的兩個步驟，使去噪和增強一次完成。在增強圖像的同時應(yīng)考慮到噪聲問題，所以對小波分解后高頻部分的處理很重要。為此，設(shè)計合適的增益函數(shù)是關(guān)

40、鍵的一步。 第九十張，PPT共一百零八頁，創(chuàng)作于2022年6月簡單的線性圖像增強的方法第一類是簡單的增益函數(shù)，即對每個方向的高頻采用不同的常數(shù) ；第二類是比較復(fù)雜的增益函數(shù)，有分段線性函數(shù)和非線性函數(shù)。 根據(jù)處理高頻增益函數(shù)的不同方法，將圖像增強分為線性增強和非線性增強 。第九十一張，PPT共一百零八頁，創(chuàng)作于2022年6月線性圖像增強方法的實現(xiàn)線性增強是對小波分解的高頻采取的常數(shù)增益，具體的算法是：a、對圖像進(jìn)行多級小波變換，得到 b、計算圖像的原始噪聲水平；C、計算圖像在各個尺度上不同方向的噪聲水平；d、去噪并增強。認(rèn)為幅值小于 的系數(shù)為噪聲，給予抑制，而幅值大于 的系數(shù)給予一定的增益。這

41、樣得到 e、對增強后的小波系數(shù)，進(jìn)行小波逆變換得到增強后的圖像。第九十二張，PPT共一百零八頁，創(chuàng)作于2022年6月基于邊緣提取的去噪和圖像增強相結(jié)合的實現(xiàn)1.對圖像進(jìn)行小波變換 2.計算小波變換的多尺度模和幅角； 3.用局部模極大值方法檢測出邊緣細(xì)節(jié)點；4.利用合適的非線性增益函數(shù)增強邊緣點，同時去除高頻噪聲點；5.對增強后的小波系數(shù)進(jìn)行小波逆變換得到增強后的圖像。 第九十三張，PPT共一百零八頁，創(chuàng)作于2022年6月基于小波變換的圖像增強方法的實驗結(jié)果對Lena圖像利用小波系數(shù)線性增強前后的效果對比第九十四張，PPT共一百零八頁，創(chuàng)作于2022年6月目 錄8.1 從傅里葉變換到小波變換的時頻分析法8.2 小波變換分類8.3 小波變換的多分辨分析特性8.4 尺度函數(shù)與小波8.5 小波變換的實現(xiàn)8.6 圖像的多分辨分解與重建8.7 小波變換在圖像邊緣檢測中的應(yīng)用8.8 小波變換在圖像去噪中的應(yīng)用8.9 小波變換在圖像融合中的應(yīng)用第九十五張，PPT共一百零八頁，創(chuàng)作于2022年6月8.9 小波變換在圖像融合中的應(yīng)用第九十六張，PPT共一百零八頁，創(chuàng)作于