请提供一个多光谱影像和全色影像融合的算法

以下是一种基于IHS变换的多光谱影像和全色影像融合算法：

1. 将多光谱影像和全色影像分别进行直方图均衡化，增强图像对比度。
2. 对多光谱影像进行IHS变换，得到三个分量：强度（I），色调（H），饱和度（S）。
3. 将全色影像转换到IHS空间，得到三个分量：强度（I_p），色调（H_p），饱和度（S_p）。
4. 将全色影像的强度分量与多光谱影像的强度分量进行融合，得到新的强度分量I_fuse。
5. 将IHS空间中的分量转换回RGB空间，得到融合后的图像。

下面是 MATLAB 代码实现：

% 读入多光谱影像和全色影像
ms_image = imread('ms_image.tif');
pan_image = imread('pan_image.tif');

% 直方图均衡化
ms_image_eq = histeq(ms_image);
pan_image_eq = histeq(pan_image);

% IHS 变换
ms_image_ihs = rgb2hsv(ms_image_eq);
pan_image_ihs = rgb2hsv(pan_image_eq);

% 融合强度分量
I_fuse = pan_image_ihs(:,:,3);
I_fuse(:,:,2) = ms_image_ihs(:,:,2);
I_fuse(:,:,1) = ms_image_ihs(:,:,1);

% 转换回 RGB 空间
fused_image = hsv2rgb(I_fuse);

% 显示和保存融合后的图像
imshow(fused_image);
imwrite(fused_image, 'fused_image.tif');

该算法的优点在于简单易实现，能够保留多光谱影像的色彩信息，同时利用全色影像的高空间分辨率，提高图像的细节表现能力。但是，该算法可能会引入椒盐噪声，同时全色影像与多光谱影像的光谱响应不同，可能会引入颜色失真。