10-12.m

来自「MATLAB数字图象处理的主要源码.包括图象去噪变换增强恢复识别压缩的各种算法」· M 代码 · 共 17 行

load sinsin
% 打开原始图像
init = 2055615866; randn('seed',init); 
x = X/18 + randn(size(X));
% 产生噪声图像
% 用wpdencmp函数除噪
[thr,sorh,keepapp,crit] = ddencmp('den','wp',x); 
% 计算缺省值
xd = wpdencmp(x,sorh,3,'sym4',crit,thr,keepapp);
% 用SURE最佳基全局阈值除噪图像
subplot(121)
imshow(x*18,map)
title('加噪图像')
subplot(122)
imshow(xd*18,map)
title('除噪图像')