tuxiangquzao2.m

来自「数字图像处理」· M 代码 · 共 25 行

kk=imread('ziji5.jpg');
X=rgb2gray(kk);
X=double(X);
figure();
axis square ;
init=2055615866
randn('seed',init);
x=X+18*randn(size(X));

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';

subplot(2,2,1);
 imshow(X,[]);
title('原始信号');
subplot(2,2,2);
imshow(x,[]);
title('含噪信号');
subplot(2,2,3);
 imshow(x,[]);
 title('消噪后的信号');