一种基于小波变换和HSI颜色空间的图像融合方法与流程

一种基于小波变换和hsi颜色空间的图像融合方法
技术领域
1.本发明涉及图像处理技术领域，具体而言，涉及一种基于小波变换和hsi颜色空间的图像融合方法。


背景技术：

2.图像融合是信息融合中可视信息部分的融合，是将多源信道所采集的关于同一目标的图像经过一定的图像处理，提取各自信道的信息，最后综合成统一图像或综合图像特征以供观察和进一步处理。就目前而言，现有的遥感图像融合方法主要采用基于空间域的算法或基于变换域的算法，基于空间域的算法存在以下缺点：当高分辨率全色图像与低分辨率多光谱图像局部差异较大时，比较容易造成融合图像的光谱扭曲；而基于变换域的算法通常需要小波变换进行处理，然而小波变换不具有平移不变性，且不能最优逼近图像的边缘纹理，容易导致图像的几何信息无法有效地捕获，容易产生模糊效果的图像。基于此，我们提供了一种基于小波变换和hsi颜色空间的图像融合方法。


技术实现要素：

3.本发明的目的在于提供一种基于小波变换和hsi颜色空间的图像融合方法，其通过融合高分辨率图像中高频细节信息和低分辨率图像中低频整体视觉信息，能够最大程度避免丢失两幅影像的信息，并且通过将小波融合后的图像替代低分辨率多光谱图像的i分量，与hs分量组合成最终的融合图像，能够保持原始图像中的色彩信息，不仅能够适应于多种传感器图像，而且对于sar、光学等图像都有很好的融合效果。
4.本发明的实施例通过以下技术方案实现：
5.一种基于小波变换和hsi颜色空间的图像融合方法，该方法的步骤包括：
6.输入第一设定图像与第二设定图像；所述第一设定图像的分辨率大于所述第二设定图像的分辨率；
7.对第二设定图像进行颜色空间转换获取i分量与hs分量，分别对第一设定图像与i分量进行小波分解，得到第一频率系数、第二频率系数、第三频率系数及第四频率系数；所述第一频率系数与第二频率系数对应第一设定图像，且所述第一频率系数的频率低于所述第二频率系数的频率，所述第三频率系数及第四频率系数对应i分量，所述第三频率系数的频率低于第四频率系数的频率；
8.基于第一融合规则融合第二频率系数与第三频率系数，并对融合后的系数进行小波重构，得到第一小波融合图像，以表征为变换后的i'分量；
9.基于第二融合规则融合变换后的i'分量与hs分量，经过颜色空间转换后，得到融合结果图像。
10.可选的，所述第一设定图像的分辨率大于所述第二设定图像的分辨率，所述第一频率系数与第二频率系数对应第一设定图像，且所述第一频率系数的频率低于所述第二频率系数的频率，所述第三频率系数及第四频率系数对应i分量，所述第三频率系数的频率低
于第四频率系数的频率。
11.可选的，所述对第二设定图像进行颜色空间转换，其具体为：将第二设定图像的颜色空间从rgb转换为his；
12.所述经过颜色空间转换后，得到融合结果图像，其具体为：将融合后的i'分量与hs分量的颜色空间从h is转换为rgb。
13.可选的，所述对第二设定图像进行颜色空间转换，其计算公式如下：
[0014][0015][0016]
其中，h为色调，s为饱和度，i为强度或亮度，r为红色，g为绿色，b为蓝色。
[0017]
可选的，所述经过颜色空间转换后，得到融合结果图像，若0
°
≤h＜120
°
，则计算公式如下：
[0018][0019]
若120
°
≤h＜240
°
，则计算公式如下：
[0020][0021]
若240
°
≤h＜360
°
，则计算公式如下：
[0022][0023]
可选的，所述第一融合规则具体为：对第二频率系数与第三频率系数进行线性加权，其计算公式如下：
[0024]
hpi＝(wa×
msi
ilp
)+(wb×
pan
ihp
)
[0025]
其中，wa和wb均为权值，且wa+wb＝1；msi
ilp
为第二设定图像i方向的低频系数，
pan
ihp
为第一设定图像i方向的高频系数。
[0026]
可选的，所述第二融合规则具体为：将i'分量与hs分量组合成his颜色空间的融合图像，其计算公式如下：
[0027]
i＝(wa×
ms
hs
)+(wb×
pani)
[0028]
其中，wa和wb均为权值，且wa+wb＝1；ms
hs
为第二设定图像i方向的低频系数；pani为第一设定图像i方向的高频系数。
[0029]
本发明实施例的技术方案至少具有如下优点和有益效果：
[0030]
本发明实施例通过融合第一设定图像中高频细节信息和第二设定图像中低频整体视觉信息，能够最大程度避免丢失两幅影像的信息，并且通过将小波融合后的图像替代低分辨率多光谱图像的i分量，与hs分量组合成最终的融合图像，能够保持原始图像中的色彩信息，不仅能够适应于多种传感器图像，而且对于sar、光学等图像都有很好的融合效果。
附图说明
[0031]
图1为本发明实施例提供的一种基于小波变换和hsi颜色空间的图像融合方法的逻辑原理示意图；
[0032]
图2为本发明实施例提供的小波分解示意图；
[0033]
图3为本发明实施例提供的小波重构示意图。
具体实施方式
[0034]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
[0035]
参照图1所示，图1为本发明实施例提供的一种基于小波变换和hsi颜色空间的图像融合方法的逻辑原理示意图。
[0036]
在一些实施方式中，一种基于小波变换和hs i颜色空间的图像融合方法，该方法的步骤包括：
[0037]
输入第一设定图像与第二设定图像；
[0038]
对第二设定图像进行颜色空间转换获取i分量与hs分量，分别对第一设定图像与i分量进行小波分解，得到第一频率系数、第二频率系数、第三频率系数及第四频率系数；
[0039]
基于第一融合规则融合第二频率系数与第三频率系数，并对融合后的系数进行小波重构，得到第一小波融合图像，以表征为变换后的i'分量；
[0040]
基于第二融合规则融合变换后的i'分量与hs分量，经过颜色空间转换后，得到融合结果图像。
[0041]
在上述实现过程中，第一设定图像为高分辨率图像，第二设定图像为低分辨率图像，在输入前需要对两者分别进行图像配准校正，在输入后，对第二设定图像进行颜色空间转换，由rgb的颜色空间变换到his的颜色空间，然后对第一设定图像和第二设定图像的i分量进行小波变换，进行两级小波分解，得到各自的高频系数和低频系数(即所述第一频率系数与第二频率系数对应第一设定图像，且所述第一频率系数(低频系数)的频率低于所述第
二频率系数(高频系数)的频率，所述第三频率系数及第四频率系数对应i分量，所述第三频率系数(低频系数)的频率低于第四频率系数(高频系数)的频率)，对第一设定图像的高频系数和第二设定图像的低频系数按照第一融合规则进行融合，融合后的小波系数进行小波重构，得到小波融合图像作为变换后的i'分量与低分辨率的hs分量按照第二融合规则进行融合，最终将颜色空间由hsi转换为rgb，得到最终的融合结果图像。
[0042]
具体的，第二设定图像得颜色空间转换具体为：将颜色空间由rgb转换为his,其转换公式如下：
[0043][0044][0045]
其中，h为色调，s为饱和度，i为强度或亮度，r为红色，g为绿色，b为蓝色，转变为hsi颜色空间图像分割成h、s、i三个单通道图像。
[0046]
参照图2所示，图2为本发明实施例提供的小波分解示意图。
[0047]
在本实施例中，对第一设定图像和第二设定图像i分量进行相同层次二维小波分解，一层二维小波分解如图2中的a所示，每进行一次小波分解，图像的低频部分会分解为四个部分，分别是低频部分ll，水平高频垂直低频部分hl，水平低频垂直高频部分lh，水平高频垂直高频(对角高频)部分hh，得到的新的低频部分又可以在下个尺度下继续分解。分解一定层数后，图像的小波系数存储如图2中的b所示。
[0048]
在一些实施方式中，本实施例在建立融合图像的每个小波系数时，必须确定哪幅源图像的小波系数对融合有利，具有更多的信息或权重，将这个信息保留在融合图像中。在一幅图像的小波变换中，高频部分被认为是体现了图像的细节特征，而低频部分认为是图像的近似特征，图像的能量集中在低频部分，绝对值较大的高频部分系数对应于边缘这些较为显著的特征，所以大部分基于小波变换的图像融合算法主要研究如何选择合成图像中的高频系数，也就是表现细节的三个高频部分，从而达到保留图像边缘的目的。而对于同一场景的不同图像，小波变换后得到的低频部分一般近似，而高频部分差别却比较大。因此，本实施例将小波分解后的高频系数(第二频率系数)和低频系数(第三频率系数)进行线性加权，并对高频系数赋予较大权值(较大权值定义为wa+wb＝1，且wa＜wb)，第一融合规则的计算公式如下所示：
[0049]
hpi＝(wa×
msi
ilp
)+(wb×
pan
ihp
)
[0050]
其中，wa和wb均为权值，且wa+wb＝1；msi
ilp
为第二设定图像i方向的低频系数，即第三频率系数；pan
ihp
为第一设定图像i方向的高频系数，即第二频率系数。
[0051]
参照图3所示，图3为本发明实施例提供的小波重构示意图。
[0052]
在一些实施方式中，经过第一融合规则融合后，得到小波分解系数，对小波分解系数进行二维小波逆变换，得到第一小波融合图像，以表征为变换后的i'分量，本实施例通过第二融合规则将小波变换后得到的融合图像i'分量和第二设定图像的hs分量组合成新的his颜色空间的融合图像，第二融合规则的计算公式如下：
[0053]
i＝(wc×
ms
hs
)+(wd×
pani)
[0054]
其中，wc和wd均为权值；ms
hs
为第二设定图像的hs分量；pani为重构后的i’。
[0055]
在本实施例中，将第二融合规则融合后图像的颜色空间由his转换为rgb，需要考虑h的取值范围0
°
≤h＜360
°
，若h的取值在rg扇区时即0
°
≤h＜120
°
，则计算公式如下：
[0056][0057]
若h的取值在gb扇区时即120
°
≤h＜240
°
，则计算公式如下：
[0058][0059]
此时公式中的h＝h-120
°
，因此计算公式为
[0060]
若h的取值在br扇区时即240
°
≤h＜360
°
，则计算公式如下：
[0061][0062]
此时公式中的h＝h-240
°
，因此计算公式为
[0063]
以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种基于小波变换和hsi颜色空间的图像融合方法，其特征在于，该方法的步骤包括：输入第一设定图像与第二设定图像；所述第一设定图像的分辨率大于所述第二设定图像的分辨率；对第二设定图像进行颜色空间转换，获取i分量与hs分量，分别对第一设定图像与i分量进行小波分解，得到第一频率系数、第二频率系数、第三频率系数及第四频率系数，所述第一频率系数与第二频率系数对应第一设定图像，且所述第一频率系数的频率低于所述第二频率系数的频率，所述第三频率系数及第四频率系数对应i分量，所述第三频率系数的频率低于第四频率系数的频率；基于第一融合规则融合第二频率系数与第三频率系数，并对融合后的系数进行小波重构，得到第一小波融合图像，以表征为变换后的i'分量；基于第二融合规则融合变换后的i'分量与hs分量，经过颜色空间转换后，得到融合结果图像。2.根据权利要求1所述的基于小波变换和hsi颜色空间的图像融合方法，其特征在于，所述对第二设定图像进行颜色空间转换，其具体为：将第二设定图像的颜色空间从rgb转换为his；所述经过颜色空间转换后，得到融合结果图像，其具体为：将融合后的i'分量与hs分量的颜色空间从his转换为rgb。3.根据权利要求2所述的基于小波变换和hsi颜色空间的图像融合方法，其特征在于，所述对第二设定图像进行颜色空间转换，其计算公式如下：所述对第二设定图像进行颜色空间转换，其计算公式如下：其中，h为色调，s为饱和度，i为强度或亮度，r为红色，g为绿色，b为蓝色。4.根据权利要求3所述的基于小波变换和hsi颜色空间的图像融合方法，其特征在于，所述经过颜色空间转换后，得到融合结果图像，若0
°
≤h＜120
°
，则计算公式如下：若120
°
≤h＜240
°
，则计算公式如下：
若240
°
≤h＜360
°
，则计算公式如下：5.根据权利要求1所述的基于小波变换和hsi颜色空间的图像融合方法，其特征在于，所述第一融合规则具体为：对第二频率系数与第三频率系数进行线性加权，其计算公式如下：hp
i
＝(w
a
×
msi
ilp
)+(w
b
×
pan
ihp
)其中，w
a
和w
b
均为权值，且w
a
+w
b
＝1；msi
ilp
为第二设定图像i方向的低频系数，pan
ihp
为第一设定图像i方向的高频系数。6.根据权利要求3所述的基于小波变换和hsi颜色空间的图像融合方法，其特征在于，所述第二融合规则具体为：将i'分量与hs分量组合成his颜色空间的融合图像，其计算公式如下：i＝(w
a
×
ms
hs
)+(w
b
×
pan
i
)其中，w
a
和w
b
均为权值，且w
a
+w
b
＝1；ms
hs
为第二设定图像i方向的低频系数；pan
i
为第一设定图像i方向的高频系数。

技术总结
本发明涉及图像处理技术领域，具体而言，涉及一种基于小波变换和HSI颜色空间的图像融合方法，该方法的步骤包括：输入第一设定图像与第二设定图像；对第二设定图像进行颜色空间转换获取I分量与HS分量，分别对第一设定图像与I分量进行小波分解，得到第一频率系数、第二频率系数、第三频率系数及第四频率系数；基于第一融合规则融合第二频率系数与第三频率系数，并对融合后的系数进行小波重构，得到第一小波融合图像，以表征为变换后的I'分量；基于第二融合规则融合变换后的I'分量与HS分量，经过颜色空间转换后，得到融合结果图像。得到融合结果图像。得到融合结果图像。


技术研发人员：陈景伟 刘建明 张丽莎 梁斯东 张宝玉
受保护的技术使用者：北京道达天际科技股份有限公司
技术研发日：2023.05.10
技术公布日：2023/7/22