matlab小波变换图像处理实现程序.doc

《matlab小波变换图像处理实现程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、（到学校的维普网上查找《基于形态学的小波变换图像去噪方法研究》论文）matlab小波变换图像处理实现程序%这个是2D-DWT的函数，是haar小波%c是图像像素矩阵steps是变换的阶数functiondwtc=dwt_haar(c,steps)%DWTC=CWT_HARR(C)-DiscreteWaveletTransformusingHaarfilter%%MDPlumbleyNov2003N=length(c)-1;%Maxindexforfilter:0..N%Ifnostepstodo,orthesequenceisasingles

2、ample,theDWTisitselfif(0==N

3、steps==0)dwtc=c;returnend%CheckthatN+1isdivisibleby2if(mod(N+1,2)~=0)disp(['Notdivisible2:'num2str(N+1)]);returnend%SettheHaaranalysisfilterh0=[1/21/2];%HaarLow-passfilterh1=[-1/21/2];%HaarHigh-passfilter%Filterthesignallowpass_c=conv(h0,c);hipas

4、s_c=conv(h1,c);%Subsamplebyfactorof2andscalec1=sqrt(2)*lowpass_c(2:2:end);d1=sqrt(2)*hipass_c(2:2:end);%Recursivelycalldwt_haaronthelow-passpart,with1fewerstepsdwtc1=dwt_haar(c1,steps-1);%ConstructtheDWTfromc1andd1dwtc=[dwtc1d1];%Donereturn--------------------------分割线-----

5、---------------------调用这个函数的例子下面的东西放在另一个文档里读入一个图像‘lena’应该是个最基础的图像了~之后分别作0阶和1阶2D-DWT的变换改变阶数可以做更高阶的clearallim=imreadreal('lena.bmp');%readimagedata%Plotfiguredwt2_step0=dwt2_haar(im,0);%2DDWTstep=0imagesc(dwt2_step0);colormapgray;axisimage;figuredwt2_step1=dwt2_haar(im,1);%2D

6、DWTstep=1imagesc(dwt2_step1);colormapgray;axisimage;---------------------分割线---------------------结果如下1阶的结果这是我的一个实验希望有所帮助