CN105719263A - 基于nsct域底层视觉特征的可见光和红外图像融合算法 - Google Patents
基于nsct域底层视觉特征的可见光和红外图像融合算法 Download PDFInfo
- Publication number
- CN105719263A CN105719263A CN201610044134.2A CN201610044134A CN105719263A CN 105719263 A CN105719263 A CN 105719263A CN 201610044134 A CN201610044134 A CN 201610044134A CN 105719263 A CN105719263 A CN 105719263A
- Authority
- CN
- China
- Prior art keywords
- theta
- image
- visible light
- sigma
- frequency
- 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
- 230000004927 fusion Effects 0.000 title claims abstract description 73
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 18
- 230000000007 visual effect Effects 0.000 title claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 36
- 230000000694 effects Effects 0.000 claims description 30
- 238000001218 confocal laser scanning microscopy Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 5
- 238000007500 overflow downdraw method Methods 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 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
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提出了一种基于非下采样轮廓波变换(NSCT)域底层视觉特征的可见光和红外图像融合算法。首先将可见光和红外图像经NSCT变换得到二者的高低频子带系数,然后联合相位一致性、邻域空间频率和邻域能量等信息,综合衡量低频子带系数的像素活跃水平,分别得到可见光和红外图像低频子带系数的融合权重,从而得到融合图像的低频子带系数;结合相位一致性、清晰度和亮度等信息衡量高频子带系数的像素活跃水平,分别得到可见光和红外图像高频子带系数的融合权重,进而得到融合图像的高频子带系数,最后利用NSCT逆变换,得到最终的融合图像。本发明能够有效地保留源图像的细节信息,同时综合可见光图像和红外图像的有用信息。
Description
技术领域
本发明涉及一种基于非下采样轮廓波变换(NSCT)域底层视觉特征的可见光和红外图像融合算法,属于图像融合数据处理技术领域。
背景技术
在图像处理领域,图像融合是一项很有发展前景的研究。图像融合技术通过综合同一场景的多传感器图像或者不同场景的同一传感器图像的有用信息来合成一幅复合图像。而合成的复合图像具有前者所有的特征信息,更适合用于后期的处理和研究。一个高效的融合方法可以根据实际需求处理多通道信息。这些优势使得图像融合在很多领域里都特别受欢迎。例如,遥感成像,计算机视觉,侦查,监控等领域。尤其是可见光和红外图像的融合,在军事领域中起到了很重要的作用。
随着图像融合研究的不断深入,很多的学者提出了很多的图像融合算法。根据对输入图像的处理方式的不同,主流的红外和可见光图像融合算法主要分为五类:基于置换的技术,基于分割技术,基于神经网络技术,基于数学形态学技术和基于多尺度变换技术。基于置换的技术在融合中容易引起光谱的变形。基于分割的技术过分依赖于分割算法,而通常一个好的图像分割算法往往是不容易得到的。基于神经网络的算法需要人工设置大量的参数。基于数学形态学的算法在图像融合中也发挥着很重要的作用,但是一些原始图像中的细节可能会在最终的融合结果中被平滑。和以上四类算法相比,多尺度变换方法更符合人眼的视觉特征。
多尺度变换对图像融合来说是一项有效的方法。近年来,多尺度变换方法已经在很多图像融合的研究中得到广泛的应用,尤其是可见光和红外图像的融合。一些常见的多尺度变换有拉普拉斯金字塔(LP)变换,金字塔变换(PT),离散小波变换(DWT),静态小波变换(SWT),轮廓波变换(CT),曲波变换(CVT)和非下采样轮廓波变换(NSCT)等。然而,拉普拉斯金字塔变换,金字塔变换和静态小波变换都无法很好地保存源图像的空间细节信息并且容易导致块效应。离散小波变换由于不能捕捉源图像中的静态特征而导致在融合图像中出现“吉布斯”现象。曲波变化也不能捕捉到源图像中的空间信息。对于轮廓波变换来说,由于在分解中存在上采样和下采样操作,缺乏平移不变性,容易导致混频现象。而非下采样轮廓波变换取消了上下采样操作,从而克服了这一缺点。
在图像融合中,除了融合工具之外,有效的融合规则也至关重要。近年来,有很多的学者提出了很多有效的基于多尺度变换的融合算法。大概可以分为三类:基于单个像素的,基于窗口的和基于区域的。基于单个像素的融合最常见的规则是“绝对值最大选择”策略;基于窗口的融合规则主要是根据当前像素的活跃水平,而活跃水平由窗口内像素之间的关系来度量,常见的活跃因子包括绝对分量和(SML)、空间频率(SF)、对比度(Contrast)等等;基于区域的融合方法的基本思想就是对不同子带系数进行分割,然后选择更清晰的区域进行融合。对于多尺度融合方法来说,清晰度因子在提高融合性能方面至关重要,通常传统的清晰度因子是基于单个图像局部特征设计的。实际上,人眼视觉系统主要是靠图像的底层特征来理解一幅图像的,也就是说由清晰度这一个底层特征所反映的视觉信息是非常有限的。
发明内容
本发明的目的在于针对现有技术的缺点和不足,提出了一种基于NSCT域底层视觉特征的可见光和红外图像融合技术方案。
本发明所采用的技术方案是一种基于NSCT域底层视觉特征的可见光和红外图像融合算法,包括以下步骤:
步骤1,输入可见光和红外图像并进行NSCT变换,分别得到可见光图像和红外图像的子带系数,所述子带系数包括低频子带系数和高频子带系数;
步骤2,根据可见光和红外图像低频子带系数的相位一致性、邻域空间频率和邻域能量信息计算融合图像的低频子带系数;
计算可见光和红外图像低频子带系数的相位一致性:
分别表示可见光和红外图像在像素点(x,y)的低频子带系数的相位一致性,L表示低频,VI和IR分别表示可见光和红外图像,n={0,1,...,J-1},k={0,1,...,K-1},J和K分别表示可见光和红外图像在傅里叶频域分解的尺度数量和方向数量,θk表示滤波器的方向角,分别表示可见光和红外图像像素点(x,y)的低频子带系数在第n尺度,θk方向角的傅里叶频域的局部能量,分别表示可见光和红外图像像素点(x,y)的低频子带系数在第n尺度,θk方向角的傅里叶频域的局部振幅信息,ε为正常量;
计算可见光和红外图像低频子带系数的空间频率:
SCMVI(x,y)、SCMIR(x,y)分别表示可见光图像和红外图像在像素点(x,y)的低频子带系数的空间频率,IL,VI(x,y)、IL,IR(x,y)分别表示可见光图像、红外图像在像素点(x,y)的低频子带系数,IL,VI(x0,y0)、IL,IR(x0,y0)分别表示可见光图像、红外图像在像素点(x0,y0)的低频子带系数,Ω0表示预设的邻域,(x0,y0)表示以像素点(x,y)为中心的邻域Ω0内任意一点;
将可见光和红外图像中所有像素的低频子带系数的空间频率构成矩阵并计算可见光和红外图像低频子带系数的邻域空间频率LSCMVI(x,y),LSCMIR(x,y);
计算可见光和红外图像低频子带系数的邻域能量:
LENVI(x,y)、LENIR(x,y)分别表示可见光图像、红外图像在像素点(x,y)的低频子带系数的邻域能量,IL,VI(x+m,y+n)、IL,IR(x+m,y+n)分别表示可见光图像、红外图像在像素点(x+m,y+n)的低频子带系数,M×N表示预设的邻域大小,(x+m,y+n)表示以(x,y)为中心的邻域M×N内任意一点;
计算可见光和红外图像低频子带系数的像素活跃水平:
NAMVI(x,y)=[PCL,VI(x,y)]α×[LSCMVI(x,y)]β×[LENVI(x,y)]γ
NAMIR(x,y)=[PCL,IR(x,y)]α×[LSCMIR(x,y)]β×[LENIR(x,y)]γ
NAMVI(x,y)、NAMIR(x,y)分别表示可见光和红外图像在像素点(x,y)的低频子带系数的活跃水平,α、β、γ分别表示相位一致性信息、邻域空间频率和邻域能量的权重;
将可见光和红外图像中所有像素的低频子带系数的活跃水平构成矩阵,设置滑动窗口Ω1,计算可见光和红外图像的低频子带系数的融合权重:
LmapVI(x,y)、LmapIR(x,y)表示可见光和红外图像在像素点(x,y)的低频子带系数的融合权重,X×Y表示滑动窗口Ω1大小,LsumVI(x,y)、LsumIR(x,y)分别表示可见光和红外图像在滑动窗口Ω1内活跃水平较大的像素数量;
表示集合的势,(x0,y0)表示滑动窗口Ω1内任意一点;
计算融合图像低频子带系数IL,F(x,y):
IL,F(x,y)=IL,VI(x,y)×LmapVI(x,y)+IL,IR(x,y)×LmapIR(x,y);
步骤3,根据可见光和红外图像高频子带系数的相位一致性信息、清晰度和图像亮度计算融合图像的高频子带系数;
计算可见光和红外图像高频子带系数的相位一致性:
分别表示可见光和红外图像在像素点(x,y)的第l尺度、第p方向高频子带的相位一致性信息,l和p分别表示高频子带系数的尺度和方向,分别表示可见光和红外图像像素点(x,y)的第l尺度、第p方向高频子带系数在第n尺度、θk方向角的傅里叶频域的局部能量,分别表示可见光和红外图像像素点(x,y)的第l尺度、第p方向高频子带系数在第n尺度、θk方向角的傅里叶频域的局部振幅,ε是正常量;
计算可见光和红外图像高频子带系数的清晰度:
分别表示可见光和红外图像在像素点(x0,y0)的第l尺度、第p方向的高频子带系数的清晰度,H表示高频子带系数,分别表示可见光和红外图像在像素点(x,y)的第l尺度,第p方向的高频子带系数,分别表示可见光和红外图像在像素点(x0,y0)的第l尺度,第p方向的高频子带系数,Ω2表示预设的领域,(x0,y0)表示以(x,y)为中心的预设领域Ω2内任意一点;
计算可见光和红外图像在像素(x,y)第l尺度,第p方向的高频子带系数的像素活跃水平:
分别表示可见光和红外图像在像素点(x,y)的第l层、第p方向的高频子带系数的活跃水平,分别表示可见光和红外图像在像素点(x,y)的第l层、第p方向的高频子带系数的亮度,α1、β1、γ1分别表示相位一致性、清晰度和图像亮度的指数权重,Ω3表示预设的领域,(x0,y0)为滑动窗口预设领域Ω3内任意一点;
将可见光和红外图像中高频子带系数的活跃水平构成矩阵,设置一个滑动窗口Ω4,计算可见光和红外图像高频子带系数的融合权重:
X1×Y1表示滑动窗口Ω4的大小,和分别表示可见光和红外图像在像素(x,y)第l尺度,第p方向的高频子带系数在滑动窗口Ω4内活跃水平较大的像素数量:
(x0,y0)为滑动窗口Ω4内任意一点;
计算融合图像第l尺度,第p方向的高频子带系数
步骤4,根据步骤2、3得到的融合图像的低频子带系数和高频子带系数,进行NSCT反变换,得到融合图像。
可见光和红外图像低频子带系数的邻域空间频率LSCMVI(x,y),LSCMIR(x,y)计算如下:
SCMVI(x+m,y+n)、SCMIR(x+m,y+n)分别表示可见光和红外图像在像素点(x+m,y+n)处的低频子带系数的空间频率,M'×N'表示预设的窗口大小,(x+m,y+n)表示以(x,y)为中心的滑动窗口内任意一点。
可见光和红外图像低频子带系数在傅里叶域的局部能量 计算如下:
和分别表示可见光和红外图像低频子带系数与log-Gabor奇偶对称滤波器的卷积结果:
分别表示log-Gabor滤波器在尺度n下的一对奇偶对称滤波器;
可见光和红外图像的低频子带系数在傅里叶域的局部振幅信息 计算如下:
可见光和红外图像第l尺度,第p方向的的高频子带系数在傅里叶域的局部能量计算如下:
和分别表示可见光和红外图像第l尺度,第p方向的高频子带系数与log-Gabor奇偶对称滤波器的卷积结果:
可见光和红外图像第l尺度,第p方向的高频子带系数在傅里叶频域的局部振幅计算如下:
log-Gabor滤波器的传递函数G(ω,θ,k)定义如下,
(log(ω),θ)表示对数极坐标,(σr,σθ)表示滤波器的半径和角度带宽,ω0表示滤波器的中心频率,θk表示滤波器中的方向角。
本发明有益效果为:本发明根据图像的高低频子带系数底层视觉特征的特点,设计了两个衡量源图像高低子带系数对于融合结果选择的重要性的活跃水平。在图像的低频子带系数的融合上,综合了相位一致性信息、邻域空间频率信息和邻域能量信息设计低频子带系数的活跃水平,在高频子带系数的融合上,综合了相位一致性信息,清晰度和亮度信息设计高频子带系数的活跃水平,更好地保留了图像的纹理和边缘等细节信息,有效综合了可见光和红外图像二者的优点;可以有效地综合可见光和红外图像的有用信息,在主观和客观上都取得良好的效果。
附图说明
图1为本发明流程图;
图2为可见光与红外源图像,其中(a)为可见光源图像,(b)为红外源图像;
图3为DWT方法处理后的融合图像;
图4为SWT方法处理后的融合图像;
图5为NSCT方法处理后的融合图像;
图6为NSCT-PCNN方法处理后的融合图像;
图7为NSCT-SR方法处理后的融合图像;
图8为NSCT-SF-PCNN方法处理后的融合图像;
图9为本发明方法处理后的融合图像。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的详细说明。
本实施例是将图2中一对已配准的可见光和红外图像进行融合。IVI表示可见光图像,IIR表示红外图像,低频、高频子带系数分别为{IL,VI,IL,IR},低频、高频子带系数的像素活跃水平分别为{NAMVI(x,y),NAMIR(x,y)},低频、高频子带系数的融合权重分别为{LmapVI(x,y),LmapIR(x,y)},对应的融合子带系数分别为IL,F和融合图像为IF。计算步骤如下:
步骤1,输入可见光和红外图像并进行NSCT变换,分别得到可见光图像和红外图像的子带系数,所述的子带系数包括低频子带系数{IL,VI(x,y),IL,IR(x,y)}和高频子带系数
步骤2,在传统的低频子带系数融合中,通常采用基于局部能量的方法,然而局部能量只能反映图像的亮度信息。因此,我们根据可见光和红外图像的低频子带系数,采用联合相位一致性特征,邻域空间频率和局部能量的策略来计算融合图像的低频子带系数;
计算可见光和红外图像低频子带系数的相位一致性:
分别表示可见光和红外图像像素点(x,y)的低频子带的相位一致性信息,L表示低频,J和K分别表示图像在傅里叶频域分解的尺度数量和方向的数量,本实施例中J和K均取4,n={0,1,...,J-1},k={0,1,...,K-1},θk表示滤波器的方向角,分别表示可见光和红外图像像素点(x,y)的低频子带系数在第n尺度,θk方向角的傅里叶频域的局部能量信息;分别表示可见光和红外图像像素点(x,y)的低频子带在第n尺度,θk方向角的傅里叶频域的局部振幅信息,ε是避免分母为0的一个很小的正的常量;
本发明实施例中涉及到的图像低频子带系数在傅里叶域的局部能量和局部振幅分别如下计算,
设(x,y)为图像某一像素点,则可见光和红外图像在该像素点的低频子带系数在傅里叶频域的局部能量信息计算如下,
其中,和定义如下,
和分别为可见光和红外图像在像素点(x,y)的低频子带系数与log-Gabor奇偶对称滤波器的卷积结果,计算如下,
分别表示log-Gabor滤波器在尺度n下的一对奇偶对称滤波器;
可见光和红外图像在该像素点的低频子带系数在傅里叶频域的局部振幅信息计算如下,
对可见光和红外图像的每个像素,分别按照以下原则求低频子带系数的空间频率信息,
IL,VI(x,y)、IL,IR(x,y)分别表示可见光图像、红外图像中像素(x,y)的低频子带系数,IL,VI(x0,y0)、IL,IR(x0,y0)分别表示可见光图像、红外图像在像素点(x0,y0)的低频子带系数,SCMVI(x,y)、SCMIR(x,y)分别表示可见光图像和红外图像在像素点(x,y)的低频子带系数的空间频率,Ω0为预设的邻域,本实施例中Ω0采用3×3,(x0,y0)为以像素点(x,y)为中心的邻域Ω0内任意一点;
对可见光和红外图像,分别按照同样的方式将图像中所有像素的低频子带系数的空间频率构成矩阵并计算可见光和红外图像低频子带系数的邻域空间频率信息:LSCMVI(x,y),LSCMIR(x,y);本发明实施例中求取分别如下,
设滑动窗口遍历图像空间频率矩阵到任一位置时的窗口中心为图像某像素点(x,y),则可见光和红外图像低频子带系数在这个窗口内的邻域空间频率定义如下,
SCMVI(x+m,y+n),SCMIR(x+m,y+n)分别表示可见光和红外图像在像素点(x+m,y+n)处的低频子带系数的空间频率,LSCMVI(x,y),LSCMIR(x,y)分别表示可见光和红外图像在像素点(x,y)的低频子带系数的邻域空间频率信息;M'×N'为预设的窗口大小,本实施例中M'×N'取值为3×3,(x+m,y+n)为以(x,y)为中心的窗口内任意一点。
对可见光和红外图像的每个像素,分别按照以下原则求低频子带系数的邻域能量信息,
IL,VI(x+m,y+n)、IL,IR(x+m,y+n)分别表示可见光图像、红外图像像素(x+m,y+n)的低频子带系数,LENVI(x,y)、LENIR(x,y)分别表示可见光图像、红外图像在像素点(x,y)的低频子带系数的邻域能量信息,M×N为预设的邻域大小,本实施例M×N取值3×3;(x+m,y+n)为以(x,y)为中心的邻域内任意一点;
设(x,y)为可见光和红外图像中的某一像素,则根据综合相位一致性信息、邻域空间频率和邻域能量的策略,可见光和红外图像在像素(x,y)的低频子带系数的活跃水平定义如下,
NAMVI(x,y)=[PCL,VI(x,y)]α×[LSCMVI(x,y)]β×[LENVI(x,y)]γ
NAMIR(x,y)=[PCL,IR(x,y)]α×[LSCMIR(x,y)]β×[LENIR(x,y)]γ
NAMVI(x,y),NAMIR(x,y)分别表示可见光、红外图像在像素(x,y)的低频子带系数的活跃水平,α,β,γ分别表示相位一致性信息、邻域空间频率和邻域能量的指数权重,本实施例中α,β,γ分别取值0.005,1,1;
对可见光和红外图像,分别按照同样的方式将图像中所有像素的低频子带系数的活跃水平构成矩阵,设置一个滑动窗口Ω1,设滑动窗口遍历图像到任一位置时的中心像素为(x,y),
根据低频子带系数的活跃水平矩阵,计算可见光和红外图像在像素点(x,y)的低频子带系数的融合权重,融合权重定义如下,
X×Y表示滑动窗口Ω1大小,LsumVI(x,y)和LsumIR(x,y)分别表示可见光和红外图像在像素(x,y)的低频子带系数在滑动窗口内占有两者中活跃水平较大像素的数量,LsumVI(x,y)和LsumIR(x,y)定义如下,
其中,为集合的势,(x0,y0)为该窗口内任意一点,本实施例中选取Ω1大小为7×7。
设(x,y)为可见光和红外图像某一像素点,根据图像低频子带系数的融合权重,最终得到融合图像中该像素点的低频子带系数IL,F(x,y)如下,
IL,F(x,y)=IL,VI(x,y)×LmapVI(x,y)+IL,IR(x,y)×LmapIR(x,y)。
步骤3,高频子带系数包含了图像中大部分的纹理信息和边缘信息,在高频子带系数中,我们引入相位一致性信息探测特征信息,引入一种新的清晰度信息衡量高频子带的质量。考虑到当前子带系数的强度和同一区域中不同子带系数之间的关系,根据可见光和红外图像的高频子带系数,最终采用基于相位一致性信息、清晰度和图像亮度的综合信息来计算融合图像的高频子带系数,实现如下:
对可见光和红外图像的每个像素,分别按照以下原则求高频子带系数的相位一致性信息,
l和p分别表示高频子带系数在第l尺度和第p方向, 分别表示可见光和红外图像像素点(x,y)的第l尺度,第p方向高频子带的相位一致性信息;分别表示可见光和红外图像像素点(x,y)的第l尺度,第p方向高频子带在第n尺度,θk方向角的傅里叶频域的局部能量信息;分别表示可见光和红外图像像素点(x,y)的第l尺度,第p方向高频子带在第n尺度,θk方向角的傅里叶频域的局部振幅信息,ε是避免分母为0的一个很小的正的常量。本发明实施例中图像的高频子带系数在傅里叶频域的局部能量和局部振幅分别如下计算,
可见光和红外图像在像素点(x,y)的第l尺度,第p方向的高频子带系数在傅里叶频域的局部能量信息计算如下,
其中,和定义如下,
和分别为可见光和红外图像在像素点(x,y)第l尺度,第p方向的高频子带系数与log-Gabor奇偶对称滤波器的卷积结果,计算如下,
可见光和红外图像在像素点(x,y)的第l尺度,第p方向的高频子带系数在傅里叶频域的局部振幅信息计算如下,
设置一个滑动窗口Ω2,设在滑动窗口遍历图像到任意位置时的窗口中心为(x,y),分别计算可见光和红外图像在窗口内的高频子带系数的清晰度信息,
分别表示可见光和红外图像在像素点(x,y)处的第l层,第p方向的高频子带系数;分别表示可见光和红外图像在像素点(x0,y0)处的第l层,第p方向的高频子带系数; 分别表示可见光和红外图像在像素点(x0,y0)处的第l尺度,第p方向的高频子带系数的清晰度信息;(x0,y0)表示以(x,y)为中心的滑动窗口内任意一点,本实施例中Ω2大小为11×11;
设(x,y)为可见光和红外图像中的某一像素,根据相位一致性、清晰度和图像亮度综合信息融合的策略,该像素点高频子带系数的融合规则如下,
设置一个滑动窗口Ω3,设滑动窗口遍历图像到任一位置时的窗口中心为(x,y),可见光和红外图像在像素(x,y)第l尺度,第p方向的高频子带系数的活跃水平定义如下,
分别表示可见光、红外图像在像素(x,y)的第l层,第p方向的高频子带系数的活跃水平;表示高频子带系数的亮度;α1,β1,γ1分别表示相位一致性信息、清晰度和图像亮度的指数权重,(x0,y0)为滑动窗口内任意一点,本实施例中,α1,β1,γ1分别取0.005,1,1,滑动窗口Ω3大小取11×11。
对可见光和红外图像,分别按照同样的方式将图像中所有像素的高频子带系数的活跃水平构成活跃水平矩阵,设置一个滑动窗口Ω4,设滑动窗口遍历图像到任一位置时的中心像素为(x,y),
根据高频子带系数的活跃水平,计算可见光和红外图像在像素点(x,y)的高频子带系数的融合权重,融合权重定义如下,
X1×Y1表示滑动窗口大小,和分别表示可见光和红外图像在像素(x,y)的第l尺度,第p方向的高频子带系数在滑动窗口分别内占有两者中活跃水平较大像素的数量,和定义如下,
(x0,y0)为滑动窗口内任意一点,本实施例中Ω1大小为7×7。
设(x,y)为可见光和红外图像某一像素点,根据图像高频子带系数的融合权重,最终得到融合图像中该像素点的第l层,第p方向的高频子带系数如下,
步骤4,根据步骤2,3所得到的融合图像的低频子带系数{IL,F(x,y)}和高频子带系数进行NSCT反变换,得到融合图像{IF(x,y)}。
本实施例中用到的log-Gabor滤波器其传递函数定义如下,
其中,(log(ω),θ)表示对数极坐标,(σr,σθ)控制滤波器的半径和角度带宽,ω0为滤波器的中心频率,θk表示在尺度k下的方向角。
实验结果如下:
图2(a)-(b)分别为可见光和红外源图像(大小均为256×256像素),图3-8分别为在经典的DWT、SWT、NSCT、NSCT-PCNN(基于非下采样轮廓波变换和脉冲耦合神经网络相结合的方法)、NSCT-SR(基于非下采样轮廓波变换和稀疏表示相结合的方法)、NSCT-SF-PCNN(非下采样轮廓波变换、空间频率和脉冲耦合神经网络相结合的方法)方法下的融合结果,图9为本发明提出方法下的融合结果。所有方法的分解尺度设置为4,所有基于NSCT的方法均采用“pkva”作为金字塔滤波器,“9-7”作为方向滤波器。DWT、SWT、NSCT方法中的高低频子带系数融合分别采用模最大选择策略和均值策略。NSCT-PCNN方法中,低频子带系数采用均值策略,高频子带系数采用PCNN最大点火时间选择策略。
为方便从主观上对比不同方法下的融合效果,对图2-9进行局部放大,并在原区域和放大区域用箭头指向细节信息。从整体上可以看出,与其他方法相比,采用DWT、SWT、NSCT和NSCT-PCNN方法的融合图像对比度较低;从局部放大区域观察可以得知,除了本发明方法的融合图像完整保留了箭头所指细节特征之外,其他方法下的融合图像都部分丢失了该细节特征。
除主观上对融合结果进行评价之外,同时采用互信息MI、边缘信息度量算子QAB/F和视觉信息保真度VIF三个客观评价指标对不同方法下的融合效果进行评价。其中MI可评价源图像有多少信息量保留到了融合图像,QAB/F衡量融合图像中包含了多少源图像的边缘信息,VIF则可通过计算源图像和融合图像之间的共享信息评价融合性能,三个评价指标的数值越大说明融合效果越好。不同融合方法的融合结果客观评价如表1所示,本发明方法客观评价的三个评价指标指均最高,说明了本发明方法的融合效果最好,从客观上证明了本发明方法的有效可行。
表1:不同融合方法的融合结果客观评价
上面结合附图对本发明的具体实施例作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (5)
1.基于NSCT域底层视觉特征的可见光和红外图像融合算法,其特征在于,包括以下步骤:
步骤1,输入可见光和红外图像并进行NSCT变换,分别得到可见光图像和红外图像的子带系数,所述子带系数包括低频子带系数和高频子带系数;
步骤2,根据可见光和红外图像低频子带系数的相位一致性、邻域空间频率和邻域能量信息计算融合图像的低频子带系数;
计算可见光和红外图像低频子带系数的相位一致性:
分别表示可见光和红外图像在像素点(x,y)的低频子带系数的相位一致性,L表示低频,VI和IR分别表示可见光和红外图像,n={0,1,...,J-1},k={0,1,...,K-1},J和K分别表示可见光和红外图像在傅里叶频域分解的尺度数量和方向数量,θk表示滤波器的方向角, 分别表示可见光和红外图像像素点(x,y)的低频子带系数在第n尺度,θk方向角的傅里叶频域的局部能量, 分别表示可见光和红外图像像素点(x,y)的低频子带系数在第n尺度,θk方向角的傅里叶频域的局部振幅信息,ε为正常量;
计算可见光和红外图像低频子带系数的邻域空间频率:
SCMVI(x,y)、SCMIR(x,y)分别表示可见光图像和红外图像在像素点(x,y)的低频子带系数的空间频率,IL,VI(x,y)、IL,IR(x,y)分别表示可见光图像、红外图像在像素点(x,y)的低频子带系数,IL,VI(x0,y0)、IL,IR(x0,y0)分别表示可见光图像、红外图像在像素点(x0,y0)的低频子带系数,Ω0表示预设的邻域,(x0,y0)表示以像素点(x,y)为中心的邻域Ω0内任意一点;
将可见光和红外图像中所有像素的低频子带系数的空间频率构成矩阵并计算可见光和红外图像低频子带系数的邻域空间频率LSCMVI(x,y),LSCMIR(x,y);
计算可见光和红外图像低频子带系数的邻域能量:
LENVI(x,y)、LENIR(x,y)分别表示可见光图像、红外图像在像素点(x,y)的低频子带系数的邻域能量,IL,VI(x+m,y+n)、IL,IR(x+m,y+n)分别表示可见光图像、红外图像在像素点(x+m,y+n)的低频子带系数,M×N表示预设的邻域大小,(x+m,y+n)表示以(x,y)为中心的邻域M×N内任意一点;
计算可见光和红外图像低频子带系数的像素活跃水平:
NAMVI(x,y)=[PCL,VI(x,y)]α×[LSCMVI(x,y)]β×[LENVI(x,y)]γ
NAMIR(x,y)=[PCL,IR(x,y)]α×[LSCMIR(x,y)]β×[LENIR(x,y)]γ
NAMVI(x,y)、NAMIR(x,y)分别表示可见光和红外图像在像素点(x,y)的低频子带系数的活跃水平,α、β、γ分别表示相位一致性信息、邻域空间频率和邻域能量的权重;
将可见光和红外图像中所有像素的低频子带系数的活跃水平构成矩阵,设置滑动窗口Ω1,计算可见光和红外图像的低频子带系数的融合权重:
LmapVI(x,y)、LmapIR(x,y)表示可见光和红外图像在像素点(x,y)的低频子带系数的融合权重,X×Y表示滑动窗口Ω1大小,LsumVI(x,y)、LsumIR(x,y)分别表示可见光和红外图像在滑动窗口Ω1内活跃水平较大的像素数量;
表示集合的势,(x0,y0)表示滑动窗口Ω1内任意一点;
计算融合图像低频子带系数IL,F(x,y):
IL,F(x,y)=IL,VI(x,y)×LmapVI(x,y)+IL,IR(x,y)×LmapIR(x,y);
步骤3,根据可见光和红外图像高频子带系数的相位一致性信息、清晰度和图像亮度计算融合图像的高频子带系数;
计算可见光和红外图像高频子带系数的相位一致性:
分别表示可见光和红外图像在像素点(x,y)的第l尺度、第p方向高频子带的相位一致性信息,l和p分别表示高频子带系数的尺度和方向, 分别表示可见光和红外图像像素点(x,y)的第l尺度、第p方向高频子带系数在第n尺度、θk方向角的傅里叶频域的局部能量, 分别表示可见光和红外图像像素点(x,y)的第l尺度、第p方向高频子带系数在第n尺度、θk方向角的傅里叶频域的局部振幅,ε是正常量;
计算可见光和红外图像高频子带系数的清晰度:
分别表示可见光和红外图像在像素点(x0,y0)的第l尺度、第p方向的高频子带系数的清晰度,H表示高频子带系数, 分别表示可见光和红外图像在像素点(x,y)的第l尺度,第p方向的高频子带系数, 分别表示可见光和红外图像在像素点(x0,y0)的第l尺度,第p方向的高频子带系数,Ω2表示预设的领域,(x0,y0)表示以(x,y)为中心的预设领域Ω2内任意一点;
计算可见光和红外图像在像素(x,y)第l尺度,第p方向的高频子带系数的像素活跃水平:
分别表示可见光和红外图像在像素点(x,y)的第l层、第p方向的高频子带系数的活跃水平, 分别表示可见光和红外图像在像素点(x,y)的第l层、第p方向的高频子带系数的亮度,α1、β1、γ1分别表示相位一致性、清晰度和图像亮度的指数权重,Ω3表示预设的领域,(x0,y0)为滑动窗口预设领域Ω3内任意一点;
将可见光和红外图像中高频子带系数的活跃水平构成矩阵,设置一个滑动窗口Ω4,计算可见光和红外图像高频子带系数的融合权重:
X1×Y1表示滑动窗口Ω4的大小,和分别表示可见光和红外图像在像素(x,y)第l尺度,第p方向的高频子带系数在滑动窗口Ω4内活跃水平较大的像素数量:
(x0,y0)为滑动窗口Ω4内任意一点;
计算融合图像第l尺度,第p方向的高频子带系数
步骤4,根据步骤2、3得到的融合图像的低频子带系数和高频子带系数,进行NSCT反变换,得到融合图像。
2.根据权利要求1所述的基于NSCT域底层视觉特征的可见光和红外图像融合算法,其特征在于:所述可见光和红外图像低频子带系数的邻域空间频率LSCMVI(x,y),LSCMIR(x,y)计算如下:
SCMVI(x+m,y+n)、SCMIR(x+m,y+n)分别表示可见光和红外图像在像素点(x+m,y+n)处的低频子带系数的空间频率,M'×N'表示预设的窗口大小,(x+m,y+n)表示以(x,y)为中心的滑动窗口内任意一点。
3.根据权利要求1所述的基于NSCT域底层视觉特征的可见光和红外图像融合算法,其特征在于:所述可见光和红外图像低频子带系数在傅里叶域的局部能量 计算如下:
其中, 和定义如下,
和分别表示可见光和红外图像低频子带系数与log-Gabor奇偶对称滤波器的卷积结果:
分别表示log-Gabor滤波器在尺度n下的一对奇偶对称滤波器;
可见光和红外图像的低频子带系数在傅里叶域的局部振幅信息 计算如下:
4.根据权利要求1所述的基于NSCT域底层视觉特征的可见光和红外图像融合算法,其特征在于:所述可见光和红外图像第l尺度,第p方向的的高频子带系数在傅里叶域的局部能量 计算如下:
其中, 和定义如下,
和分别表示可见光和红外图像第l尺度,第p方向的高频子带系数与log-Gabor奇偶对称滤波器的卷积结果:
可见光和红外图像第l尺度,第p方向的高频子带系数在傅里叶频域的局部振幅 计算如下:
5.根据权利要求1所述的基于NSCT域底层视觉特征的可见光和红外图像融合算法,其特征在于:所述log-Gabor滤波器的传递函数G(ω,θ,k)定义如下,
(log(ω),θ)表示对数极坐标,(σr,σθ)表示滤波器的半径和角度带宽,ω0表示滤波器的中心频率,θk表示滤波器中的方向角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044134.2A CN105719263B (zh) | 2016-01-22 | 2016-01-22 | 基于nsct域底层视觉特征的可见光和红外图像融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044134.2A CN105719263B (zh) | 2016-01-22 | 2016-01-22 | 基于nsct域底层视觉特征的可见光和红外图像融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105719263A true CN105719263A (zh) | 2016-06-29 |
CN105719263B CN105719263B (zh) | 2018-05-25 |
Family
ID=56153794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610044134.2A Active CN105719263B (zh) | 2016-01-22 | 2016-01-22 | 基于nsct域底层视觉特征的可见光和红外图像融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105719263B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204509A (zh) * | 2016-07-07 | 2016-12-07 | 西安电子科技大学 | 基于区域特性的红外与可见光图像融合方法 |
CN106548467A (zh) * | 2016-10-31 | 2017-03-29 | 广州飒特红外股份有限公司 | 红外图像和可见光图像融合的方法及装置 |
CN106846289A (zh) * | 2017-01-17 | 2017-06-13 | 中北大学 | 一种基于显著性迁移与细节分类的红外光强与偏振图像融合方法 |
CN106846288A (zh) * | 2017-01-17 | 2017-06-13 | 中北大学 | 一种双模态红外图像差异特征指数测度的多算法融合方法 |
CN107451984A (zh) * | 2017-07-27 | 2017-12-08 | 桂林电子科技大学 | 一种基于混合多尺度分析的红外与可见光图像融合算法 |
CN108764249A (zh) * | 2018-04-23 | 2018-11-06 | 云南民族大学 | 一种旋转不变多源影像局部特征描述方法、系统和装置 |
CN108961154A (zh) * | 2018-07-13 | 2018-12-07 | 福州大学 | 基于改进非下采样轮廓波变换的太阳电池热斑检测方法 |
CN109035189A (zh) * | 2018-07-17 | 2018-12-18 | 桂林电子科技大学 | 基于柯西模糊函数的红外与弱可见光图像融合方法 |
CN109489576A (zh) * | 2018-10-19 | 2019-03-19 | 杭州电子科技大学 | 一种基于初级视觉通路计算模型的轮廓检测方法 |
CN110211082A (zh) * | 2019-05-31 | 2019-09-06 | 浙江大华技术股份有限公司 | 一种图像融合方法、装置、电子设备及存储介质 |
CN110363731A (zh) * | 2018-04-10 | 2019-10-22 | 杭州海康威视数字技术股份有限公司 | 一种图像融合方法、装置及电子设备 |
CN110796632A (zh) * | 2019-07-30 | 2020-02-14 | 重庆小富农康农业科技服务有限公司 | 猪只计数装置 |
CN111462028A (zh) * | 2020-03-16 | 2020-07-28 | 中国地质大学(武汉) | 基于相位一致性和目标增强的红外和可见光图像融合方法 |
CN112102217A (zh) * | 2020-09-21 | 2020-12-18 | 四川轻化工大学 | 一种可见光图像与红外图像快速融合方法及系统 |
CN113947554A (zh) * | 2020-07-17 | 2022-01-18 | 四川大学 | 一种基于nsst和显著信息提取的多聚焦图像融合方法 |
CN116503454A (zh) * | 2023-06-27 | 2023-07-28 | 季华实验室 | 红外和可见光图像融合方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100239144A1 (en) * | 2009-02-20 | 2010-09-23 | Gabor Fichtinger | Marker Localization Using Intensity-Based Registration of Imaging Modalities |
CN105069768A (zh) * | 2015-08-05 | 2015-11-18 | 武汉高德红外股份有限公司 | 一种可见光图像与红外图像融合处理系统及融合方法 |
-
2016
- 2016-01-22 CN CN201610044134.2A patent/CN105719263B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100239144A1 (en) * | 2009-02-20 | 2010-09-23 | Gabor Fichtinger | Marker Localization Using Intensity-Based Registration of Imaging Modalities |
CN105069768A (zh) * | 2015-08-05 | 2015-11-18 | 武汉高德红外股份有限公司 | 一种可见光图像与红外图像融合处理系统及融合方法 |
Non-Patent Citations (3)
Title |
---|
LIN ZHANG ET AL.: "FSIM: A Feature Similarity Index for Image Quality Assessment", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
WEIWEI KONG ET AL.: "Novel fusion method for visible light and infrared images based on NSST–SF–PCNN", 《INFRARED PHYSICS & TECHNOLOGY》 * |
张蕾 等: "采用非采样Contourlet变换与区域分类的红外和可见光图像融合", 《光学精密工程》 * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204509A (zh) * | 2016-07-07 | 2016-12-07 | 西安电子科技大学 | 基于区域特性的红外与可见光图像融合方法 |
CN106204509B (zh) * | 2016-07-07 | 2019-12-17 | 西安电子科技大学 | 基于区域特性的红外与可见光图像融合方法 |
CN106548467A (zh) * | 2016-10-31 | 2017-03-29 | 广州飒特红外股份有限公司 | 红外图像和可见光图像融合的方法及装置 |
CN106846289B (zh) * | 2017-01-17 | 2019-08-23 | 中北大学 | 一种红外光强与偏振图像融合方法 |
CN106846288A (zh) * | 2017-01-17 | 2017-06-13 | 中北大学 | 一种双模态红外图像差异特征指数测度的多算法融合方法 |
CN106846288B (zh) * | 2017-01-17 | 2019-09-06 | 中北大学 | 一种双模态红外图像差异特征指数测度的多算法融合方法 |
CN106846289A (zh) * | 2017-01-17 | 2017-06-13 | 中北大学 | 一种基于显著性迁移与细节分类的红外光强与偏振图像融合方法 |
CN107451984A (zh) * | 2017-07-27 | 2017-12-08 | 桂林电子科技大学 | 一种基于混合多尺度分析的红外与可见光图像融合算法 |
CN107451984B (zh) * | 2017-07-27 | 2021-06-22 | 桂林电子科技大学 | 一种基于混合多尺度分析的红外与可见光图像融合算法 |
CN110363731B (zh) * | 2018-04-10 | 2021-09-03 | 杭州海康微影传感科技有限公司 | 一种图像融合方法、装置及电子设备 |
CN110363731A (zh) * | 2018-04-10 | 2019-10-22 | 杭州海康威视数字技术股份有限公司 | 一种图像融合方法、装置及电子设备 |
CN108764249A (zh) * | 2018-04-23 | 2018-11-06 | 云南民族大学 | 一种旋转不变多源影像局部特征描述方法、系统和装置 |
CN108764249B (zh) * | 2018-04-23 | 2021-07-09 | 云南民族大学 | 一种旋转不变多源影像局部特征描述方法、系统和装置 |
CN108961154A (zh) * | 2018-07-13 | 2018-12-07 | 福州大学 | 基于改进非下采样轮廓波变换的太阳电池热斑检测方法 |
CN108961154B (zh) * | 2018-07-13 | 2022-12-23 | 福州大学 | 基于改进非下采样轮廓波变换的太阳电池热斑检测方法 |
CN109035189A (zh) * | 2018-07-17 | 2018-12-18 | 桂林电子科技大学 | 基于柯西模糊函数的红外与弱可见光图像融合方法 |
CN109035189B (zh) * | 2018-07-17 | 2021-07-23 | 桂林电子科技大学 | 基于柯西模糊函数的红外与弱可见光图像融合方法 |
CN109489576A (zh) * | 2018-10-19 | 2019-03-19 | 杭州电子科技大学 | 一种基于初级视觉通路计算模型的轮廓检测方法 |
CN109489576B (zh) * | 2018-10-19 | 2020-03-03 | 杭州电子科技大学 | 一种基于初级视觉通路计算模型的轮廓检测方法 |
CN110211082A (zh) * | 2019-05-31 | 2019-09-06 | 浙江大华技术股份有限公司 | 一种图像融合方法、装置、电子设备及存储介质 |
CN110211082B (zh) * | 2019-05-31 | 2021-09-21 | 浙江大华技术股份有限公司 | 一种图像融合方法、装置、电子设备及存储介质 |
CN110796632A (zh) * | 2019-07-30 | 2020-02-14 | 重庆小富农康农业科技服务有限公司 | 猪只计数装置 |
CN110796632B (zh) * | 2019-07-30 | 2023-08-11 | 重庆渝通合数字科技有限公司 | 猪只计数装置 |
CN111462028A (zh) * | 2020-03-16 | 2020-07-28 | 中国地质大学(武汉) | 基于相位一致性和目标增强的红外和可见光图像融合方法 |
CN111462028B (zh) * | 2020-03-16 | 2023-06-27 | 中国地质大学(武汉) | 基于相位一致性和目标增强的红外和可见光图像融合方法 |
CN113947554A (zh) * | 2020-07-17 | 2022-01-18 | 四川大学 | 一种基于nsst和显著信息提取的多聚焦图像融合方法 |
CN113947554B (zh) * | 2020-07-17 | 2023-07-14 | 四川大学 | 一种基于nsst和显著信息提取的多聚焦图像融合方法 |
CN112102217A (zh) * | 2020-09-21 | 2020-12-18 | 四川轻化工大学 | 一种可见光图像与红外图像快速融合方法及系统 |
CN112102217B (zh) * | 2020-09-21 | 2023-05-02 | 四川轻化工大学 | 一种可见光图像与红外图像快速融合方法及系统 |
CN116503454A (zh) * | 2023-06-27 | 2023-07-28 | 季华实验室 | 红外和可见光图像融合方法、装置、电子设备及存储介质 |
CN116503454B (zh) * | 2023-06-27 | 2023-10-20 | 季华实验室 | 红外和可见光图像融合方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105719263B (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105719263B (zh) | 基于nsct域底层视觉特征的可见光和红外图像融合方法 | |
CN107194904B (zh) | 基于增补机制和pcnn的nsct域图像融合方法 | |
CN104809734B (zh) | 一种基于引导滤波的红外图像与可见光图像融合的方法 | |
Yang et al. | Wavelet u-net and the chromatic adaptation transform for single image dehazing | |
Liu et al. | Multi-modality medical image fusion based on image decomposition framework and nonsubsampled shearlet transform | |
CN108399611B (zh) | 基于梯度正则化的多聚焦图像融合方法 | |
CN113837974B (zh) | 一种基于改进beeps滤波算法的nsst域电力设备红外图像增强方法 | |
CN104616274B (zh) | 一种基于显著性区域提取的多聚焦图像融合方法 | |
CN104657965B (zh) | 基于离散连续曲波的偏振图像融合方法 | |
CN111462028B (zh) | 基于相位一致性和目标增强的红外和可见光图像融合方法 | |
CN111179208B (zh) | 基于显著图与卷积神经网络的红外-可见光图像融合方法 | |
Wan et al. | Practical remote sensing image fusion method based on guided filter and improved SML in the NSST domain | |
CN105894483A (zh) | 一种基于多尺度图像分析和块一致性验证的多聚焦图像融合方法 | |
CN102740114A (zh) | 一种视频主观质量的无参评估方法 | |
CN114120176A (zh) | 远红外与可见光视频图像融合的行为分析方法 | |
CN105184819A (zh) | 用于医学图像重建参数寻优的客观图像质量评价方法 | |
Kanimozhi et al. | Brain MR image segmentation using self organizing map | |
Duan et al. | Infrared and visible image fusion using multi-scale edge-preserving decomposition and multiple saliency features | |
CN108648174A (zh) | 一种基于自动对焦技术的多层图像融合方法及系统 | |
CN109961408A (zh) | 基于nsct和块匹配滤波的光子计数图像去噪算法 | |
CN116091917A (zh) | 一种基于光谱扩展方法的农作物目标检测方法及系统 | |
CN117575923A (zh) | 基于局部主方向多尺度的多模态内镜图像配准融合方法 | |
CN103049910B (zh) | 基于nsct与多通道互信息相结合的图像融合质量评估方法 | |
CN103778615A (zh) | 基于区域相似性的多聚焦图像融合方法 | |
CN109360194B (zh) | 一种基于离散不可分离剪切波变换和人眼视觉特性的图像质量评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210219 Address after: 650000 room 1701, 17th floor, block a, science and Technology Information Innovation Incubation Center, Chenggong District, Kunming City, Yunnan Province Patentee after: YUNNAN UNITED VISUAL TECHNOLOGY Co.,Ltd. Address before: 650093 No. 253, Xuefu Road, Wuhua District, Yunnan, Kunming Patentee before: Kunming University of Science and Technology |
|
TR01 | Transfer of patent right |