例程13-24.m

来自「小波分析Matlab实现」· M 代码 · 共 20 行

% 装载原始图像
load sinsin
nbc = size(map,1);
% 画出原始图像
subplot(221), image(wcodemat(X,nbc))
title('原始图像');
% 产生噪声图像
init = 2055615866; randn('seed',init); 
x = X/18 + randn(size(X));
% 画出噪声图像
subplot(222), image(wcodemat(x,nbc))
title('噪声图像');
% 使用wpdencmp进行图像降噪
% 寻找默认值（参见ddencmp函数）
[thr,sorh,keepapp,crit] = ddencmp('den','wp',x)
% 使用全局阈值和SURE最优基进行降噪
xd = wpdencmp(x,sorh,3,'sym4',crit,thr,keepapp);
subplot(223),image(wcodemat(xd,nbc))
title('降噪后的图像');