tuxiangquzao.m

数字图像处理

kk=imread('ziji5.jpg');
X=rgb2gray(kk);
figure;

subplot (221);
imshow(X,[]);
title('原始图像');
%colormap(map2) ;
axis square ;
init=2055615866 ;
rand('seed',init) ;
X=double(X) ;
x=X+40*randn(size(X));
subplot(222);
imshow(x,[]);
title('含噪声图像');
axis square;
[c,s]=wavedec2(x,2,'sym4');
a1=wrcoef2('a',c,s,'sym4',1);
subplot(223);
imshow(a1,[]) ;
title('第一次消噪后图像');
axis square;
a2=wrcoef2('a',c,s,'sym4',2);
subplot(224);
imshow(a2,[]);
title('第二次消噪后图像');