image_denoise.m

来自「原创：在matlab里仿真有关图像的一些数字处理技术」· M 代码 · 共 34 行

clear
%装入图像
load tire;
subplot(221);
%画出原图象
image(X);
colormap(map)
title('原始图像');
axis square
init=2055615866;
randn('seed',init)
x=X+38*randn(size(X));
subplot(222);
%画出含噪声图象
image(x);
colormap(map)
title('加噪声图像');
axis square
[c,s]=wavedec2(x,2,'sym4');
a1=wrcoef2('a',c,s,'sym4',1);
%第一次低通滤波消噪
subplot(223);
image(a1);
title('第一次消噪后图象');
axis square
%画出第一次低通滤波消噪后图象
%第二次低通滤波消噪
a2=wrcoef2('a',c,s,'sym4',2);
subplot(224);
%画出第二次低通滤波消噪后图象
image(a2);
title('第二次消噪后图象');
axis square