zao1

数字图像处理

kk=imread('ziji5.jpg');
X=rgb2gray(kk);
figure;
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('第二次消噪后图像');

kk=imread('ziji5.jpg');
X=rgb2gray(kk);
figure;
subplot (221);
imshow(X,[]);
title('原始图像');
%colormap(map2) ;
axis square ;
init=2055615866
randn('seed',init);
x=X+18*randn(size(X));
nbc=size(map,1);
wname='sym4';lev=3;
[c,s]=wavedec2(x,lev,wname);
alpha=3;m=3.5*prod(s(1,:));
[thr,nkeep]=wdcbm2(c,s,alpha,m);
sorh='s';
%colormap(pink(nbc));
figure(1);
subplot(2,2,1);
 image(wcodemat(X,[]));
title('原始信号');
subplot(2,2,2);
imshow(wcodemat(x,[]));
title('含噪信号');
subplot(2,2,3);
 imshow(wcodemat(xd,[]));
 title('消噪后的信号');