二维小波变换MATLAB实现

1、二维小波变换M ATLAB实现dwt2函数功能：二维离散小波变换格式：cA,cH,cV5cD=dwt2(X3'wname')cA,cH,cV,cD=dwt2(X丄o_D,Hi_D)说明:cAJcHJcVJcD=dwt2(X,'wname')使用指定的小波基函数'wname，对二维信号X进行二维离散 小波变幻；cA, cH,cV,cD分别为近似分量、水平 细节分量、垂直细节分量和对角细节分量；cA,cHJcVJcD=dwt2(X ,Lo_D,Hi_D)使用指定的分 解低通和高通滤波器Lo_DHiJ5分解信号X。load woman;nbcol = size

2、(mapj);%返回矩阵的行数和列数 cA1,cH1,cV1,cD1 = dwt2(X；db1');cod_x=wcodemat(X,nbcol);%返回矩阵X的编码矩阵， nbc6l为编码的最矢值cod_cA1=wcodemat(cA13nbcol);cod_cH1 =wcodemat(cH15nbcol);cod_cV1 =wcodemat(cV15nbcol);cod_cD1 =wcodemat(cD15nbcol); dec2d=cod_cA13cod_cH1 ;cod_cV13cod_cD1; subplot(1 525 1 )3imshow(cod_x,);title”量化

3、后的图像)；subplot(132,2),imshow(dec2d,);titled二维离散小波分解后的图像)；量化后的图像二维离散小浪分解后的图像'-a4 idwt2函数功能：二维离散小波反变换格式：X=idwt2(cA,cH5cV5cD5,wname,) X=idwt2(cA5cH5cV5cD5Lo_R5Hi_R) X=idwt2(cA,cH,cV,cD；wnameS) X=idwt2(cA5cH5cV5cD5Lo_R5Hi_R5S)说明：X=idwt2(cA5cH5cV5cD；wname,)由信号小波分解的 迸侦彳言号cA和细节信号cH、cH、eV、cD经小波反交強董 构原信号X

4、； X=idwt2(cA5cH5cV5cD5Lo_R5Hi R)使用指定 的重构低通和高通滤波器Lo_R和Hi_R重构原信号X； X=idwt2(cA5cH5cV5cD5,wnameS)和X=idwt2(cA,cH3cV,cD ,Lo_R,Hi_R,S)返回中心附近的 S 个 数据点。 load woman; sX=size(X); cA1 5cH1 3cV1 JcD1=dwt2(X/db4,); A0=idwt2(cA1,cH15cV1 ,cD1 ；db4',sX); subplot(1,2,1)5imshow(X5) titled原始图像)； subplot(15252)5imsh

5、ow(A05) titled由二维小波分解重构后的图像)；厲始團像由二维小滾分解重构后的图像 waved ec2 函数功能：二维信号的多层小波分解格式：C,S=wavedec2(XN'wname) C,S=wavedec2(X5N 5Lo_D5Hi_D)说明：C,S=wavedec2(X,N,*wname1)用小 菠基函数'wnameW二维信号X进行N定分 解；C,S=wavedec2(X, N 丄o_D,Hi_D)用指 定的分解低通和高通滤波器Lo_DWi_D分 解信号X。 waverec2 函数说明：二维信号的多层小波重构格式：X=waverec2(C,S5'wna

6、me') X=waverec2(C,S ,Lo_R,Hi_R)说明：X=waverec2(C,S5'wname')由多层二 维小波分解的结果C、S重构原始信号X, 'wname'为使用的小波基函数； X=waverec2(C,S ,Lo_R,Hi_R)使用重构低 通和高通滤波器Lo_R和Hi_R重构原信号。 load woman; c5s=wavedec2(X52；sym4, ); a0=waverec2(c,s,'sym4'); subplot(1,2,1 ),imshow(X,)title('原始图像)； subplot(1,

7、2,2)Jimshow(aOJ)title('重构后的图像')；原始團像重构后的團像appcoef2提取二维信号小波分解的近似分量 load woman; c5s=wavedec2(X53；db1');% 用 db1 小波进行三层分解 cal =appcoef2(c5s；db1 j;% 提取一层近似分量 ca2=appcoef2(c5s；db22);% 提取二层近似分量 imshow(X,);title('原始图橡*); figure5imshow(ca15)title ('第一层近似分量')； figure5imshow(ca25) title

8、C第二层近似分量)；detcoef2提取二维信号小波分解 的细节分量load woman; c,s=wavedec2(X,3,db1'); chd2=detcoef2(,hc,s,1); cvd2=detcoef2(V,c,s,1); cdd2=detcoef2(,d,Jc,s, 1); imshow(X,)title(源始图像f);figure熄始慘像rra *7 I :沁 VA,vZ!第一层分解的水平细节第一层分解的垂直细节第一层分解的对角细节subplot(1,3,1 ),imshow(chd2,) titled第一层分解的水平细节)； subplot(1,3,2),imshow

9、(cvd25) titled第一层分解的垂直细节)； subplot(1,3,3),imshow(cdd2,) title第一层分解的对角细节')； wrcoef2由多层小波分解重构某一层的分解信号 clear all load wbarb;image(X);colormap(map);colorbar; I 二 in d2gray(X3map);figure;l=imadjust(l5stretchlim(l),0,1);imshow(l); %j=imadjust(i,low_in;highnow_oulhigh_out)将 i $的亮度眉映射到 j 中 的新值 wname=sym

10、2' c,s=wavedec2(l,2,w name); cA1 =appcoef2(c,s,w name,1); ch1 ,cv1 ,cd1=detcoef2(,aH,c,s,1); cA2=appcoef2(c,s,w name,2); ch2,cv2,cd2=detcoef2('air,c,s,2); a1 =wrcoef2(,ac,s,wname31); h1=wrcoef2('hc5s3w name,1); v1 =wrcoef2( V,c,s3wname, 1); d1=wrcoef2(,d,c)s3wname,1); a2=wrcoef2(,a,c5s3w

11、name,2); h2=wrcoef2(,hc)s3wname,2); v2=wrcoef2('v,c)sJwname,2); d2=wrcoef2(,dc)sJwname,2); figure; subplot(2,2,1 );a1 =abs(a1 );a1 =imadjust(a13stretchlim(a1 ),0,1);imshow(a1); titlefapproximation a1') subplot(2,2,2);h1 =abs(h1 );h1 =imadjust(h1 ,stretchlim(h1 ),0,1 );imshow(h1); title('

12、horizontal detail hi') subplot(2,2,3);v1 =abs(v1 );v1 =imadjust(v1 ,stretchlim(v1 ),0,1);imshow(v1); title('vertical detail v1') subplot(2,2,4);d1 =abs(d1 );d1 =imadjust(d1 ,stretchlim(d1 ),0,1);imshow(d1); title('diagonal detail d1') figure; subplot(2,2,1);a2=abs(a2);a2=imadjust

13、(a2,stretchlim(a2),051);imshow(a2); title('approximation a2') subplot(252,2);h2=abs(h2);h2=imadjust(h23stretchlim(h2),031);imshow(h2); titlefhorizontal detail h2) subplot(232,3);v2=abs(v2);v2=imadjust(v25stretchlim(v2),0,1);imshow(v2); titlefvertical detail v2') subplot(252,4);d2=abs(d2)

14、;d2=imadjust(d2Jstretchlim(d2),031);imshow(d2);title('diagonal detail d2')approximation a1horizontal detail hivertical detail y1diagonal detail d1approximation a2horizontal detail h2vertical detail v2diagonal detail d2upcoef2由多层小波分解重构近似分量或细节分量 load woman; c,s=wavedec2(X,2；db4,); siz=s(size(s,1),:); cal=appcoef2(c,s,db4',1);重构近似信号重构水平细节信号 a1=upcoef2(,acalJ,db4,1,siz); chd1 =detcoef2(,h,c,s,1); hd1 =upcoef2(,hchd1,siz); cvd1 =detcoef2(V,c,s,1); vd1 =upcoef2('h,cvd1,siz); cdd1 =detcoef2('dc,s,1); dd1 =upcoef2(,dc