hhh.m

融合遥感图像

%%%此程序实现遥感图像的融合
%%% 输入原始图像noon_GMT1.jpg
[A1,map1]=imread('noon_GMT1.jpg');
figure;imshow(A1);colormap(map1);title('原始图像1');
%%% 输入原始图像noon_GMT2.jpg
[B1,map2]=imread('noon_GMT2.jpg');
figure;imshow(B1);colormap(map2);title('原始图像2');
[m,n]=size(A1(:,:,1));
%%%%求A1与B1的共同部分
for k=1:3
    for i=1:m
        for j=1:n
            C(i,j,k)=min(A1(i,j,k),B1(i,j,k));
        end;
    end;
end;
AA=A1-C;
BB=B1-C;
A2=A1-BB;
B2=B1-AA;
for k=1:3
    for i=1:m
        for j=1:n
            if abs(B1(i,j,k))<abs(A1(i,j,k))
                D(i,j,k)=0;
            else 
                D(i,j,k)=B1(i,j,k)-A1(i,j,k);
            end;
        end;
    end;
end;
a=0.5;b=0.4;c=0.1;
%%%根据比例因子融合
G=a*double(A2)+b*double(B2)+c*double(D);
% % 假彩色增强
% h=ones(5,5)/49;
% G=imfilter(G,h);
figure;imshow(uint8(G));colormap(map1);title('融合后的图像');