在MATLAB中，能否提供详细的代码来说明如何使用StarFM算法融合TM遥感影像和两幅MODIS遥感影像的过程？

在MATLAB中，星形因子模型（StarFM）通常用于光谱数据融合，比如将TM（Thematic Mapper）遥感影像和MODIS（Moderate Resolution Imaging Spectroradiometer）遥感影像结合。以下是使用StarFM算法的一个简化版示例代码：

% 首先，假设我们有四个矩阵，A是TM影像的数据，B和C是MODIS影像的数据（这里假设A和B对应同一时间点）

% 对于每个波段，计算星形因子
SF_TM = (A - mean(A)) ./ std(A);
SF_MODIS_B = (B - mean(B)) ./ std(B);
SF_MODIS_C = (C - mean(C)) ./ std(C);

% 创建融合权重矩阵，可以根据需要调整权重
weights = [0.5; 0.25; 0.25]; % 假设TM权重0.5，MODIS B和C各占0.25

% 使用星形因子和权重融合图像
fused_image = weights(1) * SF_TM + weights(2) * SF_MODIS_B + weights(3) * SF_MODIS_C;

% 融合后的数据需转换回原始尺度
fused_image = fused_image .* std(A) + mean(A); % 回归到原数据范围

% 显示结果
imagesc(fused_image);
colormap(gray); % 可视化灰度图像
title('Fused Image using StarFM Algorithm');

%