MATLAB绘制加噪和去噪图像2

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('除噪图像')

MATLAB图像处理—孙兆林编著—清华大学出版社