example7_15.m

小波分析教材算例的源程序

%装载原始图像（woman.mat和wbarb.mat）并显示
load wbarb;
X1=X;map1=map;
subplot(2,2,1);
image(X1);colormap(map1);
title('图像wbarb');
load woman;
X2=X;map2=map;
subplot(2,2,2);
image(X2);colormap(map2);
title('图像woman');
%=====================================
%对上述二图像进行分解
[c1,l1]=wavedec2(X1,2,'sym4');
[c2,l2]=wavedec2(X2,2,'sym4');
%对分解系数进行融合
c=c1+c2;
%=====================================
%应用融合系数进行图像重构并显示
XX=waverec2(c,l1,'sym4');
subplot(2,2,3);
image(XX);
title('融合图像1');
Csize1=size(c1);
%=====================================
%对图像进行增强处理
for i=1:Csize1(2)
    c1(i)=1.2*c1(i);
end
Csize2=size(c2);
for j=1:Csize2(2)
    c2(j)=0.8*c2(j);
end
%=====================================
%通过减小融合系数以减小图像的亮度
c=0.5*(c1+c2);
%=====================================
%对融合系数进行图像重构
XXX=waverec2(c,l2,'sym4');
%=====================================
%显示重构结果
subplot(2,2,4);
image(XXX);
title('融合图像2');