一、获取代码方式

获取代码方式1:
通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码。

获取代码方式2:
通过紫极神光博客主页开通CSDN会员,凭支付凭证,私信博主,可获得此代码。

获取代码方式3:
完整代码已上传我的资源:【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】

备注:开通CSDN会员,仅只能免费获得1份代码(有效期为开通日起,三天内有效);
订阅紫极神光博客付费专栏,可免费获得1份代码(有效期为订阅日起,三天内有效);

二、IHS简介

IHS图像融合方法千千万,我的只是其中一种!小白刚上手,记录每一分钟!希望高人指点!
颜色模型通常采用的红、绿、蓝(RGB)三原色模型
另外一种广泛应用的颜色模型是亮度、色调、饱和度(IHS) 颜色模型
亮度表示光谱的整体亮度大小,对应图像的空间分辨率,与地物表面粗糙度相对应, 代表地物的空间几何特征。(光源或物体明亮程度)
色调描述纯色的属性,决定与光谱的主波长。(物体什么颜色)
饱和度表征光谱的主波长在强度中的比例,色调和饱和度代表图像的光谱分辨率。(物体颜色的浓淡)
下图为IHS 图像融合的步骤
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_matlab
下图为IHS 图像此方法的变换矩阵
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_matlab_02
那么以上就是我们需要的知识。下面是图像

MS(多光谱) 图像
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_计算机视觉_03
PAN(全色)图像
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_matlab_04

三、部分源代码

clear
g_R=0;                                        %r清晰度描述
g_G=0;                                        %g清晰度描述
g_B=0;                                        %b清晰度描述
h_R=0;                                        %熵的描述 
h_G=0;
h_B=0;
u1=zeros(3,1);
u2=zeros(3,1);
v11=zeros(3,1);               %高分辨率RGB值
v12=zeros(3,1);               %低分辨率RGB值
AA=zeros(3,3);
BB=zeros(3,3);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
AA=[1/sqrt(3),1/sqrt(3),1/sqrt(3);       
    1/sqrt(6),1/sqrt(6),-2/sqrt(6);                    %RGB->IHS正变换矩阵
    1/sqrt(2),-1/sqrt(2),0];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
BB=[1/sqrt(3),1/sqrt(6),1/sqrt(2);
    1/sqrt(3),1/sqrt(6),-1/sqrt(2);                    %IHS->RGB逆变换矩阵
    1/sqrt(3),-2/sqrt(6),0];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
RGB_up=imread('high.jpg');         %读图像
RGB_low=imread('low.jpg');         
figure(1)
imshow(RGB_up);   
title('IHS-RGB表示的高分辨率图像');
figure(2)                              
imshow(RGB_low); 
title('IHS-RGB表示的低分辨率图像');

[M,N,color]=size(RGB_up);

R_up=RGB_up(:,:,1);
G_up=RGB_up(:,:,2);
B_up=RGB_up(:,:,3);
[hang,lie,color]=size(RGB_up);
for i=1:hang
    for j=1:lie
       v11(1)=double(RGB_up(i,j,1));             %把函数类型unit8为doublev11(2)=double(RGB_up(i,j,2));
       v11(3)=double(RGB_up(i,j,3));  
       v12(1)=double(RGB_low(i,j,1));            %把函数类型unit8为doublev12(2)=double(RGB_low(i,j,2));
       v12(3)=double(RGB_low(i,j,3)); 
       u1=AA*v11;                               %正变换
       u2=AA*v12;
       u2(1)=u1(1);                             %交换空间明度
                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%空间信息明度定义????????
       v12=BB*u2;
       RGB(i,j,1)=v12(1);                       %逆变换
       RGB(i,j,2)=v12(2);
       RGB(i,j,3)=v12(3);   
   end
end
%subplot(2,3,5);
r=RGB(:,:,1);
g=RGB(:,:,2);
b=RGB(:,:,3);
RGB=uint8(round(RGB)); 
figure(3)
 
title('IHS-RGB表示的转化图像');

四、运行结果

【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_图像融合_05
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_开发语言_06
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_matlab图像处理_07
【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】_matlab_08

五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社,2020.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社,2013.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社,2013.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社,2015.