一种基于多方向梯度滤波的多模态遥感图像融合方法

1.本发明属于图像处理领域，具体涉及一种基于多方向梯度滤波的多模态图像融合方法。


背景技术：

2.近年来，多传感器技术和计算机视觉技术发展如火如荼，在此背景下，多传感器图像融合技术得到广泛应用。例如，在遥感领域，大量遥感图像的融合将具有互补性的源图像进行合并，获得清晰度高、信息量大且纹理细节丰富的融合图像，实现对探测场景的精准描述，为更方便更全面地认识环境和自然资源提供了技术途径。在医学领域，多模态图像融合技术极大地提高了计算机辅助诊疗的能力，在很大程度上节约了医学检查时间。在各个领域的需求牵引下，各国学者对多传感器图像融合技术的研究也越来越重视。
3.多模态图像融合是指将不同模态的传感器所获得的同一地理坐标下的景物图像或同一传感器在不同时刻所获得的同一地理坐标下的景物图像经过滤波去噪、时间配准、空间对齐以及重采样之后，通过图像融合技术进行融合；融合后的图像突破了单源图像在空间分辨率、物理属性、信息量以及清晰度等方面的局限性，具有较高的图像质量，从而有利于对物理现象和事件进行定位、识别和精确解译。
4.根据信息表征层次的不同、多模态图像融合可分为三个等级，分别为像素级、特征级以及决策级。与特征级和决策级图像融合相比，像素级图像融合更具有普适性，其处理对象为图像像素，利用某种图像融合策略提取源图像中的有用信息剔除冗余信息，更大程度上保留源图像的特征和背景信息，所获得的融合结果表现出高度的全面性、精准性以及可靠性。
5.现有的像素级多模态图像融合方法存在多种。例如，有nsst(nonsubsampled shearlet transform，非下采样剪切波变换)方法以及latlrr(latent low-rank representation，潜在低秩表示)方法等。这类现有方法的大致过程包括：通过对多模态图像进行图层分解，获取图像不同性质的表征信息，通过构建相应的融合策略进行图层融合，最后经过逆重构操作得到相应的多模融合图像。
6.然而，利用现有方法所得到的多模融合图像存在对比度低、边缘模糊等问题，图像质量有待进一步提高。


技术实现要素：

7.为了解决现有技术中所存在的上述问题，本发明提供了一种基于多方向梯度滤波的多模态遥感图像融合方法。
8.本发明要解决的技术问题通过以下技术方案实现：
9.一种基于多方向梯度滤波的多模态图像融合方法，包括：
10.获取待进行融合的多张不同模态的源图像，并对每张所述源图像均进行拉普拉斯分解，得到多个拉普拉斯金字塔序列；其中，每个所述拉普拉斯金字塔序列均包括一个低频
层和多个高频层；
11.将每个所述低频层均代入至一滤波优化问题中进行迭代求解，得到该低频层对应的滤波结果，并以该滤波结果为该低频层对应的显著层，计算该低频层对应的残差层；其中，所述滤波优化问题是以滤波效果为优化目标的优化问题，所述优化目标至少包括滤波前后图像的亮度分布相似性以及滤波后的八方向总梯度；
12.基于目标源图像的显著层计算显著层融合因子，并基于各残差层计算残差层融合因子；其中，所述目标源图像为所述多张不同模态的源图像中的红外图像或sar图像；
13.基于所述显著层融合因子和所述残差层融合因子对不同模态的源图像分别进行显著层融合以及残差层融合，并利用显著层融合结果和残差层融合结果得到低频层融合结果；
14.对不同模态的源图像进行高频层融合，并基于高频层融合结果以及所述低频层融合结果采用拉普拉斯逆变换法进行图像重构，得到多模融合图像。
15.可选地，所述滤波优化问题为：
[0016][0017]
其中，用o表示迭代过程中的滤波图像，oc表示所述滤波图像的任一像素点c的像素值，该像素点c的坐标为(x,y)，ic表示低频层i中坐标同为(x,y)的像素点的像素值；||
·
||2表示l2范数；范数；分别表示左、右、上、下、右上、右下、左上、左下共八个方向的梯度算子；p
c,l
、p
c,r
、p
c,u
、p
c,d
、p
c,ru
、p
c,rd
、p
c,lu
、p
c,ld
分别为为对应的惩罚系数；每个所述惩罚系数均是根据对应的梯度算子进行计算得到的；o
*
为迭代求解得到的最优解，即所述滤波结果。
[0018]
可选地，所述惩罚系数的计算方式如下：
[0019][0020]
其中，direct表示方向，分别为左、右、上、下、右上、右下、左上、左下；α为预设的梯度阈值，β为预设的斜率控制因子，p
max
为预设的最大惩罚系数。
[0021]
可选地，基于目标源图像的显著层计算显著层融合因子的方式，包括：
[0022]
m＝mean2(a1)，
[0023]
λ＝10*m,
[0024]
s＝std2(a1)，
[0025]
r＝|a1|，
[0026][0027][0028]
其中，粗体的a1表示目标源图像的显著层，未加粗的a1表示该显著层的任一像素点，r为a1的模值，r
max
为该显著层的所有像素点中的最大模值；mean2(
·
)表示求平均值，std2(
·
)表示求标准差；w为对应于所述目标源图像的一显著层融合因子，对应于非目标源图像的显著层融合因子根据w确定。
[0029]
可选地，所述基于各残差层计算残差层融合因子，包括：
[0030]
采用sobel算子计算每个残差层的sobel梯度图；
[0031]
利用高斯滤波器对每个所述sobel梯度图进行扩散，得到对应的扩散梯度图；
[0032]
基于各扩散梯度图分别计算各残差层对应的残差层融合因子。
[0033]
可选地，所述基于各扩散梯度图分别计算各残差层对应的残差层融合因子，包括：
[0034][0035]
其中，表示第k∈[1，k]个残差层对应的扩散梯度图，sr表示各扩散梯度图的总和，k为所述多张不同模态的源图像的总张数，k≥2。
[0036]
可选地，所述高斯滤波器的高斯滤波模板以及标准差的值均取5。
[0037]
可选地，所述对不同模态的源图像进行高频层融合，包括：
[0038]
利用绝对值取大策略对不同模态的源图像进行高频层融合。
[0039]
可选地，所述多张不同模态的源图像具有不可忽略的视差；
[0040]
在对每张所述源图像均进行拉普拉斯分解之前，所述多模态图像融合方法还包括：对所获取的多张不同模态的源图像进行配准对齐。
[0041]
可选地，所述多张不同模态的源图像中包含有可见光图像；
[0042]
所述多模态图像融合方法还包括：
[0043]
在对所获取的多张不同模态的源图像进行配准对齐之前，对所述可见光图像进行ihs变换，并提取变换结果中的i分量形成新的图像，来替换所述可见光图像作为源图像；
[0044]
以及，在得到所述多模融合图像之后，对所述多模融合图像进行ihs逆变换。
[0045]
本发明提供的基于多方向梯度滤波的多模态遥感图像融合方法中，利用拉普拉斯金字塔的多尺度性以及多方向梯度滤波的多方向性提取了源图像中不同尺度的基础信息和不同方向的纹理信息，由此可剔除源图像中的噪声边缘和细小的纹理信息，更好地保留了图像的整体亮度和边缘信息，从而使多模融合图像具有较高的对比度、且边缘清晰，提高了多模融合图像的图像质量。
[0046]
本发明中，基于目标源图像的显著层来动态地计算显著层融合因子，即采用了一种自适应模糊耦合机制实现显著层融合；并且，本发明中在基于各残差层计算残差层融合因子的过程中，采用了一种梯度显著性映射策略来实现残差层的融合，具体是利用高斯滤波器对残差层的sobel梯度图进行扩散，并基于各扩散梯度图分别计算各残差层对应的残
差层融合因子。由此，在利用显著层融合结果和残差层融合结果得到低频层融合结果后，可以使低频层融合结果能够尽可能保留红外图像或sar图像中的显著信息以及可见光图像中的纹理细节信息，进一步地提高了多模融合图像的质量。
[0047]
此外，本发明中还利用了绝对值取大策略对不同模态的源图像进行高频层融合，使得融合后的高频层能最大程度地保留源图像的边缘信息，由此可更进一步提高多模融合图像的质量。
[0048]
以下将结合附图及对本发明做进一步详细说明。
附图说明
[0049]
图1是本发明实施例提供的一种基于多方向梯度滤波的多模态图像融合方法的流程图；
[0050]
图2是拉普拉斯金字塔分解的过程示意图；
[0051]
图3是根据图1所示方法对包含有可见光图像的2张多模图像进行融合时的流程示意图；
[0052]
图4(a)示出了一组待进行融合的可见光图像和sar图像；
[0053]
图4(b)中分别示出了利用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法对图4(a)中两张图像进行融合得到的多模融合图像；
[0054]
图5(a)示出了另一组待进行融合的可见光图像和sar图像；
[0055]
图5(b)中分别示出了利用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法对图5(a)中两张图像进行融合得到的多模融合图像；
[0056]
图6(a)为一组待进行融合的可见光图像和红外图像；
[0057]
图6(b)中分别示出了利用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法对图6(a)中两张图像进行融合得到的多模融合图像；
[0058]
图7(a)为一组待进行融合的可见光图像和红外图像；
[0059]
图7(b)中分别示出了利用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法对图7(a)中两张图像进行融合得到的多模融合图像。
具体实施方式
[0060]
下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。
[0061]
多模态图像融合技术也称异源图像融合处理技术，其可对同一地物目标所对应的多传感器图像信息进行多方面、深层次的组合筛选，在整合优化互补信息的同时剔除掉冗余信息。
[0062]
异源图像的成像机理以及图像表征形式差别明显，能分别获取到目标在不同电磁波谱段的信息，具有明显的互补性和融合价值。但是不容忽视的是，异源图像间的灰度差异以及因不可忽略的视差导致的几何扭曲使得一般的图像融合结果出现光谱失真以及对比度低等现象，图像质量较差。
[0063]
因此，针对以上问题，为了提高多模态图像融合质量，本发明实施例提供了一种基于多方向梯度滤波的多模态图像融合方法，来实现异源图像的高质量融合，从而获取信息
度丰富的融合结果。本发明实施例提供的多模态图像融合方法是一种像素级的多模态图像融合方法。
[0064]
参见图1所示，本发明实施例提供的该多模态图像融合方法包括以下步骤：
[0065]
s1：获取待进行融合的多张不同模态的源图像，并对每张源图像均进行拉普拉斯分解，得到多个拉普拉斯金字塔序列；其中，每个拉普拉斯金字塔序列均包括一个低频层和多个高频层。
[0066]
这里，待进行融合的多张不同模态的源图像可以包括可见光图像、红外图像以及sar(合成孔径雷达)图像中的任意两种或全部的三种。其中，当源图像包括可见光图像时，可以对可见光图像进行ihs(也称作hsi)变换，并提取变换结果中的i分量形成新的图像，来替换可见光图像作为源图像。ihs是一种数字图像模型，它以色调(hue，简写为h)、饱和度(saturation，简写为s)和亮度(intensity，简写为i)三种基本特征分量来感知颜色。
[0067]
其中，ihs变换的过程如下：
[0068][0069][0070][0071]
其中，r、g、b分别是可见光图像的r分量、g分量以及b分量；v1和v2是中间变量；h、s、i分别是his变换结果中的h分量、s分量、i分量。
[0072]
金字塔分解是多尺度融合领域的一种主要实现手段，也是研究最广泛的图像分解算法。拉普拉斯金字塔是对高斯金字塔的一种改进，也被称之为预测残差金字塔，其优势在于可以对图像进行最大程度的还原且分解速度较快。拉普拉斯金字塔分解的过程可以参见图2所示，图2表示的是在一维空间中由低层到高层的分解过程，实际的分解过程是在二维空间进行的。
[0073]
具体而言，将原始的源图像作为高斯金字塔的最底层图像，该图像可以看作r行c列的数组i0，每一个像素值代表该位置的强度。金字塔第一层的图像是底层图像的低通滤波和下采样之后的结果，然后金字塔第二层的图像是在第一层图像基础上进行低通滤波和下采样所得到的，以此类推构造高斯金字塔。分解过程可用公式进行表示如下：
[0074][0075]
1≤l≤n,0≤j≤r
l
,0≤j≤c
l
；
[0076]
其中，i
l
(i,j)表示第l层高斯金字塔图像，n表示金字塔总层数。r
l
和c
l
表示第l层
高斯金字塔图像的行数和列数，ω(m,n)为一个二维的窗函数，其定义为：ω(m,n)＝h(m)
×
h(n)；其中，h(
·
)为满足以下约束条件的高斯密度分布函数：
[0077][0078]
h(i)＝h(-i)，
[0079]
h(0)+h(-2)+h(2)＝h(-1)+h(1)。
[0080]
根据窗函数的定义，通过计算可得实际分解过程中使用的窗函数ω(m,n)为：
[0081][0082]
通过以上公式，对每张源图像均进行拉普拉斯分解，可以得到该张源图像的高斯金字塔子图像序列i0,i1,
…
,in。
[0083]
构造完所有的高斯金字塔子图像序列之后，将相邻层之间两两相减即可得到对应的拉普拉斯金字塔序列。在这一过程中，由于相邻层的高斯金字塔子图的尺寸各不相同，不满足矩阵差分运算的要求；因此，在每次做矩阵差分运算之前对上层相邻的高斯金字塔子图进行升采样，使得上层和下层的高斯金字塔子图的尺寸一致，从而满足矩阵差分运算的要求。
[0084]
具体的，将i
l
进行插值膨胀之后变为此时的尺寸大小和i
l-1
一致，用公式表示如下：
[0085][0086]
1≤l≤n,0≤j≤r
l
,0≤j≤c
l
；
[0087]
其中，
[0088][0089]
经过插值膨胀之后，可得到一系列的插值序列然后用该插值序列实现相邻层之间两两相减，即可得到拉普拉斯金字塔序列lp1,lp2,
…
,lpn。其中，第l层拉普拉斯金字塔图像的定义如下：
[0090][0091]
其中，经历升采样的次数最多、图像最为模糊的lpn为低频层，剩余的均为高频图
层。
[0092]
s2：将每个低频层均代入至一滤波优化问题中进行迭代求解，得到该低频层对应的滤波结果，并以该滤波结果为该低频层对应的显著层，计算该低频层对应的残差层；其中，滤波优化问题是以滤波效果为优化目标的优化问题，优化目标至少包括滤波前后图像的亮度分布相似性以及滤波后的八方向总梯度。
[0093]
在一种实现方式中，上述滤波优化问题可以为：
[0094][0095][0096][0097]
其中，为了保留源图像的亮度分布，对滤波图像的亮度分布进行了限制，使其像素的亮度分布与源图像尽可能相似，因此定义了误差e1；并且，为了尽可能去除较小的无用梯度变化，比如噪声边缘的梯度变化，对滤波图像在八个方向的总梯度e2进行了限制。
[0098]
该滤波优化问题的表达式中，o表示迭代过程中的滤波图像，oc表示所述滤波图像的任一像素点c的像素值，该像素点c的坐标为(x,y)，ic表示低频层i中坐标同为(x,y)的像素点的像素值；||
·
||2表示l2范数；范数；分别表示左、右、上、下、右上、右下、左上、左下共八个方向的梯度算子；o
*
为迭代求解得到的最优解，即迭代求解得到的滤波结果。
[0099]
可以理解的是，基于该滤波优化问题所得到的滤波结果，既可以在亮度分布方面与源图像保持相似，又可以对细小纹理以及噪声等无用的小梯度值进行去除，保留更多有用的梯度。
[0100]
在另一种实现方式中，为了在剔除细小纹理或噪声等同时更多地保留大梯度的边缘，可以通过引入梯度阈值来区分边缘、细小纹理和噪声。同时，基于该梯度阈值设计惩罚系数，对于梯度较大的边缘给予较小的惩罚系数，而对于梯度较小的纹理和噪声应当给予较大的惩罚系数；由此，上述滤波优化问题可以进一步优化为：
[0101][0102]
其中，p
c,l
、p
c,r
、p
c,u
、p
c,d
、p
c,ru
、p
c,rd
、p
c,lu
、p
c,ld
分别为分别为对应的惩罚系数；每个惩罚系数均是根据对应的梯度算子进行计算得到的，计算方式如下：
[0103][0104]
其中，direct表示方向，分别为左、右、上、下、右上、右下、左上、左下；α为预设的梯度阈值，β为预设的斜率控制因子，p
max
为预设的最大惩罚系数。
[0105]
利用滤波优化问题得到低频层对应的滤波结果后，该滤波结果即为低频层对应的显著层。然后，用低频层减去显著层，即得到低频层对应的残差层。
[0106]
s3：基于目标源图像的显著层计算显著层融合因子，并基于各残差层计算残差层融合因子；其中，目标源图像为多张不同模态的源图像中的红外图像或sar图像。
[0107]
具体的，当待进行融合的多张不同模态的源图像包括可见光图像以及红外图像时，目标源图像即为其中的红外图像。当待进行融合的多张不同模态的源图像包括可见光图像以及sar图像时，目标源图像即为其中的sar图像。当待进行融合的多张不同模态的源图像包括红外图像以及sar图像时，可以在这两者之间选取目标源图像，具体选取原则可根据源图像的成像质量来确定，优先选取成像质量更高的一种模态下的图像作为目标源图像。同理，当待进行融合的多张不同模态的源图像包括可见光图像、红外图像以及sar图像三种图像时，选取目标源图像的原则不变。
[0108]
其中，基于目标源图像的显著层计算显著层融合因子的方式，可以基于一种自适应模糊耦合函数来计算，其定义如下：
[0109]
m＝mean2(a1)，
[0110]
λ＝10*m,
[0111]
s＝std2(a1)，
[0112]
r＝|a1|，
[0113][0114][0115]
其中，粗体的a1表示目标源图像的显著层，是一个矩阵图像；未加粗的a1表示该矩阵图像的任一像素点，r为a1的模值，r
max
为该矩阵图像的所有像素点中的最大模值；mean2(
·
)表示求平均值，std2(
·
)表示求标准差；p∈[0,1]；w为对应于目标源图像的一显著层融合因子，具体是对应于目标源图像的显著层中像素点a1的因子；对应于非目标源图像的显著层融合因子根据w确定。具体而言，若非目标源图像仅包括可见光图像，则对应于可见光图像的一显著层融合因子为1-w；若非目标源图像除了包括可见光图像之外，还包括一红外图像或sar图像，则对应于该红外图像或sar图像的显著层融合因子为w'，0＜w'＜w，对应于可见光图像的显著层融合因子为1-w-w'。
[0116]
本发明实施例中，w反映了目标源图像的显著层的明暗特性分布，w越大代表明暗特性越明显。在实际的融合过程中，为避免融合图像的对比度的损失，在构建自适应模糊耦
合函数时，应使w尽量大。因此，在其他计算显著层融合因子的方式中，λ也可以取10*m以外的其他值，如8*m或者12*m等，以此来试出更大的w。或者，也可以采用其他方式来对r进行归一化得到p，这都是合理且可实现的；任何基于目标源图像的显著层动态自适应地计算显著层融合因子的方式，均可以应用到本发明实施例中。
[0117]
残差层主要反应低频图层中的小梯度变化，包括细小纹理和噪声，因此可以根据图像的梯度构造残差层的梯度显著性映射矩阵，从而基于该梯度显著性映射矩阵实现残差层的融合，使得融合后的残差层尽可能保留可见光图像中的纹理信息。
[0118]
具体的，该步骤s3中基于各残差层计算残差层融合因子的方式存在多种。
[0119]
示例性的，在一种计算方式中，基于各残差层计算残差层融合因子，可以包括：
[0120]
(1)采用sobel算子计算每个残差层的sobel梯度图；
[0121]
这里，sobel梯度图即是所构造的梯度显著性映射矩阵；所采用的sobel算子用公式表示为：
[0122]
sobel(s(x,y))＝|(s
(x+1,y-1)
+2s
(x+1,y)
+s
(x+1,y+1)-s
(x-1,y-1)
+2s
(x-1,y)
+s
(x-1,y+1)
)|+|(s
(x-1,y+1)
+2s
(x,y+1)
+s
(x+1,y+1)-s
(x-1,y-1)
+2s
(x,y-1)
+s
(x+1,y-1)
)|；
[0123]
其中，sobel(
·
)表示sobel算子，s(x,y)表示残差层中坐标为(x,y)的像素点，s
(x,y-1)
、s
(x,y+1)
、s
(x-1,y)
、s
(x+1,y)
、s
(x-1,y-1)
、s
(x+1,y-1)
、s
(x-1,y+1)
、s
(x+1,y+1)
分别为该像素点的上、下、左、右、左上、右上、左下、右下的像素点。
[0124]
假设用sk表示残差层图像，gk表示残差层图像对应的sobel梯度图，则计算残差层的sobel梯度图的过程可以用如下公式来表示：
[0125]gk
＝sobel(sk)；
[0126]
其中，k∈[1，k]，k为待进行融合的多张不同模态的源图像的总张数，k≥2。
[0127]
(2)基于各sobel梯度图分别计算各残差层对应的残差层融合因子。
[0128]
具体的计算方式如下：
[0129][0130]
其中，sk表示第k个残差层对应的sobel梯度图，s表示各sobel梯度图的总和。
[0131]
在另一种计算方式中，为了增大梯度对残差层融合因子的影响，基于各残差层计算残差层融合因子，可以包括：
[0132]
(1)采用sobel算子计算每个残差层的sobel梯度图。
[0133]
该步骤(1)的计算方式与上一种计算方式中相同。
[0134]
(2)利用高斯滤波器对每个sobel梯度图进行扩散，得到对应的扩散梯度图。
[0135]
具体的，该步骤可以用公式进行表示如下：
[0136][0137]
其中，gaussian(
·
)表示高斯滤波算子，gk表示残差层图像对应的sobel梯度图，δ表示其标准差，r表示其高斯滤波模板，δ和r的具体取值可预先试出；例如，δ和r可以均取5，并不局限于此。
[0138]
可以理解的试，此时扩散梯度图即是所构造的梯度显著性映射矩阵。
[0139]
(3)基于各扩散梯度图分别计算各残差层对应的残差层融合因子。
[0140]
具体的，计算方式如下：
[0141][0142]
其中，表示第k∈[1，k]个残差层对应的扩散梯度图，sr表示各扩散梯度图的总和。
[0143]
s4：基于显著层融合因子和残差层融合因子对不同模态的源图像分别进行显著层融合以及残差层融合，并利用显著层融合结果和残差层融合结果得到低频层融合结果。
[0144]
其中，显著层融合的过程为求取对应于各模态源图像的显著层的加权平均。
[0145]
举例而言，若各模态源图像包括可见光图像和红外图像，则显著层融合的过程可表示为fa＝w*a1+(1-w)*a2；其中，a1此时表示的是红外图像的显著层，w为其对应的显著层融合因子，a2此时表示的是可见光图像的显著层；fa表示显著层融合结果。
[0146]
同理，残差层层融合的过程也是求取对应于各模态源图像的残差层的加权平均。
[0147]
举例而言，若各模态源图像包括可见光图像和sar图像，则显著层融合的方式可表示为其中，s1代表sar图像对应的残差层，为其对应的残差层融合因子；s2代表可见光图像对应的残差层，为其对应的残差层融合因子，fs表示残差层融合结果。
[0148]
得到显著层融合结果和残差层融合结果后，便可以得到低频层融合结果fl＝fa+fs。
[0149]
s5：对不同模态的源图像进行高频层融合，并基于高频层融合结果以及低频层融合结果采用拉普拉斯逆变换法进行图像重构，得到多模融合图像。
[0150]
优选地，可以采用绝对值取大策略对不同模态的源图像进行高频层融合。可以理解的是，图像的高频部分通常包含图像的边缘轮廓信息，并且这部分信息可以表示图像相应位置的信息丰富度。在图像分解过程中，绝对值较大的高频系数对应于亮度突然变化的点，即图像中对比度变化较大的边缘特征，如边界、亮线和区域轮廓。高频部分绝对值越大对应位置的信息越详细，因此利用绝对值取大策略实现高频层的融合可以得到边缘信息丰富的高频融合图像。
[0151]
举例而言，假设不同模态的源图像包括两张，则利用绝对值取大策略进行高频层融合的过程可以用公式表示如下：
[0152][0153]
其中，0≤l《n，和为分别位于两张不同模态源图像的高频层中相同位置处的两个像素点，fh
l
为这两个像素点的高频层融合值，位于目标源图像中。
[0154]
得到高频层融合结果后，结合低频层融合结果，采用拉普拉斯逆变换进行图像重构，具体过程如下：
[0155]
[0156]
其中，将低频层融合结果fl作为第l＝n层的高低频融合结果fn，逐层对第l∈[0，n-1]层的高频层融合结果fh
l
以及第l+1层的高低频融合结果f
l+1
进行融合，最终得到多模融合图像f0。融合过程中，为了便于进行矩阵求和，对f
l+1
进行插值膨胀，得到
[0157]
在一个实施例中，参加图3所示，若待进行配准的多张不同模态的源图像中包含有可见光图像，则在得到多模融合图像之后，再对该多模融合图像进行ihs逆变换，从而恢复颜色信息，得到彩色的多模融合图像。
[0158]
其中，ihs逆变换的过程如下所示：
[0159]v1
＝s cos(h)，
[0160]v2
＝s sin(h)，
[0161][0162]
其中，分量i即是待进行ihs变换的多模融合图像，其余各参数的释义可参加上文对ihs变换的描述，此处不再赘述。
[0163]
本发明实施例提供的基于多方向梯度滤波的多模态遥感图像融合方法中，利用拉普拉斯金字塔的多尺度性以及多方向梯度滤波的多方向性提取了源图像中不同尺度的基础信息和不同方向的纹理信息，由此可剔除源图像中的噪声边缘和细小的纹理信息，更好地保留了图像的整体亮度和边缘信息，从而使多模融合图像具有较高的对比度、且边缘清晰，提高了融合图像的图像质量。
[0164]
本发明实施例中，基于目标源图像的显著层来动态地计算显著层融合因子，即采用了一种自适应模糊耦合机制实现显著层融合；并且，本发明中采用了一种梯度显著性映射策略来实现残差层的融合，具体是在基于各残差层计算残差层融合因子的过程中，利用高斯滤波器对残差层的sobel梯度图进行扩散，并基于各扩散梯度图分别计算各残差层对应的残差层融合因子。由此，在利用显著层融合结果和残差层融合结果得到低频层融合结果后，可以使低频层融合结果能够尽可能保留红外图像或sar图像中的显著信息以及可见光图像中的纹理细节信息，进一步地提高了多模融合图像的质量。
[0165]
此外，本发明实施例中还可利用绝对值取大策略对不同模态的源图像进行高频层融合，使得融合后的高频层能最大程度地保留源图像的边缘信息，由此也可更进一步提高多模融合图像的质量。
[0166]
可选地，在一种实现方式中，若多张不同模态的源图像具有不可忽略的视差，则在对每张源图像均进行拉普拉斯分解之前，还可以先行对所获取的多张不同模态的源图像进行配准对齐。
[0167]
这里，对图像进行配准对齐可以采用现有的sift(scale-invariant feature transform，尺度不变特征变换)方法实现，本发明实施例不再进行赘述。
[0168]
下面采用试验结果对本发明实施例的有益效果进行进一步验证和说明。
[0169]
图4(a)为一组待进行融合的可见光图像和sar图像；其中，左侧为可见光图像，右侧为sar图像；图4(b)中分别示出了采用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法，对图4(a)中两张图像进行融合得到的融合图像。
[0170]
图5(a)为另二组待进行融合的可见光图像和sar图像，其中，左侧为可见光图像，右侧为sar图像；图5(b)中分别示出了采用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法，对图5(a)中两张图像进行融合得到的多模融合图像。
[0171]
图4(b)和图5(b)在实际中均为彩色图像，其中利用本发明实施例所提供方法得到的多模融合图像具有更好的视觉效果，色彩失真小，色彩信息保持得较好。
[0172]
基于以上试验结果，表1给出了利用三种不同方法所得到的多模融合图像的性能分析对比结果。
[0173]
表1
[0174][0175]
该表1中，ag表示平均梯度，e表示信息熵，psnr表示峰值信噪比，qabf表示融合质量，rmse表示均方根误差，sf表示空间频率，ssim表示结构相似性，vif为基于视觉信息保真度的指标。
[0176]
表1中，加粗的数值为同一指标中最好的结果。从表1可以看出，利用本发明实施例所提供方法得到多模融合图像锁具有的高性能指标更多，这说明了本发明实施例所提供方法具有更好的优势。
[0177]
图6(a)为一组待进行融合的可见光图像和红外图像，其中，左侧为可见光图像，右侧为红外图像；图6(b)中分别示出了采用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法，对图6(a)中两张图像进行融合得到的多模融合图像。
[0178]
图7(a)为另一组待进行融合的可见光图像和红外图像，其中左侧为可见光图像，右侧为红外图像；图7(b)中分别示出了采用现有nsst方法、现有latlrr方法以及本发明实施例所提供方法，对图7(a)中两张图像进行融合得到的多模融合图像。
[0179]
图6(b)和图7(b)在实际中均为彩色图像，其中利用本发明实施例所提供方法得到的多模融合图像具有更好的视觉效果，色彩信息保持得较好；并且，利用本发明实施例所提供方法得到的多模融合图像中，对红外图像的显著性区域保持得较好；例如，在图6(b)的各
子图的左下角处，本发明更多的保留了屋顶以及树冠的显著性特征，而同一区域在另外两种对比方法中的显著性则不够强。而在图7(b)的各子图中，本发明更多的保留了远处行人的显著性特征，而同一区域在另外两种对比方法中的显著性则不够强。
[0180]
基于以上试验结果，表2给出了利用三种不同方法所得到的多模融合图像的性能分析对比结果。
[0181]
表2
[0182][0183][0184]
同样从表2可以看出，本发明实施例所提供方法在相比之下具有较好的优势。
[0185]
本发明实施例中，首先将拉普拉斯金字塔分解与多方向梯度滤波进行结合提取源图像不同尺度的基础信息和不同方向的细节信息，以此来剔除源图像中的噪声边缘和细小的纹理信息，更好地保留图像的整体亮度和边缘信息。其次，本发明实施例设计来自适应模糊耦合机制和梯度显著性映射策略实现显著层和残差层的融合，再对融合后的显著层和残差层进行相加重构作为融合后的低频层，使得低频层融合结果能够尽可能保留红外图像或sar图像中的显著信息以及可见光图像中的纹理细节信息。最后设计绝对值取大策略实现高频层的融合，使得融合后的高频层能最大程度地保留源图像的边缘信息，将融合后的低频层和高频层进行重构，得到的多模融合图像具有较高的图像质量。
[0186]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。此外，本领域的技术人员可以将本说明书中描述的不同实施例或示例进行接合和组合。
[0187]
尽管在此结合各实施例对本技术进行了描述，然而，在实施所要求保护的本技术过程中，本领域技术人员通过查看所述附图、公开内容、以及所附权利要求书，可理解并实现所述公开实施例的其他变化。
[0188]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定
本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。