Infrared and Visible Image Fusion Based on Saliency Detection and Latent Low-Rank Representation
-
摘要: 针对红外与可见光图像融合过程中细节信息的缺失、融合结果对比度较低等问题，提出一种基于显著性检测与多层潜在低秩表示的红外与可见光图像融合方法。首先，使用基于显著性检测的方法对红外与可见光图像进行预融合；然后，使用多层潜在低秩表示方法依次将红外图像、可见光图像和预融合图像分解为低秩层和细节层；其中细节层采用结构相似性和L2范数相结合的方法进行融合，低秩层使用基于能量属性的方法进行融合；最后，将低秩层和细节层的融合结果重构便得到最终的融合图像。文中将该方法与11种具有代表性的图像融合方法进行了评估比较，通过对比多组融合图像的主客观评价，其结果表明，相较于对比方法，本方法能够保留红外与可见光图像融合过程中源图像的有效细节，且融合结果具有较高的对比度，更符合人们的视觉理解。Abstract: To address the problems of missing detail and low contrast in the fusion of infrared and visible images, this study proposes a fusion method based on saliency detection and latent low-rank representation. First, a pre-fusion image is obtained by saliency detection for the infrared and visible images. Then, the infrared, visible, and pre-fused images are decomposed into low-rank and detail layers by the multilevel latent low-rank representation method. The detail layer is fused by combining the hyperspherical L2 norm and structural similarities, while the low-rank layer is fused using an approach based on the energy property. The final fused image is obtained by adding the fusion results of the low-rank and detail layers. The proposed method is compared with 11 representative image fusion methods by conducting subjective and objective evaluations of multiple groups of fused images. The results show that the image fusion method enhances the effective detail information and improves the image contrast, yielding a fusion result that is more in line with people's visual understanding.
-
Keywords:
- image fusion /
- saliency detection /
- latent low-rank representation /
- infrared image /
- visible image
-
0. 引言
随着导弹拦截技术的不断发展，拦截弹具有更快的飞行速度、更高的制导精度和更强的机动响应能力，给进攻方弹道目标带来了严重的威胁[1]。同时，拦截弹在攻方防区外发射，攻方阵地探测设备难以实现对拦截弹的远距离全方位探测，目标信息获取难度大[2]。因此，利用弹载平台对拦截弹进行探测跟踪，从而进行规避、突防，是一种必要可行的办法。而红外探测器体积小、重量轻、功耗小且不受战场复杂电磁环境影响，能在保持较大作用距离的同时还保证了平台的隐蔽性，不论是在性能仿真还是在实测应用中都取得了较好的效果，是弹载平台首选的探测手段[3-4]。为了验证这种手段的有效性，需要根据拦截弹的红外辐射特点，给出弹载红外探测器性能参数，在此基础上分析其对拦截弹的探测能力。
目前，大部分文献的研究主要集中在防御方采用地基、空基平台对进攻方飞行器、巡航导弹等武器的红外辐射建模以及作用距离分析上[5-9]，或者是对进攻方弹道导弹中段突防策略的研究上[10-11]。而有关拦截弹探测的研究文献较少，文献[12]分析了大气层外动能拦截弹的特点，给出了弹载红外告警装置对大气层外动能杀伤拦截器（Exoatmospheric Kill Vehicle, EKV）探测所需的探测谱段、探测视场、探测距离等指标，但文中只讨论了对中段防御拦截弹的探测需求，没有对末段防御拦截弹的探测分析。同时，前述文献中涉及到的红外探测器性能参数大多是直接给出，探测波段覆盖整个中波红外（3~5 μm）或长波红外（8~12 μm），而实际上红外探测器的响应波段并不能完全覆盖，与制备的红外探测器成品参数不符。此外，前述文献对于探测能力的分析主要体现在红外探测器对一定辐射强度目标的极限探测距离上，即作用距离；而实际上探测能力还体现在红外探测器对一定距离、一定红外辐射强度目标的检测能力上，即探测概率[13]。
因此，在已有文献的研究基础上，本文根据拦截弹的红外辐射特点，结合国内外红外探测器的实际研制进展，推测了弹载红外探测器的性能参数。同时，给出了一定虚警概率下的探测概率和输入信噪比的拟合公式，并根据基于辐射差的信噪比模型推导了探测概率模型和作用距离模型。最后，调研了末段、中段防御拦截弹的性能指标参数，分析了相应的探测需求，并通过仿真计算了探测概率及作用距离，分析了弹载红外探测器的探测能力水平。
1. 红外探测器性能参数估计
为满足弹载红外探测器对拦截弹的探测，通过分析拦截弹的红外辐射特性，并参考国内外红外探测器研制现状以及导引头、空天探测平台等典型红外探测器参数，对适用于拦截弹探测的红外探测器参数进行分析。
1）探测波段选择
拦截弹的红外辐射主要来自弹体蒙皮、尾喷口和尾焰及废气柱这3个方面。在助推段拦截弹的红外辐射主要来自尾喷口及尾焰，根据维恩位移定律可知，红外辐射集中在短波红外、中波红外；在中、末制导段的红外辐射主要来自弹体及姿轨控发动机，红外辐射集中在中波红外、长波红外。由于中波红外可以兼顾助推段、制导段的探测，同时为了保证对深空背景下拦截弹弹头的稳定探测，因此弹载红外探测器采用中波红外和长波红外的探测波段。
2）红外探测器类型及尺寸
根据国内外的研制进展和实际应用情况，考虑前述探测波段需求，可以采用中/长波双色红外探测器，常用的探测器材料为碲镉汞或二类超晶格。目前，主流中、长波红外探测器面阵大小为320×256、640×512、1280×512，像元尺寸为10~30 μm。考虑红外探测器成本，同时兼顾探测视场及冷屏效率，凝视红外焦平面阵列可以采用640×512像元。
3）红外探测器参数
目前，国外Raytheon、SELEX、Lynred等公司对中/长波双色探测器的设计、生产及应用已十分成熟，国内昆明物理研究所、中电科光科技有限公司、高德红外公司等也已取得丰硕的成果，其中华北光电技术研究所已成功研制出前述分辨率的制冷型红外探测器阵列[14]，采用二类超晶格材料，工作温度为70 K，像元尺寸为20 μm×20 μm，响应波段为3.5~4.8 μm和7.5~9.5 μm，采用边积分边读取的工作方式，本文后续对拦截弹探测的分析将基于该探测器进行。帧频取为100 Hz[15]，积分时间为8 ms。根据中/长波双色探测器研制情况，3.5~4.8 μm波段的比探测率可以达到3×109 m⋅Hz1/2⋅W-1，7.5~9.5 μm波段的比探测率可以达到7×108 m⋅Hz1/2⋅W-1。
4）红外探测口径及焦距
由于弹载平台多为突防弹头，空间尺寸有限，参考国内外主要型号弹道导弹，弹头一般长度为1.5~3 m，直径为0.5~1 m。而在像元尺寸d确定的情况下，焦距f（Focal length）越大，则瞬时视场角ω（Instantaneous Field Of View, IFOV）越小，空间分辨率越高。因此，在满足弹载平台尺寸限制及光学分辨率的前提下，焦距应取一个较为合适的值。当像元尺寸为20 μm时，瞬时视场角随焦距变化的曲线如图 1所示。可以看出，当焦距大于0.6 m时，焦距的变化对瞬时视场角的变化影响很小，故可以选取焦距的范围为0.4~0.6 m，本文选取f=0.5 m。
现有的制冷型红外探测器F数一般为2或4。为了保证100%的冷光阑效率，光学系统的F数要与探测器的冷屏F数相匹配。为了获取更多目标辐射能量，从而增加探测距离，应采用较小的F数，即F=2。根据F数的定义可知，光学口径D=f/F=0.25 m。综上，弹载红外探测器的性能参数估计如表 1所示。
表 1 弹载红外探测器参数Table 1. Parameters of missile-borne infrared detectorParameters MWIR LWIR Operating band /μm 3.5~4.8 7.5~9.5 Resolution 640×512 640×512 Pixel size/μm 20 20 Integration time/ms 8 8 Specific detectivity D*/(m·Hz1/2·W-1) 3×109 7×108 Focal length/m 0.5 0.5 Optical aperture/m 0.25 0.25 IFOV/μrad 40 40 2. 红外探测能力分析
2.1 红外探测器探测概率
红外探测器的探测概率主要是由输入信噪比和阈值信噪比决定的。通常阈值信噪比是一个固定值，红外探测器的探测概率Pd与输入信噪比fSNR有如下关系：
$$ {P_{\text{d}}} = \frac{1}{{\sqrt {{\text{2}\pi}} }}\int_{{f_{{\text{TNR}}}} - {f_{{\text{SNR}}}}}^{ + \infty } {\exp \left( { - \frac{{{x^2}}}{2}} \right){\text{d}}x} $$ (1) 式中：fTNR为探测概率为50%时的阈值信噪比，在虚警概率为Pfa=10-5时取值为4.26。由此计算不同输入信噪比条件下的探测概率如表 2所示。
表 2 红外探测器探测概率与信噪比关系Table 2. The relationship between detection probability and SNR of infrared detectorPd 1.0 0.95 0.9 0.8 0.7 0.6 fSNR 7.56 5.91 5.55 5.1 4.79 4.52 Pd 0.5 0.4 0.3 0.2 0.1 0.01 fSNR 4.26 4.01 3.74 3.42 2.98 1.94 利用MATLAB对表中数据进行拟合，得到探测概率Pd与输入信噪比fSNR的函数表达式：
$$ \begin{aligned} P_{\mathrm{d}}= & 0.1578 \sin \left(1.237 f_{\mathrm{SNR}}-5.274\right)+ \\ & 8.428 \sin \left(-0.0227 f_{\mathrm{SNR}}+3.179\right) \end{aligned} $$ (2) 通常红外探测器对远距离目标进行探测时，探测器上所成的像不足或仅占一个像元。但由于光学系统衍射、大气抖动造成的弥散效应，使得目标成像能量分散，实际成像会占据多个像元。假设目标成像占N个像元且能量分布均匀，则基于辐射差的信噪比模型为[16]：
$$ \begin{gathered} {f_{{\text{SNR}}}} = \delta \left( {\left| {{{\left( {{L_{\text{t}}} - {L_{\text{b}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{L_{\text{t}}} - {L_{\text{b}}}} \right)} N}} \right. } N}} \right|{A_{\text{t}}}} \right){\tau _{\text{a}}}{\tau _0}{A_0}\frac{{{D^ * }}}{{\sqrt {{A_{\text{d}}}\Delta f} {R^2}}} \hfill \\ {\text{ = }}\frac{{\delta {A_0}{\tau _{\text{a}}}{\tau _0}{D^ * }}}{{\sqrt {{A_{\text{d}}}\Delta f} }} \cdot \frac{{\Delta I}}{{{R^2}}} \hfill \\ \end{gathered} $$ (3) 式中：δ为信号提取因子；Lt为目标的辐射亮度；Lb为目标所在背景的辐射亮度；N为目标所占像元数目；At为目标的有效辐射面积；τa为大气透过率；τ0为光学系统透过率；A0=πD2/4为光学系统入射孔径面积；D为光学系统的孔径；D*为红外探测器比探测率；Ad为探测器像元面积；Δf为等效噪声带宽；R为探测距离；ΔI=|(Lt-Lb)/N|⋅At，为单个像元对应的红外辐射强度。
通过式(3)、式(2)可知探测概率Pd、探测距离R和单个像元对应红外辐射强度ΔI之间的函数关系为：
$$ \begin{gathered} {P_{\text{d}}} = 0.1578 \cdot {\text{sin}}\left( {1.237\frac{{\delta {A_0}{\tau _{\text{a}}}{\tau _0}{D^ * }}}{{\sqrt {{A_{\text{d}}}\Delta f} }} \cdot \frac{{\Delta I}}{{{R^2}}} - 5.274} \right) + \hfill \\ \quad \quad 8.428 \cdot {\text{sin}}\left( { - 0.0227\frac{{\delta {A_0}{\tau _{\text{a}}}{\tau _0}{D^ * }}}{{\sqrt {{A_{\text{d}}}\Delta f} }} \cdot \frac{{\Delta I}}{{{R^2}}} + 3.179} \right) \hfill \\ \end{gathered} $$ (4) 2.2 红外探测器作用距离
当红外探测器距离目标一定距离时，所接收到的目标红外辐射恰好满足信号提取的要求，此时刚好能把目标从背景中提取出来，该距离称为红外探测器的作用距离。由式(3)可以得到作用距离公式如下：
$$ R = {\left[ {\delta \cdot \left( {\left| {{{\left( {{L_{\text{t}}} - {L_{\text{b}}}} \right)} \mathord{\left/ {\vphantom {{\left( {{L_{\text{t}}} - {L_{\text{b}}}} \right)} N}} \right. } N}} \right| \cdot {A_{\text{t}}} \cdot {\tau _{\text{a}}}} \right) \cdot \left( {{A_0}{\tau _0}} \right) \cdot {D^*} \cdot \frac{1}{{{\text{SNR}}\sqrt {\Delta f \cdot {A_{\text{d}}}} }}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} $$ (5) 此时，SNR为满足系统探测的最低信噪比。
3. 仿真分析
为了分析弹载红外传感器的探测能力，根据前文估计的弹载红外探测器参数，从探测概率、探测距离对典型拦截弹的探测能力进行分析，相关导弹防御系统的性能指标如表 3所示[17]。
表 3 导弹防御系统性能指标Table 3. Performance indicators of missile defense systemParameters Terminal defense system Midcourse defense system Missile airframe size/m 6×0.37 16.6×1.3 Warhead size/m 2.3×0.37 1.4×0.6 Capture distance/km 50 200 Intercept distance/km 200 5000 Intercept speed/(km/s) 2.5~2.8 2.7~3.2 Intercept height/km 40~150 2000 3.1 探测概率分析
3.1.1 参数设置
对末段防御拦截弹探测的大气透过率采用PCModWin软件计算，对中段防御拦截弹探测的大气透过率τa=1，光学系统透过率τ0=0.9，其余参数由表 1中取值或经计算得到。对于目标所占像元数目N，为保证95%以上的探测概率，根据Johnson判别准则[18]，像元数目至少要大于3×3，故N=9。对于等效噪声带宽Δf，计算公式为：
$$ \Delta f = \frac{1}{{2{\tau _{{int} }}}} $$ (6) 式中：τint为探测器积分时间，表 1中取8 ms，则Δf=62.5 Hz。
3.1.2 对末段防御拦截弹探测概率分析
根据表 3中拦截系统的性能指标，对于相对距离较近的末段防御系统，希望弹载红外探测器尽可能在拦截弹发射后便能够探测发现目标，因此探测距离要能达到拦截弹的最大拦截距离，即200 km。经计算，3.5~4.8 μm波段的大气透过率为0.73，7.5~9.5 μm波段的大气透过率为0.72，故在探测距离需求的约束下，两个探测波段内单个像元对应的不同红外辐射强度与探测概率计算结果如图 2所示。
![]()
图 2 末段防御拦截弹的探测概率Figure 2. Detection probability of terminal defense interception missile由计算结果可知，在3.5~4.8 μm波段内，对拦截弹的探测概率随单个像元对应红外辐射强度的增加迅速变大，当辐射强度达到0.70 W/sr时，探测概率达到100%。在7.5~9.5 μm波段内，拦截弹探测概率的增加速度弱于3.5~4.8 μm波段，当单个像元对应的红外辐射强度达到3.04 W/sr时，探测概率达到100%。而末段防御拦截弹在助推段起始时，尾焰及废气柱是主要的红外辐射源，温度可以达到900~1000 K，有效辐射面积变化范围为0.9~13 m2。按照尾焰及废气柱的平均发射率为0.5计算，3.5~4.8 μm波段的红外辐射强度为1220~17623.5 W/sr，单个像元对应的红外辐射强度为136~1958 W/sr；7.5~9.5 μm波段的红外辐射强度为443.7~6409.5 W/sr，单个像元对应的红外辐射强度为49.3~712.2 W/sr。在两个探测波段，单个像元对应的红外辐射强度远大于探测概率为100%时所需的最小辐射强度，可以保证在末段防御拦截弹的最远拦截距离处实现发射后即被探测。
3.1.3 对中段防御拦截弹探测概率分析
对于相对距离较远的中段防御系统，在外部信息的支援引导下，希望弹载红外探测器尽可能在被拦截弹导引头发现前探测到来袭的目标，从而进行机动规避，因此探测距离至少要能达到拦截弹导引头的捕获距离，即200 km。在探测距离需求的约束下，两个波段内不同辐射强度的拦截弹探测概率计算结果如图 3所示。
![]()
图 3 中段防御拦截弹的探测概率Figure 3. Detection probability of intermediate defense interception missile由计算结果可知，在3.5~4.8 μm波段，当单个像元对应的红外辐射强度达到0.51 W/sr时，探测概率达到100%。在7.5~9.5 μm波段，当单个像元对应的红外辐射强度达到2.19 W/sr时，探测概率达到100%。由于中段防御的拦截弹飞行在大气层外，弹头蒙皮是主要的红外辐射源，在日照区时，弹头表面温度受太阳直接辐射的影响较大，而在日影区时，由于弹头所在的位置被地球遮挡而无法被太阳直接辐射，因此在进行分析时，应该区分拦截弹弹头是位于日影区还是位于日照区。通常拦截弹弹头在大气层外飞行弹头会很快到达热平衡，在日照区约为400 K，在日影区约为200 K[19]。根据表 3中中段防御拦截弹弹头外型参数计算有效辐射面积为0.28~0.81 m2，蒙皮发射率取0.8，则在日照区时，3.5~4.8 μm波段的红外辐射强度为4.87~14.1 W/sr，单个像元对应的红外辐射强度为0.54~1.57 W/sr，由图 3可知探测概率为100%；7.5~9.5 μm波段的红外辐射强度为17.61~50.94 W/sr，单个像元对应的红外辐射强度为2.0~5.6 W/sr，探测概率大于98%。因此，在日照区能够保证在拦截弹导引头发现前探测到拦截弹。
在日影区时，3.5~4.8 μm波段的红外辐射强度为0.0012~0.0036 W/sr，7.5~9.5 μm波段的红外辐射强度为0.25~0.73 W/sr，与日照区的红外辐射相差较大。根据图 3可知，两个波段内的探测概率接近于0。考虑日影区弹头红外辐射强度变化不大，红外辐射强度取有效辐射面积中值计算，探测距离范围取20~200 km，经计算，在3.5~4.8 μm波段，探测概率恒为0，即中红外波段在日影区基本无探测作用。在7.5~9.5 μm波段的探测距离与探测概率的关系如图 4所示。
![]()
图 4 7.5~9.5 μm波段探测距离与探测概率的关系Figure 4. The relationship between detection distance and detection probability in the 7.5-9.5 μm band由图 4可以看出，探测概率在100 km处随相对距离的逐渐减小而呈指数级增加。当相对距离为35 km时，探测概率达到90%。根据文献[12]，弹载平台实施突防时的相对距离不小于100 km，此时的最短反应时间约为12.3 s。因此，在日影区时，弹载平台无法保证探测距离满足突防需求。为了增加在日影区的探测距离，可以采取牺牲一定的探测概率和虚警概率，从而提高作用距离，相关内容将在3.2.3节进行分析。
3.2 作用距离分析
3.2.1 参数设置
采用式(5)进行作用距离分析，相关参数参考3.1.1节。对于最小可探测信噪比SNR，主要由虚警概率和探测概率确定，表 4给出了一定虚警概率及探测概率下对应的信噪比。为了保证对目标探测的可靠性，假设要求虚警概率不大于10-4，从表中可以看出，当SNR=6时，探测概率已大于95%，可以满足对目标的探测要求，故取最小可探测信噪比为6。
表 4 不同探测概率与虚警概率下的最小可探测信噪比Table 4. Minimum detectable SNR under different detection probabilities and false alarm probabilitiesSNR 5 6 7 Pfa 10-3 10-4 10-5 10-3 10-4 10-5 10-3 10-4 10-5 Pd 0.972 0.900 0.769 0.998 0.989 0.958 0.999 0.999 0.997 3.2.2 对末段防御拦截弹探测距离分析
对于末段防御的拦截弹，重点分析3个探测阶段，一是在拦截弹助推段起始时，需要分析此时对拦截弹的尽早尽远发现能力；二是在拦截弹助推段结束后，弹头与弹体分离，有效辐射面积突然减小，需要分析此时对弹头的探测能力；三是拦截弹进入制导段高速飞向目标，需要分析此时对拦截弹的持续探测能力。
1）第一个探测阶段
拦截弹助推段起始时，飞行高度较低，背景辐射主要来自地面背景、大气背景。设拦截弹发射5 s时，高度约为1.3 km，速度约为520 m/s，采用文献[20-21]中的方法对拦截弹及背景的红外辐射特性进行建模，并计算红外辐射结果如表 5所示。
表 5 拦截弹助推段起始时仿真参数Table 5. Simulation parameters at the beginning of the boost phase of interceptor missileParameters Temperature/K Radiating area/m2 Infrared radiation/(W/m2·sr) 3.5~4.8μm 7.5~9.5μm Envelope 390 0.11(head-on)
2.2(side-looking)14.02 56.32 Tail nozzle 1800 - 20009.53 3202.99 Tail flame and exhaust column 990 0.86(head-on)
12.98(side-looking)1944.42 605.09 Ground background 298 - 1.17 17.31 Atmospheric background 290 - 0.90 15.50 从表 5可以看出，尾喷口在两个探测波段中的红外辐射亮度最大，但辐射面积小，且只能在追尾时才能探测到，因此蒙皮、尾焰及废气柱是主要红外辐射源。尤其是尾焰及废气柱的辐射面积大，红外辐射亮度明显，且在3.7~4.8 μm波段占据主导。在此阶段，对拦截弹的探测视角通常为一定角度下的俯视探测，以迎头时为0°，侧视时为90°，3.5~4.8 μm波段的大气透过率为0.73，7.5~9.5 μm波段的大气透过率为0.72，计算作用距离如图 5所示。
![]()
图 5 助推起始时对末段防御拦截弹的作用距离Figure 5. The operating distance of the terminal defense interception missile at the beginning of the booster thrust可以看出，3.7~4.8 μm波段的作用距离远大于7.5~9.5 μm波段，这是由于尾焰及废气柱的红外辐射主要集中在中波红外。两个探测波段在迎头时的作用距离最小，均远大于拦截弹的最远拦截距离，能够做到尽早尽远发现目标，这也与3.1.2节的分析结果相符。
2）第二个探测阶段
拦截弹助推段结束时，助推器停止工作，尾焰及废气柱消失，同时弹头与弹体分离，此时弹头成为主要的红外辐射源，有效辐射面积突然减小。设助推段结束时弹头高度约为11 km，速度达到最大，约为7Ma左右，采用前述方法计算弹头及背景的红外辐射结果如表 6所示。
表 6 拦截弹助推段结束时仿真参数Table 6. Simulation parameters at the end of the boost phase of interceptor missileParameters Temperature /K Radiating area /m2 Infrared radiation/(W/m2·sr) 3.5~4.8μm 7.5~9.5μm Envelope 1958.77 0.11(head-on)
0.84(side-looking)21494.21 3240.26 Ground background 298 - 1.17 17.31 Atmospheric background 230 - 0.05 3.39 从表 6中可以看出，高速弹头因气动加热效应温度达到1958 K，在3.7~4.8 μm的红外辐射亮度最大。在此阶段，对拦截弹的探测视角仍为俯视探测，3.5~4.8 μm波段的大气透过率为0.82，7.5~9.5 μm波段的大气透过率为0.91，作用距离如图 6所示。
![]()
图 6 助推结束时对末段防御拦截弹的作用距离Figure 6. The operating distance of the terminal defense interception missile at the end of the boost可以看出，助推结束时虽然有效辐射面积减小，但弹头蒙皮气动加热明显，红外辐射强度明显大于助推起始时，此时中波红外辐射仍占据主导，3.7~4.8 μm波段的作用距离远大于7.5~9.5 μm波段。两个波段的作用距离均能满足对此时弹头的探测。
3）第三个探测阶段
进入制导段后，弹头在姿轨控发动机的控制下以4~7Ma的速度飞向目标，速度逐渐减小，取平均速度5Ma，飞行高度为25 km，此时弹头蒙皮气动加热效应仍十分明显，不考虑姿轨控发动机脉冲点火的影响，同理计算弹头及背景的红外辐射结果如表 7所示。
表 7 拦截弹制导段仿真参数Table 7. Simulation parameters at the guidance phase of interception missileParameters Temperature/K Radiating area /m2 Infrared radiation/(W/m2·sr) 3.5~4.8μm 7.5~9.5μm Envelope 1130 0.11(head-on)
0.84(side-looking)4948.77 1269.52 Ground background 298 - 1.17 17.31 Atmospheric background 225 - 0.03 2.88 从表 7中可以看出，在制导段弹头温度约为1130 K，红外辐射强度较上一阶段明显减弱。在此阶段，对拦截弹的探测视角为迎头探测或俯视探测，3.5~4.8 μm波段的大气透过率为0.85，7.5~9.5 μm波段的大气透过率为0.97，作用距离如图 7所示。
![]()
图 7 对制导段末段防御拦截弹的作用距离Figure 7. The operating distance of the terminal defense interception missile at the guidance phase可以看出，在拦截弹制导段时，两个波段的最小作用距离分别达到6684 km和1747 km，远大于弹载平台与拦截弹的相对距离，因此可以保证对拦截弹的持续探测。
3.2.3 对中段防御拦截弹作用距离分析
中段拦截弹的弹头一般飞行在80 km以上的空间中，在助推段结束时，其表面温度最高，在进入制导段后，弹头温度快速下降并达到热平衡，弹头及背景的红外辐射结果如表 8所示。
表 8 拦截弹制导段仿真参数Table 8. Simulation parameters at the guidance phase of interception missileParameters Temperature/K Radiating area /m2 Infrared radiation/(W/m2·sr) 3.5~4.8μm 7.5~9.5μm Envelope 400(sunshine)
200(sun shadow)0.28(head-on)
0.81(side-looking)17.41
0.004562.89
0.91Deep space background 4 — 0 $ 9.978 \times {10^{ - 164}} $ 从表 8中可以看出，深空背景辐射可以忽略不计，蒙皮的红外辐射在7.5~9.5 μm波段占据优势，尤其是当弹头处于日影区时，3.5~4.8 μm波段的红外辐射低于7.5~9.5 μm波段两个数量级。在此阶段，对拦截弹的探测视角为迎头探测或俯视探测，作用距离如图 8所示。
![]()
图 8 对中段防御拦截弹的作用距离Figure 8. The operating distance for intermediate defense interception missile可以看出，在日照区，两个探测波段的作用距离在各个探测视角下均大于200 km，可以保证被拦截弹导引头发现前探测到来袭的目标。而在日影区，3.5~4.8 μm波段的作用距离几乎可以忽略不计，而7.5~9.5 μm波段的作用距离小于前文提到的最短突防距离。
为了保证弹载平台突防的探测距离需求，可以牺牲一定的虚警概率和探测概率，减小最小可探测信噪比SNR和所需像元数目N，从而提高作用距离。由于日影区弹头目标的红外辐射峰值波长位于长波波段，因此只讨论7.5~9.5 μm波段的作用距离。
由图 9(a)和图 9(b)对比可以看出，随着虚警概率的增加，相同探测概率对应的最小可探测信噪比减小。而随着探测概率的减小，所需像元数目也相应减小，作用距离随之增加。对于虚警概率为10-5，探测概率为0.5时，探测视角大于42°所对应的作用距离大于100 km；探测概率小于0.4时，探测视角大于22°所对应的作用距离大于100 km。对于虚警概率为10-4，探测概率为0.5时，探测视角大于19°所对应的作用距离大于100 km；探测概率为0.4时，探测视角大于9°所对应的作用距离大于100 km。综上，在一定的虚警概率下，只有探测视角在一定范围时，才能够保证在100 km以上的距离有40%~50%的概率发现拦截弹头。
![]()
图 9 7.5~9.5 μm波段作用距离分析Figure 9. Analysis of the operating range in the 7.5-9.5 μm band4. 结束语
本文推测了用于拦截弹探测的弹载红外探测器参数，并对其探测能力进行了仿真分析，研究结果对红外探测器能力分析、装备论证及作战应用具有一定的借鉴意义。但是对弹载红外探测器能力的分析是一个复杂且精细的问题，本文在目标及背景红外辐射特性、大气透过率等作了一定的简化计算，下一步若能对此进行更精确的分析和建模，将进一步提高本文方法的准确性。
-
表 1 数据集‘Camp’的融合结果客观评价指标
Table 1 Objective evaluation indicators of fusion results on dataset 'Camp'
Methods AG CE EI EN MI QCV SD SF SSIM VIF Time/s CBF 4.862 0.464 51.008 6.548 1.062 510.614 27.418 11.644 1.260 0.284 7.716 CNN 4.569 0.239 47.485 7.102 1.394 349.232 37.269 11.506 1.414 0.469 23.71 GFCE 6.747 1.569 70.022 7.202 1.295 603.566 38.121 15.649 1.123 0.510 0.386 GFF 4.045 0.478 42.519 6.524 1.163 444.798 27.393 10.554 1.433 0.272 0.086 HMSD_GF 5.204 0.379 53.686 6.936 1.356 416.013 36.358 13.164 1.412 0.540 0.363 IVFusion 8.720 0.859 91.261 7.355 0.948 796.213 42.972 20.731 1.054 1.030 14.22 LP 5.371 0.294 56.580 7.289 1.167 601.295 45.316 12.937 1.389 0.703 0.339 MDLatLRR 8.398 0.250 88.341 7.192 1.094 653.114 42.825 20.818 1.252 1.161 52.55 NSST_PCNN 4.135 0.601 43.284 6.649 1.189 295.207 29.113 10.722 1.471 0.408 62.50 TIF 4.394 0.672 46.186 6.639 1.307 323.957 29.249 11.043 1.513 0.557 0.025 WLS 4.628 0.543 47.058 6.603 1.091 400.285 28.251 11.782 1.461 0.452 0.572 Proposed 9.549 0.372 100.363 7.489 1.143 183.030 54.758 22.880 1.019 0.977 12.42 
下载: 导出CSV
表 2 数据集‘Kaptein’的融合结果客观评价指标
Table 2 Objective evaluation indicators of fusion results on dataset 'Kaptein '
Methods AG CE EI EN MI QCV SD SF SSIM VIF Time/s CBF 5.599 1.672 58.013 6.878 1.295 319.046 36.413 13.255 1.079 0.327 16.24 CNN 4.493 1.073 45.282 7.312 1.710 145.364 57.149 11.534 1.365 0.533 55.98 GFCE 8.015 1.734 76.406 7.584 1.418 553.108 54.349 19.889 0.888 0.559 2.327 GFF 4.133 1.103 41.579 7.129 3.341 374.493 52.447 10.947 1.353 0.381 0.523 HMSD_GF 5.147 1.250 51.318 7.253 1.619 129.287 53.568 13.133 1.359 0.602 0.916 IVFusion 8.774 0.843 87.343 7.258 0.744 243.631 43.409 22.083 1.004 0.894 70.26 LP 5.113 1.012 51.605 7.369 1.663 504.230 58.894 12.689 1.309 0.697 0.533 MDLatLRR 8.798 1.145 89.001 7.196 0.885 353.961 50.387 22.812 1.064 0.839 145.5 NSST_PCNN 4.256 2.007 42.916 6.739 1.427 124.723 43.423 11.163 1.403 0.452 124.1 TIF 4.218 1.964 43.059 6.613 1.058 152.787 33.914 10.801 1.392 0.489 0.044 WLS 4.728 1.865 46.134 6.743 1.398 157.646 42.077 12.158 1.395 0.469 3.834 Proposed 9.630 1.186 95.924 7.171 1.169 528.788 70.799 25.193 0.927 0.913 63.67
下载: 导出CSV
表 3 数据集‘Marne’的融合结果客观评价指标对比
Table 3 Objective evaluation indicators of fusion results on dataset 'Marne'
Methods AG CE EI EN MI QCV SD SF SSIM VIF Time/s CBF 5.510 1.172 57.215 6.957 0.741 859.46 31.876 12.906 0.918 0.210 16.31 CNN 3.013 0.640 29.321 7.411 1.362 1066.3 45.254 7.576 1.357 0.606 53.88 GFCE 5.366 1.011 47.239 7.444 0.939 1296.6 42.973 14.087 1.062 0.659 2.414 GFF 2.513 0.582 22.849 7.113 3.002 443.47 40.391 7.010 1.362 0.119 0.570 HMSD_GF 3.429 0.712 33.045 7.450 1.687 1372.6 45.971 8.704 1.345 1.065 2.442 IVFusion 6.547 0.639 59.878 7.782 1.547 1707.7 59.069 16.922 0.985 1.144 39.02 LP 3.346 0.404 32.674 7.368 1.519 405.41 45.814 8.056 1.244 0.519 0.485 MDLatLRR 5.746 0.725 56.777 7.254 0.857 899.58 38.720 14.517 1.146 0.962 149.2 NSST_PCNN 2.936 0.799 28.508 7.230 1.673 808.54 38.255 7.431 1.394 0.567 126.4 TIF 2.684 1.162 27.618 6.839 1.004 1293.7 28.530 6.678 1.381 0.475 0.036 WLS 3.861 0.750 36.552 7.182 1.271 882.73 39.416 10.166 1.339 0.691 3.642 Proposed 7.360 0.590 65.360 7.805 1.437 459.2 65.821 20.057 0.898 0.702 42.86
下载: 导出CSV
表 4 数据集‘Airplane in Trees’的融合结果客观评价指标对比
Table 4 Objective evaluation indicators of fusion results on dataset 'Airplane in Trees'
Methods AG CE EI EN MI QCV SD SF SSIM VIF Time/s CBF 3.934 1.298 42.149 6.522 1.367 1431.60 40.857 9.671 1.362 0.332 13.62 CNN 2.484 1.443 26.972 6.654 2.391 181.686 59.796 6.373 1.654 0.546 40.85 GFCE 5.283 1.094 57.008 7.387 1.893 302.870 61.272 11.670 1.187 1.142 1.601 GFF 2.431 1.502 26.516 6.494 2.257 964.541 50.078 6.179 1.652 0.429 0.402 HMSD_GF 2.732 1.137 29.404 6.661 2.076 198.739 59.171 7.196 1.652 0.587 1.669 IVFusion 5.419 1.733 58.824 7.115 0.824 739.613 35.871 12.672 1.246 1.304 36.75 LP 2.799 1.255 30.534 6.987 1.689 381.239 62.178 6.667 1.626 0.776 0.176 MDLatLRR 4.843 1.032 52.793 7.129 1.233 402.472 45.285 12.417 1.401 1.047 173.6 NSST_PCNN 2.379 1.224 25.960 6.683 1.934 189.736 48.505 6.016 1.676 0.451 150.2 TIF 2.379 1.430 25.784 6.505 1.442 366.907 34.257 6.189 1.670 0.462 0.048 WLS 2.404 1.221 25.689 6.611 1.885 224.499 46.689 5.920 1.687 0.489 2.484 Proposed 5.897 2.373 64.020 7.162 1.501 260.734 81.195 15.109 1.039 1.453 32.75
下载: 导出CSV
-
[1] MA J, TANG L, XU M, et al. STD FusionNet: An infrared and visible image fusion network based on salient target detection[J]. IEEE Transactions on Instrumentation and Measurement, 2021, 70: 1-13.
[2] MA J, MA Y, LI C. Infrared and visible image fusion methods and applications: A survey[J]. Information Fusion, 2019, 45: 153-178. DOI: 10.1016/j.inffus.2018.02.004
[3] YAN H, ZHANG J X, ZHANG X. Injected infrared and visible image fusion via l1 decomposition model and guided filtering[J]. IEEE Transactions on Computational Imaging, 2022, 8: 162-173. DOI: 10.1109/TCI.2022.3151472
[4] 唐霖峰, 张浩, 徐涵, 等. 基于深度学习的图像融合方法综述[J]. 中国图象图形学报, 2023, 28(1): 3-36. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTB202301002.htm TANG L F, ZHANG H, XU H, et al. Deep learning-based image fusion: a survey[J]. Journal of Image and Graphics, 2023, 28(1): 3-36. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTB202301002.htm
[5] ZHANG X. Deep learning-based multi-focus image fusion: A survey and a comparative study[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2021, 44(9): 4819-4838.
[6] REN L, PAN Z, CAO J, et al. Infrared and visible image fusion based on weighted variance guided filter and image contrast enhancement[J]. Infrared Physics & Technology, 2021, 114: 103662.
[7] LI G, LIN Y, QU X. An infrared and visible image fusion method based on multi-scale transformation and norm optimization[J]. Information Fusion, 2021, 71: 109-129. DOI: 10.1016/j.inffus.2021.02.008
[8] TAN X, GUO L. Visible and infrared image fusion based on visual saliency detection[C]//2020 19th International Symposium on Distributed Computing and Applications for Business Engineering and Science (DCABES). IEEE, 2020: 134-137.
[9] Bavirisetti D P, Dhuli R. Two-scale image fusion of visible and infrared images using Sali ency detection[J]. Infrared Physics & Technology, 2016, 76: 52-64.
[10] 霍星, 邹韵, 陈影, 等. 双尺度分解和显著性分析相结合的红外与可见光图像融合[J]. 中国图象图形学报, 2021, 26(12): 2813-2825. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTB202112005.htm HUO X, ZOU Y, CHEN Y, et al. Dual-scale decomposition and saliency analysis based infrared and visible image fusion[J]. Journal of Image and Graphics, 2021, 26(12): 2813-2825 https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTB202112005.htm
[11] LIU G, YAN S. Latent low-rank representation for subspace segmentation and feature extraction[C]//2011 International Conference on Computer Vision. IEEE, 2011: 1615-1622.
[12] 孙彬, 诸葛吴为, 高云翔, 等. 基于潜在低秩表示的红外和可见光图像融合[J]. 红外技术, 2022, 44(8): 853-862. http://hwjs.nvir.cn/article/id/7fc3a60d-61bb-454f-ad00-e925eeb54576 SUN B, ZHUGE W W, GAO Y X, et al. Infrared and visible image fusion based on latent low-rank representation[J]. Infrared Technology, 2022, 44(8): 853-862. http://hwjs.nvir.cn/article/id/7fc3a60d-61bb-454f-ad00-e925eeb54576
[13] LI H, WU X J, Kittler J. MDLatLRR: A novel decomposition method for infrared and visible image fusion[J]. IEEE Transactions on Image Processing, 2020, 29: 4733-4746.
[14] TAN W, Tiwari P, Pandey H M, et al. Multimodal medical image fusion algorithm in the era of big data[J/OL]. Neural Computing and Applications, 2020, https://doi.org/10.1007/s00521-020-05173-2.
[15] WANG Z, Bovik A C, Sheikh H R, et al. Image quality assessment: from error visibility to structural similarity[J]. IEEE Transactions on Image Processing, 2004, 13(4): 600–612.
[16] 张蕾, 金龙旭, 韩双丽, 等. 采用非采样Contourlet变换与区域分类的红外和可见光图像融合[J]. 光学精密工程, 2015, 23(3): 810-818. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201503027.htm Zhang L, Jing L X, Han S L, et al. Fusion of infrared and visual images based on non-sampled Contourlet transform and region classification[J]. Optics and Precision Engineering, 2015, 23(3): 810-818. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201503027.htm
[17] Alexander Toet. TNO Image Fusion Dataset[EB/OL]. 2014. https://figshare.com/articles/dataset/TNO_Image_Fusion_Dataset/1008029.
[18] Shreyamsha Kumar B K. Image fusion based on pixel significance using cross bilateral filter[J]. Signal, Image and Video Processing, 2015, 9(5): 1193-1204.
[19] LIU Y, CHEN X, CHENG J, et al. Infrared and visible image fusion with convolutional neural networks[J]. International Journal of Wavelets, Multiresolution and Information Processing, 2018, 16(3): 1850018.
[20] ZHOU Z, DONG M, XIE X, et al. Fusion of infrared and visible images for night-vision context enhancement[J]. Applied Optics, 2016, 55(23): 6480-6490.
[21] LI S, KANG X, HU J. Image fusion with guided filtering[J]. IEEE Transactions on Image Processing, 2013, 22(7): 2864-2875.
[22] Burt P J, Adelson E H. The laplacian pyramid as a compact image code[J]. IEEE Transactions on Communications, 1983, 31(4): 532-540.
[23] MA J, ZHOU Z, WANG B, et al. Infrared and visible image fusion based on visual saliency map and weighted least square optimization[J]. Infrared Physics & Technology, 2017, 82: 8-17.
[24] ZHANG X, YE P, XIAO G. VIFB: A visible and infrared image fusion benchmark[C]//Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops, 2020: 104-105.
计量
- 文章访问数: 182
- HTML全文浏览量: 40
- PDF下载量: 40
