本发明涉及图像融合技术领域，具体涉及一种红外图像与可见光图像的融合方法、系统及相关设备，目的在于克服现有技术方法在图像融合时易出现的模态混叠现象，以及融合结果中高频细节信息丢失、图像对比度降低的问题。本发明的融合方法包括：基于W‑BEMD分解算法分别将两种源图像均分解成一系列高频细节子图像和一个低频轮廓子图像；基于局部区域方差选择与加权的融合规则，将两种源图像分解得到的高频细节子图像进行融合；基于区域能量选择与加权的融合规则，将两种源图像分解得到的低频轮廓子图像进行融合；将融合后的高频细节子图像与融合后的低频轮廓子图像进行W‑BEMD逆变换，得到融合图像。本发明算法冗余度低，而且运算速度快。
红外图像与可见光图像的融合方法、系统及相关设备
技术领域
[0001] 本发明涉及图像融合技术领域，具体涉及一种红外图像与可见光图像的融合方法、系统及相关设备。
背景技术
[0002] 由于成像方式的不同，红外传感器和可见光传感器获取的关于同一场景的图像存在很大的差异。红外传感器根据目标物体与环境的整体热辐射差异进行成像，不易受到照明条件和恶劣天气的影响，具有较强的抗干扰能力，所得图像能够很好地反映复杂背景中隐藏的目标及其位置信息，但是对比度较低，细节信息不明显。而可见光传感器的成像则利用了目标与场景的反射光信息，所得图像较红外图像具有更清晰的背景信息，但不能有效地显示有遮挡的目标，且抗干扰能力弱，在照明、天气等自然条件恶劣的情况下，成像质量易受影响。因此，同一场景的红外图像和可见光图像是两类差异大、互补信息和冗余信息丰富的图像。充分利用这两类图像的特征，将其融合生成一副既目标突出，又对比度高、细节清晰的高质量图像具有重要的意义。
[0003] 随着图像融合研究的深入，新的算法被不断提出。其中，基于多分辨分析理论的图像融合方法逐渐引起国内外学者们的广泛关注。该类方法的基本原理是：首先利用已有的多分辨分析工具对源图像进行多尺度多方向分解，然后选用有效的融合规则对源图像的各对应部分分别进行融合。由于多分辨分析工具能较好地捕获图像的细节信息，所以基于多分辨分析理论的图像融合方法通常能获得较满意的图像融合结果。然而，以图像金字塔、小波变换以及包括Contourlet和非下采样Contourlet变换在内的多尺度几何变换为代表的多分辨分析工具对图像的分解并非完全基于图像自身的特性，而是与滤波器的选取有关，自适应能力较差。而经验模态分解(EMD)是一种自适应的信号分析工具，其分解过程完全从数据本身出发，不依赖于任何预先指定的滤波器或基函数，是一种比小波变换性能更优的时频局部化分析方法。
[0004] 文献“郑有志,覃征.基于二维经验模态分解的医学图像融合算法[J].软件学报,2009,20(5):1096-1105”将BEMD应用于医学图像融合，所提出的算法具有一定的自适应性，无需预先给定滤波器或小波函数，具有比传统融合算法更好的性能。文献“胡钢,郑皎月,秦新强,结合局部邻域特性和C-BEMD的图像融合方法[J],计算机辅助设计与图形学学报,29(4),2017:607-615”基于三角剖分和协调操作，提出了一种可协调二维经验模态分解方法(C-BEMD)，并将其应用于多聚焦图像、医学图像、以及红外与可见光图像融合中。文献“MU Ahmed,DP Mandic,Image fusion based on Fast and Adaptive Bidimensional
Empirical Mode Decomposition[J],Information Fusion,2011:1-6.”构建了一种新的BEMD实现方法，该方法为每个源图像生成相同数量的BIMF，加快了图像融合的速度，并获得了较好的图像融合效果。针对红外与可见光图像的融合，文献“朱攀,黄战华,基于二维经验模态分解和高斯模糊逻辑的红外与可见光图像融合[J],光电子·激光,28(10),2017:
1157-1162”提出了基于BEMD和高斯模糊逻辑的融合方法。文献“P.Zhu,Z.H.Huang,H.Lei,Fusion of infrared and visible images based on BEMD and NSDFB[J],Infrared Physics&Technology,77(2016)82–93.”提出了基于BEMD和非下采样方向滤波器组(NSDFB)的融合方法。
[0005] 现有文献(特别是以上文献)公开的方法，主要是通过构建或选择不同的BEMD工具和不同的融合规则，甚至结合其他的数学分析工具来改进不同类型图像的融合结果。但是，各种BEMD方法中都或多或少存在模态混叠现象，另外，由于对经验模态分解结果的后续处理未能较好地捕获高频成分中的奇异点，从而导致融合结果中高频细节信息丢失、图像对比度降低。
发明内容
[0006] 为了解决现有技术中的上述问题，即为了解决现有技术中在图像融合时易出现的模态混叠现象，以及融合结果中高频细节信息丢失、图像对比度降低等问题，本发明提出了一种红外图像与可见光图像的融合方法、系统及相关设备，提高了运算速度，并且提高了融合后图像的清晰度。
[0007] 本发明的一方面，提出一种红外图像与可见光图像的融合方法，所述方法包括：
[0008] 步骤S10，基于W-BEMD分解算法，将红外图像分解成预设数量的第一高频细节子图像和一个第一低频轮廓子图像，并将可见光图像分解成所述预设数量的第二高频细节子图像和一个第二低频轮廓子图像；
[0009] 步骤S20，基于局部区域方差选择与加权的融合规则，将所述预设数量的所述第一高频细节子图像与所述预设数量的所述第二高频细节子图像进行融合，得到所述预设数量的第三高频细节子图像；
[0010] 步骤S30，基于区域能量选择与加权的融合规则，将所述第一低频轮廓子图像和所述第二低频轮廓子图像进行融合，得到第三低频轮廓子图像；
[0011] 步骤S40，将所述预设数量的所述第三高频细节子图像与所述第三低频轮廓子图像进行W-BEMD逆变换，得到融合图像。
[0012] 优选地，步骤S10中“基于W-BEMD分解算法，将红外图像分解成预设数量的第一高频细节子图像和一个第一低频轮廓子图像，并将可见光图像分解成所述预设数量的第二高频细节子图像和一个第二低频轮廓子图像”的步骤包括：
[0013] 步骤S11，将所述红外图像和所述可见光图像的残差分量分别设置为RIR,j-1＝XIR、RVI,j-1＝XVI；
[0014] 其中，RIR,j-1和RVI,j-1分别为所述红外图像的第j-1个残差分量和所述可见光图像的残差分量，j-1表示残差分量的序号，且j＝1；为残差分量的序号，XIR和XVI分别为预设的所述红外图像的残差分量的初始值和预设的所述可见光图像的残差分量的初始值；
[0015] 步骤S12，根据所述红外图像的残差分量RIR,j-1对所述红外图像进行BEMD分解得到所述红外图像的内蕴模函数IMFIR,j和残差分量RIR,j；根据所述可见光图像的残差分量RVI,j-1对所述可见光图像进行BEMD分解得到所述可见光图像的内蕴模函数IMFVI,j和残差分量RVI,j；
[0016] 步骤S13：分别对所述红外图像的残差分量RIR,j和所述可见光图像的残差分量RVI,j进行W变换，并对变换所得低频系数进行W逆变换，得到残差分量四分之一大小的低频轮廓图像 和
[0017] 步骤S14，根据所述红外图像的残差分量RIR,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFIR,j求和，计算出第j个所述第一高频细节子图像，并更新残差分量RIR,j；根据所述可见光图像的残差分量RVI,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFVI,j求和计算出第j个所述第
二高频细节子图像，并更新残差分量RVI,j；
[0018] 步骤S15，j＝j+1；若j≤L，转至步骤S13；否则，将更新后的残差分量RIR,L与RVI,L分别作为所述第一低频轮廓子图像WRIR,L与所述第二低频轮廓子图像WRVI,L；其中，L为所述预设数量。
[0019] 优选地，步骤S14中“根据所述红外图像的残差分量RIR,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFIR,j求和，计算出第j个所述第一高频细节子图像，并更新残差分量RIR,j；根据所述可见光图像的残差分量RVI,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFVI,j求和，计算出
第j个所述第二高频细节子图像，并更新残差分量RVI,j”的步骤包括：
[0020] 步骤S141，对所述红外图像的低频轮廓图像 进行插值膨胀，得到与RIR,j具有相同大小的图像
[0021] 步骤S142，按照下式计算残差图像RIR,j的高频细节图像：
[0022]
[0023] 步骤S143，按照下式计算第j个所述第一高频细节子图像：
[0024]
[0025] 并更新残差图像
[0026] 步骤S144，对所述可见光图像的低频轮廓图像 进行插值膨胀，得到与RVI,j具有相同大小的图像
[0027] 步骤S145，按照下式计算残差图像RVI,j的高频细节图像：
[0028]
[0029] 步骤S146，按照下式计算第j个所述第二高频细节子图像：
[0030]
[0031] 并更新残差图像
[0032] 优选地，步骤S20中“基于局部区域方差选择与加权的融合规则，将所述预设数量的所述第一高频细节子图像与所述预设数量的所述第二高频细节子图像进行融合，得到所述预设数量的第三高频细节子图像”的步骤包括：
[0033] 步骤S21，分别计算所述第一高频细节子图像与所述第二高频细节子图像在各像素点(x,y)处的局部区域方差：
[0034]
[0035]
[0036] 其中，WIMFIR,j和WIMFVI,j分别为第j个所述第一高频细节子图像和第j个所述第二高频细节子图像；j＝1,2,...,L，L为所述预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽； 和 分别代表图像WIMFIR,j和WIMFVI,j在以(x,y)为中心的M×N邻域内的灰度平均值；ω(m,n)代表预设的第一权值矩阵；
[0037] 步骤S22，计算所述第一高频细节子图像WIMFIR,j和所述第二高频细节子图像WIMFVI,j在各像素点(x,y)处的局部区域方差匹配度：
[0038]
[0039] 其中，
[0040]
[0041]
[0042] 为WIMFIR,j和WIMFVI,j在(x,y)处的相关局部区域方差；
[0043] 步骤S23，若 则根据下式计算所述第三高频细节子图像：
[0044]
[0045] 否则，根据下式计算所述第三高频细节子图像：
[0046]
[0047] 其中，
[0048]
[0049] β(x,y)＝1-α(x,y)
[0050] α(x,y)和β(x,y)均为加权系数，T1为预设的第一匹配度阈值。
[0051] 优选地，步骤S30中“基于区域能量选择与加权的融合规则，将所述第一低频轮廓子图像和所述第二低频轮廓子图像进行融合，得到所述第三低频轮廓子图像”的步骤包括：
[0052] 步骤S31，分别计算所述第一低频轮廓子图像与所述第二低频轮廓子图像在各像素点(x,y)处的局部区域能量：
[0053]
[0054]
[0055] 其中，WRIR,L和WRVI,L分别为所述第一低频轮廓子图像和所述第二低频轮廓子图像；L为所述预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽；ω′(m,n)代表预设的第二权值矩阵；
[0056] 步骤S32，计算所述第一低频轮廓子图像WRIR,L和所述第二低频轮廓子图像WRVI,L在各像素点(x,y)处的局部区域能量匹配度：
[0057]
[0058] 其中，
[0059]
[0060]
[0061] 为WRIR,L和WRVI,L在(x,y)处的相关局部区域能量；
[0062] 和 分别代表图像WRIR,L和WRVI,L在以(x,y)为中心的M×N邻域内的灰度平均值；
[0063] 步骤S33，若 则根据下式计算所述第三低频轮廓子图像：
[0064]
[0065] 否则，根据下式计算所述第三低频轮廓子图像：
[0066]
[0067] 其中，
[0068]
[0069] β′(x,y)＝1-α′(x,y)
[0070] α′(x,y)和β′(x,y)均为加权系数，T2为预设的第二匹配度阈值。
[0071] 本发明的第二方面，提出一种红外图像与可见光图像的融合系统，所述系统包括：
[0072] 待融合图像分解模块，配置为：基于W-BEMD分解算法，将红外图像分解成预设数量的第一高频细节子图像和一个第一低频轮廓子图像，并将可见光图像分解成所述预设数量的第二高频细节子图像和一个第二低频轮廓子图像；
[0073] 高频细节子图像融合模块，配置为：基于局部区域方差选择与加权的融合规则，将所述预设数量的所述第一高频细节子图像与所述预设数量的所述第二高频细节子图像进行融合，得到所述预设数量的第三高频细节子图像；
[0074] 低频轮廓子图像融合模块，配置为：基于区域能量选择与加权的融合规则，将所述第一低频轮廓子图像和所述第二低频轮廓子图像进行融合，得到第三低频轮廓子图像；
[0075] 融合图像生成模块，配置为：将所述预设数量的所述第三高频细节子图像与所述第三低频轮廓子图像进行W-BEMD逆变换，得到融合图像。
[0076] 优选地，所述待融合图像分解模块包括：
[0077] 残差分量初始化单元，配置为：将所述红外图像和所述可见光图像的残差分量分别设置为RIR,j-1＝XIR、RVI,j-1＝XVI；
[0078] 其中，RIR,j-1和RVI,j-1分别为所述红外图像的第j-1个残差分量和所述可见光图像的残差分量，j-1表示残差分量的序号，且j＝1；为残差分量的序号，XIR和XVI分别为预设的所述红外图像的残差分量的初始值和预设的所述可见光图像的残差分量的初始值；
[0079] 分解单元，配置为：根据所述红外图像的残差分量RIR,j-1对所述红外图像进行BEMD分解得到所述红外图像的内蕴模函数IMFIR,j和残差分量RIR,j；根据所述可见光图像的残差分量RVI,j-1对所述可见光图像进行BEMD分解得到所述可见光图像的内蕴模函数IMFVI,j和残差分量RVI,j；
[0080] 低频轮廓图像计算单元：配置为：分别对所述红外图像的残差分量RIR,j和所述可见光图像的残差分量RVI,j进行W变换，并对变换所得低频系数进行W逆变换，得到残差分量四分之一大小的低频轮廓图像 和
[0081] 高频细节子图像计算单元，配置为：根据所述红外图像的残差分量RIR,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFIR,j求和，计算出第j个所述第一高频细节子图像，并更新残差分量RIR,j；根据所述可见光图像的残差分量RVI,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFVI,j求和
计算出第j个所述第二高频细节子图像，并更新残差分量RVI,j；
[0082] 控制单元，配置为：计算j＝j+1；若j≤L，则继续调用所述低频轮廓图像计算单元与所述高频细节子图像计算单元；否则，将更新后的残差分量RIR,L与RVI,L分别作为所述第一低频轮廓子图像WRIR,L与所述第二低频轮廓子图像WRVI,L；其中，L为所述预设数量。
[0083] 优选地，所述高频细节子图像计算单元包括：
[0084] 第一插值膨胀子单元，配置为：对所述红外图像的低频轮廓图像 进行插值膨胀，得到与RIR,j具有相同大小的图像
[0085] 第一高频细节子图像计算子单元，配置为：按照下式计算残差图像RIR,j的高频细节图像：
[0086]
[0087] 按照下式计算第j个所述第一高频细节子图像：
[0088]
[0089] 并更新残差图像
[0090] 第二插值膨胀子单元，配置为：对所述可见光图像的低频轮廓图像 进行插值膨胀，得到与RVI,j具有相同大小的图像
[0091] 第二高频细节子图像计算子单元，配置为：按照下式计算残差图像RVI,j的高频细节图像：
[0092]
[0093] 按照下式计算第j个所述第二高频细节子图像：
[0094]
[0095] 并更新残差图像
[0096] 优选地，所述高频细节子图像融合模块包括：
[0097] 局部区域方差计算单元，配置为：分别计算所述第一高频细节子图像与所述第二高频细节子图像在各像素点(x,y)处的局部区域方差：
[0098]
[0099]
[0100] 其中，WIMFIR,j和WIMFVI,j分别为第j个所述第一高频细节子图像和第j个所述第二高频细节子图像；j＝1,2,...,L，L为所述预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽； 和 分别代表图像WIMFIR,j和WIMFVI,j在以(x,y)为中心的M×N邻域内的灰度平均值；ω(m,n)代表预设的第一权值矩阵；
[0101] 局部区域方差匹配度计算单元，配置为：计算所述第一高频细节子图像WIMFIR,j和所述第二高频细节子图像WIMFVI,j在各像素点(x,y)处的局部区域方差匹配度：
[0102]
[0103] 其中，
[0104]
[0105]
[0106] 为WIMFIR,j和WIMFVI,j在(x,y)处的相关局部区域方差；
[0107] 第三高频细节子图像计算单元，配置为：若 则根据下式计算所述第三高频细节子图像：
[0108]
[0109] 否则，根据下式计算所述第三高频细节子图像：
[0110]
[0111] 其中，
[0112]
[0113] β(x,y)＝1-α(x,y)
[0114] α(x,y)和β(x,y)均为加权系数，T1为预设的第一匹配度阈值。
[0115] 优选地，所述低频轮廓子图像融合模块包括：
[0116] 局部区域能量计算单元，配置为：分别计算所述第一低频轮廓子图像与所述第二低频轮廓子图像在各像素点(x,y)处的局部区域能量：
[0117]
[0118]
[0119] 其中，WRIR,L和WRVI,L分别为所述第一低频轮廓子图像和所述第二低频轮廓子图像；L为所述预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽；ω′(m,n)代表预设的第二权值矩阵；
[0120] 局部区域能量匹配度计算单元，配置为：计算所述第一低频轮廓子图像WRIR,L和所述第二低频轮廓子图像WRVI,L在各像素点(x,y)处的局部区域能量匹配度：
[0121]
[0122] 其中，
[0123]
[0124]
[0125] 为WRIR,L和WRVI,L在(x,y)处的相关局部区域能量；
[0126] 和 分别代表图像WRIR,L和WRVI,L在以(x,y)为中心的M×N邻域内的灰度平均值；
[0127] 第三低频轮廓子图像计算单元，配置为：若 则根据下式计算所述第三低频轮廓子图像：
[0128]
[0129] 否则，根据下式计算所述第三低频轮廓子图像：
[0130]
[0131] 其中，
[0132]
[0133] β′(x,y)＝1-α′(x,y)
[0134] α′(x,y)和β′(x,y)均为加权系数，T2为预设的第二匹配度阈值。
[0135] 与最接近的现有技术相比，本发明具有如下有益效果：
[0136] 本发明提出的红外图像与可见光图像的融合方法，基于BEMD分解和W变换两种技术进行图像融合，该方法通过BEMD分解和W变换对待融合的源图像进行多尺度分解，然后根据高低频各分量的特性选择合适的融合规则进行融合。该方法保持了W变换和经验模态分解的优点，具有算法冗余度低，运算速度快的优点，且保留了经验模态分解的自适应分析能力和W变换的捕获高频奇异信息的能力，使得采用该方法所得融合结果图像更加清晰、与源图像相关程度较高。解决了现有技术中在图像融合时易出现的模态混叠现象，以及融合结果中高频细节信息丢失、图像对比度降低等问题。
附图说明
[0137] 图1是本发明的红外图像与可见光图像的融合方法实施例的主要步骤示意图；
[0138] 图2是本发明红外图像与可见光图像的融合系统实施例的主要构成示意图。
具体实施方式
[0139] 下面参照附图来描述本发明的优选实施方式。本领域技术人员应当理解的是，这些实施方式仅用于解释本发明的技术原理，并非旨在限制本发明的保护范围。
[0140] 基于BEMD优良的自适应分析能力和W变换捕获高频奇异信息的能力，本发明首先对待融合的源图像(红外图像和可见光图像)分别进行一次BEMD分解，得到图像的内蕴模函数分量和残差分量；然后，采用W变换对所得残差分量进行低频提取，得到原残差分量四分之一大小的低频轮廓信息图像，将其作为新的残差分量(低频成分，记做WR)。将原残差分量与新的残差分量WR做差，得到残差剩余细节图像。考虑到内蕴模函数分量与残差剩余细节图像都具有高频性质，将二者的和作为新的IMF分量(高频成分，记做WIMF)，至此，完成了图像的一次W-BEMD分解；对新的残差分量WR重复上述操作，可得到所需的多层WIMFs分量和最终残差分量；接着，对源图像的各对应WIMFs分量分别采用基于局部区域方差及其匹配度的融合规则进行融合，而对WR分量则采用基于局部区域能量选择和加权的融合规则，得到融合图像的BEMD分解；最后，通过W-BEMD逆变换得到最终融合图像。
[0141] 本实施例中，由于所得WIMFs分量和WR分量具有不同的图像特征，故针对W-BEMD分解下的WIMFs分量和WR分量分别设计不同的融合规则，使最终图像融合的效果达到最优。
[0142] 图1是本发明的红外图像与可见光图像的融合方法实施例的主要步骤示意图。如图1所示，本实施例的融合方法包括步骤S10-S40：
[0143] 步骤S10，基于W-BEMD分解算法，将红外图像分解成预设数量的第一高频细节子图像和一个第一低频轮廓子图像，并将可见光图像分解成预设数量的第二高频细节子图像和一个第二低频轮廓子图像。
[0144] 具体包括步骤S11-S15：
[0145] 步骤S11，将红外图像和可见光图像的残差分量分别设置为RIR,j-1＝XIR、RVI,j-1＝XVI。
[0146] 其中，RIR,j-1和RVI,j-1分别为红外图像的第j-1个残差分量和可见光图像的残差分量，j-1表示残差分量的序号，且j＝1；为残差分量的序号，XIR和XVI分别为预设的红外图像的残差分量的初始值和预设的可见光图像的残差分量的初始值。
[0147] 步骤S12，根据红外图像的残差分量RIR,j-1对红外图像进行BEMD分解得到红外图像的内蕴模函数IMFIR,j和残差分量RIR,j；根据可见光图像的残差分量RVI,j-1对可见光图像进行BEMD分解得到可见光图像的内蕴模函数IMFVI,j和残差分量RVI,j。
[0148] 步骤S13：分别对红外图像的残差分量RIR,j和可见光图像的残差分量RVI,j进行W变换，并对变换所得低频系数进行W逆变换，得到残差分量四分之一大小的低频轮廓图像和
[0149] 步骤S14，根据红外图像的残差分量RIR,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFIR,j求和作为新的内蕴模函数分量WIMFIR,j，即计算出第j个第一高频细节子图像，并更新残差分量RIR,j；根据可见光图像的残差分量RVI,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFVI,j求和作为新的内蕴模函数分量WIMFVI,j，即计算出第j个第二高频细节子图像，并更新残差分量RVI,j。可以具体包括步骤S141-S146：
[0150] 步骤S141，对红外图像的低频轮廓图像RIloR,wj进行插值膨胀，得到与RIR,j具有相同大小的图像
[0151] 步骤S142，按照公式(1)计算残差图像RIR,j的高频细节图像：
[0152]
[0153] 步骤S143，按照公式(2)计算第j个第一高频细节子图像：
[0154]
[0155] 并更新残差图像
[0156] 步骤S144，对可见光图像的低频轮廓图像RlVoI,wj进行插值膨胀，得到与RVI,j具有相同大小的图像
[0157] 步骤S145，按照公式(3)计算残差图像RVI,j的高频细节图像：
[0158]
[0159] 步骤S146，按照公式(4)计算第j个第二高频细节子图像：
[0160]
[0161] 并更新残差图像
[0162] 步骤S15，j＝j+1；若j≤L，转至步骤S13；否则，将更新后的残差分量RIR,L与RVI,L分别作为第一低频轮廓子图像WRIR,L与第二低频轮廓子图像WRVI,L；其中，L为预设数量。
[0163] 步骤S20，基于局部区域方差选择与加权的融合规则，将预设数量的第一高频细节子图像与预设数量的第二高频细节子图像进行融合，得到预设数量的第三高频细节子图像。
[0164] 图像的WIMF分量往往保留了源图像的高频细节信息，如边缘、纹理、线条等，因此合理选择高频融合规则对保持源图像的细节特征具有重要的意义。区域方差作为一种像素点邻域特征的描述指标，表征了区域内像素点之间的偏离程度，该值越大。说明区域内像素点周围变化程度越大，区域内所含信息量越多。据此，所提图像融合方法采用基于局部区域方差选择与加权的融合规则来获得融合系数，从而对WIMF分量进行融合。这种方法符合单像素点变化不易被察觉、而局部区域变化更易被人眼感知的特性。
[0165] 步骤S20可以具体包括步骤S21-S23：
[0166] 步骤S21，分别根据公式(5)、(6)计算第一高频细节子图像与第二高频细节子图像在各像素点(x,y)处的局部区域方差：
[0167]
[0168]
[0169] 其中，WIMFIR,j和WIMFVI,j分别为第j个第一高频细节子图像和第j个第二高频细节子图像；j＝1,2,...,L，L为预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽； 和 分别代表图像WIMFIR,j和WIMFVI,j在以(x,y)为中心的M×N邻域内的灰度平均值；ω(m,n)代表预设的第一权值矩阵，本实施例中选取高斯分布所生成的矩阵作为权值矩阵。
[0170] 步骤S22，计算第一高频细节子图像WIMFIR,j和第二高频细节子图像WIMFVI,j在各像素点(x,y)处的局部区域方差匹配度，如公式(7)所示：
[0171]
[0172] 其中，
[0173]
[0174]
[0175] 为WIMFIR,j和WIMFVI,j在(x,y)处的相关局部区域方差；
[0176] 步骤S23，若 则根据公式(10)计算第三高频细节子图像：
[0177]
[0178] 否则，根据公式(11)计算第三高频细节子图像：
[0179]
[0180] 其中，
[0181]
[0182] β(x,y)＝1-α(x,y)(13)
[0183] α(x,y)和β(x,y)均为加权系数，T1为预设的第一匹配度阈值。
[0184] 步骤S30，基于区域能量选择与加权的融合规则，将第一低频轮廓子图像和第二低频轮廓子图像进行融合，得到第三低频轮廓子图像。
[0185] 残差分量WR是源图像的一种近似描述。为了在融合结果中更多地保留红外图像的目标信息和可见光图像丰富的场景信息，本文对WR分量采用基于区域能量选择与加权的融合规则。
[0186] 步骤S30可以具体包括步骤S31-S33：
[0187] 步骤S31，分别根据公式(14)、(15)计算第一低频轮廓子图像与第二低频轮廓子图像在各像素点(x,y)处的局部区域能量：
[0188]
[0189]
[0190] 其中，WRIR,L和WRVI,L分别为第一低频轮廓子图像和第二低频轮廓子图像；L为预设数量；M和N分别为以点(x,y)为中心的局部邻域窗口的长和宽；ω′(m,n)代表预设的第二权值矩阵；
[0191] 步骤S32，计算第一低频轮廓子图像WRIR,L和第二低频轮廓子图像WRVI,L在各像素点(x,y)处的局部区域能量匹配度，如公式(16)所示：
[0192]
[0193] 其中，
[0194]
[0195]
[0196] 为WRIR,L和WRVI,L在(x,y)处的相关局部区域能量；和 分别代表图像WRIR,L和WRVI,L在以(x,y)为中心的M×N邻域内的灰度
平均值；
[0197] 步骤S33，若 则根据公式(19)计算第三低频轮廓子图像：
[0198]
[0199] 否则，根据公式(20)计算第三低频轮廓子图像：
[0200]
[0201] 其中，
[0202]
[0203] β′(x,y)＝1-α′(x,y) (22)
[0204] α′(x,y)和β′(x,y)均为加权系数，T2为预设的第二匹配度阈值。
[0205] 步骤S40，将预设数量的第三高频细节子图像与第三低频轮廓子图像进行W-BEMD逆变换，得到融合图像。
[0206] 上述实施例中虽然将各个步骤按照上述先后次序的方式进行了描述，但是本领域技术人员可以理解，为了实现本实施例的效果，不同的步骤之间不必按照这样的次序执行，其可以同时(并行)执行或以颠倒的次序执行，这些简单的变化都在本发明的保护范围之内。
[0207] 基于与方法实施例相同的技术构思，本发明还提供了一种红外图像与可见光图像的融合系统，下面进行具体说明。
[0208] 图2是本发明红外图像与可见光图像的融合系统实施例的主要构成示意图。如图2所示，本实施例的红外图像与可见光图像的融合系统1包括：待融合图像分解模块10、高频细节子图像融合模块20、低频轮廓子图像融合模块30和融合图像生成模块40。
[0209] 其中，待融合图像分解模块10配置为：基于W-BEMD分解算法，将红外图像分解成预设数量的第一高频细节子图像和一个第一低频轮廓子图像，并将可见光图像分解成预设数量的第二高频细节子图像和一个第二低频轮廓子图像；高频细节子图像融合模块20配置为：基于局部区域方差选择与加权的融合规则，将预设数量的第一高频细节子图像与预设数量的第二高频细节子图像进行融合，得到预设数量的第三高频细节子图像；低频轮廓子图像融合模块30配置为：基于区域能量选择与加权的融合规则，将第一低频轮廓子图像和第二低频轮廓子图像进行融合，得到第三低频轮廓子图像；融合图像生成模块40配置为：将预设数量的第三高频细节子图像与第三低频轮廓子图像进行W-BEMD逆变换，得到融合图像。
[0210] 本实施例中，待融合图像分解模块10包括：残差分量初始化单元11、分解单元12、低频轮廓图像计算单元13、高频细节子图像计算单元14和控制单元15。
[0211] 其中，残差分量初始化单元11配置为：将红外图像和可见光图像的残差分量分别设置为RIR,j-1＝XIR、RVI,j-1＝XVI；其中，RIR,j-1和RVI,j-1分别为红外图像的第j-1个残差分量和可见光图像的残差分量，j-1表示残差分量的序号，且j＝1；为残差分量的序号，XIR和XVI分别为预设的红外图像的残差分量的初始值和预设的可见光图像的残差分量的初始值；分解单元12配置为：根据红外图像的残差分量RIR,j-1对红外图像进行BEMD分解得到红外图像的内蕴模函数IMFIR,j和残差分量RIR,j；根据可见光图像的残差分量RVI,j-1对可见光图像进行BEMD分解得到可见光图像的内蕴模函数IMFVI,j和残差分量RVI,j；低频轮廓图像计算单元13配置为：分别对红外图像的残差分量RIR,j和可见光图像的残差分量RVI,j进行W变换，并对变换所得低频系数进行W逆变换，得到残差分量四分之一大小的低频轮廓图像 和 高频细节子图像计算单元14配置为：根据红外图像的残差分量RIR,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFIR,j求和，计算出第j个第一高频细节子图像，并更新残差分量RIR,j；根据可见光图像的残差分量RVI,j、低频轮廓图像 计算出残差分量的高频细节图像 将 与内蕴模函数IMFVI,j求和计算出第j个第二高频细节子图像，并更新残差分量RVI,j；控制单元15配置为：计算j＝j+1；若j≤L，则继续调用低频轮廓图像计算单元与高频细节子图像计算单元；否则，将更新后的残差分量RIR,L与RVI,L分别作为第一低频轮廓子图像WRIR,L与第二低频轮廓子图像WRVI,L；其中，L为预设数量。
[0212] 本实施例中，高频细节子图像计算单元14包括：第一插值膨胀子单元、第一高频细节子图像计算子单元、第二插值膨胀子单元和第二高频细节子图像计算子单元。
[0213] 其中，第一插值膨胀子单元配置为：对红外图像的低频轮廓图像 进行插值膨胀，得到与RIR,j具有相同大小的图像 第一高频细节子图像计算子单元配置为：按照公式(1)计算残差图像RIR,j的高频细节图像，按照公式(2)计算第j个第一高频细节子图像，并更新残差图像 第二插值膨胀子单元配置为：对可见光图像的低频轮廓图像进行插值膨胀，得到与RVI,j具有相同大小的图像 第二高频细节子图像计算子单元配置为：按照公式(3)计算残差图像RVI,j的高频细节图像，按照公式(4)计算第j个第二高频细节子图像，并更新残差图像
[0214] 本实施例中，高频细节子图像融合模块20包括：局部区域方差计算单元21、局部区域方差匹配度计算单元22和第三高频细节子图像计算单元23。
[0215] 其中，局部区域方差计算单元21配置为：分别根据公式(5)、(6)计算第一高频细节子图像与第二高频细节子图像在各像素点(x,y)处的局部区域方差；局部区域方差匹配度计算单元22配置为：根据公式(7)-(9)计算第一高频细节子图像WIMFIR,j和第二高频细节子图像WIMFVI,j在各像素点(x,y)处的局部区域方差匹配度；第三高频细节子图像计算单元23配置为：若 则根据公式(10)计算第三高频细节子图像：
[0216] 否则，根据公式(11)-(13)计算第三高频细节子图像；
[0217] 本实施例中，低频轮廓子图像融合模块30包括：局部区域能量计算单元31、局部区域能量匹配度计算单元32和第三低频轮廓子图像计算单元33。
[0218] 其中，局部区域能量计算单元31配置为：分别根据公式(14)、(15)计算第一低频轮廓子图像与第二低频轮廓子图像在各像素点(x,y)处的局部区域能量；局部区域能量匹配度计算单元32配置为：根据公式(16)-(18)计算第一低频轮廓子图像WRIR,L和第二低频轮廓子图像WRVI,L在各像素点(x,y)处的局部区域能量匹配度；第三低频轮廓子图像计算单元33配置为：若 则根据公式(19)计算第三低频轮廓子图像，否则，根据公式(20)-(22)计算第三低频轮廓子图像。
[0219] 本领域技术人员应该能够意识到，结合本文中所公开的实施例描述的各示例的方法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明电子硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以电子硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。本领域技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
[0220] 至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征做出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。
