program_16_15.m

小波分析理论与matlab2007配套程序

%装载原始图像
load sinsin;
subplot(2,2,1);image(X);colormap(map);title('原始图像');
%对图像加噪声并显示出含噪图像
init=2788605800;
randn('seed',init);
X=X+10*randn(size(X));
subplot(2,2,2);image(X);colormap(map);title('含噪图像');
%应用中值滤波进行图像平滑处理
[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);
            end
        end
        Xtemp=Xtemp/9;
        X1(i,j)=Xtemp;
    end
end
subplot(2,2,3);image(X1);colormap(map);
title('平滑滤波后图像');
%应用小波变换进行消噪
wname = 'coif2'; lev = 3;
[c,s] = wavedec2(X,lev,wname);
det1 = detcoef2('compact',c,s,1);
sigma = median(abs(det1))/0.6745;
alpha = 1.2;
thr = wbmpen(c,s,sigma,alpha);
keepapp = 1;
xd = wdencmp('gbl',c,s,wname,lev,thr,'s',keepapp);
subplot(2,2,4); image(xd);colormap(map);title(' 小波消噪后图像');