Infrared and Visible Image Fusion Based on NSCT Combined with Saliency Map and Region Energy
-
摘要: 针对传统的红外与可见光图像融合出现的清晰度和对比度偏低,目标不够突出的问题,本文提出了一种基于Non-subsampled Contourlet(NSCT)变换结合显著图与区域能量的融合方法。首先,使用改进的频率调谐(Frequency-tuned, FT)方法求出红外图像显著图并归一化得到显著图权重,单尺度Retinex(Single-scale Retinex, SSR)处理可见光图像。其次,使用NSCT分解红外与可见光图像,并基于归一化显著图与区域能量设计新的融和权重来指导低频系数融合,解决了区域能量自适应加权容易引入噪声的问题;采用改进的“加权拉普拉斯能量和”指导高频系数融合。最后,通过逆NSCT变换求出融合图像。本文方法与7种经典方法在6组图像中进行对比实验,在信息熵、互信息、平均梯度和标准差指标中最优,在空间频率中第一组图像为次优,其余图像均为最优结果。融合图像信息量丰富、清晰度高、对比度高并且亮度适中易于人眼观察,验证了本文方法的有效性。
-
关键词:
- 图像融合 /
- Non-subsampled Contourlet变换 /
- 区域能量自适应加权 /
- 拉普拉斯能量和
Abstract: To address the problems of low clarity and contrast of indistinct targets in traditional infrared and visible image-fusion algorithms, this study proposes a fusion method based on non-subsampled contourlet transform (NSCT) combined with a saliency map and region energy. First, an improved frequency-tuning (FT) method is used to obtain the infrared image saliency map, which is subsequently normalized to obtain the saliency map weight. A single-scale retinex (SSR) algorithm is then used to enhance the visible image. Second, NSCT is used to decompose the infrared and visible images, and a new fusion weight is designed based on the normalized saliency map and region energy to guide low-frequency coefficient fusion. This solves the problem of region-energy adaptive weighting being prone to introducing noise, and the improved "weighted Laplace energy sum" is used to guide the fusion of high-frequency coefficients. Finally, the fused image is obtained by inverse NSCT. Six groups of images were used to compare the proposed method with seven classical methods. The proposed method outperformed others in terms of information entropy, mutual information, average gradient, and standard deviation. Regarding spatial frequency, the first group of images was second best, and the remaining images exhibited the best results. The fusion images displayed rich information, high resolution, high contrast, and moderate brightness, demonstrating suitability for human observation, which verifies the effectiveness of this method. -
0. 引言
光电探测电路是激光模拟实兵对抗系统的重要组成部分,基于光电效应原理,实现光电转换与信号处理[1-2]。在设计中,高灵敏度、低噪声、低失真和高增益需要重点考虑,其接收灵敏度与信噪比对系统的性能具有重要影响[3]。而电路的噪声是影响信噪比的主要因素,探测器在进行光电信号转换过程中会伴随着噪声电流(电压)的存在,经放大电路放大,噪声也随之放大,从而严重影响光电探测器的效率[4-6]。
当射击距离较长,激光脉冲信号到达光电探测器时变得微弱,甚至被噪声完全埋没,从而使得探测器对光信号的探测本领受到限制。因此有效抑制环境光产生的噪声干扰成为光电探测电路设计中一个非常重要的挑战[7-8]。系统保持接收信号强度稳定是实现高通信速率、低误码率的必要条件[9-10]。为了确保输出信号稳定,方便后续处理,需要信号处理电路能够相应的改变放大器增益,提高信噪比与对微弱信号检测的能力[11-14]。
Li等人[15]利用对数放大器和积分放大器的功能,有效地提取信号,设计了一种微弱光信号检测系统。He等人[16]基于LabVIEW设计了一种微弱光电信号采集与处理系统。Wu等人[17]提出了一种锁相放大技术来检测微弱信号,具有一定应用价值。
实兵交战系统中的光电探测属于微弱信号探测领域,所需放大倍数高,对电路低噪声要求高。为解决环境光以及噪声对系统探测性能的影响,本研究基于信号差分和自动增益控制原理设计了一种性能优良的光电探测电路。采用对称跨阻放大回路结构消除暗电流噪声与环境光引起的直流偏移;通过自动增益控制电路对信号进行自适应放大,使信号幅值在正常范围内保持稳定,提高信号采集效率;采用多次滤波的方式去除信号的直流分量并能够较好地抑制背景噪声;借助主控芯片的ADC模块对信号进行采集,提取信号中的信息。
1. 方案设计
在激光模拟实兵对抗系统中,光电探测器采用光伏特性PIN光电二极管,为准确获取激光传递的信息,需设计与之相匹配的检测电路,达到信号高效检测的目的[18-20]。整体结构如图 1所示,首先,光电探测器接收激光信号,并将其转换为电流信号,通过跨阻放大器将光电探测器输出的电流信号转换为电压信号同时进行放大;其次,为了扩大系统增益,在跨阻放大电路后增加一级放大电路对信号进行放大;随后将放大后的电压信号接入带通滤波电路进行降噪;最后,经单片机将带通滤波电路输出的模拟信号转换为数字信号做采样处理。
为消除光电探测器的暗电流噪声与环境光引起的直流偏移,系统采用两路对称的跨阻放大器,一路接光电探测器正极,另一路接光电探测器负极,将两路信号通过差分放大电路进行叠加得到脉冲电压信号,从而实现对光电探测器噪声与跨阻放大器噪声抑制;采用自动增益控制(automatic gain control,AGC)电路解决由于射击距离不同导致接收到信号强弱不一致的问题,在输入信号幅值变化很大的情况下,使输出信号仅在较小范围内变化,保证了信号输出的稳定性;采用带通滤波器可以滤除所需频率范围之外的干扰,提高信噪比。
2. 光电探测器选型
光电探测器的探测模式主要有两种,光电导模式和光伏模式。光电导式需外加偏置电压,适合于高速光电转换的应用。光电压式存在入射光时才会产生光电流,适合精确测量[21-22]。本文采用光伏式,可有效降低噪声的影响,从而实现对微弱信号的精确探测。光电探测器是根据光电转换定律将光辐射量转化为光电流量。其光电转化模型如下:
$$ i(t)=\frac{e \eta}{h \nu} P(t) $$ (1) 式中;i(t)为产生的光电流;P(t)为光功率;hν为光子能量;η为探测器量子效率;e是电子电荷;可以看出,光电流主要受光功率与探测器的量子效率的影响,当探测器确定时,探测器的量子效率是一个常值,故光电流只与光功率有关。
不同光电探测器性能不同,因此选择一个合适的光电探测器是实现光电检测的前提[23]。光电探测器选型要充分考虑灵敏度、光谱利用率系数、频率等,同时还应考虑无规则噪声的影响。这里选用深圳市诚强光电数码有限公司生产的CHQ0038L型PIN光电探测器接收激光信号并实现光电转化。该光电探测器具有接收角度宽、距离长、可靠性高、抗干扰能力强、光功率转化效率显著的特点。对噪声或干扰信号引起的杂散输出脉冲具有一定的抑制作用。其性能参数如表 1所示。
表 1 光电探测器性能参数Table 1. Performance parameters of photoelectric detectorType CHQ0038L Spectral response/nm 500-1100 Central wavelength λp/nm 960 Short circuit current/μA 63 Dark current/nA 1 Operating voltage/V 2.7-5.5 Sensitivity@980 nm /(A/W) 0.59 Angle of arrival/° ±45 3. 光电探测电路设计
由于激光在传输过程中受环境天气等各种因素影响,能量在大气中的衰减严重,此时光电探测器接收到的光功率微弱,容易受其他信号的干扰[24-26]。因此,设计合理的信号处理电路,选择匹配的电气参数是提高电路抗干扰能力与稳定性的关键[27-29]。
3.1 前置放大电路设计
为使光电探测器的输出信号易于检测,首先需要考虑如何将光电探测器输出的电流信号转换为易于检测的电压信号。这里采用TI公司生产的OPA380型运算放大器实现光电信号的I-V转换,该放大器具有低噪声、低偏置电流和高速特性,是宽带跨阻抗和电压放大器应用的理想解决方案,其电路如图 2所示。
该前置放大电路由3部分构成,分别是前面两个对称的跨阻放大电路和将两个跨阻放大电路信号进行叠加的差值电路。两个跨阻放大电路的输出通过积分器反馈给正输入端,可有效补偿由于太阳光或其他外部光源引起的直流输出偏移。采用差分配置更好地抑制由于耦合干扰引起的噪声,抵消共模干扰信号,达到提高电路信噪比与探测精度的目的。
为了补偿光电二极管的内部结电容带来的振荡问题以及RC滤波后导致的相位迟滞问题,电路采用补偿电容C3,C18与反馈电阻R并联用于滤除信号中的高频噪声以防止电路自激振荡。电路中并联两个相反方向的肖特基二极管D1, D2改变跨阻放大器不同方向的反馈量,起到施密特效果。为了提高电路的灵敏度,降低噪声,反馈电阻采用T型网络布局,通过减小寄生电容来实现更大的跨阻增益。
3.2 自动增益控制电路
系统为了扩展入射光强的范围,在前置放大电路后加入自动增益控制电路,使电路增益随信号强弱自适应调整,放大电路输出信号经滤波器滤除截止频率为33 kHz的低频调制分量和噪声后,产生用以控制增益受控放大器的电压。使得在射击距离近,激光能量强时,放大电路增益减小;射击距离远,激光能量弱时,放大电路增益增大。始终将信号保持在0.5~5 V之间,便于后续采集与处理。STC15F2K60S2主控芯片通过脉冲峰值的大小来判断放大电路的增益是否处于合适档位,继而控制放大电路改变增益,从而实现自增益控制。
电路选用TI公司的一款可变增益芯片VCA810,这是一款高增益,可调范围宽带的压控放大器,具有很好的带宽和噪声特性,是精密AGC应用的理想选择,电路如图 3所示。AGC的环路部分采用OPA820作为电压比较器,AD8307用于幅值检测构成对数检波环路。比较器将VCA810输出与DAC输入的电压进行比较,输出线性电压,经过电容滤除直流分量再输入至AD8307检波环路,对VCA810增益进行自动控制。电压比较模块中电阻R7避免高频时产生的自激震荡,电容C15作为相位补偿构成AGC反馈回路。当输入信号小的时候提高VCA810的增益,在信号大的时候降低增益,但不改变响应的基本对数性质。其增益与控制电压的关系可表示为:
$$ G_{\mathrm{dB}}=-40 \times\left(V_{\mathrm{c}}+1\right) $$ (2) 式中:Vc为检波器AD8307的输出电压,通过单片机控制电压Vc从-2~0 V变化,从而改变放大器增益GdB在-40~40 dB范围内线性变化。当VCA810的输出大于DAC输入引脚的电压时,OPA820输出正向摆动,电压保持电容上的电压随之上升,VCA810的增益控制引脚电压随之上升,从而降低VCA810的增益。当VCA810的输出小于DAC输入引脚的电压时,OPA820输出负向摆动,VCA810的增益控制引脚电压随之下降,从而升高VCA810的增益。由此种负反馈机制的电路,完成了自动增益控制功能,以满足信号处理对距离衰减效应的补偿的要求,保持很好的线性输出。
3.3 选频滤波电路设计
系统信号中心频率为38 kHz,低于或高于这个频率都会对信号产生影响,为了提高信噪比消除噪声的干扰,这里对信号进行选频滤波使得信号的质量更佳。其电路如图 4所示,采用双运放LM324芯片组成多路反馈拓扑结构的四阶巴特沃斯带通滤波电路滤除通频带以外的噪声信号,方便后续对信号进行采集处理。
确定选频滤波电路中心频率为38 kHz,增益为0 dB,经仿真软件TINA-TI对信号与噪声进行分析可知通频带Δf越小,电路的输出噪声越小。但减小通频带,会影响信号的波形。综合考虑这两方面因素,最终确定滤波电路-3 dB通频带带宽为10 kHz,取C1=C2=C3=C4=C=10 nF。根据文献[30]计算得到元器件R1, R2, R3, R4, R5, R6的值。确定好各元器件参数后仿真得到带通滤波电路的频率响应曲线如图 5所示,结果显示该滤波器在38 kHz附近增益曲线平坦,滤波效果优良,性能稳定,能够很好地保持信号波形,抑制噪声。
3.4 控制电路设计
整个激光模拟对抗光电探测系统采用STC公司的STC15F2K60S2型号单片机作为控制单元,具有高速、高可靠、低功耗、超强抗干扰等优点。通过单片机内置的A/D转换模块把检测到的模拟信号转换为计算机能够识别的数字量,并通过高低电平时间的长短对信号进行解码获得激光传输的信息,然后通过D/A转换变为电压或电流信号,驱动执行元件。
3.5 软件设计
主控芯片STC15F2K60S2的控制程序采用keil uVision5软件编写实现,系统程序采用C语言编写。单片机控制部分主要包括系统初始化、信号采集、DA输出、与上位机通信等,主要流程如图 6所示。
系统在运行前初始化,比较输入信号幅值的强弱来确定DA输出的电压值,进一步对自动增益控制环路的增益进行控制,保证输入单片机的信号幅值在规定范围内,便于对信号进行采集。将采集到的信号经过A/D转换之后进行解码获得信号传输的信息,将信息经D/A转化为电压信号驱动执行元件响应。
4. 实验结果与分析
在完成各个子模块的设计后,为提高实验结果的准确性与说服力,信号处理模块采用PCB来实现,所设计的PCB如图 7(a)所示,将设计好的PCB进行制板焊接搭建实验平台,如图 7(c)所示。图 7(b)为激光发射端实物图,所采用的激光波长为980 nm。
为了验证激光接收电路对接收信号的处理效果,对电路性能进行实验分析,图 8为电路的波形图,其中图 8(a)为本研究输入输出波形,图 8(b)为传统系统输入输出波形。由图 8(a)看到输入输出波形具有很好的保持,说明系统性能稳定可靠。对比图 8(a), (b)电路优化前后系统噪声抑制效果明显,输出信号稳定,误码率降低30%,信号的传输效果得到改善。
理想情况下,增益自适应控制电路输出信号的幅度应在不同输入信号强度下保持不变。为了验证电路增益的自适应性能,通过改变激光管的驱动电压(或是改变发射与接收端之间的距离)观察输出波形,并测量光电接收电路输出信号幅值,与无自动增益光电检测电路进行对比如图 9所示。从图上可以看到移动距离在1~10 m范围内,采用固定增益,不同移动距离下输出信号的幅度变化很大,且存在距离过小输出信号饱和,距离过大输出信号弱,从而导致后续信号处理时高低电平出现误判的现象。而采用增益自适应控制电路输出信号幅度基本保持恒定,可有效解决固定增益存在的信号误判问题。
利用TINA-TI仿真测试光电探测电路在增益最大时的频率特性曲线如图 10所示,可以看出电路在38 kHz频率附近曲线平坦,说明在此范围内电路性能稳定。幅频特性曲线的通频带增益大小与输入信号强度有关,但是其曲线形状、上限截止频率却与信号强弱无关,输入信号越弱,幅频特性曲线增益越大。
由于整个光电探测电路中光电探测器和前置放大电路的噪声影响比较大,为了验证前置放大电路的降噪效果,采用TINA-TI软件对前置放大电路进行仿真分析,得到该电路与传统前置放大电路的输出噪声曲线与信噪比曲线如图 11、12所示,可以看到在频率为38 kHz处电路噪声为3.46 μV,信噪比为109 dB,相较于传统前置放大电路输出噪声降低效果显著,系统噪声得到明显抑制,纹波干扰得到消除,电路信噪比得到提高,方便信号的进一步处理,对微弱信号的探测性能得到提升。
为验证电路的选频滤波效果,通过TINA-TI软件的扫频功能改变输入信号的频率,测得电路在不同频率下的输出特性曲线如图 13所示,从图中可以看到频率越接近中心频率38 kHz,其波形保持的越完整,证明该电路具有很好的选频滤波效果,可有效滤除低频与高频噪声,保证信号的输出质量。
5. 结论
本文依据信号叠加原理与自动增益控制理论设计并实现了一种增益自适应光电探测电路,采用两个对称的跨阻放大电路叠加的结构,有效抑制了系统噪声与暗电流,采用自动增益控制电路对信号进行放大,避免因距离太近导致信号饱和,实现对光电信号的高增益、低噪声放大。通过实验分析与TINA-TI仿真结果表明:系统可实现对激光脉冲信号进行采集、转换和放大,设计能够保证输出信号幅值基本稳定在2 V左右,提高探测系统的信噪比,具有良好的稳定性和噪声抑制效果,电路的最低噪声为3.46 μV,信噪比为109 dB。探测能力得到大幅提升,满足实际应用需求,为光电检测技术的研究提供了参考。下一步设想在激光模拟实兵对抗系统中得到应用。
-
表 1 Nato_camp客观评价结果
Table 1 Objective evaluation results of Nato_camp
Methods IE SF MI AG SD LP 6.6022 13.9059 1.8726 5.6680 28.3882 RP 6.5073 14.7712 1.7414 5.8658 27.3435 DTCWT 6.4292 12.8851 1.8237 5.1054 25.6576 CVT 6.3256 11.8559 1.7603 4.7871 24.0154 NSCT 6.6369 13.5828 1.8622 5.6284 28.6006 DCHWT 6.3231 11.0041 1.7605 4.4255 25.8427 Hybrid_MSD 6.7020 15.1560 2.0825 6.0049 29.1284 Proposed 7.0235 14.8216 2.6280 6.2154 37.1931 表 2 Tree评价结果
Table 2 Objective evaluation results of Tree
Methods IE SF MI AG SD LP 5.8484 7.7347 1.6911 3.0154 14.7225 RP 5.8788 9.0936 1.5948 3.2999 15.1486 DTCWT 5.7484 7.1703 1.6836 2.7127 13.5645 CVT 5.7172 6.8121 1.6988 2.5939 13.0307 NSCT 5.8498 7.4456 1.6643 2.9577 14.6721 DCHWT 6.0325 6.1470 1.5695 2.4196 16.3405 Hybrid_MSD 6.2954 8.5220 1.8939 3.2627 20.2859 Proposed 6.8678 11.3701 3.1874 4.8995 31.3029 表 3 Duine客观评价结果
Table 3 Objective evaluation results of Duine
Methods IE SF MI AG SD LP 5.9780 7.4565 1.5175 3.3497 15.4549 RP 5.8923 6.9660 1.5174 3.0050 14.5825 DTCWT 5.8808 6.7566 1.5059 3.0241 14.4632 CVT 5.8322 6.5233 1.4829 2.9103 14.0198 NSCT 5.9806 7.0812 1.5181 3.2809 15.5069 DCHWT 5.8015 5.7144 1.5240 2.5810 13.6918 Hybrid_MSD 5.9472 8.2933 1.5553 3.6047 15.2660 Proposed 7.1126 13.2301 3.1155 6.1093 35.0350 表 4 APC_4客观评价结果
Table 4 Objective evaluation results of APC_4
Methods IE SF MI AG SD LP 5.8574 12.7384 0.9065 5.5624 14.5848 RP 5.6555 12.4937 0.7502 5.0838 12.9853 DTCWT 5.6766 11.7014 0.8259 5.0618 12.8683 CVT 5.5498 11.3298 0.7902 4.8656 11.6973 NSCT 5.8743 12.1900 0.8733 5.4444 14.7461 DCHWT 5.5262 9.7761 0.8671 4.2252 11.4534 Hybrid_MSD 5.9514 14.0328 1.0230 5.9508 15.5852 Proposed 6.6951 16.5562 2.8247 7.5079 25.6268 表 5 Kaptein_1654客观评价结果
Table 5 Objective evaluation results of Kaptein_1654
Methods IE SF MI AG SD LP 6.6557 19.0519 2.2753 6.8009 36.8878 RP 6.7122 19.7509 2.2324 6.9286 34.3867 DTCWT 6.4858 17.9294 2.2036 6.3067 31.2040 CVT 6.3880 16.3414 2.2711 5.6601 28.7014 NSCT 6.6451 18.7667 2.2189 6.7824 35.7430 DCHWT 6.7437 15.4424 2.2980 5.5617 36.3196 Hybrid_MSD 6.8692 20.4999 2.1659 7.1597 37.6306 Proposed 7.0479 20.6661 3.0611 7.7455 53.3902 表 6 Movie_18客观评价结果
Table 6 Objective evaluation results of Movie_18
Methods IE SF MI AG SD LP 5.9545 8.3646 1.7206 3.0967 17.9892 RP 5.7520 9.0084 1.5003 3.1423 16.3032 DTCWT 5.6640 7.8037 1.6294 2.8594 14.5863 CVT 5.5991 7.6734 1.6439 2.8605 13.9172 NSCT 5.8977 8.2244 1.6706 3.0956 17.1809 DCHWT 5.5755 6.6402 1.6596 2.4455 15.0237 Hybrid_MSD 6.2333 9.2791 2.3197 3.3676 21.1429 Proposed 6.8108 10.0577 3.3102 3.7343 47.7950 -
[1] 杨孙运, 奚峥皓, 王汉东, 等. 基于NSCT和最小化-局部平均梯度的图像融合[J]. 红外技术, 2021, 43(1): 13-20. http://hwjs.nvir.cn/article/id/144252d1-978c-4c1e-85ad-e0b8d5e03bf6 YANG Sunyun, XI Zhenghao, WANG Handong, et al. Image fusion based on nsct and minimum-local mean gradient[J]. Infrared Technology, 2021, 43(1): 13-20. http://hwjs.nvir.cn/article/id/144252d1-978c-4c1e-85ad-e0b8d5e03bf6
[2] 肖中杰. 基于NSCT红外与可见光图像融合算法优化研究[J]. 红外技术, 2017, 39(12): 1127-1130. http://hwjs.nvir.cn/article/id/hwjs201712010 XIAO Zhongjie. Improved infrared and visible light image fusion algorithm based on NSCT[J]. Infrared Technology, 2017, 39(12): 1127-1130. http://hwjs.nvir.cn/article/id/hwjs201712010
[3] ZHANG B, LU X, PEI H, et al. A fusion algorithm for infrared and visible images based on saliency analysis and non-subsampled Shearlet transform[J]. Infrared Physics & Technology, 2015, 73: 286-297.
[4] ZHOU Z, WANG B, LI S, et al. Perceptual fusion of infrared and visible images through a hybrid multi-scale decomposition with Gaussian and bilateral filters[J]. Information Fusion, 2016, 30: 15-26. DOI: 10.1016/j.inffus.2015.11.003
[5] MA J, YU W, LIANG P, et al. Fusion GAN: A generative adversarial network for infrared and visible image fusion[J]. Information Fusion, 2019, 48: 11-26. DOI: 10.1016/j.inffus.2018.09.004
[6] ZHANG Y, ZHANG L, BAI X, et al. Infrared and visual image fusion through infrared feature extraction and visual information preservation[J]. Infrared Physics & Technology, 2017, 83: 227-237.
[7] Mertens T, Kautz J, Reeth F V. Exposure fusion[C]//Proceedings of Pacific Conference on Computer Graphics and Applications, 2007: 382-390.
[8] ZHANG Z, Blum R S. A categorization of multiscale-decomposition-based image fusion schemes with a performance study for a digital camera application[C]//Proc. IEEE, 1999, 87(8): 1315-1326.
[9] Lewis J J, Callaghan R J O, Nikolov S G, et al. Pixel- and region-based image fusion with complex wavelets[J]. Inf. Fus., 2007, 8(2): 119-130. DOI: 10.1016/j.inffus.2005.09.006
[10] Nencini F, Garzelli A, Baronti S, et al. Remote sensing image fusion using the curvelet transform[J]. Inf. Fus., 2007, 8(2): 143-156. DOI: 10.1016/j.inffus.2006.02.001
[11] YANG S, WANG M, JIAO L, et al. Image fusion based on a new contourlet packet[J]. Inf. Fus., 2010, 11(2): 78-84. DOI: 10.1016/j.inffus.2009.05.001
[12] Cunha A L, ZHOU J P, DO M N. The nonsubsampled Contourlet transform: theory, design, and applications [J]. IEEE Transactions on Image Processing, 2006, 15(10): 3089-3101. DOI: 10.1109/TIP.2006.877507
[13] 郭明, 符拯, 奚晓梁. 基于局部能量的NSCT域红外与可见光图像融合算法[J]. 红外与激光工程, 2012, 41(8): 2229-2235. GUO Ming, FU Zheng, XI Xiaoliang. Novel fusion algorithm for infrared and visible images based on local energy in NSCT domain[J]. Infrared and Laser Engineering, 2012, 41(8): 2229-2235.
[14] Cands E J, Donoho D L. Curvelets and curvilinear integrals[J]. J. Approximation Theor., 2001, 113(1): 59-90. DOI: 10.1006/jath.2001.3624
[15] DO M N, Vetterli M. Contourlets: a directional multiresolution image representation[C]//Proceedings of IEEE International Conference on Image Processing, 2002, 1: I-357-I-360.
[16] ZHAI Y, SHAH M. Visual attention detection in video se-quencesusing spatiotemporal cues[C]//Proceedings of the 14th ACM Conference on Multimedia, 2006: 815-824.
[17] Achanta R, Hemami S, Estrada F, et al. Frequency-tuned salient region detection[C]// IEEE Conference on Computer Vision & Pattern Recognition, 2009: 1597-1604.
[18] 叶坤涛, 李文, 舒蕾蕾, 等. 结合改进显著性检测与NSST的红外与可见光图像融合方法[J]. 红外技术, 2021, 43(12): 1212-1221. http://hwjs.nvir.cn/article/id/bfd9f932-e0bd-4669-b698-b02d42e31805 YE Kuntao, LI Wen, SHU Leilei, LI Sheng. Infrared and visible image fusion method based on improved saliency detection and non-subsampled shearlet transform[J]. Infrared Technology, 2021, 43(12): 1212-1221. http://hwjs.nvir.cn/article/id/bfd9f932-e0bd-4669-b698-b02d42e31805
[19] 唐中剑, 毛春. 基于显著导向的可见光与红外图像融合算法[J]. 太赫兹科学与电子信息学报, 2021, 19(1): 125-131. TANG Zhongjian, MAO Chun. Visible and infrared image fusion algorithm based on saliency guidance[J]. Journal of Terahertz Science and Electronic Information Technology, 2021, 19(1): 125-131.
[20] 王惠琴, 吕佳芸, 张伟. 基于双边滤波-BM3D算法的GPR图像去噪[J]. 兰州理工大学学报, 2022, 48(1): 91-97. WANG Huiqin, LYU Jiayun, ZHANG Wei. GPR image denoising based on bilateral filtering BM3D algorithm[J]. Journal of Lanzhou University of Technology, 2022, 48(1): 91-97.
[21] Edwin H Land. The retinex theory of color vision [J]. Scientific American, 1977, 237(6): 108-129. DOI: 10.1038/scientificamerican1277-108
[22] 程铁栋, 卢晓亮, 易其文, 等. 一种结合单尺度Retinex与引导滤波的红外图像增强方法[J]. 红外技术, 2021, 43(11): 1081-1088. http://hwjs.nvir.cn/article/id/b49a0a09-e295-40e6-9736-24a58971206e CHENG Tiedong, LU Xiaoliang, YI Qiwen, et al. Research on infrared image enhancement method combined with single-scale retinex and guided image filter[J]. Infrared Technology, 2021, 43(11): 1081-1088. http://hwjs.nvir.cn/article/id/b49a0a09-e295-40e6-9736-24a58971206e
[23] 翟海祥, 何嘉奇, 王正家, 等. 改进Retinex与多图像融合算法用于低照度图像增强[J]. 红外技术, 2021, 43(10): 987-993. http://hwjs.nvir.cn/article/id/23500140-4bab-40b3-9bef-282f14dc105e ZHAI Haixiang, HE Jiaqi, WANG Zhengjia, et al. Improved Retinex and multi-image fusion algorithm for low illumination image enhancemen[J]. Infrared Technology, 2021, 43(10): 987-993. http://hwjs.nvir.cn/article/id/23500140-4bab-40b3-9bef-282f14dc105e
[24] LIU Y, LIU S, WANG Z. A general framework for image fusion based on multi-scale transform and sparse representation[J]. Information Fusion, 2015, 24: 147-167.
[25] Kumar B. Multifocus and multispectral image fusion based on pixel significance using discrete cosine harmonic wavelet transform[J]. Signal, Image & Video Processing, 2013, 7: 1125-1143. DOI: 10.1007/s11760-012-0361-x
[26] ZHOU Z, BO W, SUN L, et al. Perceptual fusion of infrared and visible images through a hybrid multi-scale decomposition with Gaussian and bilateral filters[J]. Information Fusion, 2016, 30: 15-26.
[27] Toet Alexander. TNO image fusion dataset[EB/OL]. 2014, https://doi.org/10.6084/m9.figshare.1008029.v1.
[28] 张小利, 李雄飞, 李军. 融合图像质量评价指标的相关性分析及性能评估[J]. 自动化学报, 2014, 40(2): 306-315. ZHANG Xiaoli, LI Xiongfei, LI Jun. Validation and correlation analysis of metrics for evaluating performance of image fusion[J]. Acta Automatica Sinica, 2014, 40(2): 306-315.