二維小波變換與圖像處理的MATLAB仿真課件

1、二維小波變換與圖像處理二維信號(hào)也稱圖像信號(hào)。為了避免引進(jìn)第二維之后問題的復(fù)雜性，我們可以把圖像信號(hào)分解成沿行和列的一維問題來處理。二維小波變換圖像的自身的特點(diǎn)決定了我們?cè)趯⑿〔ㄗ儞Q應(yīng)用到圖像處理中時(shí)，必須把小波變換從一維推廣到二維。二維連續(xù)小波定義令 表示一個(gè)二維信號(hào)，x1、x2分別是其橫坐標(biāo)和縱坐標(biāo)。 表示二維基本小波，二維連續(xù)小波定義：二維連續(xù)小波定義則二維連續(xù)小波變換為：式中因子 是為了保證小波伸縮前后其能量不變而引入的歸一因子。圖像的小波變換實(shí)現(xiàn)圖像小波分解的正變換可以依據(jù)二維小波變換按如下方式擴(kuò)展，在變換的每一層次，圖像都被分解為4個(gè)四分之一大小的圖像。Mallat二維多分辨率分解與

2、重構(gòu)合成重構(gòu)工具idwt2單尺度逆二維離散小波變換waverec2多尺度二維小波重構(gòu)wrcoef2對(duì)二維小波系數(shù)進(jìn)行單支重構(gòu)upcoef2對(duì)二維小波分解的直接重構(gòu)分解結(jié)構(gòu)工具detcoef2提取二維小波分解高頻系數(shù)appcoef2提取二維小波分解低頻系數(shù)upwlev2二維小波分解的單尺度重構(gòu)二維離散平穩(wěn)小波變換swt2二維離散平穩(wěn)小波變換iswt2二維離散平穩(wěn)小波逆變換圖像之所以能夠壓縮，數(shù)學(xué)機(jī)理主要有以下兩點(diǎn)（1）原始圖像數(shù)據(jù)往往存在各種信息的冗余，數(shù)據(jù)之間存在相關(guān)性，臨近像素的灰度，高度往往是相關(guān)的。（2）在多媒體應(yīng)用領(lǐng)域中，人眼作為圖像信息的接受端，其視覺對(duì)邊緣急劇變化不敏感，對(duì)圖像亮度

3、信息敏感，對(duì)顏色分辨率弱，因此在高壓縮比的情況下，解壓縮后的圖像信號(hào)仍有滿意的主觀質(zhì)量。所謂圖像壓縮就是去掉冗余，保留主要信息。小波變化通過多分辨分析過程，將一副圖像分成近似和細(xì)節(jié)部分，細(xì)節(jié)對(duì)應(yīng)的是小尺度的瞬間，它在本尺度內(nèi)很穩(wěn)定。因此將細(xì)節(jié)存儲(chǔ)起來，對(duì)近似部分在下一個(gè)尺度進(jìn)行分解，重復(fù)該過程即可。近似與細(xì)節(jié)在正交鏡像濾波器算法中分別對(duì)應(yīng)于高通和低通濾波器。ca1=wcodemat(cA1,440,mat,0);ca1=0.5*ca1;subplot(2,2,3);image(ca1);colormap(pink);title(第一次壓縮圖像);axis square;disp(第一次壓縮圖像

4、的大?。?;whos(ca1);cA2=appcoef2(c,s,bior3.7,2);ca2=wcodemat(cA2,440,mat,0);ca2=0.5*ca2;subplot(2,2,4);image(ca2);colormap(pink);title(第二次壓縮圖像);disp(第二次壓縮圖像的大小：);whos(ca2)壓縮前圖像的大?。?Name Size Bytes Class Attributes X 256x256 65536 uint8 壓縮前圖像的大?。?Name Size Bytes Class Attributes X 256x256 65536 uint8 第一次

5、壓縮圖像的大小 Name Size Bytes Class Attributes ca1 135x135 145800 double 第二次壓縮圖像的大小 Name Size Bytes Class Attributes ca2 75x75 45000 double 小波消噪的基本原理常用的圖像消噪方法是小波閾值消噪方法，就是對(duì)小波分解后的各層系數(shù)模大于和小于某閾值的系數(shù)分別進(jìn)行處理，然后利用處理后的小波系數(shù)重構(gòu)出消噪后的圖像。小波消噪X=imread(C:UsershmPictures1234.jpg);X=rgb2gray(X);figuresubplot(2,2,1);imshow(X)

6、;colormap(pink);title(原始圖像);B=imnoise(X,gaussian,0,0.01);subplot(2,2,2);image(B);colormap(pink);title(含噪圖像);c,l=wavedec2(B,2,sym5);a1=wrcoef2(a,c,l,sym5,1);a2=wrcoef2(a,c,l,sym5,2);subplot(2,2,3);image(a1);colormap(pink);title(第1層重構(gòu)圖像);subplot(2,2,4);image(a2);colormap(pink);title(第2層重構(gòu)圖像);小波平滑p,q=size(X);for i=2:p-1 for j=2:q-1 Xtemp=0; for m=1:3 for n=1:3 Xtemp=Xtemp+X(i+m-2,j+n-2)