资源描述:
《二维小波算法程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、appcoef2函数%当前延拓模式是补零%装载原始图像loadsinsin;%绘制原始图像subplot(2,2,1);image(X);colormcip(mtip);titleC原始图像');%X包含装载的图像%使用dbl对X进行尺度为2的分解[c,s]=wavcdcc2(X,2,'dbl');sizex=size(X)sizee=size(c)val_s=s%提取尺度为2时的低频系数ca2=appcocf2(c,s,'dbl',2);sizeca2=size(ca2)%绘制尺度为2时的低频图像subplot(2,2
2、,2);image(ca2):colormap(map);titleC尺度为2时的图像');%计算尺度为1时的低频系数cal=appcocf2(c,s,'dbl',1);sizecal=size(cal)%绘制尺度为1时的低频图像subplot(2,2、3);colormap(map);titleC尺度为1时的图像');nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnndetcoef2函数%当前延扌石模式是补零%装载原始图像1oadwoman;%X包含原始图像%绘制原始图像subplot(3,3,1);i
3、mage(X);colormap(map);titleC原始图像');%使用dbl对X进行尺度为2的分解[c,s]=wavedec2(X,2,'dbl');sizex=size(X)sizec=size(c)val_s=s%对丁•各个方向,从小波分解结构[c,s]中提取尺度为2时的高频系数[chd2,cvd2,cdd2]=detcoef2('all',c,s,2);sizecd2=size(chd2)%绘制尺度为2时高频图像subplot(3,3,4);imcige(chd2);titleC尺度为2时高频图像的水平部分'
4、);subplot(3,3,5);titleC尺度为2时高频图像的垂直部分');subplot(3,3,6);image(cdd2);titleC尺度为2时高频图像的对角部分');%对于各个方向,从小波分解结构[c,s]中提取尺度为1时的高频系数[chdl,cvdl,cddl]=detcoef2(,all*,c,s,1);sizecdl=size(chdl)%绘制尺度为1时高频图像subplot(3,3,7);image(chdl);titleC尺度为1时高频图像的水平部分');subplot(3,3,8);image(
5、cvdl);titleC尺度为1时高频图像的垂直部分');subplot(3,3,9);image(cddl);titleC尺度为1时高频图像的对角部分');#################################idwt2函数%当前延拓模式是补零%装碱原始图像loadsinsin;%X包含装载的图像%绘制原始图像subplot(1,2,1);image(X);colormap(map);titleC原始图像');sX=size(X);%使用db4对X进行单尺度分解[cAl,cHl,cVl,cDl]=dvt2
6、(X,'db4');%利用尺度为1时的系数直接重构AO=idwt2(cAl,cHl,cVl,cDl,'db4',sX);%绘制重构图像subplot(1,2,2);image(AO);titleC重构图像');%检查重构效果max(max(abs(X-AO)))iswt2函数%载入原始图像locidsinsin;%X包含装载的图像%绘制原始图像subplot(2,2,1);image(X);colormap(map);titleC原始图像');%使用symd对X在第3层进行SWT分解swc=swt2(X,3,'sym4'
7、);%第二种用法[ca,chd,cvd,cdd]=swt2(X,3,'sym4,);%从小波分解结构swc,重构saO=isvt2(swc,'sym4,);%绘制重构图像1subplot(2,2,3);colormap(map);titleC重构图像1J;%第二种用法aO二iswt2(ca,chd,cvd,cdd,Jsym4,);%绘制璽构图像2subplot(2,2,4):image(aO);colormap(map);titleC重构图像2');%检查璽构的效果err=max(max(abs(X-aO)))errb
8、is=max(max(abs(X~aObis)))nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnwavedec2窗数%当前延扌石模式是补零%装载原始图像1oadwoman;%X包含装载的图像%使用dbl对X进行尺度为2时的小波分解[c,s]=wavedec2(X,2,'dbl');%分解结构siz