一、简介
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:hang
for j=1:lie
v11(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表示的转化图像');
三、运行结果
四、备注
版本：2014a
