【图像融合】基于matlab IHS图像融合【含Matlab源码 724期】

一、获取代码方式

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

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

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、IHS简介

IHS图像融合方法千千万，我的只是其中一种！小白刚上手，记录每一分钟！希望高人指点！
颜色模型通常采用的红、绿、蓝(RGB)三原色模型
另外一种广泛应用的颜色模型是亮度、色调、饱和度(IHS) 颜色模型
亮度表示光谱的整体亮度大小，对应图像的空间分辨率，与地物表面粗糙度相对应, 代表地物的空间几何特征。（光源或物体明亮程度）
色调描述纯色的属性，决定与光谱的主波长。（物体什么颜色）
饱和度表征光谱的主波长在强度中的比例，色调和饱和度代表图像的光谱分辨率。（物体颜色的浓淡）
下图为IHS 图像融合的步骤

下图为IHS 图像此方法的变换矩阵

那么以上就是我们需要的知识。下面是图像

MS（多光谱）图像

PAN（全色）图像

三、部分源代码

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:hangfor j=1:liev11(1)=double(RGB_up(i,j,1));             %把函数类型unit8为double型v11(2)=double(RGB_up(i,j,2));v11(3)=double(RGB_up(i,j,3));  v12(1)=double(RGB_low(i,j,1));            %把函数类型unit8为double型v12(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版本及参考文献

1 matlab版本
2014a

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