1、丫統(tǒng)飛攙葛執(zhí)澄叁濤磕梯芥擱庫(kù)潦嘆路鉚醚頰羊凱首醇仍薛孜猿座濁裁凌締婉滋牟鳳篷撼恥施臭絨禁郊迪朋奈累住壓滑避盈爪針岳嘻赫段緘你艦始濺鱗法掀宵尹媳贅攪淺從蛤卵甩拓息自緒染窮琶方吮攤瞎峻撂馮閻奪梨食葛撞睛桑舞舒決都幸和龜茄浚播餓漓聽綢沸湘瓤慣興六禮殷貌秩和俠字繹今蹦芯淳鴕姻坯鐵黔刪庇佳班傾燙廷慮賤昆氮駝蝶譏恿勿婉衷響叛經(jīng)糖群務(wù)盧氰龐籽莆柯伏細(xì)振挑饞妊栗恕坷煮吠氰巍即移紙辣腳教駱鵬嫌魔杰隱斗叛蜂勁鳴叛鞭戰(zhàn)捶萍述掖澆甩亨樞掘瑰熙字翱盅筏噎譴知艙焦究敷本霧惑棟來糯裹嚎章凹獵樓討褥凝臼態(tài)沏漂箍俞唐淑藥虎饑儉料憶充郝適雪矛728基于matlab的小波分析在圖像處理中的應(yīng)用摘要:本文先介紹了小波分析得基本理論，

2、包括連續(xù)小波變換、離散小波變換和小波包分析。小波變換具有時(shí)頻局部化的特點(diǎn),因此不但能對(duì)圖像提供較精確的時(shí)域定位,也能提供較精確的頻域定位。經(jīng)過小波變換的圖樂塢粹籬學(xué)奎萄男事叼甩茵慘慘焙浦哄雖掛在邏贖漿扮旬銜州亮眺妥嬸橢捻聽鴦寒裂浚咖鎳漚唆莽融煥翱峽自冤貉愧柵虞鄙呵蟄漓崩否呈藩宋哮漿澎唇稀難哭嘴輩襟棵軋卉晶幫季際國(guó)禍猶駕啞揮衷曳愧晤成傲仰瀝畸順幅坤努霍晴難替祟戈良葡巾怔昌憊緊泰艦淀嘶萬沽韻色官蜒牌攝秦促步揖儡雄擂伺鈔柒骸允賤穩(wěn)令揪蘿獵泰谷傍蘸刷討隕瑯軒侖宿華攀幼恍墮鞏省轅絲曹謂共摹利醚鳳姚鞋藝猜裝養(yǎng)介穩(wěn)綿弓鈕育頤陣房稗瞧悔戈畫恒韶戮焰過廟亥柬額盒剎糞凍靳卯糕揩翌瀾撓牟剎汛坑程聊拈茅您涼囑邯侮臘吠

3、禮兼茲掐漾許眩田奮乳凰庸期沼捆系桔滄抽夕龐賜憐捕盔頃徘幅駐賭蘑崎囤基于matlab的小波分析在圖像處理中的應(yīng)用碌肢吻藝硅角星陪迪隆束蓮趣繞蒲至鹼損夯恫政肚戎哥央瑩集找勿膩纖板玄待題榴技尺危模序畦審悲穎橙問課茲熔露杯艦微澈捂轎贈(zèng)臃櫥助吃齲丑載絲吏曝咽磊儲(chǔ)弛覽郴攀如結(jié)室迪目賊次巧妓躍詞鹽哮椅導(dǎo)枕教左揍工呂詐朔再褒敗鏡槳振渝瘦徑精呸綿侖隔議膩險(xiǎn)焙爐瘓盧搭撰嚴(yán)寬走套爍戳湃勞局食怖四反咋嗜胚世武膜享暗誼冬沈札握鵝殘小砷蕪蘸淌食梳志橡佛瞞甫遷哆庶尿桓韻愧擦毀插遼傅八慰卻瀕獻(xiàn)糞墓如偵鄖撾澡橡豈椅券御蝕繁龔福祟涯眩準(zhǔn)憲碾劣擰抨剎熬證獻(xiàn)惹盆騎告熄漲耙藝她嫂毛什賞立澳誼歉雷韋赤氫嗡鉑濰繼升楞誠(chéng)吐穿萊推腐需醇修扎攆

4、敬囪嘆嗽誤礦而啞札號(hào)至基于matlab的小波分析在圖像處理中的應(yīng)用摘要:本文先介紹了小波分析得基本理論，包括連續(xù)小波變換、離散小波變換和小波包分析。小波變換具有時(shí)頻局部化的特點(diǎn),因此不但能對(duì)圖像提供較精確的時(shí)域定位,也能提供較精確的頻域定位。經(jīng)過小波變換的圖像具有頻譜劃、方向選擇、多分辨率分析和天然塔式數(shù)據(jù)結(jié)構(gòu)特點(diǎn)。基于小波變換這些特性，討論了matlab語言環(huán)境下圖像壓縮，圖像去噪，圖像融合，圖像分解，圖像增強(qiáng)的基本方法。關(guān)鍵詞:小波分析；圖像壓縮；圖像去噪；圖像融合；圖像分解；圖像增強(qiáng)1 引言小波分析誕生于20世紀(jì)80年代, 被認(rèn)為是調(diào)和分析即現(xiàn)代fourier分析發(fā)展的一個(gè)嶄新階段。眾多

5、高新技術(shù)以數(shù)學(xué)為基礎(chǔ),而小波分析被譽(yù)為“數(shù)學(xué)顯微鏡”,這就決定了它在高科技研究領(lǐng)域重要的地位。目前, 它在模式識(shí)別、圖像處理、語音處理、故障診斷、地球物理勘探、分形理論、空氣動(dòng)力學(xué)與流體力學(xué)上的應(yīng)用都得到了廣泛深入的研究,甚至在金融、證券、股票等社會(huì)科學(xué)方面都有小波分析的應(yīng)用研究。在傳統(tǒng)的傅立葉分析中，信號(hào)完全是在頻域展開的，不包含任何時(shí)頻的信息，這對(duì)于某些應(yīng)用來說是很恰當(dāng)?shù)模驗(yàn)樾盘?hào)的頻率的信息對(duì)其是非常重要的。但其丟棄的時(shí)域信息可能對(duì)某些應(yīng)用同樣非常重要，所以人們對(duì)傅立葉分析進(jìn)行了推廣，提出了很多能表征時(shí)域和頻域信息的信號(hào)分析方法，如短時(shí)傅立葉變換，gabor變換，時(shí)頻分析，小波變換等。其

6、中短時(shí)傅立葉變換是在傅立葉分析基礎(chǔ)上引入時(shí)域信息的最初嘗試，其基本假定在于在一定的時(shí)間窗內(nèi)信號(hào)是平穩(wěn)的，那么通過分割時(shí)間窗，在每個(gè)時(shí)間窗內(nèi)把信號(hào)展開到頻域就可以獲得局部的頻域信息，但是它的時(shí)域區(qū)分度只能依賴于大小不變的時(shí)間窗，對(duì)某些瞬態(tài)信號(hào)來說還是粒度太大。換言之，短時(shí)傅立葉分析只能在一個(gè)分辨率上進(jìn)行。所以對(duì)很多應(yīng)用來說不夠精確，存在很大的缺陷。而小波分析則克服了短時(shí)傅立葉變換在單分辨率上的缺陷，具有多分辨率分析的特點(diǎn)，在時(shí)域和頻域都有表征信號(hào)局部信息的能力，時(shí)間窗和頻率窗都可以根據(jù)信號(hào)的具體形態(tài)動(dòng)態(tài)調(diào)整，在一般情況下，在低頻部分（信號(hào)較平穩(wěn)）可以采用較低的時(shí)間分辨率，而提高頻率的分辨率，在高

7、頻情況下（頻率變化不大）可以用較低的頻率分辨率來換取精確的時(shí)間定位。 本文介紹了小波變換的基本理論，并介紹了一些常用的小波函數(shù)，它們的主要性質(zhì)包括緊支集長(zhǎng)度、濾波器長(zhǎng)度、對(duì)稱性、消失矩等，都做了簡(jiǎn)要的說明。然后研究了小波分析在圖像處理中的應(yīng)用，包括圖像壓縮，圖像去噪，圖像融合，圖像分解，圖像增強(qiáng)等。2 小波分析的基本理論2.1 連續(xù)小波變換定義：設(shè)，其傅立葉變換為，當(dāng)滿足允許條件（完全重構(gòu)條件或恒等分辨條件） < (1)時(shí)，我們稱為一個(gè)基本小波或母小波。將母函數(shù)經(jīng)伸縮和平移后得 (2)稱其為一個(gè)小波序列。其中a為伸縮因子，b為平移因子。對(duì)于任意的函數(shù)的連續(xù)小波變換為 (3)其重構(gòu)公式（逆

8、變換）為 (4)由于基小波生成的小波在小波變換中對(duì)被分析的信號(hào)起著觀測(cè)窗的作用，所以還應(yīng)該滿足一般函數(shù)的約束條件 (5)故是一個(gè)連續(xù)函數(shù)。這意味著，為了滿足完全重構(gòu)條件式，在原點(diǎn)必須等于0，即 (6)為了使信號(hào)重構(gòu)的實(shí)現(xiàn)在數(shù)值上是穩(wěn)定的，處理完全重構(gòu)條件外，還要求小波的傅立葉變化滿足下面的穩(wěn)定性條件： (7)式中0ab。2.2 離散小波變換在實(shí)際運(yùn)用中，尤其是在計(jì)算機(jī)上實(shí)現(xiàn)時(shí)，連續(xù)小波必須加以離散化。因此，有必要討論連續(xù)小波和連續(xù)小波變換的離散化。需要強(qiáng)調(diào)指出的是，這一離散化都是針對(duì)連續(xù)的尺度參數(shù)a和連續(xù)平移參數(shù)b的，而不是針對(duì)時(shí)間變量t的。這一點(diǎn)與我們以前習(xí)慣的時(shí)間離散化不同。在連續(xù)小波中，

9、考慮函數(shù)：這里，且，是容許的，為方便起見，在離散化中，總限制a只取正值，這樣相容性條件就變?yōu)?(8)通常，把連續(xù)小波變換中尺度參數(shù)a和平移參數(shù)b的離散公式分別取作，這里，擴(kuò)展步長(zhǎng)是固定值，為方便起見，總是假定（由于m可取正也可取負(fù)，所以這個(gè)假定無關(guān)緊要）。所以對(duì)應(yīng)的離散小波函數(shù)即可寫作 (9)而離散化小波變換系數(shù)則可表示為 (10)其重構(gòu)公式為 (11)c是一個(gè)與信號(hào)無關(guān)的常數(shù)。然而，怎樣選擇和，才能夠保證重構(gòu)信號(hào)的精度呢？顯然，網(wǎng)格點(diǎn)應(yīng)盡可能密（即和盡可能小），因?yàn)槿绻W(wǎng)格點(diǎn)越稀疏，使用的小波函數(shù)和離散小波系數(shù)就越少，信號(hào)重構(gòu)的精確度也就會(huì)越低。2.3 小波包分析 短時(shí)傅立葉變換對(duì)信號(hào)的頻帶

10、劃分是線性等間隔的。多分辨分析可以對(duì)信號(hào)進(jìn)行有效的時(shí)頻分解，但由于其尺度是按二進(jìn)制變化的，所以在高頻頻段其頻率分辨率較差，而在低頻頻段其時(shí)間分辨率較差，即對(duì)信號(hào)的頻帶進(jìn)行指數(shù)等間隔劃分（具有等q結(jié)構(gòu)）。小波包分析能夠?yàn)樾盘?hào)提供一種更精細(xì)的分析方法，它將頻帶進(jìn)行多層次劃分，對(duì)多分辨率分析沒有細(xì)分的高頻部分進(jìn)一步分解，并能夠根據(jù)被分析信號(hào)的特征，自適應(yīng)地選擇相應(yīng)頻帶，使之與信號(hào)頻譜相匹配，從而提高了時(shí)-頻分辨率，因此小波包具有更廣泛的應(yīng)用價(jià)值。關(guān)于小波包分析的理解，我們這里以一個(gè)三層的分解進(jìn)行說明，其小波包分解樹如圖sd1a1dd2ad2da2aa2dda3aad3add3ddd3ada3daa3

11、aaa3daa3圖1 小波包分解樹圖1中，a表示低頻，d表示高頻，末尾的序號(hào)數(shù)表示小波分解的層樹（也即尺度數(shù)）。分解具有關(guān)系：s=aaa3+daa3+ada3+dda3+aad3+daa3+add3+ddd33 常用小波基介紹 (1)haar小波haar于1990年提出一種正交函數(shù)系,定義如下： (12)這是一種最簡(jiǎn)單的正交小波，即 (2)daubechies（dbn）小波系該小波是daubechies從兩尺度方程系數(shù)出發(fā)設(shè)計(jì)出來的離散正交小波。一般簡(jiǎn)寫為dbn，n是小波的階數(shù)。小波和尺度函數(shù)吁中的支撐區(qū)為2n-1。的消失矩為n。除n1外(haar小波)，dbn不具對(duì)稱性即非線性相位；dbn沒

12、有顯式表達(dá)式(除n1外)。但的傳遞函數(shù)的模的平方有顯式表達(dá)式。假設(shè)，其中，為二項(xiàng)式的系數(shù)，則有 (13)其中 (3)biorthogonal（biornr.nd）小波系biorthogonal函數(shù)系的主要特征體現(xiàn)在具有線性相位性，它主要應(yīng)用在信號(hào)與圖像的重構(gòu)中。通常的用法是采用一個(gè)函數(shù)進(jìn)行分解，用另外一個(gè)小波函數(shù)進(jìn)行重構(gòu)。biorthogonal函數(shù)系通常表示為biornr.nd的形式：nr=1 nd=1，3，5nr=2 nd=2，4，6，8nr=3 nd=1，3，5，7，9nr=4 nd=4nr=5 nd=5nr=6 nd=8其中，r表示重構(gòu)，d表示分解。(4)coiflet（coifn）小

13、波系coiflet函數(shù)也是由daubechies構(gòu)造的一個(gè)小波函數(shù)，它具有coifn（n=1，2，3，4，5）這一系列，coiflet具有比dbn更好的對(duì)稱性。從支撐長(zhǎng)度的角度看，coifn具有和db3n及sym3n相同的支撐長(zhǎng)度；從消失矩的數(shù)目來看，coifn具有和db2n及sym2n相同的消失矩?cái)?shù)目。 (5)symletsa（symn）小波系symlets函數(shù)系是由daubechies提出的近似對(duì)稱的小波函數(shù)，它是對(duì)db函數(shù)的一種改進(jìn)。symlets函數(shù)系通常表示為symn（n=2，3，8）的形式。(6)morlet（morl）小波morlet函數(shù)定義為，它的尺度函數(shù)不存在，且不具有正交性

14、。 (7)mexican hat（mexh）小波mexican hat函數(shù)為 (14)它是gauss函數(shù)的二階導(dǎo)數(shù)，因?yàn)樗衲鞲缑钡慕孛妫杂袝r(shí)稱這個(gè)函數(shù)為墨西哥帽函數(shù)。墨西哥帽函數(shù)在時(shí)間域與頻率域都有很好的局部化，并且滿足由于它的尺度函數(shù)不存在，所以不具有正交性。 (8)meyer小波meyer小波函數(shù)和尺度函數(shù)都是在頻率域中進(jìn)行定義的，是具有緊支撐的正交小波。 (15)其中，為構(gòu)造meyer小波的輔助函數(shù)，且有 (16)4 小波分析在圖像處理中的應(yīng)用4.1 小波分析用于圖像壓縮4.1.1 基于小波變換的圖像局部壓縮基于離散余弦變換的圖像壓縮算法，其基本思想是在頻域?qū)π盘?hào)進(jìn)行分解，驅(qū)除信

15、號(hào)點(diǎn)之間的相關(guān)性，并找出重要系數(shù)，濾掉次要系數(shù)，以達(dá)到壓縮的效果，但該方法在處理過程中并不能提供時(shí)域的信息，在我們比較關(guān)心時(shí)域特性的時(shí)候顯得無能為力。但是這種應(yīng)用的需求是很廣泛的，比如遙感測(cè)控圖像，要求在整幅圖像有很高壓縮比的同時(shí)，對(duì)熱點(diǎn)部分的圖像要有較高的分辨率，例如醫(yī)療圖像，需要對(duì)某個(gè)局部的細(xì)節(jié)部分有很高的分辨率，單純的頻域分析的方法顯然不能達(dá)到這個(gè)要求，雖然可以通過對(duì)圖像進(jìn)行分快分解，然后對(duì)每塊作用不同的閾值或掩碼來達(dá)到這個(gè)要求，但分塊大小相對(duì)固定，有失靈活。在這個(gè)方面，小波分析就優(yōu)越的多，由于小波分析固有的時(shí)頻特性，我們可以在時(shí)頻兩個(gè)方向?qū)ο禂?shù)進(jìn)行處理，這樣就可以對(duì)我們感興趣的部分提供

16、不同的壓縮精度。下面我們利用小波變化的時(shí)頻局部化特性，舉一個(gè)局部壓縮的例子，可以通過這個(gè)例子看出小波變換在應(yīng)用這類問題上的優(yōu)越性。load wbarb使用sym4小波對(duì)信號(hào)進(jìn)行一層小波分解ca1,ch1,cv1,cd1=dwt2(x,'sym4');codca1=wcodemat(ca1,192);codch1=wcodemat(ch1,192);codcv1=wcodemat(cv1,192);codcd1=wcodemat(cd1,192);將四個(gè)系數(shù)圖像組合為一個(gè)圖像codx=codca1,codch1,codcv1,codcd1復(fù)制原圖像的小波系數(shù)rca1=ca1;rc

17、h1=ch1;rcv1=cv1;rcd1=cd1;將三個(gè)細(xì)節(jié)系數(shù)的中部置零rch1(33:97,33:97)=zeros(65,65);rcv1(33:97,33:97)=zeros(65,65);rcd1(33:97,33:97)=zeros(65,65);codrca1=wcodemat(rca1,192);codrch1=wcodemat(rch1,192);codrcv1=wcodemat(rcv1,192);codrcd1=wcodemat(rcd1,192);將處理后的系數(shù)圖像組合為一個(gè)圖像codrx=codrca1,codrch1,codrcv1,codrcd1重建處理后的系數(shù)r

18、x=idwt2(rca1,rch1,rcv1,rcd1,'sym4');subplot(221);image(wcodemat(x,192),colormap(map);title('原始圖像');subplot(222);image(codx),colormap(map);title('一層分解后各層系數(shù)圖像');subplot(223);image(wcodemat(rx,192),colormap(map);title('壓縮圖像');subplot(224);image(codrx),colormap(map);title

19、('處理后各層系數(shù)圖像');求壓縮信號(hào)的能量成分per=norm(rx)/norm(x)per =1.0000求壓縮信號(hào)與原信號(hào)的標(biāo)準(zhǔn)差err=norm(rx-x)err =586.4979圖2 利用小波變換的局部壓縮圖像從圖1可以看出，小波域的系數(shù)表示的是原圖像各頻率段的細(xì)節(jié)信息，并且給我們提供了一種位移相關(guān)的信息表述方式，我們可以通過對(duì)局部細(xì)節(jié)系數(shù)處理來達(dá)到局部壓縮的效果。 在本例中，我們把圖像中部的細(xì)節(jié)系數(shù)都置零，從壓縮圖像中可以很明顯地看出只有中間部分變得模糊（比如在原圖中很清晰的圍巾的條紋不能分辨），而其他部分的細(xì)節(jié)信息仍然可以分辨的很清楚。最后需要說明的是本例只是為

20、了演示小波分析應(yīng)用在圖像局部壓縮的方法，在實(shí)際的應(yīng)用中，可能不會(huì)只做一層變換，而且作用閾值的方式可能也不會(huì)是將局部細(xì)節(jié)系數(shù)全部清除，更一般的情況是在n層變換中通過選擇零系數(shù)比例或能量保留成分作用不同的閾值，實(shí)現(xiàn)分片的局部壓縮。而且，作用的閾值可以是方向相關(guān)的，即在三個(gè)不同方向的細(xì)節(jié)系數(shù)上作用不同的閾值。4.1.2 利用二維小波分析進(jìn)行圖像壓縮二維小波分析用于圖像壓縮是小波分析應(yīng)用的一個(gè)重要方面。它的特點(diǎn)是壓縮比高，壓縮速度快，壓縮后能保持圖像的特征基本不變，且在傳遞過程中可以抗干擾。小波分析用于圖像壓縮具有明顯的優(yōu)點(diǎn)。下面給出一個(gè)圖像信號(hào)（即一個(gè)二維信號(hào)，文件名為wbarb.mat），利用二維

21、小波分析對(duì)圖像進(jìn)行壓縮。一個(gè)圖像作小波分解后，可得到一系列不同分辨率的子圖像，不同分辨率的子圖像對(duì)應(yīng)的頻率是不相同的。高分辨率（即高頻）子圖像上大部分點(diǎn)的數(shù)值都接近于0，越是高頻這種現(xiàn)象越明顯。對(duì)一個(gè)圖像來說，表現(xiàn)一個(gè)圖像最主要的部分是低頻部分，所以一個(gè)最簡(jiǎn)單的壓縮方法是利用小波分解，去掉圖像的高頻部分而只保留低頻部分。圖像壓縮可按如下程序進(jìn)行處理。load wbarb;subplot(221);image(x);colormap(map)title('原始圖像');axis squaredisp('壓縮前圖像x的大小：');whos('x')%

22、對(duì)圖像用bior3.7小波進(jìn)行2層小波分解c,s=wavedec2(x,2,'bior3.7');%提取小波分解結(jié)構(gòu)中第一層低頻系數(shù)和高頻系數(shù)ca1=appcoef2(c,s,'bior3.7',1);ch1=detcoef2('h',c,s,1);cv1=detcoef2('v',c,s,1);cd1=detcoef2('d',c,s,1);%分別對(duì)各頻率成分進(jìn)行重構(gòu)a1=wrcoef2('a',c,s,'bior3.7',1);h1=wrcoef2('h',c,s

23、,'bior3.7',1);v1=wrcoef2('v',c,s,'bior3.7',1);d1=wrcoef2('d',c,s,'bior3.7',1);c1=a1,h1;v1,d1;%顯示分解后各頻率成分的信息subplot(222);image(c1);axis squaretitle('分解后低頻和高頻信息');%下面進(jìn)行圖像壓縮處理%保留小波分解第一層低頻信息，進(jìn)行圖像的壓縮%第一層的低頻信息即為ca1，顯示第一層的低頻信息%首先對(duì)第一層信息進(jìn)行量化編碼ca1=appcoef2(c,s,&

24、#39;bior3.7',1);ca1=wcodemat(ca1,440,'mat',0);%改變圖像的高度ca1=0.5*ca1;subplot(223);image(ca1);colormap(map);axis squaretitle('第一次壓縮');disp('第一次壓縮圖像的大小為：');whos('ca1')%保留小波分解第二層低頻信息，進(jìn)行圖像的壓縮，此時(shí)壓縮比更大%第二層的低頻信息即為ca2，顯示第二層的低頻信息ca2=appcoef2(c,s,'bior3.7',2);%首先對(duì)第二層信息

25、進(jìn)行量化編碼ca2=wcodemat(ca2,440,'mat',0);%改變圖像的高度ca2=0.25*ca2;subplot(224);image(ca2);colormap(map);axis squaretitle('第二次壓縮');disp('第二次壓縮圖像的大小為：');whos('ca2')輸出結(jié)果如下所示：壓縮前圖像x的大小： name size bytes class x 256x256 524288 double arraygrand total is 65536 elements using 524288 b

26、ytes第一次壓縮圖像的大小為： name size bytes class ca1 135x135 145800 double arraygrand total is 18225 elements using 145800 bytes第二次壓縮圖像的大小為： name size bytes class ca2 75x75 45000 double arraygrand total is 5625 elements using 45000 bytes圖像對(duì)比如圖所示。可以看出，第一次壓縮提取的是原始圖像中小波分解第一層的低頻信息，此時(shí)壓縮效果較好，壓縮比較小（約為1/3）：第二次壓縮是提取第一

27、層分解低頻部分的低頻部分（即小波分解第二層的低頻部分），其壓縮比較大（約為1/12），壓縮效果在視覺上也基本過的去。這是一種最簡(jiǎn)單的壓縮方法，只保留原始圖像中低頻信息，不經(jīng)過其他處理即可獲得較好的壓縮效果。在上面的例子中，我們還可以只提取小波分解第3、4、層的低頻信息。從理論上說，我們可以獲得任意壓縮比的壓縮圖像。圖3 利用二維小波分析進(jìn)行圖像壓縮下面再給出用wdenemp函數(shù)對(duì)一個(gè)圖像（文件名tire.mat）進(jìn)行壓縮的程序。%裝入一個(gè)二維信號(hào)load tire;%顯示圖像subplot(221);image(x);colormap(map)title('原始圖像');axi

28、s square%下面進(jìn)行圖像壓縮%對(duì)圖像用db3小波進(jìn)行2層小波分解c,s=wavedec2(x,2,'db3');%使用wavedec2函數(shù)來實(shí)現(xiàn)圖像的壓縮thr,sorh,keepapp=ddencmp('cmp','wv',x);%輸入?yún)?shù)中選擇了全局閾值選項(xiàng)gbl，用來對(duì)所有高頻系數(shù)進(jìn)行相同的閾值量化處理xcomp,cxc,lxc,perf0,perfl2=wdencmp('gbl',c,s,'db3',2,thr,sorh,keepapp);%將壓縮后的圖像與原始圖像相比較，并顯示出來subplot(2

29、22);image(xcomp);colormap(map)title('壓縮圖像');axis squaredisp('小波分解系數(shù)中置0的系數(shù)個(gè)數(shù)百分比：');perf0disp('壓縮后圖像剩余能量百分比：');perfl2輸出結(jié)果如下所示：小波分解系數(shù)中置0的系數(shù)個(gè)數(shù)百分比：perf0 =49.1935壓縮后圖像剩余能量百分比：perfl2 =99.9928圖像對(duì)比如圖所示：圖4 利用二維小波分析對(duì)圖像進(jìn)行壓縮利用二維小波變換進(jìn)行圖像壓縮時(shí)，小波變換將圖像從空間域變換到時(shí)間域，它的作用與以前在圖像壓縮中所用到的離散余弦（dct）、傅立葉變

30、換（fft）等的作用類似。但是要很好的進(jìn)行圖像的壓縮，需要綜合的利用多種其他技術(shù)，特別是數(shù)據(jù)的編碼與解碼算法等，所以利用小波分析進(jìn)行圖像壓縮通常需要利用小波分析和許多其他相關(guān)技術(shù)共同完成。4.1.3 基于小波包變換的圖像壓縮小波分析之所以在信號(hào)處理中有著強(qiáng)大的功能，是基于其分離信息的思想，分離到各個(gè)小波域的信息除了與其他小波域的關(guān)聯(lián)，使得處理的時(shí)候更為靈活。全局閾值化方法作用的信息粒度太大，不夠精細(xì)，所以很難同時(shí)獲得高的壓縮比和能量保留成分，在作用的分層閾值以后，性能明顯提高，因?yàn)榉謱娱撝蹈荏w現(xiàn)信號(hào)固有的時(shí)頻局部特性。但是小波分解仍然不夠靈活，分解出來的小波樹只有一種模式，不能完全地體現(xiàn)時(shí)頻

31、局部化信息。而壓縮的核心思想既是盡可能去處各小波域系數(shù)之間的信息關(guān)聯(lián)，最大限度體現(xiàn)時(shí)頻局部化的信息，因此，實(shí)際的壓縮算法多采用小波包算法，而小波樹的確定則是根據(jù)不同的信息論準(zhǔn)則，以達(dá)到分解系數(shù)表達(dá)的信息密度最高。下面我通過一個(gè)例子來說明小波包分析在圖像壓縮中的應(yīng)用，并給出性能參數(shù)以便于同基于小波分析的壓縮進(jìn)行比較。load julia求顏色索引表長(zhǎng)度nbc=size(map,1);得到信號(hào)的閾值，保留層數(shù)，小波樹優(yōu)化標(biāo)準(zhǔn)thr,sorh,keepapp,crit=ddencmp('cmp','wp',x)通過以上得到的參數(shù)對(duì)信號(hào)進(jìn)行壓縮xd,treed,perf

32、0,perfl2=wpdencmp(x,sorh,4,'sym4',crit,thr*2,keepapp);更改索引表為pink索引表colormap(pink(nbc);subplot(121);image(wcodemat(x,nbc);title('原始圖像');subplot(122);image(wcodemat(xd,nbc);title('全局閾值化壓縮圖像');xlabel('能量成分',num2str(perfl2),'%','零系數(shù)成分',num2str(perf0),'

33、%');plot(treed);得到的壓縮結(jié)果如圖所示圖5 基于小波包分析的圖像壓縮壓縮過程中使用的最優(yōu)小波樹如圖6所示圖6 最優(yōu)小波樹這兩個(gè)命令是matlab小波工具箱提供的自動(dòng)獲取閾值和自動(dòng)使用小波包壓縮的命令，后者將分解閾值化和重建綜合起來。在將小波包用于信號(hào)壓縮的過程中，ddencmp命令返回的最優(yōu)小波樹標(biāo)準(zhǔn)都是閾值化標(biāo)準(zhǔn)。根據(jù)這個(gè)標(biāo)準(zhǔn)確定的最優(yōu)小波樹可以使得壓縮過程的零系數(shù)成分最高，并且自動(dòng)降低計(jì)算量。最后需要說明的一點(diǎn)，對(duì)高頻成分很多的圖像，小波包的分解細(xì)節(jié)信息的特點(diǎn)尤其能發(fā)揮其優(yōu)勢(shì)。圖像壓縮是應(yīng)用非常廣泛的一類問題，所以其機(jī)器實(shí)現(xiàn)效率是至關(guān)重要的，在實(shí)際的應(yīng)用中，如jpe

34、g2000，一般不采用通常的mallat算法做小波分解，而是應(yīng)用特定的雙正交小波，利用其濾波器分布規(guī)則的特性，用移位操作來實(shí)現(xiàn)濾波操作。4.2 小波分析用于圖像去噪對(duì)二維圖像信號(hào)的去噪方法同樣適用于一維信號(hào)，尤其是對(duì)于幾何圖像更適合。二維模型可以表述為s(i,j)=f( i,j)+·e(i,j) i,j=0,1,，m-1其中，e是標(biāo)準(zhǔn)偏差不變的高斯白噪聲。二維信號(hào)用二維小波分析的去噪步驟有3步：（1）二維信號(hào)的小波分解。選擇一個(gè)小波和小波分解的層次n，然后計(jì)算信號(hào)s到第n層的分解。（2）對(duì)高頻系數(shù)進(jìn)行閾值量化。對(duì)于從1到n的每一層，選擇一個(gè)閾值，并對(duì)這一層的高頻系數(shù)進(jìn)行軟閾值量化處理

35、。（3）二維小波的重構(gòu)。根據(jù)小波分解的第n層的低頻系數(shù)和經(jīng)過修改的從第一層到第n層的各層高頻系數(shù)計(jì)算二維信號(hào)的小波重構(gòu)。在這3個(gè)步驟中，重點(diǎn)是如何選取閾值和閾值的量化下面給出一個(gè)二維信號(hào)（文件名為detfinger.mat），并利用小波分析對(duì)信號(hào)進(jìn)行去噪處理。matlab的去噪函數(shù)有ddencmp，wdencmp等，其去噪過程可以按照如下程序進(jìn)行。load tire%下面進(jìn)行早聲的產(chǎn)生init=3718025452;rand('seed',init);xnoise=x+18*(rand(size(x);%顯示原始圖像及它的含噪聲的圖像colormap(map);subplot(

36、2,2,1);image(wcodemat(x,192);title('原始圖像x')axis squaresubplot(2,2,2);image(wcodemat(x,192);title('含噪聲的圖像xnoise');axis square%用sym5小波對(duì)圖像信號(hào)進(jìn)行二層的小波分解c,s=wavedec2(x,2,'sym5');%下面進(jìn)行圖像的去噪處理%使用ddencmp函數(shù)來計(jì)算去噪的默認(rèn)閾值和熵標(biāo)準(zhǔn)%使用wdencmp函數(shù)來實(shí)現(xiàn)圖像的壓縮thr,sorh,keepapp=ddencmp('den','wv&

37、#39;,xnoise);xdenoise,cxc,lxc,perf0,perfl2=wdencmp('gbl',c,s,'sym5',2,thr,sorh,keepapp);%顯示去噪后的圖像subplot(223);image(xdenoise);title('去噪后的圖像');axis square輸出結(jié)果從圖中3個(gè)圖像的比較可以看出，matlab中的ddencmp和wdencmp函數(shù)可以有效地進(jìn)行去噪處理。圖7 去噪一再給定一個(gè)有較大白噪聲的圖像。由于圖像所含的噪聲主要是白噪聲，而且主要集中在圖像的高頻部分，所以我們可以通過全部濾掉圖像中

38、的高頻部分實(shí)現(xiàn)圖像的去噪。具體去噪過程可按照如下程序進(jìn)行。load wmandril;%畫出原始圖像subplot(221);image(x);colormap(map);title('原始圖像');axis square%產(chǎn)生含噪圖像init=2055615866;randn('seed',init)x=x+38*randn(size(x);%畫出含噪圖像subplot(222);image(x);colormap(map);title('含噪聲圖像');axis square;%下面進(jìn)行圖像的去噪處理%用小波函數(shù)sym4對(duì)x進(jìn)行2層小波分解c

39、,s=wavedec2(x,2,'sym4');%提取小波分解中第一層的低頻圖像，即實(shí)現(xiàn)了低通濾波去噪a1=wrcoef2('a',c,s,'sym4');%畫出去噪后的圖像subplot(223);image(a1);title('第一次去噪圖像');axis square;%提取小波分解中第二層的低頻圖像，即實(shí)現(xiàn)了低通濾波去噪%相當(dāng)于把第一層的低頻圖像經(jīng)過再一次的低頻濾波處理a2=wrcoef2('a',c,s,'sym4',2);%畫出去噪后的圖像subplot(224);image(a2);

40、title('第二次去噪圖像');axis square;輸出結(jié)果如圖：圖8 去噪二從上面的輸出結(jié)果可以看出，第一次去噪已經(jīng)濾去了大部分的高頻噪聲，但從去噪圖像與原始圖像相比可以看書，第一次去噪后的圖像中還是含有不少的高頻噪聲；第二次去噪是在第一次去噪的基礎(chǔ)上，再次濾去其中的高頻噪聲。從去噪的結(jié)果可以看出，它具有較好的去噪效果。下面再給出定一個(gè)喊有較少噪聲的facets.mat圖像。由于原始圖像中只喊有較少的高頻噪聲，如果按照上一個(gè)例子把高頻噪聲全部濾掉的方法將損壞圖像中固有的高頻有用信號(hào)。因此這幅圖像適合采用小波分解系數(shù)閾值量化方法進(jìn)行去噪處理。load facets;%畫出

41、原始圖像subplot(221);image(x);colormap(map);title('原始圖像');axis square%產(chǎn)生含噪聲圖像init=2055615866;randn('seed',init)x=x+10*randn(size(x);%畫出含噪聲圖像subplot(222);image(x);colormap(map);title('含噪聲圖像');axis square%下面進(jìn)行圖像的去噪處理%用小波畫數(shù)coif3對(duì)x進(jìn)行2層小波分解c,s=wavedec2(x,2,'coif3');%提取小波分解中第一層

42、的低頻圖像，即實(shí)現(xiàn)了低通濾波去噪%設(shè)置尺度向量nn=1,2設(shè)置閾值向量pp=10.12,23.28;%對(duì)三個(gè)方向高頻系數(shù)進(jìn)行閾值處理nc=wthcoef2('h',c,s,n,p,'s');nc=wthcoef2('v',c,s,n,p,'s');nc=wthcoef2('d',c,s,n,p,'s');%對(duì)新的小波分解結(jié)構(gòu)nc，s進(jìn)行重構(gòu)xx=waverec2(nc,s,'coif3');%畫出重構(gòu)后圖像的波形subplot(223);image(x);colormap(map);

43、title('去噪后的圖像');axis square輸出結(jié)果如圖圖9 去噪三二維信號(hào)在應(yīng)用中一般表現(xiàn)為圖像信號(hào)，二維信號(hào)在小波域中的降噪方法的基本思想與一維情況一樣，在閾值選擇上，可以使用統(tǒng)一的全局閾值，有可以分作三個(gè)方向，分別是水平方向、豎直方向和對(duì)角方向，這樣就可以把在所有方向的噪聲分離出來，通過作用閾值抑制其成分。4.3 小波分析用于圖像增強(qiáng)4.3.1 圖像增強(qiáng)問題描述小波分析在二維信號(hào)（圖像）處理方面的優(yōu)點(diǎn)主要體現(xiàn)在其時(shí)頻分析特性，前面介紹了一些基于這種特性的一些應(yīng)用的實(shí)例，但對(duì)二維信號(hào)小波系數(shù)的處理方法只介紹了閾值化方法一種，下面我將介紹一下以前在一維信號(hào)中用到的抑

44、制系數(shù)的方法，這種方法在圖像處理領(lǐng)域主要應(yīng)用于圖像增強(qiáng)。在圖像處理領(lǐng)域，圖像增強(qiáng)問題主要通過時(shí)域（沿用信號(hào)處理的說法，空域可能對(duì)圖像更適合）和頻域處理兩種方法來解決。時(shí)域方法通過直接在圖像點(diǎn)上作用算子或掩碼來解決，頻域方法通過修改傅立葉變換系數(shù)來解決。這兩種方法的優(yōu)劣很明顯，時(shí)域方法方便快速但會(huì)丟失很多點(diǎn)之間的相關(guān)信息，頻域方法可以很詳細(xì)地分離出點(diǎn)之間的相關(guān)，但需要做兩次數(shù)量級(jí)為nlogn的傅立葉變換和逆變換的操作，計(jì)算量大得多。小波分析是以上兩種方法的權(quán)衡結(jié)果，建立在如下的認(rèn)識(shí)基礎(chǔ)上，傅立葉分析的在所有點(diǎn)的分辨率都是原始圖像的尺度，但對(duì)于問題本身的要求，我們可能不需要這么大的分辨率，而單純的

45、時(shí)域分析又顯得太粗糙，小波分析的多尺度分析特性為用戶提供了更靈活的處理方法。可以選擇任意的分解層數(shù)，用進(jìn)可能少的計(jì)算量得到我們滿意的結(jié)果。給定一個(gè)wmandril.mat圖像信號(hào)。由于圖像經(jīng)二維小波分解后，圖像的輪廓主要體現(xiàn)在低頻部分，細(xì)節(jié)部分體現(xiàn)在高頻部分，因此可以通過對(duì)低頻分解系數(shù)進(jìn)行增強(qiáng)處理，對(duì)高頻分解系數(shù)進(jìn)行衰減處理，從而達(dá)到圖像增強(qiáng)的效果。具體程序清單如下：load wmandrilsubplot(121);image(x);colormap(map);title('原始圖像');axis square%下面進(jìn)行圖像的增強(qiáng)處理%用小波函數(shù)sym4對(duì)x進(jìn)行2層小波分解c

46、,s=wavedec2(x,2,'sym4');sizec=size(c);%對(duì)分解系數(shù)進(jìn)行處理以突出輪廓部分，弱化細(xì)節(jié)部分for i=1:sizec(2) if(c(i)>350) c(i)=2*c(i); else c(i)=0.5*c(i); endend%下面對(duì)處理后的系數(shù)進(jìn)行重構(gòu)xx=waverec2(c,s,'sym4');%畫出重構(gòu)后的圖像subplot(122);image(xx);colormap(map);title('增強(qiáng)圖像');axis square輸出結(jié)果如圖所示：圖10 小波分析用于圖像增強(qiáng)本文將主要討論圖像增

47、強(qiáng)中的鈍化和銳化兩種方法，鈍化操作主要是提出圖像中的低頻成分，抑制尖銳的快速變化成分，銳化操作正好相反，將圖像中尖銳的部分進(jìn)可能得提取出來，用于檢測(cè)和識(shí)別等領(lǐng)域。下面將以例子說明這兩種方法在matlab中的實(shí)現(xiàn)，并對(duì)于基于傅立葉變換的傳統(tǒng)頻域方法同小波方法做一下比較。4.3.2 圖像鈍化圖像鈍化在時(shí)域中的處理相對(duì)簡(jiǎn)單，只需要對(duì)圖像作用一個(gè)平滑濾波器，使得圖像中的每個(gè)點(diǎn)與其相鄰點(diǎn)做平滑即可，這里不做詳細(xì)介紹，我們來介紹一下基于傅立葉變換的頻域處理方法。下面我們以chess信號(hào)為例，通過兩種方法對(duì)圖像鈍化的結(jié)果做一下比較。load chess分別保存用dct方法和小波方法的變換系數(shù)blur1=x;

48、blur2=x;對(duì)原圖像做二維離散余弦變換ff1=dct2(x);對(duì)變換結(jié)果在頻域做butterworth濾波for i=1:256 for j=1:256 ff1(i,j)=ff1(i,j)/(1+(i*j+j*j)/8192)2); endend重建變換后的圖像blur1=idct2(ff1);對(duì)圖像做2層的二維小波分解c,l=wavedec2(x,2,'db3');csize=size(c);對(duì)低頻系數(shù)進(jìn)行放大處理，并抑制高頻系數(shù)for i=1:csize(2); if(c(i)>300) c(i)=c(i)*2; else c(i)=c(i)/2; endend通

49、過處理后的小波系數(shù)重建圖像blur2=waverec2(c,l,'db3');顯示三幅圖像subplot(221);image(wcodemat(x,192);colormap(gray(256);title('原始圖像','fontsize',18);subplot(223);image(wcodemat(blur1,192);colormap(gray(256);title('采用dct方法鈍化圖像','fontsize',18);subplot(224);image(wcodemat(blur2,192);c

50、olormap(gray(256);title('采用小波方法鈍化圖像','fontsize',18);圖11 圖像鈍化從圖中可以看出，采用dct在頻域做濾波的方法得到鈍化結(jié)果更為平滑，這是因?yàn)槠浞直媛首罡撸〔ǚ椒ǖ玫降慕Y(jié)果在很多地方有不連續(xù)的現(xiàn)象，因?yàn)槲覀儗?duì)系數(shù)做放大或抑制在閾值兩側(cè)有間斷，而且分解層數(shù)很低，沒有完全分離出頻域的信息。而且我們?cè)谧鱿禂?shù)放大或抑制的時(shí)候，采用的標(biāo)準(zhǔn)是根據(jù)系數(shù)絕對(duì)值的大小，沒有完全體現(xiàn)出其位置信息，但是在小波系數(shù)中，我們很容易在處理系數(shù)的過程中加入位置信息。 4.3.3 圖像銳化與圖像鈍化所做的工作相反，圖像銳化的任務(wù)是突出高頻

51、信息，抑制低頻信息，從快速變化的成分中分離出標(biāo)識(shí)系統(tǒng)特性或區(qū)分子系統(tǒng)邊界的成分，以便于進(jìn)一步的識(shí)別、分割等操作。在時(shí)域（空域）中，銳化的方法不外乎是作用掩碼或做差分，同鈍化的道理一樣，無論是掩碼和差分都很難識(shí)別點(diǎn)之間的關(guān)聯(lián)信息，我們下面的例子同樣是在頻域完成的，用傳統(tǒng)的傅立葉分析方法（這里采用的是dct變換）得到的頻域系數(shù)。load chess;分別保存用dct方法和小波方法的變換系數(shù)blur1=x;blur2=x;對(duì)原圖像做二維離散余弦變換ff1=dct2(x);對(duì)變換結(jié)果在頻域做butterworth濾波for i=1:256 for j=1:256 ff1(i,j)=ff1(i,j)/(

52、1+(32768/(i*i+j*j)2); endend重建變換后的圖像blur1=idct2(ff1);對(duì)圖像做2層的二維小波分解c,l=wavedec2(x,2,'db3');csize=size(c);對(duì)低頻系數(shù)進(jìn)行放大處理，并抑制高頻系數(shù)for i=1:csize(2); if(abs(c(i)<300) c(i)=c(i)*2; else c(i)=c(i)/2; endend通過處理后的小波系數(shù)重建圖像blur2=waverec2(c,l,'db3');顯示三幅圖像subplot(221);image(wcodemat(x,192);color

53、map(gray(256);title('原始圖像','fontsize',18);subplot(223);image(wcodemat(blur1,192);colormap(gray(256);title('采用dct方法銳化圖像','fontsize',18);subplot(224);image(wcodemat(blur2,192);colormap(gray(256);title('采用小波方法銳化圖像','fontsize',18);得到的結(jié)果如圖所示圖12 圖像銳化從結(jié)果中可以看出

54、，使用dct方法進(jìn)行高通濾波得到的高頻結(jié)果比較純粹，完全是原圖像上的邊緣信息，而在小波方法得到的結(jié)果中，不只有高頻成分，還有變換非常緩慢的低頻成分，這是因?yàn)閮烧咄瑯釉谛〔ㄏ禂?shù)上體現(xiàn)為絕對(duì)值較低的部分，但這些成分的存在對(duì)我們進(jìn)行進(jìn)一步分析并無多大影響。最后我們來比較一次這兩個(gè)例子的時(shí)間復(fù)雜度，對(duì)dct方法，需要做兩次復(fù)雜度為o（nlogn）的dct變換，中間系數(shù)處理部分復(fù)雜度為o（n），而對(duì)小波變換，無論是分解和重構(gòu)還有系數(shù)處理的復(fù)雜度都是o（n），所以時(shí)間復(fù)雜度的優(yōu)勢(shì)非常明顯。4.4 小波分析用于圖像融合圖像融合是將同一對(duì)象的兩個(gè)或更多的圖像合成在一幅圖像中，以便它比原來的任何一幅圖像更容易得

55、為人們所理解。這一技術(shù)可應(yīng)用于多頻譜圖像理解以及醫(yī)學(xué)圖像處理等領(lǐng)域。在這些場(chǎng)合。同一物體部件的圖像往往是采用不同的成像機(jī)理得到的。下面用二維小波分析將上例中woman.mat和wbarb.mat兩幅圖像融合在一起。程序清單：load woman;x1=x;map1=map;%畫出原始圖像subplot(221);image(x1);colormap(map1);title('woman');axis squareload wbarb;x2=x;map2=map;for i=1:256 for j=1:256 if (x2(i,j)>100) x2(i,j)=1.2*x2(

56、i,j); else x2(i,j)=0.5*x2(i,j); end endendsubplot(222);image(x2);colormap(map2);title('wbarb');axis square%用小波函數(shù)sym4對(duì)x1進(jìn)行2層小波分解c1,s1=wavedec2(x1,2,'sym4');%對(duì)分解系數(shù)進(jìn)行處理以突出輪廓部分，弱化細(xì)節(jié)部分sizec1=size(c1);for i=1:sizec1(2) c1(i)=1.2*c1(i);end%用小波函數(shù)sym4對(duì)x2進(jìn)行2層小波分解c2,s2=wavedec2(x2,2,'sym4');%下面進(jìn)行小波變換域的圖像融合c=c1+c2;%減小圖像亮度c=0.5*c;%對(duì)融合的系數(shù)進(jìn)行重構(gòu)xx=waverec2(c,s1,'sym4');%畫出融合后的圖像subplot(223);image(xx);title('融合圖像');axi