%圖象融合
load woman; %裝入原圖像
X1=X;map1=map;
subplot(221);
image(X1);
colormap(map1);
title('woman');
axis square %畫出woman 圖像
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 %畫出wbarb 圖像
%對上述二圖像進行分解
[c1,s1]=wavedec2(X1,2,'sym4');
[c2,s2]=wavedec2(X2,2,'sym4');
%對分解的圖像進行融合
c=c1+c2;
%應用融合系數進行圖像重構并顯示
xx=waverec2(c,s1,'sym4');
subplot(223);image(xx);
title('融合圖象1');
sizec1=size(c1);
for I=1:sizec1(2)
 c1( I )=1.2*c1( I );
end
sizec2=size(c2);
for j=1:sizec2(2)
 c2( j)=0.8*c2( j );
end
%通過減小融合系數以減小圖象寬度
c=0.5*(c1+c2);
%對融合系數進行圖象重構
xxx=waverec2(c,s2,'sym4');
%顯示重構結果
subplot(224);image(xxx);
title('融合圖象2');