【图像融合】基于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: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表示的转化图像');


  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
  
15
  
16
  
17
  
18
  
19
  
20
  
21
  
22
  
23
  
24
  
25
  
26
  
27
  
28
  
29
  
30
  
31
  
32
  
33
  
34
  
35
  
36
  
37
  
38
  
39
  
40
  
41
  
42
  
43
  
44
  
45
  
46
  
47
  
48
  
49
  
50
  
51
  
52
  
53
  
54
  
55
  
56
  
57
  
58
  
59
  
60
  
61
  
62
  
63
  
64
 

四、运行结果

五、matlab版本及参考文献

1 matlab版本
2014a

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

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/115481813