基于增补机制和PCNN的NSCT域图像融合方法与流程

本发明涉及一种nsct域图像融合方法，特别涉及一种基于增补机制和pcnn的nsct域图像融合方法。

背景技术：

由于成像机理和技术限制，单一成像传感器获取的图像并不能反映被观测对象的所有特征，因此需要对不同传感器图像中的有用信息进行提取，并融合成一副具备更完整信息，并有助于人类观察和处理的图像。

图像融合技术可以综合两种不同类型图像的互补信息，它不是简单地对不同类型图像进行叠加处理，通过融合得到的图像具备两种图像各自的优点，克服了单一传感器在应用环境、使用范围以及特定目标获取上存在的局限性，同时还能够提高图像的空间分辨率和清晰度，便于图像的理解和识别，有效地提高图像数据的利用率。图像融合技术不仅能够去除多传感器图像数据间的冗余信息，减少传输量，减轻观察者的观测压力，还能够获得关于同一场景或目标全面而精确的地理空间信息，因此在军事或民用领域的应用中，图像融合具有非常重要的研究价值。

近年来，针对图像的融合，已经有大量的图像融合算法被相继提出，其中基于多尺度分析的图像融合算法从多尺度金字塔变换开始，经由小波变换，ridgelet变换，curvelet变换等一系列算法发展，已取得了良好的融合效果。2006年，minh.n.do等人提出非下采样轮廓波变换(nonsubsampledcontourlet，nsct)，nsct是在contourlet变换的基础上引入非下采样思想，保证了分解过程中的平移不变性，解决了传统算法在图像分解时的信号频谱混叠问题，获得了渐进最优的图像分析结果。在此基础上，国内外学者们相继提出了不少行之有效的新算法，但大多数融合图像纹理细节缺失、失真和目标显著性等问题难点仍然没有得到完全解决。

文献“基于补偿机制的nsct域红外与可见光图像融合，仪器仪表学报,2016,37(4):861-870”公开了一种基于补偿机制的nsct域红外与可见光图像融合方法，该方法采用nsct变换对红外与可见光图像分别进行分解，得到低频子带和高频子带，对低频子带进行小波变换，利用合适的融合规则对对小波变换的子带进行融合，将其变换得到的子带重构作为nsct重构的低频子带，接着将nsct变换的高频子带的最高层利用绝对值取大方法，其它层利用pcnn方法进行融合，然后将得到的低频子带和高频子带进行nsct重构获得融合后的图像。文献所述方法得到融合后图像对比度较低，可将光图像的背景细节信息丢失，无法体现源图像的纹理信息，并且融合图像部分区域存在一定的扭曲与失真现象。

技术实现要素：

为了克服现有nsct域图像融合方法存在图像失真现象的不足，本发明提供一种基于增补机制和pcnn的nsct域图像融合方法。该方法采用增补的小波变换对nsct分解的低频子带进行融合处理，以便尽可能多的保留图像背景的细节信息；利用改进的高斯加权sml的方法进行融合来增强图像细节；利用边缘梯度信息激励pcnn的方法进行融合来增强图像边缘信息；实验证明本发明方法与现有的图像融合方法相比具有更好的融合效果，在对克服图像的纹理细节缺失、失真不足的同时，目标显著性得到很大的提高，进一步提高了图像质量。

本发明解决其技术问题所采用的技术方案：一种基于增补机制和pcnn的nsct域图像融合方法，其特点是包括以下步骤：

步骤一、对已配准的待融合图像a和b分别进行nsct分解，得到各自的nsct系数和其中为a图像第j尺度下第k个方向的nsct高频子带系数，l^a为a图像的nsct低频系数，为b图像第j尺度下第k个方向的nsct高频子带系数，l^b为b图像的nsct低频系数，j表示nsct分解尺度，j＝1，2，λ，j，其中j为最大分解尺度，k表示分解尺度j下的分解方向数，满足k＝2^j。

步骤二、按融合规则对步骤一得到的图像a和b各自的nsct系数和进行融合，得到融合后的nsct系数，融合规则为：

a.nsct低频子带系数融合策略。

对nsct分解的低频子带进行增补机制的小波变换分解的低频部分采用局部加权区域能量进行融合。计算以像素点(i,j)为中心，大小为3×3窗口区域局部加权能量定义如下式：

式中sl(i,j)表示低频图像点(i,j)的灰度，ωl(x,y)为局部区域窗口，表示卷积，l表示待融合图像，步骤一中a＝vi,b＝ir,vi是可见光图像，ir是红外图像。图像中的目标物体在红外图像中的局部邻域能量比其在可见光图像中的局部邻域能量大，低频子带融合规则如下式所示：

式中，

式中，分别是第l层第k个方向对应低频子带融合系数、红外图像和可见光图像的低频子带系数，分别是第l层第k个方向对应红外图像和可见光图像的局部区域加权能量，分别是第l层第k个方向对应红外图像和可见光图像的局部区域加权能量的平均值。ω1为权重系数，为图像信息熵，m为图像总的灰度级数，pi＝mi/m为图像中像素灰度值为i的概率，k1，k2为图像权重系数，取k1＝5，k2＝1.6。

由于拉普拉斯分量绝对和sml能够较好的反映图像边缘、纹理等细节信息并且能表示图像的锐化度和清晰度，利用改进的高斯加权sml对增补小波变换分解后的三个方向子带系数分别进行融合。

设经过二维小波变换分解后第l层第k个方向的方向子带系数为h^l,k，sml定义如下：

ml^l,k(i,j)＝-h^l,k(i-1,j-1)-4h^l,k(i-1,j)-h^l,k(i-1,j+1)

-4h^l,k(i,j-1)+20h^l,k(i,j)-4h^l,k(i,j+1)

-h^l,k(i+1,j-1)-4h^l,k(i+1,j)-h^l,k(i+1,j+1)(8)

式中采用3×3窗口，采用式(7)计算的sml认为窗口内所有点对sml的贡献一致，无法突出中心的的重要性。采用高斯函数对ml加权求和，高斯函数满足强调中心像素以及距离中心像素越近的像素权值越大的要求，归一化后权值之和为1。高斯函数计算公式如下式所示：

其中σ表示方差。

改进的高斯加权sml如下所示：

其中，权值函数是按式(9)取值，表示卷积。

引入图像区域锐化度比率，定义如下：

其中,

式中：j＝vi,ir，m×n表示图像的区域大小，rrs^vi,ir表示待融合图像vi和ir的区域锐化度比率，mean^j表示图像的均值，std^j表示图像的方差。

结合不同谱段源图像的方向子带细节信息，利用加权组合的方式融合方向子带系数。设源图像的方向子带系数为则融合后的系数为及权值分别如式(14)、式(15)

b.nsct高频子带系数融合策略。

对nsct高频子带最高层融合：

在图像最高层尺度n上采用改进的高斯加权sml。

对nsct高频子带其它层融合：

对nsct高频子带其它层采用改进的脉冲神经网络pcnn的图像融合规则。设改进的pcnn的迭代次数为n，随着迭代次数的增加，各pcnn中外部刺激输入较大的神经元点火次数将会明显增加。完成n次迭代后统计各nsct变换分解子图像每个像素的点火次数，根据点火次数的大小来选择融合图像系数，然后经过nsct逆变换得到融合图像。而人眼对边缘信息更敏感，边缘梯度信息能更好地代表图像特征。边缘梯度信息为

式中，d是以位置(i,j)为中心的周围邻域窗为3×3，对每个子带以滑动窗口为单位计算边缘梯度系数。c(i,j)是像素经过nsct变换后对应点(i,j)的系数。

在pcnn中，其数学表达式如下：

fij(n)＝sij(17)

uij(n)＝fij(n)(1+βlij(n))(19)

tij＝tij(n-1)+yij(n)(22)

式中：fij表示第(i,j)个神经元的反馈输入，lij表示第(i,j)个神经元的连接输入项，β为突触之间的连接强度，θij为第(i,j)个神经元的动态门限阈值，uij为第(i,j)个神经元的内部活动项。yij为第(i,j)个神经元的脉冲输出，tij为迭代之后第(i,j)个神经元的点火总次数，yij,pq为第(i,j)个神经元的邻域神经元(p,q)的脉冲输出，n为迭代次数，sij为第(i,j)个神经元的外部输入刺激信号，wij,pq为第(i,j)个神经元的连接矩阵中第p行q列元素，αl、αθ分别为连接输入衰减系数和阈值衰减系数，vl、vθ分别为连接幅值和阈值幅值。

由式(17)～(22)得出，改进型pcnn模型有4个参数β、αθ、vθ、w。其中w表示的是两个像素间距离平方的导数，神经元ij和神经元pq的连接权重为：

改进的pcnn图像融合过程如下：

[1].利用nsct变换对红外和可见光图像进行分解，得到尺度n以外的其它层带通子带的系数和

[2].用式(16)计算边缘梯度信息，得到图像系数的边缘梯度信息特征系数和用和激励pcnn，根据式(17)～(22)流程分别得到点火次数矩阵和并采取以下融合规则对系数和进行融合处理。

步骤三、图像重构。

对融合图像的各子带系数进行nsct逆变换，得到最终的融合图像。

本发明的有益效果是：该方法采用增补的小波变换对nsct分解的低频子带进行融合处理，以便尽可能多的保留图像背景的细节信息；利用改进的高斯加权sml的方法进行融合来增强图像细节；利用边缘梯度信息激励pcnn的方法进行融合来增强图像边缘信息；实验证明本发明方法与现有的图像融合方法相比具有更好的融合效果，在对克服图像的纹理细节缺失、失真不足的同时，目标显著性得到很大的提高，进一步提高了图像质量。

下面结合附图和具体实施方式对本发明作详细说明。

附图说明

图1是本发明基于增补机制和pcnn的nsct域图像融合方法的流程图。

图2是本发明方法构建的pcnn神经元模型图。

图3是本发明方法实施例中实验图像数据集图。

图4是红外光与可见光图像集融合后图像的对比图。

具体实施方式

参照图1-4。本发明基于增补机制和pcnn的nsct域图像融合方法具体步骤如下：

用于实施的硬件环境是：实验环境为cpuintelcorei5-5200u@2.20ghz，内存为4gb，采用matlabr2014a编程。本发明采用“uncamp”图像集(320×240)、“octec”图像集(640×480)、“quad”图像集(256×256)和“seascape”图像集(256×256)4组红外与可见光图像集。

步骤一：对已配准的待融合图像a和b分别进行nsct分解，得到各自的nsct系数和其中为a图像第j尺度下第k个方向的nsct高频子带系数，l^a为a图像的nsct低频系数，为b图像第j尺度下第k个方向的nsct高频子带系数，l^b为b图像的nsct低频系数，所述的j表示nsct分解尺度，j＝1，2，λ，j，j为最大分解尺度，k表示分解尺度j下的分解方向数，满足k＝2^j。

步骤二：按融合规则对步骤1得到的图像a和b各自的nsct系数和进行融合，得到融合后的nsct系数。

融合规则为：

a.nsct低频子带系数融合策略。

图像经nsct分解后得到的低频子带是源图像的近似描述，包含了图像中的大部分能力特征。本发明采用基于增补小波的nsct低频融合规则，在一定程度上能更好的表征图像的边缘纹理等细节信息。

对nsct分解的低频子带进行增补机制的小波变换分解，得到1个低频子带和3个方向子带。对低频子带采用局部加权区域能量进行融合，对3个方向子带分别采用改进的高斯加权sml对三个方向子带系数进行融合，然后进行小波重构得到用于nsct重构的低频子带。

对增补机制的小波变换分解的低频部分采用局部加权区域能量进行融合。计算以像素点(i,j)为中心，大小为3×3窗口区域局部加权能量定义如下式(1)：

式中sl(i,j)表示低频图像点(i,j)的灰度，ωl(x,y)为局部区域窗口，表示卷积，l表示待融合图像，以红外和可见光为例，步骤1中a＝vi,b＝ir,vi是可见光图像，ir是红外图像。由于图像中的目标物体在红外图像中的局部邻域能量比其在可见光图像中的局部邻域能量大，因此为了更好地突显出红外图像中的目标物体与可见光图像中的背景信息，本发明的低频子带融合规则如下式所示：

式中：

式中：分别是第l层第k个方向对应低频子带融合系数、红外图像和可见光图像的低频子带系数，分别是第l层第k个方向对应红外图像和可见光图像的局部区域加权能量，分别是第l层第k个方向对应红外图像和可见光图像的局部区域加权能量的平均值。ω1为权重系数，为图像信息熵，m为图像总的灰度级数，pi＝mi/m为图像中像素灰度值为i的概率，k1，k2为图像权重系数，为了突出融合后图像的质量效果，本发明中k1＝5，k2＝1.6。经过二维小波变换分解后的三个方向子带系数表示源图像的细节信息。

由于拉普拉斯分量绝对和sml能够较好的反映图像边缘、纹理等细节信息并且能表示图像的锐化度和清晰度，本发明采用改进的高斯加权sml对三个方向子带系数分别进行融合。

设经过二维小波变换分解后第l层第k个方向的方向子带系数为h^l,k，sml定义如下

采用3×3窗口，采用式(7)计算的sml认为窗口内所有点对sml的贡献一致，无法突出中心的的重要性。本发明采用高斯函数对ml加权求和，高斯函数满足强调中心像素以及距离中心像素越近的像素权值越大的要求，归一化后权值之和为1。高斯函数计算公式如下式所示：

其中σ表示方差。

改进的高斯加权sml如下所示：

其中，权值函数是按式(9)取值，表示卷积。

为了尽可能多的突出图像的细节信息和图像的清晰度，使融合后图像的质量效果达到最佳，本发明引入图像区域锐化度比率，定义如下：

其中：

式中：j＝vi,ir，m×n表示图像的区域大小，rrs^vi,ir表示待融合图像vi和ir的区域锐化度比率，mean^j表示图像的均值，std^j表示图像的方差。

结合不同谱段源图像的方向子带细节信息，本发明采用加权组合的方式融合方向子带系数。设源图像的方向子带系数为则融合后的系数为及权值分别如式(14)、(15)

b.nsct高频子带系数融合策略

在nsct分解所得的高频子带中，高频系数的值表示图像细节的变化情况，直接反映图像中边缘信息和纹理细节分布。为了更多的反映图像边缘、纹理等细节信息表示的图像锐化度和清晰度，对nsct分解所得的高频子带分为其它层和最高层融合。

对nsct高频子带最高层融合。

在图像最高层尺度n上采用改进的高斯加权sml，这种方法与低频部分的小波变换分解的三个方向子带的融合方法相同，这里不再叙述。

对nsct高频子带其它层融合

为了更好的提取待融合图像的细节等信息，提高融合图像的质量，本发明在除最高层n以外的其它层带通子带，采用改进的脉冲神经网络pcnn的图像融合规则。在改进的pcnn图像融合算法中，点火次数较大的神经元所对应的像素点对应于图像中的显著特征，因而可以根据神经元的点火次数来选择融合图像的nsct系数。设pcnn的迭代次数为n，随着迭代次数的增加，各pcnn中外部刺激输入较大的神经元点火次数将会明显增加。完成n次迭代后统计各nsct变换分解子图像每个像素的点火次数，根据点火次数的大小来选择融合图像系数，然后经过nsct逆变换得到融合图像。而人眼对边缘信息更敏感，边缘梯度信息能更好地代表图像特征。边缘梯度信息为

式中，d是以位置(i,j)为中心的周围邻域窗为3×3，对每个子带以滑动窗口为单位计算边缘梯度系数。c(i,j)是像素经过nsct变换后对应点(i,j)的系数。

在pcnn中，每个神经元均由接收部分、调制部分和脉冲产生部分3个部分组成，其数学表达式如下：

fij(n)＝sij(17)

uij(n)＝fij(n)(1+βlij(n))(19)

tij＝tij(n-1)+yij(n)(22)

式中：fij表示第(i,j)个神经元的反馈输入，lij表示第(i,j)个神经元的连接输入项，β为突触之间的连接强度，θij为第(i,j)个神经元的动态门限阈值，uij为第(i,j)个神经元的内部活动项。yij为第(i,j)个神经元的脉冲输出，tij为迭代之后第(i,j)个神经元的点火总次数，yij,pq为第(i,j)个神经元的邻域神经元(p,q)的脉冲输出，n为迭代次数，sij为第(i,j)个神经元的外部输入刺激信号，wij,pq为第(i,j)个神经元的连接矩阵中第p行q列元素，αl、αθ分别为连接输入衰减系数和阈值衰减系数，vl、vθ分别为连接幅值和阈值幅值。

由式(17)～(22)可以看出，改进型pcnn模型主要由4个参数β、αθ、vθ、w。其中w表示的是两个像素间距离平方的导数，神经元ij和神经元pq的连接权重为：

改进的pcnn图像融合过程如下：

[1].利用nsct变换对红外和可见光图像进行分解，得到尺度n以外的其它n-1层带通子带的系数和

[2].用式(16)计算边缘梯度信息，得到图像系数的边缘梯度信息特征系数和用和激励pcnn，根据式(17)～(22)流程分别得到点火次数矩阵和并采取以下融合规则对系数和进行融合处理。

步骤三：图像重构。

对融合图像的各子带系数进行nsct逆变换从而得到最终的融合图像。

下面结合附图3和附图4对本发明的效果作进一步的描述。

附图3是4组红外与可见光图像集。其中附图3(a)是“uncamp”图像集(320×240)，附图3(b)是“octec”图像集(640×480)，附图3(c)是“quad”图像集(256×256)，附图3(d)是“seascape”图像集(256×256)。

1.实验条件。

实验环境为cpuintelcorei5-5200u@2.20ghz，内存为4gb，采用matlabr2014a编程。本发明采用“uncamp”图像集(320×240)、“octec”图像集(640×480)、“quad”图像集(256×256)和“seascape”图像集(256×256)4组红外与可见光图像集。

2.实验内容。

附图4是4组红外与可见光图像集融合后图像的对比图。

实验1，用本发明的方法和现有的四种融合方法对附图3中4组红外与可见光图像集进行融合实验，附图4(a)融合结果从左及右依次是，文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753.文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84.文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227.文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.和本发明的图像融合结果图。

从附图4(a)可得，本发明的融合方法与现有的四种融合方法相比，在灰度范围和整体亮度上与源图像较为匹配，视觉效果更高，目标和场景更清晰。文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753、文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84、文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227和文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.的融合方法的融合结果和本发明方法相比对比度较低，可将光图像的背景细节信息丢失，无法体现源图像的纹理信息，并且融合图像部分区域存在一定的扭曲与失真现象。

实验2，用本发明的方法和现有的四种融合方法对附图3中4组红外与可见光图像集进行融合实验，附图4(b)融合结果从左及右依次是，文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753.文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84.文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227.文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.和本发明的图像融合结果图。

从附图4(b)可得，本发明的融合方法与现有的四种融合方法相比，在灰度范围和整体亮度上与源图像较为匹配，视觉效果更高，目标和场景更清晰。文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753、文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84、文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227和文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.的融合方法的融合结果和本发明方法相比对比度较低，可将光图像的背景细节信息丢失，无法体现源图像的纹理信息，并且融合图像部分区域存在一定的扭曲与失真现象。

实验3，用本发明的方法和现有的四种融合方法对附图3中4组红外与可见光图像集进行融合实验，附图4(c)融合结果从左及右依次是，文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753.文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84.文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227.文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.和本发明的图像融合结果图。

从附图4(c)可得，本发明的融合方法与现有的四种融合方法相比，在灰度范围和整体亮度上与源图像较为匹配，视觉效果更高，目标和场景更清晰。文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753、文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84、文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227和文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.的融合方法的融合结果和本发明方法相比对比度较低，可将光图像的背景细节信息丢失，无法体现源图像的纹理信息，并且融合图像部分区域存在一定的扭曲与失真现象。

实验4，用本发明的方法和现有的四种融合方法对附图3中4组红外与可见光图像集进行融合实验，附图4(d)融合结果从左及右依次是，文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753.文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84.文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227.文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.和本发明的图像融合结果图。

从附图4(d)可得，本发明的融合方法与现有的四种融合方法相比，在灰度范围和整体亮度上与源图像较为匹配，视觉效果更高，目标和场景更清晰。文章《一种基于小波变换的图像融合算法》，电子学报，2004,32(5):750-753、文章《imagefusionbasedonanewcontourletpacket》，informationfusion,2010,11(2):78-84、文章《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227和文章《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.的融合方法的融合结果和本发明方法相比对比度较低，可将光图像的背景细节信息丢失，无法体现源图像的纹理信息，并且融合图像部分区域存在一定的扭曲与失真现象。

将本发明的融合方法与文章1《一种基于小波变换的图像融合算法》，电子学报，2004，32(5)：750-753.文章2《imagefusionbasedonanewcontourletpacket》，informationfusion，2010，11(2)：78-84.文章3《high-qualityfusionforvisibleandinfraredimagesbasedonthedoublensct》，ieee7thinternationalcongressonimageandsignalprocessing,2014:223-227.文章4《基于补偿机制的nsct域红外与可见光图像融合》，仪器仪表学报,2016,37(4):861-870.的融合方法在4种图像质量评价指标上进行比较，来客观评价本发明的效果。5种融合方法的融合客观评价指标如下表：

表1“uncamp”图像集融合质量评价指标对比

表2“octec”图像集融合质量评价指标对比

表3“quad”图像集融合质量评价指标对比

表4“seascape”图像集融合质量评价指标对比

ie表示信息熵，ie最大，表明本发明方法融合结果所包含的信息量最丰富；ssim表示信息结构相似度，ssim最大，表明本发明方法融合图像与源图像的相似程度最大，所包含源图像的纹理等细节信息最丰富；sd表示标准差，sd最大，则融合图像灰度分布离散程度最大，对比度更高，且视觉效果更好；psnr表示峰值信噪比，psnr最大，反映了融合图像失真程度最小，更加可以凸显红外目标，保留源图像的更多纹理细节信息。

由以上表1～4分析可得，本发明方法融合图像对比度高，目标信息突出，边缘纹理等细节信息更丰富，从主客观评价看，具有更好的融合性能，融合图像质量达到最佳，符合人眼视觉效果，表明了本发明所提方法的有效性。