r21.m

该程序基于小波变换

%装载并显示原始图象
A=imread('m.jpg');
%=============================================
%转换为索引图
A=double(A);
Xrgb=0.2990*A(:,:,1)+0.5870*A(:,:,2)+0.1140*A(:,:,3);
Nbcolors=255;
X=wcodemat(Xrgb,Nbcolors);
map=pink(Nbcolors);
figure(1)
subplot(1,3,1);
image(X);
colormap(map);
title('原始图象');
axis square;
%=============================================
%加噪并显示图象
init=2055615866;
rand('seed',init);
x=X+10*randn(size(X));
subplot(1,3,2);
image(x);
colormap(map);
title('含噪图象');
axis square;
%===============================================
%对图象消噪
wnane='sym4',lev=3;
[c,s]=wavedec2(x,lev,wnane);
[thr3,nkeep3]=wdcbm2(c,s,1.5,2.7*prod(s(1,:)));
[xd,cxd,sxd,perf0,perf12]=wdencmp('lvd',c,s,'sym4',3,thr3,'s');
%==============================================
%显示消噪后的结果
subplot(1,3,3);
image(xd);
colormap(map);
title('消噪图象');
axis square;
figure(2)
image(X);
colormap(map);
title('原始图象');
axis square;
figure(3)
image(x);
colormap(map);
title('含噪图象');
axis square;
figure(4)
image(xd);
colormap(map);
title('消噪图象');
axis square;