myfusion.m

介绍了一中简单实用的图像数据融合程序,该程序运行快,融合结果优!

clear
[imA,map1] = imread('A.tif');
%imA = double(imA) / 256;
figure,imshow(imA);
[imB,map2] = imread('B.tif');
%imB = double(imB) / 256;
figure,imshow(imB);

X1=double(imA); 
X11=uint8(X1);
figure
subplot(2,2,1);imshow(X11);
% colormap(map1);
colormap(gray);
title('input1');
axis square


X2=double(imB);
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

X22=uint8(X2);
subplot(2,2,2);imshow(X22);
% colormap(map2)
colormap(gray);
title('input2');
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;
xx=waverec2(c,s1,'sym4');
subplot(223);imshow(xx,[]);
colormap(gray);
title('融合图像');
axis square;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%