matlab IHS图像融合代码

以下是一份基于Matlab实现的IHS图像融合代码，仅供参考：

% IHS图像融合
% 作者：XXX
% 时间：XXXX-XX-XX

% 读入待融合的三幅图像
I_H = imread('high_resolution_image.jpg');
I_S = imread('panchromatic_image.jpg');
I_MS = imread('multispectral_image.jpg');

% 将三幅图像转换到IHS空间
I_H_IHS = rgb2hsv(I_H);
I_S_IHS = rgb2hsv(I_S);
I_MS_IHS = rgb2hsv(I_MS);

% 分离IHS空间通道
I_H_I = I_H_IHS(:, :, 3);
I_H_H = I_H_IHS(:, :, 1);
I_H_S = I_H_IHS(:, :, 2);

I_S_I = I_S_IHS(:, :, 3);
I_S_H = I_S_IHS(:, :, 1);
I_S_S = I_S_IHS(:, :, 2);

I_MS_I = I_MS_IHS(:, :, 3);
I_MS_H = I_MS_IHS(:, :, 1);
I_MS_S = I_MS_IHS(:, :, 2);

% 对I和S通道进行直方图匹配，使其在直方图上更加接近
I_H_I_matched = histeq(I_H_I, histeq(I_S_I));
I_H_S_matched = histeq(I_H_S, histeq(I_S_S));
I_MS_I_matched = histeq(I_MS_I, histeq(I_S_I));
I_MS_S_matched = histeq(I_MS_S, histeq(I_S_S));

% 将匹配后的I和S通道和原始的H通道合并到一起
I_H_IHS_matched = cat(3, I_H_H, I_H_S_matched, I_H_I_matched);
I_MS_IHS_matched = cat(3, I_MS_H, I_MS_S_matched, I_MS_I_matched);

% 将匹配后的IHS空间图像转换回RGB空间
I_H_matched = hsv2rgb(I_H_IHS_matched);
I_MS_matched = hsv2rgb(I_MS_IHS_matched);

% 显示结果
subplot(1, 3, 1), imshow(I_H), title('High Resolution Image');
subplot(1, 3, 2), imshow(I_MS), title('Multispectral Image');
subplot(1, 3, 3), imshow(I_MS_matched), title('Fused Image');

需要注意的是，这份代码中使用的是rgb2hsv函数将图像转换到IHS空间，而不是真正的IHS变换。此外，在直方图匹配时，我们只对I和S通道进行匹配，而保持H通道不变。最后，注意将匹配后的IHS空间图像转换回RGB空间以得到最终的融合图像。