l3.22(融合).m

简单的图象融合

load woman;
X1=X;map1=map;
subplot(221);image(X1);colormap(map1);
title('woman');
axis square
load wbarb;
X2=X;map2=map;
for i=1:256
    for j=1:256
        if(X2(i,j)>100)
            X2(i,j)=1.2* X2(i,j);
        else
            X2(i,j)=0.5* X2(i,j);
        end
    end
end
subplot(222);image(X2);colormap(map2);
title('wbarb');
axis square

[c1,s1]=wavedec2(X1,2,'sym4');
sizec1=size(c1);
for i=1:sizec1(2)
    c1(i)=1.2*c1(i);
end
[c2,s2]=wavedec2(X2,2,'sym4');
c=c1+c2;
c=0.2*c;
xx=waverec2(c,s1,'sym4');
subplot(223);image(xx);
title('融合图象');
axis square