f11_2.m

《Matlab时频分析技术及其应用》电子书。

%装载原始图像信号
load wgatlin;
subplot(2,2,1);
image(X);
colormap(map);
title('原始图像');

%生成含噪图像
init=2055615866;
rand('seed',init);
XX=X+12*randn(size(X));
subplot(2,2,2);
image(XX);
colormap(map);
title('含噪图像');

%下面对图像消噪
%首先用sym5小波函数对图像进行2层分解
[c,l]=wavedec2(XX,2,'sym5');
%图像第一层的重构逼近信号
a1=wrcoef2('a',c,l,'sym5',1);
%图像第二层的重构逼近信号
a2=wrcoef2('a',c,l,'sym5',2);

%图示消噪处理后的结果
subplot(2,2,3);
image(a1);
colormap(map);
title('第一层重构图像');

subplot(2,2,4);
image(a2);
colormap(map);
title('第二层重构图像');