yyqz.asv

来自「数字图像处理」· ASV 代码 · 共 31 行

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';
xd=wdencmp('lvd',c,s,wname,lev,thr,sorh);

%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('消噪后的信号');