figure16.m

来自「是一本教科书的各个章节的源程序代码,适合初学数字图像的人」· M 代码 · 共 21 行

clear
load facets;
subplot(221); image(X); colormap(map);
title('原始图像');
axis square
init=1992625743; randn('seed', init)
x=X+10*randn(size(X));
subplot(222); image(x); colormap(map);
title('含噪声图像');
axis square
[c, s]=wavedec2(x, 2, 'coif3');
n=[1, 2]
p=[11.20, 22.81];
nc=wthcoef2('h', c, s, n, p, 's');
nc=wthcoef2('v', c, s, n, p, 's');
nc=wthcoef2('d', c, s, n, p, 's');
xx=waverec2(nc, s, 'coif3');
subplot(223); image(xx); colormap(map);
title('去噪后图像');
axis square;