小波變換的特殊例子

1、小波變換與圖像處理作者:雪漫一：引言本文從二維小波理論出發(fā)，對其在圖像處理的應用上進行了一些分析和處理，力圖反映出小波分析在圖像處理方面有著其獨特的特點。本文就以下幾點進行闡述：小波基本概念圖像壓縮圖像消噪圖象增強圖象平滑處理二：小波基本概念小波定義：設,其傅立葉變換為,當滿足允許條件，即完全重構條件或恒等分辨條件.時，我們稱為一個基本小波或母小波，將母函數經伸縮和平移后，得。我們稱其為一個小波序列。其中a為伸縮因子，b為平移因子。小波變換是一種信號的時間尺度分析方法，他具有多分辨率分析的特點，而且在時頻兩域都具有表征信號局部特征的能力，是一種窗口大小固定不變但其形狀可變，時間窗和頻率窗都可變

2、的時頻局部化分析方法。即再低頻部分具有較高的頻率分辨率和時間分辨率，在高頻部分具有較高的時間分辨率和較低的頻率分辨率，很適合探測正常信號中夾帶的瞬態(tài)反?，F象并展示其成分，所以被譽為分析信號的顯微鏡。波分析是把信號分解成低頻al和高頻dl兩部分，在分解中，低頻al中失去的信息由高頻dl捕獲。在下一層的分解中，又將al分解成低頻a2和高頻d2兩部分，低頻a2中失去的信息由高頻d2捕獲，如此類推下去，可以進行更深層次的分解。二維小波函數是通過一維小波函數經過張量積變換得到的，二維小波函數分解是把尺度j的低頻部分分解成四部分：尺度j+1的低頻部分和三個方向（水平、垂直、斜線）的高頻部分。三：圖像壓縮對

3、于圖像來說，如果需要進行快速或實時傳輸以及大量存儲，就需要對圖像數據進行壓縮。在同樣的通信容量下，如果圖像數據壓縮后在傳輸，就可以傳輸更多的圖像信息。例如，用普通的電話線傳輸圖像信息。圖像壓縮研究的就是尋找高壓縮比的方法且壓縮后的圖像要有合適的信噪比,在壓縮傳輸后還要恢復原信號，斌且在壓縮、傳輸、恢復的過程中，還要求圖像的失真度小。這就是圖像壓縮的研究問題。圖像數據往往存在各種信息的冗余、如空間冗余、信息熵冗余、視覺冗余和結構冗余等等。所謂壓縮就是去掉各種冗余，保留對我們有用的信息。圖像壓縮的過程常稱為編碼。相對的，圖像的恢復當然就是解碼了。圖像壓縮的方法通?？煞譃橛惺д婢幋a和無失真編碼兩大類

4、：無失真編碼方法如改進的霍夫曼編碼。有失真編碼方法的還原圖像較之原始圖像存在著一些誤差，但視覺效果是可以接受的。常見的方法有預測編碼、變換編碼、量化編碼、信息熵編碼、分頻帶編碼和結構編碼等等。而將小波分析引入圖像壓縮的范疇也是一個重要的手段，并且有著它自己的特點。它的特點在于壓縮比高、壓縮速度快，壓縮后能保持信號與圖像的特征基本不變，且在傳遞過程中可以抗干擾等等。下面我們就舉一個粒子來說明怎樣用小波分析進行圖像壓縮。例如現在有一個二維圖像（文件名為），我們利用二維小波分析來進行圖像壓縮。由原理可知，一個圖像作小波分解后，可得到一系列不同分辨率的子圖像，不同分辨率的子圖像對應的頻率是不相同的。高

5、分辨率（高頻）子圖像上大部分點的數值都接近于0，越是高就越是明顯。而對于一個圖像來說，表現一個圖像的最主要的部分是低頻部分，所以最簡單的壓縮方法是利用小波分解去掉圖像的高頻部分而只保留低頻部分。程序大致如下：clear%裝入圖像loadwbarb;%顯示圖像symsX;subplot(221);image(coast);colormap(map)title('原始圖像');axissquaredisp('壓縮前圖像X的大小');whos('coast')%對圖像用小波進行層小波分解c,s=wavedec2(X,2,'bior3.7'

6、;);%提取小波分解結構中的一層的低頻系數和高頻系數cal=appcoef2(c,s,'bior3.7',1);%水平方向ch1=detcoef2('h',c,s,1);%垂直方向cv1=detcoef2('v',c,s,1);%斜線方向cd1=detcoef2('d',c,s,1);%各頻率成份重構a1=wrcoef2('a',c,s,'bior3.7',1);h1=wrcoef2('h',c,s,'bior3.7',1);v1=wrcoef2('v'

7、;,c,s,'bior3.7',1);d1=wrcoef2('d',c,s,'bior3.7',1);c1=a1,h1;v1,d1;%顯示分頻信息subplot(222);image(c1);axissquare;title('分解后低頻和高頻信息');%進行圖像壓縮%保留小波分解第一層低頻信息%首先對第一層信息進行量化編碼ca1=appcoef(c,s,'bior3.7',1);ca1=wcodemat(ca1,440,'mat',0);%改變圖像高度并顯示ca1=0.5*ca1;subplot(

8、223);image(ca1);colormap(map);axissquare;title('第一次壓縮圖像');disp('第一次壓縮圖像的大小為：');whos('ca1')%保留小波分解第二層低頻信息進行壓縮ca2=appcoef2(c,s,'bior3.7',2);%首先對第二層信息進行量化編碼ca2=wcodemat(ca2,440,'mat',0);%改變圖像高度并顯示ca2=0.25*ca2;subplot(224);image(ca2);colormap(map);axissquare;title

9、('第二次壓縮圖像');disp('第二次壓縮圖像的大小為：');whos('ca2')輸出結果如圖：NameSizeBytesclass壓縮前圖像X256x256524288Doublearray第一次壓縮圖像Ca1135x135145800Doublearray第二次壓縮圖像Ca275x7545000Doublearray在這里可以看出，第一次壓縮我們是提取原始圖像中小波分解第一層的低頻信息，此時壓縮效果較好，壓縮比較小(約為1/3大小)。第二次壓縮實提取第一層分解低頻部分的低頻部分(即第二層的低頻部分),其壓縮比較大(約為1/12)，壓縮效

10、果在視覺上也基本過得去。上面的保留原始圖像中低頻信息的壓縮辦法只是一種最簡單的壓縮辦法。它不需經過其他處理即可獲得較好的壓縮效果。當然，對于上面的例子我們還可以只提取小波分解的第三、第四層的低頻信息。從理論上說，我們可以獲得任意壓縮比的壓縮圖像。只不過在對壓縮比和圖像質量都有較高要求時，它就不如其他編碼方法了。下面我們在舉一個例子，這一次用中函數來對上圖進行壓縮。Clear;%裝入圖形信號loadwbarb;%顯示圖像subplot(221);image(X);colormap(map);title('原始圖像');disp('壓縮前圖像的大小');whos(&

11、#39;X');axissquare;%對圖像進行壓縮%對圖像用db3小波進行二層小波分解c,s=wavedec2(X,5,'db3');thr,sorh,keepapp=ddencmp('cmp','wv',X);Xcomp,cxc,lxc,perf0,perfl2=wdencmp('gbl',c,s,'db3',5,thr,sorh,keepapp);%將壓縮后的圖像于原始圖像相比較subplot(222);image(Xcomp);colormap(map);title('壓縮后的圖像'

12、;);disp('壓縮后圖像的大小');whos('Xcomp')%顯示有關參數disp('小波分解系數中值為0的系數個數百分比');disp(perf0);disp('壓縮后剩余能量百分比');disp(perfl2);輸出結果如下：小波分解系數中值為0的系數個數百分比：49.8088壓縮后剩余能量百分比：99.9754總之，是事無絕對。一種壓縮圖像的方法不可能盡善盡美。要想很好的進行圖像的壓縮，就需要綜合的利用多種其他技術，特別是數據編碼和解碼算法。四：圖像消噪圖像消噪方法的一般說明對二維圖像信號的消噪方法同樣適用于一維信號,

13、尤其是對于幾何圖像更適合。二維模型可以表述為其中，e是標準偏差不變得高斯白噪聲。二維信號的消噪步驟與一維信號的消噪步驟完全相同，也有三步，只是用二維小波分析工具代替了一維小波分析工具。如果用固定的閥值形式，測選擇的閥值用mH代替了一維信號中的n。著三步是：(1) 二維信號的小波分解。選擇一個小波和小波分解的層次N,然后計算信號s到第N層的分解。(2) 對高頻系數進行閥值量化。對于從一到N的每一層，選擇一個閥值，斌對著一層的高頻系數進行軟閥值化處理。(3) 二維小波的重構。根據小波分解的第N層的低頻系數和經過修改的從第1層到第N層的各層高頻系數，來計算二維信號的小波重構。在這三個步驟中，重點內容

14、就是如何選取閥值和如何進行閥值的量化。請注意，了一維信號自動消噪的情況，對于其他的情況，一維信號的消噪和壓縮用的是wdencmp,這對于二維信號也是一樣的。編程給定一個有較大白噪聲的圖象，利用二維小波分析進行信號消噪處理。分析：由于圖象所含的噪聲主要是白噪聲，且集中于高部分，故用第通實現消去噪聲。程序如下。loadtire;subplot(221);image(X);colormap(map);title('原圖');axissquare;畫出原圖象init=2055615866;randn('seed',init)x=X+38*randn(size(X);su

15、bplot(222);image(x);colormap(map);title('含噪聲圖象');axissquare;畫出含噪聲圖象c,s=wavedec2(x,2,'sym4');a1=wrcoef2('a',c,s,'sym4',1);第一次低通濾波消噪subplot(223);image(a1);title('第一次消噪后圖象');axissquare;畫出第一次低通濾波消噪后圖象a2=wrcoef2('a',c,s,'sym4',2);第二次低通濾波消噪subplot(22

16、4);image(a2);title('第二次消噪后圖象');axissquare;畫出第二次低通濾波消噪后圖象分析：第一次消噪濾去了大部分高頻噪聲，但與原圖比較，依然有不少高頻噪聲，第二次消噪在第一次消噪基礎上，再次濾去高頻噪聲，消噪效果較好，但圖像質量比原圖稍差。五：圖象增強說明小波變換將一幅圖象分解為大小、位置和方向都不同的分量。在做逆變換之前可以改變小波變換域中某些系數的大小，這樣就能夠洋選擇的放大所感興趣的分量而減小不需要的分量編程：給定一個圖象信號，用二維小波分析對圖象進行增強處理。分析由于圖象經二維小波分解后，圖象的輪廓主要體現在低頻部分，而細節(jié)部分則體現子高頻部

17、分，因此，可以通過對低頻分解系數進行增強處理，對高頻分解系數進行衰減處理，即可以達到圖象增強的作用。具體處理過程如下程序：loadwoman;subplot(121);image(X);colormap(map);title(原始圖象)；axissquare;畫出原圖象*c,s+=wavedec2(X,2,'sym4');%進行二層小波分解sizec=size(c);處理分解系數，突出輪廓，弱化細節(jié)forI=1:sizec(2)if(c(I)>350)c(I)=2*c(I);elsec(I)=0.5*c(I);end分解系數重構畫出增強圖像endxx=waverec2(c

18、,s,'sym4');subplot(122);image(xx);title(增強圖象)axissquare;結果分析：打到了圖像增強的效果，試圖像對比更加明顯，但由于細節(jié)上的弱化，卻使得圖像產生模糊的感覺。就那婦女托著下腮的手來說，增強后的圖像幾乎就不能辨認。手指更是消失了。六：圖象融合說明圖象融合是將同一對象的兩個或更多的圖象合成在一幅圖象中，以便他比原來的任何一幅更能容易的為人們所理解。真一技術可應用于多頻譜圖象理解以及醫(yī)學圖象處理等領域，再這些場合，同一物體部件的圖象往往是采用不同的成象機理得到的。編程：用二維小波分析將兩幅圖象融合在一起。處理過程如下：loadwom

19、an;裝入原圖像X1=X;map1=map;subplot(221);image(X1);colormap(map1);title(woman');axissquare%畫出woman圖像loadwbarb;裝入原圖像X2=X;map2=map;forI=1:256forj=1:256if(X2(I,j)>100)X2(I,j)=1.2*X2(I,j);elseX2(I,j)=0.5*X2(I,j);endendendsubplot(222);image(X2);colormap(map2);title(wbarb');axissquare%畫出wbarb圖像*cl,sl

20、+=wavedec2(X1,2,'sym4');sizec1=size(c1);forI=1:sizec1(2)c1(I)=1.2*c1(I);endc2,s2+=wavedec2(X2,2,'sym4');c=c1+c2;c=0.5*c;xx=waverec2(c,s,'sym4');subplot(223);image(xx);title(融合圖象)；axissquare%畫出融合后的圖結果分析：一幅圖像和他某一部分放大后的圖像融合，融合后的圖像給人一種朦朦朧朧夢幻般的感覺，對較深的背景部分則做了淡化處理。七：圖象平滑處理說明圖像平滑的主要目

21、的是為了減少噪聲，一般情況下，在空間域內可以用于平均來減少噪聲。在頻率域，因為噪聲浦多在高頻段，因此可以曹用各種形式的低通濾波的辦法來減少噪聲。編程給定一個含噪聲的圖象，用二維小波分析和圖象的中值濾波進行圖象的平滑。分析這是一個圖象平滑處理問題。首先，對圖象在頻域內進行增強，然后在空域內加入較大的白噪聲。通過對含噪圖象進行平滑處理，即可以使含噪圖象具有較好的平滑效果。具體處理過程如下：loadwoman;裝入原圖X1=X;map1=map;subplot(221);image(X1);colormap(map1);title('woman');axissquare畫出原圖二層分解小波信號頻域里增強圖像系數重構加入噪聲c,s+=wavedec2(X,2,'sym4');sizec=size(c);forI=1:sizec(2)if(c(I)>35