cwq3.m

图象的消噪,图象的压缩,图象的融合,图象的增强

%基于小波的图象消噪
load wifs;
subplot(221);
image(X);
colormap(map);
title('原图 ');
axis square; %画出原图象
%生成含噪图象并显示
init=2055615866;
randn('seed',init)
XX=X+randn(size(X));
subplot(222);
image(XX);
colormap(map);
title('含噪声图象 ');
axis square; %画出含噪声图象
%下面对图象消噪
%首先用sym4小波函数对x进行两层分解
[c,s]=wavedec2(XX,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; %画出第二次低通滤波消噪后图象