%裝載原圖像1
load bust
X1=X;
map1=map;
%畫出woman圖像
subplot(221);image(X1);
colormap(map1);title('原圖像1');
axis square
%裝載原圖像2
load mask
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
%畫出wbarb圖像
subplot(222);
image(X2);colormap(map2);title('原圖像2');
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.5*c;
s=s1+s2;
s=0.5*s;
xx=waverec2(c,s,'sym4');
%畫出融合后的圖像
subplot(223);image(xx);title('融合圖像');
axis square