load tartan;
X1=X;map1=map;
load sinsin;
X2=X;map2=map;
% 打開圖像
subplot(121)
image(X1),colormap(map1);
title('圖像tartan')
subplot(122)
image(X2),colormap(map2);
title('圖像sinsin')

%%%    圖像直接融合
figure,subplot(131)
image((X1+X2)/2),colormap(map2);
% 在空域內圖像直接融合
title('兩圖像直接相加融合')

%%%    圖像直接融合
figure,subplot(131)
image((X1+X2)/2),colormap(map2);
% 在空域內圖像直接融合
title('兩圖像直接相加融合')

%%%    圖像小波變換融合
[C1,L1]=wavedec2(X1,2,'sym4');
[C2,L2]=wavedec2(X2,2,'sym4');
% 分別對上述兩圖像進行分解
C=C1+C2;
% 對分解系數進行融合
X=waverec2(C,L1,'sym4');
% 融合后的信號進行圖像重構
subplot(133)
image(X/2),colormap(map2);
% 顯示小波變換的圖像融合
title('兩圖像小波變換融合')