1. 引言
随着科技不断发展,传感器技术不断进步。可获得大量不同种类的信息与数据。同一类型的传感器只能从一个方面获取信息,不能提供下游计算需要的所有信息。因此,在现代应用技术中,将多传感器信息融合为一体发挥着越来越重要的作用。信息融合技术是将不同传感器获得的信息融合为一体的技术。融合后的信息具有多个源信息的互补特征。
可见光图像一般由普通相机获取,和人类视觉几乎一致。获得的信息是具有高清细节、高分辨率的图像。但在下雨、夜晚、雾霾等情况下获得的图像几乎不能使用。红外图像是指红外传感器捕获信息源发出的热辐射能而形成的图像。其特点通常是分辨率低、细节纹理较模糊等。但不会受到来自恶劣天气、普通遮挡的影响。
由于可见光和红外图像优缺点具有互补性,互相结合后具有更大的利用价值。于是可见光和红外图像融合技术顺应而生。可见光与红外图像融合是指将可见光图像与红外图像融合为一张图像,融合后的图像具有两种类型图像互补信息,表达力更加丰富。
融合后的信息广泛应用于计算机视觉领域 [1] ,如:目标监控、目标跟踪、目标识别。农业领域,如水果采摘、水果质量评估、农作物特征检测。医学领域、遥感领域等 [2] 。
融合算法始于上个世纪,距今已有几十年的历史。根据研究者提出的融合算法架构,大致可分为多尺度分解法、稀疏表示法、神经网络法、其他混合法等。接下来章节分别对其种类进行概述,最后一小节对可见光与红外图像融合算法进行总结与展望。
2. 多尺度分解法
本多尺度分解法是深度学习算法提出之前最火的方法之一。其主要思想是首先将源图像分解为若干层,然后利用相应的融合规则对每一层进行融合,最后对每层融合后的图像进行逆变换,最终的得到融合图像。根据其融合框架可得知,分解算法与融合规则决定最终融合质量。
2.1. 金字塔分解
最早基于多尺度分解的可见光和红外图像融合算法是基于金字塔分解法。文献 [3] 产生拉普拉斯金字塔的多尺度变换包括四个步骤:低通滤波、子采样、插值和差分,所有这些都是在不同的尺度上完成的。图像的融合是通过组合源图像生成的两个拉普拉斯金字塔的差图像,然后应用逆变换得到最终融合图像。还有学者对传统金字塔进行改进,提出可操纵金字塔,具有多尺度、多方向和自反转的图像分解优点。
活动水平测量和权重分配是图像融合的两个关键部分。Li [4] 提出了一种新的基于主成分分析网络(PCANet)和图像金字塔的红外和可见光融合方法。首先,使用轻量级的学习网络PCANet来获得源图像的活动水平测量和权重分配。然后,通过图像金字塔将权重和源图像分解为多个尺度,并在每个尺度上应用加权平均融合规则。最后,通过重构得到融合后的图像。
2.2. 小波变换
第二种是基于离散小波变换法。文献 [5] 通过离散小波变换获得源图像的稀疏系数。然后通过改进的熵加权融合规则和基于绝对值最大的融合规则对红外和可见光图像的低系数和高系数进行融合。在使用随机高斯矩阵的局部线性投影之后,通过压缩采样匹配追踪算法重建融合图像。取得了不错的融合效果。
普通离散小波本身具有其局限性,如不具有平移不变性、提取多方向信息不丰富等缺点。为了克服这些缺点,学者们提出双树小波小波、四元数小波、提升平稳小波等,有效的克服普通小波的缺点。
还有其小波算法的变种,如Chen提出了一种基于深度小波密集网络(WT DenseNet) [6] 的红外和可见光图像融合方法。WT DenseNet包括三个网络层,即混合特征提取层、融合层和图像重建层。最后,通过图像重建层输出融合图像。
2.3. 多尺度几何分析
第三种多尺度分解法是多尺度几何分析(MGA)法。其方法主要思想来源也是克服小波表示图像的缺点。比如小波在二维空间中捕获方向信息的能力有限。大多数小波变换都无法从信号稀疏性的角度得到图像的最优表示等。多尺度几何分析法应用到可见光和红外图像融合中最初是轮廓波法。
Shen提出了一种基于非下采样轮廓波变换(NSCT)和改进的双边滤波器的红外和可见光图像融合算法 [7] ,该算法利用NSCT将图像分解为低频分量和高频分量。采用改进的双边滤波方法对两幅图像的高频分量进行处理,对图像的噪声进行滤波,计算红外图像高频分量的图像细节。对于低频系数的融合规则,采用了局部标准方差系数法。最后,根据NSCT的逆变换,得到融合图像。
多尺度几何分析除了轮廓波还有曲波变换法、剪切波法等等。
3. 稀疏表示法
由于多尺度分解法的融合规则几乎都是高频系数取绝对值最大,低频系数取平均。在一定程度上会造成融合结果模糊。于是学者将稀疏表示理论引入图像融合领域中。
稀疏表示的理论指任何一组信号都可以基于一个过完备字典去线性表示。即使用尽量少的字典中原子线性的去描述一组信号称为稀疏表示。基于稀疏表示 [8] 的可见光和红外图像融合算法基本步骤主要有四步。首先,使用滑动窗口策略将每个源图像分解为步长为1的重叠的小块。第二,学习过完备字典,过完备字典一般有离线学习和在线学习两种学习方式。并对每个块进行稀疏编码,利用学习到的过完备字典获取稀疏表示系数。第三,对稀疏表示系数进行融合。最后,利用过完备字典与融合后的系数重构出最终图像。过完备字典、稀疏编码、融合规则是融合质量的关键。
3.1. 离线字典学习
过完备字典一般有离线学习和在线学习两种学习方式。离线学习字典主要是采用固定基构造或通过自然图像训练得到。文献 [9] 采用了三种离线字典方式进行可见光和红外图像融合。第一种对不同频率的余弦波进行采样,构造过完备离散余弦变换(DCT)字典。第二种是混合字典,它由DCT基、小波“db1”基、Gabor基和脊波基组成。第三种是使用迭代K-SVD算法从自然样本中获得训练字典。训练数据从50幅自然图像的数据库中随机挑选,然后分割成50,000个8 × 8个图像块进行训练。
3.2. 在线字典学习
过完备字典在线学习法是指采用融合前的源图像进行学习构造字典。文献 [10] 提出了一种基于联合块聚类的高效字典学习方法。根据结构相似性,聚集来自不同传感器模态源图像的所有块,并且只选择几个有效描述每个联合块的主要部分,以形成过完备的字典。构造信息丰富且紧凑字典是研究目标。为了在字典构建中获得足够的稀疏表示信息,文献 [11] 提出一种新的几何字典构建方法。首先基于形态学相似性将源图像中的图像块分类为不同的组。然后根据随机坐标编码提取相应的图像块信息,最后构建成字典。
对于使用稀疏表示的图像融合方法,自适应字典和融合规则对多模态图像融合有很大影响,最大L1范数融合规则可能会造成融合结果的灰色不一致的现象。为了解决这个问题,Wang根据联合补丁聚类的自适应字典与稀疏表示相关理论,提出了一种改进的图像融合方法 [12] 。首先,使用高斯滤波器来分离高频和低频信息。第二,采用局部能量加权策略来完成低频融合。第三,使用联合补丁聚类算法重建过完备自适应学习字典。根据设计混合融合规则,融合高频部分。最后,通过将频域转换为空间域获得融合结果。
4. 神经网络法
近年来,深度学习算法受到学者们青睐。正是由于神经网络通常由许多神经元组成,可以模仿人脑的感知行为机制。神经网络还具备强大的特征提取和重构能力。深度学习的发展也促进了图像融合的巨大进步,使得可见光和红外图像融合很有前景 [13] 。基于神经网络的可见光和红外图像融合算法一般分为基于预训练网络法、基于自编码法、基于端对端网络法3类。
4.1. 基于预训练网络法
基于预训练网络法主要思想是将提前训练好的神经网络引入到可见光和红外图像融合。算法框架是首先利用通过固定网络获取图像深度特征,然后采用合适策略进行特征融合。最后用恰当方法来重建融合图像。如文献Xu [14] 提出将解耦表示应用于可见光和红外图像融合。首先根据不同源图像中的信息源进行分解表示。然后应用不同的策略来融合这些不同类型的表示。最后将融合的表示输入到预训练的生成器中,以生成融合结果。
4.2. 基于自编码网络法
自编码网络可见光与红外图像融合算法一般由编码器、融合层和解码器三部分组成。其中编码器从源图像中提取充分的深度特征,融合层根据相应融合策略用来融合深层特征。解码器根据编码的特征自适应地重建原始数据。通过解码器网络对融合后的图像进行重构。文献 [15] 的编码器包含C1与密集块(DenseBlock)两个部分。第一部分C1包含3 × 3的滤波器来提取粗略特征,第二部分DenseBlock包含三个卷积层(每个层的输出为下一层的输入),其中滤波器大小也是3 × 3。融合层选择了加法策略和L1范数策略两种融合策略来组合编码器获得的显著特征图。融合层的输出将作为解码器的输入。解码器包含四个卷积层,每层也是3 × 3大小的滤波器。此方法取得了不错的融合效果。但此方法在深层特征方向上比较冗余。
Ren [16] 提出了一种自编码网络法并命名为基于变分自动编码器法(VAE),该所提出的网络分为图像融合网络和红外特征补偿网络。首先,在图像融合网络中,创建图像融合网络的编码器,从源图像中生成隐藏层中的潜在向量。然后根据高斯概率密度的乘积,把两个不同的潜在向量合并为一个。随着损失函数值的下降,解码器开始重建融合图像。为了提高了网络训练的效率,在网络中加入了残差块和对称跳跃连接方法。并且为了克服融合网络中损失函数设置的缺陷,还设计了一个网络来补偿红外图像的关键辐射特征。此方法在一定程度上改善了特征中方向信息比较冗余问题,使融合后的特征更精细。
4.3. 基于端对端网络法
端对端网络法结构比较简单,输入源图像与输出最终融合结果之间就只有一个设计好的融合网络。关键优势在于,直接通过训练学习得到源图像和融合图像之间的映射,因此不需要过多的处理步骤。避免了像传统方法那样手动设计复杂的活动水平测量和融合规则。
为了克服手动且复杂设计融合规则等缺点,Ma等人提出一种生成对抗性端到端网络FusionGAN算法,并设计了一个专门用于红外和可见光图像融合的损失函数。可以从输入源图像自动生成融合图像,而无需手动设计活动水平测量或融合规则。该方法在生成器和鉴别器之间建立了一个对抗性游戏,其中生成器主要生成的融合图像的红外强度信息与可见光图像梯度信息,鉴别器旨在使融合图像具有可见光图像中更多的细节信息。相对来说,融合结果不太清晰,细节较丰富。
更大程度为了使融合结果目标更清晰,细节更丰富。Ma [17] 继续改进FusionGAN提出了一种基于细节保留对抗性学习的端到端模型。它能够克服传统融合方法中活动水平测量和融合规则手工设计和复杂设计的局限性。考虑到红外和可见光图像的特定信息,设计了两个损失函数。包括细节损失和目标边缘增强损失,以提高细节信息的质量,并在生成对抗性网络的框架下锐化红外目标的边缘。此方法使融合图像能够保留热辐射的同时,还可以锐化红外图像中的目标边界和可见光图像中丰富的纹理细节。其融合结果取得了优异成绩。
为了进一步避免手工设计的复杂且不准确的融合规则,与充分利用红外和可见光图像中的信息。Li提出了一种实用的融合方法RCGAN [18] 。将耦合生成对抗性网络开创性引入到图像融合领域且提出了使用预融合图像作为标签的策略。实验表明RCGAN法融合结果优秀,能够有效地保持可见光图像中丰富的纹理和红外图像中的热辐射信息。
值的注意的是,一个良好的模型应该应用相对广泛。Xu提出了一种统一且无监督的端到端图像融合网络U2Fusion [19] ,它根据特征提取和信息测量,能够自动估计被融合图像的重要性。能对多模态、多曝光和多聚焦图像进行融合。通过实验,验证了U2Fusion在三个典型图像融合任务上的有效性和通用性。
5. 其他方法
除了上述总结的三类方法以外,学者还提出了基于其他理论法以及混合方法等。比如刘等人提出多尺度变换与稀疏表示结合的方法 [20] 。其融合框架中,首先对每个源图像执行多尺度变换,获得它们的低频系数和高频系数。然后,使用基于稀疏表示法合并低频系数,高频系数的合并采用绝对值最大法。最后对合并后的系数执行多尺度逆变换,最终获得融合后的图像。还有基于主成分分析理论(PCA)法、独立成分分析法等。
为了解决红外和可见光图像融合中亮度和细节信息丢失的问题,Li [21] 提出了一种利用滚动引导滤波和梯度显著性图进行红外和可见图像融合的有效方法。滚动引导滤波用于将输入图像分解为近似层和残差层;能量属性融合模型用于融合近似层;通过重构融合的近似层子图像和残差层子图像来生成融合图像。
然而,现有的大多数融合方法都是用于灰度图像融合任务,不能应用于真实的融合成像系统。为了解决这些问题,Mo [22] 提出了一种属性引导的红外和可见光图像融合算法,首先,对红外图像进行属性滤波,得到显著物体的初始权重图。然后,利用边缘保持滤波器来细化初始权重图。最后,采用拉普拉斯金字塔融合策略将源图像与所获得的权重图进行融合。与以往的方法相比,该方法利用了红外图像的固有特性,设计了一种基于属性滤波器的特征提取规则,无需训练步骤或先验信息,即可有效地提取红外图像中的突出目标。
6. 总结与展望
本文主要对红外与可见光图像融合算法研究进展做梳理。使新学者快速掌握此研究方向脉络,并为本领域做出贡献。对主流算法与结果的主要总结如下:
1) 融合算法大致可分为多尺度变换、稀疏表示、神经网络以及其他方法混合法等几大类。并对每个大类进行细分类与探讨。
2) 神经网络深度学习算法目前最火,但其融合结果并不是完全碾压传统算法 [23] 。
3) 主观上看,多尺度分解法融合结果比较稳定,但细节保留不足,尤其是对红外目标淡化。
4) 在纹理复杂的图像中,稀疏表示法有时细节模糊,可能原因是字典的表达能力不足,且其最大的缺点就是计算代价太大,不适用于实时任务。
近几十年的发展,红外和可见光图像融合方面已经取得了长足的进步,但仍有一些问题需要进一步研究解决。对本课题展望如下:
1) 大多数基于多尺度变换法的融合规则都比较简单(如:高频绝对值最大,低频取平均),如何设计有价值融合规则值得进一步探讨。
2) 一般稀疏表示法采用滑动窗口技术对图像分块,然后融合每个块,最终融合图像由融合后的块拼接而成。其并未考虑块与块之间的关联性,且此算法计算开销很大,不适合实时处理的系统。如何考虑关联性,值得关注。
3) 神经网络算法中端到端算法较火,但设计损失函数非常主观。如何有效的考虑融合与后续应用之间的相关性来设计设计损失函数。是一个有趣的研究点。
4) 由于深度学习算法不是完全碾压传统算法,学者大可不必对传统算法保持悲观。开发有效的传统特征提取算法或者传统算法与深度学习算法结合是一个突破点。如文献 [24] 提出的显著性检测与多层潜在低秩表示的融合方法。
5) 研究理论的目标是最终落实到实时应用。然而应用到真实的成像系统较少,开发计算代价小与实时处理的平台也是一个重要的发展点。
基金项目
重庆对外经贸学院科学研究项目(KYKJ202002)。