遥感图像融合用Matlab代码

遥感图像融合是将多个遥感图像融合成一幅高质量的图像。Matlab提供了许多图像处理工具箱，可以用于遥感图像融合。以下是一些可能有用的Matlab代码示例：

1. 像素级融合

function fusedImage = pixelFusion(image1,image2)
%像素级融合函数，将两幅图像像素逐个相加并除以二
%输入：image1，image2，待融合的两幅图像
%输出：fusedImage，融合后的图像

fusedImage = (double(image1) + double(image2)) / 2;
fusedImage = uint8(fusedImage);
end

2. 小波变换融合

function fusedImage = waveletFusion(image1,image2)
%小波变换融合函数，将两幅图像分别进行小波变换，融合后进行反变换
%输入：image1，image2，待融合的两幅图像
%输出：fusedImage，融合后的图像

[cA1, cH1, cV1, cD1] = dwt2(image1,'haar');
[cA2, cH2, cV2, cD2] = dwt2(image2,'haar');

cA3 = (cA1 + cA2) / 2;
cH3 = (cH1 + cH2) / 2;
cV3 = (cV1 + cV2) / 2;
cD3 = (cD1 + cD2) / 2;

fusedImage = idwt2(cA3,cH3,cV3,cD3,'haar');
fusedImage = uint8(fusedImage);
end

3. PCA融合

function fusedImage = PCAFusion(image1,image2)
%PCA融合函数，将两幅图像进行主成分分析，提取主成分相加后还原
%输入：image1，image2，待融合的两幅图像
%输出：fusedImage，融合后的图像

data = double([image1(:),image2(:)]);

[coeff,score,latent] = pca(data);

fusedScore = score(:,1) + score(:,2);

fusedImage = reshape(coeff * fusedScore',size(image1));
fusedImage = uint8(fusedImage);
end

以上是三个常用的遥感图像融合方法的Matlab代码示例，你可以根据自己的需求进行调整。