zao2

数字图像处理

kk=imread('ziji5.jpg');
X=rgb2gray(kk);
X=double(X);
imshow(X,[]);
%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);
 imshow(X,[]);
title('原始信号');
subplot(2,2,2);
imshow(x,[]);
title('含噪信号');
subplot(2,2,3);
 imshow(x,[]);
 title('消噪后的信号');