CN109345496A - 一种总变分和结构张量的图像融合方法及装置 - Google Patents
一种总变分和结构张量的图像融合方法及装置 Download PDFInfo
- Publication number
- CN109345496A CN109345496A CN201811055297.6A CN201811055297A CN109345496A CN 109345496 A CN109345496 A CN 109345496A CN 201811055297 A CN201811055297 A CN 201811055297A CN 109345496 A CN109345496 A CN 109345496A
- Authority
- CN
- China
- Prior art keywords
- image
- total variation
- gradient
- follows
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000002156 mixing Methods 0.000 claims abstract description 69
- 230000004927 fusion Effects 0.000 claims abstract description 43
- 230000006870 function Effects 0.000 claims description 37
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000000203 mixture Substances 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 8
- 230000008602 contraction Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 230000009021 linear effect Effects 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims 2
- 239000000155 melt Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 15
- 230000007704 transition Effects 0.000 abstract description 5
- 230000008447 perception Effects 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 108010057516 2-hydroxymuconic semialdehyde dehydrogenase Proteins 0.000 description 2
- 102100033983 Serpin-like protein HMSD Human genes 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007500 overflow downdraw method Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 208000027534 Emotional disease Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000005267 amalgamation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 125000001475 halogen functional group Chemical group 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000001757 thermogravimetry curve Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本发明提出一种总变分和结构张量的图像融合方法及装置，总变分TV模型利用图像的结构张量获取融合图像的目标梯度场，在目标梯度场内部利用Soft‑Max函数柔化梯度场的权重，将融合梯度场放入到TV模型的正则项内，保障融合图像的梯度场与目标梯度场接近，利用分裂Bregman算法求解全新的TV模型，获取最终的融合图像，通过采用本方案可以更好地表达图像的特征，将微光的纹理信息与红外图像的显著性特征完美结合，取得较好的融合效果，可以较完整地保留了丰富的纹理信息和细节信息，并且图像过渡自然。
Description
技术领域
本发明涉及图像融合领域，特别涉及一种总变分和结构张量的图像融合方法及装置。
背景技术
图像融合作为信息融合的一门分支，是一门包含信号处理、图像处理、计算机技术、图像探测器技术等多种领域的综合高新技术。图像融合对某一场景获取的多幅图像信息加以提取，最终将两幅或多幅包含多样、冗余信息的图像合成为一幅含有更丰富信息的新图像。融合图像克服了单幅图像在空间分辨率、物理属性、时间分辨率以及光谱信息等方面的局限性，能更好地为重要目标的定位、识别和理解提供支撑。
红外(IR)与可见光(VI)传感器之间的信息融合是目前应用最多，同时也是最具有现实意义的融合方式。通过将可见光的纹理细节信息与红外热辐射目标的显著信息相结合，有利于在隐藏、伪装和迷惑的情况下更快、更精确地探测目标。针对IR与VI图像之间的融合，最为常见的就是加权平均法。该方法只是将源图像之间的亮度进行空间叠加，容易导致融合图像解析度下降，并且丢失很多纹理细节信息。为了能够获得更优良的融合效果，使融合图像符合人眼视觉特性，人们相继提出了新的融合方法。总变分法(TV)是一种新颖的融合方法，其独特之处在于将融合问题转化为图像优化和重建问题。由于其内部参数较少，并且能够完整地保持图像局部特征，因此目前受到了学者们的广泛关注。TV模型在解决融合问题时，往往在正则项内约束融合图像与给定VI图像具有相似的梯度特征，然而包含在IR图像的纹理信息却被忽略，导致融合图像出边缘细节缺失的问题。为此，我们希望可以直接获取融合图像的梯度场，将其放入到正则项内，从而保证源图像内部的纹理梯度信息较为全面的转移到融合图像中。
作为另一个重要的分支，基于梯度域的图像融合方法也得到了的广泛应用。人类视觉对图像局部细节的变化非常敏感，而图像局部细节实际上就是图像梯度。在梯度域内对图像进行处理，可以使融合图像保留输入图像的重要特征。文献《G.Piella,Imagefusion for enhanced visualization:a variational approach[J],Int.J.Comput.Vision 83(1)(2009)1-11.》等人采用结构张量来获取融合图像的目标梯度场，并从目标梯度场中重建融合图像。该方法为本文的研究提供了一定的思路，我们可以将获取的目标梯度场放入到TV模型的正则项内，从而取代原有的VI图像梯度场。通过全新的约束关系，使融合图像的梯度场最大程度接近于目标梯度场，抑制了纹理细节的丢失。除此外，鉴于红外热像图与可见光图像的较大光谱差异性，单一地使用加权法来设计目标梯度场的光谱梯度权重，会导致融合结果出现光晕伪影现象。除此外，针对TV模型的求解问题，学者们一般采用梯度下降法。然而，利用该方法获取的欧拉-拉格朗日方程往往具有非线性和不可微性，并且由于其最小化速度较慢，会影响求解速度。
发明内容
本发明实施例提供了一种总变分和结构张量的图像融合方法及装置，可以弥补两种源图像较大的光谱差异性，同时较完整地保留了丰富的纹理信息和细节信息，并且图像过渡自然。
第一方面，本发明提供一种总变分和结构张量的图像融合方法，所述方法包括：
获取来自同一场景并进行图像配准的红外图像和可见光图像；
利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场；
利用Soft-Max函数柔化所述目标梯度场的权重；
将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近；
利用分裂Bregman算法求解TV模型得到融合图像。
可选地，所述获取来自同一场景并进行图像配准的红外图像和可见光图像，包括：
分别获取同一场景的红外图像和可见光图像；
对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
可选地，所述利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场，包括：
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
可选地，所述利用Soft-Max函数柔化所述目标梯度场的权重，包括：
利用Soft-Max函数构成目标梯度场的梯度权重，表达式如下：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
可选地，所述将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近，包括：
利用结构张量获取的目标梯度场放入到总变分TV模型的正则项内，约束融合图像的融合梯度场与目标梯度场Vf接近，表达式如下：
在数据保真项内加入权重函数，用于表达源图像对融合图像的贡献程度，表达式如下：
其中f代表融合图像，i代表红外图像，v代表可见光图像，α与β是两个非负权重函数，定义权重函数的表达式为：
其中F(x,y)表示两种源图像中像素(x,y)的强度值，p表示像素强度，Mp表示其强度等于p的像素的数量，并且L是灰度级的数量，S(x,y)代表源图像的显著特征矩阵。
可选地，所述利用分裂Bregman算法求解TV模型得到融合图像，包括：
总变分模型的建立
将TV模型的数据保真项与正则项相结合，得到如下表达式：
将TV模型用离散化的l2范数与l1范数再次简化，表示为：
其中λ为正则参数，并且：
总变分模型的求解
采用分裂Bregman算法实现TV模型的最优求解，其求解过程如下：
引入变量d，将公式(8)变为以下形式：
对fk+1可以按照如下方式求解：
对f的最小化方程由以下最优性条件给出：
(f-αi-βv)-Δf+γdiv(Vf+dk-bk)＝0； (12)
根据公式(12)是线性的性质，把公式(12)重写为如下形式：
f-Δf＝αi+βv-γdiv(Vf+dk-bk)； (13)
对公式(13)的进行快速傅立叶FFT形式为：
F(K)F(fk+1)＝F(rhsk)； (14)
K＝I-Δ； (15)
rhk＝αi+βv-γdiv(Vf+dk-bk)； (16)
其中I表示单位矩阵，F表示快速傅里叶变换；
公式(14)中fk+1的闭合形式解为：
dk+1可以按照如下方式求解：
公式(18)的解使用如下的软阈值形式获得：
由Shrink算子得到dk+1的解为：
对每一步迭代中对其做出如下约束：
fk+1＝max[min(fk+1,1),0] (21)。
第二方面，本发明提供一种总变分和结构张量的图像融合装置，所述装置包括：
图像配准单元，用于获取来自同一场景并进行图像配准的红外图像和可见光图像；
目标梯度场获取单元，用于利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场；
权重获取单元，用于利用Soft-Max函数柔化所述目标梯度场的权重；
处理单元，用于将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近；
融合单元，用于利用分裂Bregman算法求解TV模型得到融合图像。
可选地，所述图像配准单元具体用于：
分别获取同一场景的红外图像和可见光图像；
对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
可选地，所述目标梯度场获取单元具体用于：
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
可选地，所述处理单元具体用于：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
从以上技术方案可以看出，本发明实施例具有以下优点：
本发明提出一种总变分和结构张量的图像融合方法及装置，总变分TV模型利用图像的结构张量获取融合图像的目标梯度场，在目标梯度场内部利用Soft-Max函数柔化梯度场的权重，将融合梯度场放入到TV模型的正则项内，保障融合图像的梯度场与目标梯度场接近，利用分裂Bregman算法求解全新的TV模型，获取最终的融合图像，通过采用本方案可以更好地表达图像的特征，将微光的纹理信息与红外图像的显著性特征完美结合，取得较好的融合效果，可以较完整地保留了丰富的纹理信息和细节信息，并且图像过渡自然。
附图说明
图1是本发明实施例中的总变分和结构张量的图像融合方法中一种实施例的流程图；
图2是本发明实施例中的总变分和结构张量的图像融合方法中一种实施例的图像融合示意图；
图3是本发明实施例中的总变分和结构张量的图像融合方法中一种实施例的Soft-Max函数示意图；
图4是本发明实施例中的总变分和结构张量的图像融合方法中一种实施例的融合图像对比示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
结合图1和2所示，本发明提供一种总变分和结构张量的图像融合方法，所述方法包括：
S101、获取来自同一场景并进行图像配准的红外图像和可见光图像。
选取来自同一场景并且经过严格配准的红外图像与可见光图像，对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
S102、利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场。
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
S103、利用Soft-Max函数柔化所述目标梯度场的权重。
结合图3所示，利用Soft-Max函数构成目标梯度场的梯度权重，表达式如下：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
S104、将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近。
利用结构张量获取的目标梯度场放入到总变分TV模型的正则项内，约束融合图像的融合梯度场与目标梯度场Vf接近，表达式如下：
在数据保真项内加入权重函数，用于表达源图像对融合图像的贡献程度，表达式如下：
其中f代表融合图像，i代表红外图像，v代表可见光图像，α与β是两个非负权重函数，定义权重函数的表达式为：
其中F(x,y)表示两种源图像中像素(x,y)的强度值，p表示像素强度，Mp表示其强度等于p的像素的数量，并且L是灰度级的数量，S(x,y)代表源图像的显著特征矩阵。
S105、利用分裂Bregman算法求解TV模型得到融合图像。
总变分模型的建立
将TV模型的数据保真项与正则项相结合，得到如下表达式：
将TV模型用离散化的l2范数与l1范数再次简化，表示为：
其中λ为正则参数，并且：
总变分模型的求解
采用分裂Bregman算法实现TV模型的最优求解，其求解过程如下：
引入变量d，将公式(8)变为以下形式：
对fk+1可以按照如下方式求解：
对f的最小化方程由以下最优性条件给出：
(f-αi-βv)-Δf+γdiv(Vf+dk-bk)＝0； (12)
根据公式(12)是线性的性质，把公式(12)重写为如下形式：
f-Δf＝αi+βv-γdiv(Vf+dk-bk)； (13)
对公式(13)的进行快速傅立叶FFT形式为：
F(K)F(fk+1)＝F(rhsk)； (14)
K＝I-Δ； (15)
rhk＝αi+βv-γdiv(Vf+dk-bk)； (16)
其中I表示单位矩阵，F表示快速傅里叶变换；
公式(14)中fk+1的闭合形式解为：
dk+1可以按照如下方式求解：
公式(18)的解使用如下的软阈值形式获得：
由Shrink算子得到dk+1的解为：
对每一步迭代中对其做出如下约束：
fk+1＝max[min(fk+1,1),0]。 (21)
本发明提出一种总变分和结构张量的图像融合方法，总变分TV模型利用图像的结构张量获取融合图像的目标梯度场，在目标梯度场内部利用Soft-Max函数柔化梯度场的权重，将融合梯度场放入到TV模型的正则项内，保障融合图像的梯度场与目标梯度场接近，利用分裂Bregman算法求解全新的TV模型，获取最终的融合图像，通过采用本方案可以更好地表达图像的特征，将微光的纹理信息与红外图像的显著性特征完美结合，取得较好的融合效果，可以较完整地保留了丰富的纹理信息和细节信息，并且图像过渡自然。
具体地，结合图2和3所示，本发明的总变分和结构张量的图像融合方法还提供了一种实施例，包括以下步骤：
步骤1：本发明选取来自同一场景并且经过严格配准红外与可见光图像。
步骤2：利用图像的结构张量和Soft-Max函数获取融合图像的加权目标梯度：
1.对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
2.考虑到合成梯度方向与原各波段梯度方向一般都存在一定的一致性，我们可以使合成梯度的方向与各图像的平均梯度方向在同一侧，于是可定义目标梯度为：
3.由于异源图像的光谱差异较大，融合的图像容易产生光晕和伪影等现象，为此本文利用Soft-Max函数构成目标梯度场的梯度权重，其表达式如下所示：
其中a为收缩因子，它控制着Soft-Max函数曲线的形状。
步骤3：本文利用结构张量获取的目标梯度场放入到总变分模型的正则项内，从而约束融合图像的梯度场与目标梯度场Vf接近。其表达式如下：
步骤4：总变分模型的数据保真项是使融合图像与两幅源图像相差不大，因此它决定了最终的主观融合视觉效果。如果只约束融合源图像与其中一种源图像接近，这显然是不妥当的。为此，我们在数据保真项内加入权重函数，用以表达源图像对融合图像的贡献程度，其表达式如下：
其中f代表最终的融合图像，i代表红外图像，v代表可见光图像，α与β是两个非负权重函数，它可以根据源图像灰度值的变化而变化，从而平衡源图像与融合图像之间的像素强度配比。由于数据保真项还需要兼容两种源图像的光谱特性，因此其非负权重函数的选择决定了最终的融合视觉特性。为此能充分挖掘源图像内部的显著性特征，本文采用如下方式定义权重函数：
其中F(x,y)表示两种源图像中像素(x,y)的强度值，p表示像素强度，Mp表示其强度等于p的像素的数量，并且L是灰度级的数量(本文中为256)，S(x,y)代表源图像的显著特征矩阵。
步骤5：总变分模型的建立
1.我们将TV模型的数据保真项与正则项相结合，得到如下表达式：
2.将上述融合模型用离散化的l2范数与l1范数再次简化，将其表示为：
其中λ为正则参数，并且：
步骤6：总变分模型的求解
鉴于分裂Bregman算法是求解带有l1正则项的快速算法，在本节中我们采用该方法实现全新模型的最优求解，其求解过程如下：
第一步：引入变量d，将公式(8)变为以下形式：
第二步：1.fk+1可以按照如下方式求解：
2.f的最小化方程由以下最优性条件给出：
(f-αi-βv)-Δf+γdiv(Vf+dk-bk)＝0； (12)
3.由于公式(12)是线性的，它可以使用快速傅立叶变换实现求解。为此，把公式(12)重写为如下形式：
f-Δf＝αi+βv-γdiv(Vf+dk-bk)； (13)
4.公式(13)的快速傅立叶(FFT)形式为：
F(K)F(fk+1)＝F(rhsk)； (14)
K＝I-Δ； (15)
rhk＝αi+βv-γdiv(Vf+dk-bk)； (16)
其中I表示单位矩阵，F表示快速傅里叶变换。至此，公式(14)中fk+1的闭合形式解为：
第三步：1.dk+1可以按照如下方式求解：
2.公式(18)的解可以使用如下的软阈值形式获得：
3.由Shrink算子得到dk+1的解为：
第四步：为防止fk+1超出范围，在每一步迭代中对其做出如下约束：
fk+1＝max[min(fk+1,1),0]； (21)
图像融合结果的性能可以通过主观和客观的方式进行评价。在大多数情况下，融合结果之间的差异不大，则难以用主观方式正确地评价融合结果。因此，在客观质量评价的基础上对融合效果进行评价也是非常必要的，本文选取了以下5个客观质量指标作为评判标准：1.平均梯度(AVG)2.信息熵(IE)3.边缘信息保留量(QAB/F)4.空间频率(SF)5.标准差(SD)。详细地说，AVG用来体现图像中的微小细节反差和纹理变化特征；IE可以直接用来衡量图像信息的丰富程度；SF用来反映图像在空间域的总体活跃程度；SD用来反应像素灰度值的分布状况；QAB/F评价从源图像传递到融合图像的边缘信息量。对于这5个指标，数值越大代表性能越好。
结合图4所示，下面通过采用多种算法进行仿真对比各种算法的图像融合效果，本仿真采用基于WLS的融合算法、基于PCNN的融合算法、基于GFF的融合算法、基于HMSD的融合算法、基于MWGF的融合算法、基于GTF的融合算法、基于l1-TV的融合算法等作为对比，其中，图4(a)为红外图像，图4(b)为可见光图像，图4(c)是基于WLS的融合效果图，图4(d)是基于PCNN的融合效果图，图4(e)是基于CFF的融合效果图，图4(f)是基于HMSD的融合效果图，图4(g)是基于MWGF的融合效果图，图4(h)是基于GTF的融合效果图，图4(i)是基于l1-TV的融合效果图，图4(j)是基于本发明算法的融合效果图。
表1表示不同算法的评价指标对比。
从图4(c)-(j)可以看出，图4(d)方法的融合图像依旧还有伪影噪点。图4(e)方法的融合图像较为接近VI图像，因此丢失了IR图像的显著性特征。与图4(e)的观感相反，图4(g)-(i)的观感依旧接近IR图像，因此丢了很多包含在VI图像的纹理细节信息。图4(c)和图4(f)方法的融合图像与图4(j)较为接近，但是本文算法的人物更加清晰，同时背景的梯度特征也保留的最多。同时从表一也可以发现，除了SF值低于VSM-WLS算法外，剩下的评价参数值均高于其它对比算法，这一定程度与主观观感相似，说明本文算法能保证在人眼视觉基础的可见光场景下，极大的融合了红外目标与背景信息，同时在融合过程中没有遗失边缘细节信息，因此观感最佳。由于全新的优化模型在数据保真项加入了显著性权重函数，因此融合图像在主体观感上兼容了源图像的光谱特性，凸显了各自的显著性信息，更加符合人眼的视觉观感。除此外，正则项利用融合的梯度场还原了更多的梯度特征，使融合图像细节更加丰富。
相应地，本发明提供一种总变分和结构张量的图像融合装置，所述装置包括：
图像配准单元，用于获取来自同一场景并进行图像配准的红外图像和可见光图像；
目标梯度场获取单元，用于利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场；
权重获取单元，用于利用Soft-Max函数柔化所述目标梯度场的权重；
处理单元，用于将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近；
融合单元，用于利用分裂Bregman算法求解TV模型得到融合图像。
可选地，所述图像配准单元具体用于：
分别获取同一场景的红外图像和可见光图像；
对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
可选地，所述目标梯度场获取单元具体用于：
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
可选地，所述处理单元具体用于：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
本发明提出一种总变分和结构张量的图像融合装置，总变分TV模型利用图像的结构张量获取融合图像的目标梯度场，在目标梯度场内部利用Soft-Max函数柔化梯度场的权重，将融合梯度场放入到TV模型的正则项内，保障融合图像的梯度场与目标梯度场接近，利用分裂Bregman算法求解全新的TV模型，获取最终的融合图像，通过采用本方案可以更好地表达图像的特征，将微光的纹理信息与红外图像的显著性特征完美结合，取得较好的融合效果，可以较完整地保留了丰富的纹理信息和细节信息，并且图像过渡自然。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器(ROM，Read Only Memory)、随机存取存储器(RAM，RandomAccess Memory)、磁盘或光盘等。
以上对本发明所提供的一种总变分和结构张量的图像融合方法及装置进行了详细介绍，对于本领域的一般技术人员，依据本发明实施例的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种总变分和结构张量的图像融合方法，其特征在于，所述方法包括：
获取来自同一场景并进行图像配准的红外图像和可见光图像；
利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场；
利用Soft-Max函数柔化所述目标梯度场的权重；
将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近；
利用分裂Bregman算法求解TV模型得到融合图像。
2.根据权利要求1所述的总变分和结构张量的图像融合方法，其特征在于，所述获取来自同一场景并进行图像配准的红外图像和可见光图像，包括：
分别获取同一场景的红外图像和可见光图像；
对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
3.根据权利要求1所述的总变分和结构张量的图像融合方法，其特征在于，所述利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场，包括：
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
4.根据权利要求3所述的总变分和结构张量的图像融合方法，其特征在于，所述利用Soft-Max函数柔化所述目标梯度场的权重，包括：
利用Soft-Max函数构成目标梯度场的梯度权重，表达式如下：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
5.根据权利要求4所述的总变分和结构张量的图像融合方法，其特征在于，所述将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近，包括：
利用结构张量获取的目标梯度场放入到总变分TV模型的正则项内，约束融合图像的融合梯度场与目标梯度场Vf接近，表达式如下：
在数据保真项内加入权重函数，用于表达源图像对融合图像的贡献程度，表达式如下：
其中f代表融合图像，i代表红外图像，v代表可见光图像，α与β是两个非负权重函数，定义权重函数的表达式为：
其中F(x,y)表示两种源图像中像素(x,y)的强度值，p表示像素强度，Mp表示其强度等于p的像素的数量，并且L是灰度级的数量，S(x,y)代表源图像的显著特征矩阵。
6.根据权利要求5所述的总变分和结构张量的图像融合方法，其特征在于，所述利用分裂Bregman算法求解TV模型得到融合图像，包括：
总变分模型的建立
将TV模型的数据保真项与正则项相结合，得到如下表达式：
将TV模型用离散化的l2范数与l1范数再次简化，表示为：
其中λ为正则参数，并且：
总变分模型的求解
采用分裂Bregman算法实现TV模型的最优求解，其求解过程如下：
引入变量d，将公式(8)变为以下形式：
对fk+1可以按照如下方式求解：
对f的最小化方程由以下最优性条件给出：
(f-αi-βv)-Δf+γdiv(Vf+dk-bk)＝0； (12)
根据公式(12)是线性的性质，把公式(12)重写为如下形式：
f-Δf＝αi+βv-γdiv(Vf+dk-bk)； (13)
对公式(13)的进行快速傅立叶FFT形式为：
F(K)F(fk+1)＝F(rhsk)； (14)
K＝I-Δ； (15)
rhk＝αi+βv-γdiv(Vf+dk-bk)； (16)
其中I表示单位矩阵，F表示快速傅里叶变换；
公式(14)中fk+1的闭合形式解为：
dk+1可以按照如下方式求解：
公式(18)的解使用如下的软阈值形式获得：
由Shrink算子得到dk+1的解为：
对每一步迭代中对其做出如下约束：
fk+1＝max[min(fk+1,1),0] (21)。
7.一种总变分和结构张量的图像融合装置，其特征在于，所述装置包括：
图像配准单元，用于获取来自同一场景并进行图像配准的红外图像和可见光图像；
目标梯度场获取单元，用于利用红外图像和可见光图像的结构张量获取融合图像的目标梯度场；
权重获取单元，用于利用Soft-Max函数柔化所述目标梯度场的权重；
处理单元，用于将融合梯度场放入总变分TV模型的正则项内，以使得融合图像的融合梯度场和目标梯度场接近；
融合单元，用于利用分裂Bregman算法求解TV模型得到融合图像。
8.根据权利要求7所述的总变分和结构张量的图像融合装置，其特征在于，所述图像配准单元具体用于：
分别获取同一场景的红外图像和可见光图像；
对红外图像和可见光图像进行特征提取得到特征点，通过进行相似性度量找到匹配的特征点对，通过匹配的特征点对得到图像空间坐标变换参数，由坐标变换参数进行红外图像和可见光图像的图像配准。
9.根据权利要求7所述的总变分和结构张量的图像融合装置，其特征在于，所述目标梯度场获取单元具体用于：
对于每一幅源图像In(n＝1,2,····,N)及每一个像素元(x,y)，定义一个归一化权重Sn(x,y)，表示在图像In点(x,y)附近的显著性指标，可得加权的G矩阵为：
利用合成梯度方向与原各波段梯度方向存在一致性使得合成梯度的方向与各图像的平均梯度方向在同一侧，定义目标梯度场VF为：
10.根据权利要求9所述的总变分和结构张量的图像融合装置，其特征在于，所述处理单元具体用于：
其中a为收缩因子，用于控制Soft-Max函数曲线的形状。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811055297.6A CN109345496B (zh) | 2018-09-11 | 2018-09-11 | 一种总变分和结构张量的图像融合方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811055297.6A CN109345496B (zh) | 2018-09-11 | 2018-09-11 | 一种总变分和结构张量的图像融合方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109345496A true CN109345496A (zh) | 2019-02-15 |
| CN109345496B CN109345496B (zh) | 2021-05-14 |
Family
ID=65305323
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811055297.6A Active CN109345496B (zh) | 2018-09-11 | 2018-09-11 | 一种总变分和结构张量的图像融合方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109345496B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110084748A (zh) * | 2019-03-26 | 2019-08-02 | 温州晶彩光电有限公司 | 一种基于整体变分的红外和可见光图像融合方法 |
| CN110853058A (zh) * | 2019-11-12 | 2020-02-28 | 电子科技大学 | 基于视觉显著性检测的高分辨率遥感影像道路提取方法 |
| CN112596104A (zh) * | 2020-12-09 | 2021-04-02 | 成都理工大学 | 一种结合张量分解和全变分的地震资料去噪方法 |
| CN113126149A (zh) * | 2018-12-28 | 2021-07-16 | 中国石油化工股份有限公司 | 增强地质构造保真性的地震图像处理的方法和系统 |
| CN113362362A (zh) * | 2021-06-17 | 2021-09-07 | 易普森智慧健康科技（深圳）有限公司 | 基于总变分区域选取的明视野显微镜全景图像对齐算法 |
| CN115908221A (zh) * | 2023-03-08 | 2023-04-04 | 荣耀终端有限公司 | 图像处理方法、电子设备及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504670A (zh) * | 2014-12-11 | 2015-04-08 | 上海理工大学 | 多尺度梯度域图像融合算法 |
| CN105608679A (zh) * | 2016-01-28 | 2016-05-25 | 重庆邮电大学 | 一种融合结构张量与非局域全变分的图像去噪方法 |
| CN107679242A (zh) * | 2017-10-30 | 2018-02-09 | 河海大学 | 融合多信息源耦合张量分解的标签推荐方法 |
| CN107945145A (zh) * | 2017-11-17 | 2018-04-20 | 西安电子科技大学 | 基于梯度置信变分模型的红外图像融合增强方法 |
| CN108399611A (zh) * | 2018-01-31 | 2018-08-14 | 西北工业大学 | 基于梯度正则化的多聚焦图像融合方法 |
-
2018
- 2018-09-11 CN CN201811055297.6A patent/CN109345496B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504670A (zh) * | 2014-12-11 | 2015-04-08 | 上海理工大学 | 多尺度梯度域图像融合算法 |
| CN105608679A (zh) * | 2016-01-28 | 2016-05-25 | 重庆邮电大学 | 一种融合结构张量与非局域全变分的图像去噪方法 |
| CN107679242A (zh) * | 2017-10-30 | 2018-02-09 | 河海大学 | 融合多信息源耦合张量分解的标签推荐方法 |
| CN107945145A (zh) * | 2017-11-17 | 2018-04-20 | 西安电子科技大学 | 基于梯度置信变分模型的红外图像融合增强方法 |
| CN108399611A (zh) * | 2018-01-31 | 2018-08-14 | 西北工业大学 | 基于梯度正则化的多聚焦图像融合方法 |
Non-Patent Citations (2)
| Title |
|---|
| AROOJ NAWAZ: ""Fusion of Color and Infrared Images Using Gradient Transfer and Total Variation Minimization"", 《IEEE》 * |
| SANJAY VISWANATH: ""MULTI-RESOLUTION BASED SPATIALLY ADAPTIVE MULTI-ORDER TOTAL VARIATION"", 《2018 IEEE 15TH INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING (ISBI 2018)》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113126149A (zh) * | 2018-12-28 | 2021-07-16 | 中国石油化工股份有限公司 | 增强地质构造保真性的地震图像处理的方法和系统 |
| CN113126149B (zh) * | 2018-12-28 | 2024-04-09 | 中国石油化工股份有限公司 | 增强地质构造保真性的地震图像处理的方法和系统 |
| CN110084748A (zh) * | 2019-03-26 | 2019-08-02 | 温州晶彩光电有限公司 | 一种基于整体变分的红外和可见光图像融合方法 |
| CN110084748B (zh) * | 2019-03-26 | 2023-04-25 | 温州晶彩光电有限公司 | 一种基于整体变分的红外和可见光图像融合方法 |
| CN110853058A (zh) * | 2019-11-12 | 2020-02-28 | 电子科技大学 | 基于视觉显著性检测的高分辨率遥感影像道路提取方法 |
| CN112596104A (zh) * | 2020-12-09 | 2021-04-02 | 成都理工大学 | 一种结合张量分解和全变分的地震资料去噪方法 |
| CN112596104B (zh) * | 2020-12-09 | 2022-04-01 | 成都理工大学 | 一种结合张量分解和全变分的地震资料去噪方法 |
| CN113362362A (zh) * | 2021-06-17 | 2021-09-07 | 易普森智慧健康科技（深圳）有限公司 | 基于总变分区域选取的明视野显微镜全景图像对齐算法 |
| CN115908221A (zh) * | 2023-03-08 | 2023-04-04 | 荣耀终端有限公司 | 图像处理方法、电子设备及存储介质 |
| CN115908221B (zh) * | 2023-03-08 | 2023-12-08 | 荣耀终端有限公司 | 图像处理方法、电子设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109345496B (zh) | 2021-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109345496A (zh) | 一种总变分和结构张量的图像融合方法及装置 | |
| Bavirisetti et al. | Multi-scale guided image and video fusion: A fast and efficient approach | |
| CN106201535B (zh) | 工具栏背景颜色随着图片的主颜色而变换的方法和装置 | |
| CA2963861C (en) | Methods, systems and computer program products for visualizing anatomical structures and blood flow and perfusion physiology using imaging techniques | |
| JP6837551B2 (ja) | Ｈｍｄｓに基づく医学画像形成装置 | |
| CN109472270A (zh) | 图像风格转换方法、装置及设备 | |
| CN110120047A (zh) | 图像分割模型训练方法、图像分割方法、装置、设备及介质 | |
| CN109345494A (zh) | 基于潜在低秩表示和结构张量的图像融合方法及装置 | |
| US10032297B2 (en) | Simulation system, simulation device, and product explanation assistance method | |
| CN110414607A (zh) | 胶囊内窥镜图像的分类方法、装置、设备及介质 | |
| CN106650215A (zh) | 一种基于云平台的肤质检测及个性化评价系统及方法 | |
| CN109102485A (zh) | 基于nsst和自适应双通道pcnn的图像融合方法及装置 | |
| CN107809582A (zh) | 图像处理方法、电子装置和计算机可读存储介质 | |
| CN106339201A (zh) | 贴图处理方法、装置和电子设备 | |
| CN102663741B (zh) | 对彩色数字图像进行视觉立体感知增强的方法及系统 | |
| CN106462986A (zh) | 用于数字图像的颜色处理的方法和系统 | |
| WO1998015226A1 (fr) | Procede et appareil de formation et d'affichage d'une image a partir d'une pluralite d'images partielles | |
| CN110009656B (zh) | 目标对象的确定方法、装置、存储介质及电子装置 | |
| CN110363296B (zh) | 任务模型获取方法和装置、存储介质及电子装置 | |
| CN109242812A (zh) | 基于显著性检测与奇异值分解的图像融合方法及装置 | |
| KR100993875B1 (ko) | 퍼지 기반의 퍼스널 컬러 진단 시스템 및 그 방법 | |
| CN109242813A (zh) | 基于潜在低秩表示与nsst的图像融合方法及装置 | |
| KR102405613B1 (ko) | 공간 이미지 분석 기반 배색 추천 장치 및 방법 | |
| CN109891519A (zh) | 信息处理装置、信息处理方法和程序 | |
| EP2790153A2 (en) | Environment evaluation apparatus, method and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |