Research on Infrared Target Characteristics of the "Ali Burke" Class Destroyer
-
摘要:
在海战对抗中，红外导弹是打击驱逐舰目标的重要武器。为了能有效论证红外导弹对驱逐舰目标的探测能力，首先需要获取典型驱逐舰红外目标特性，为此本文开展“阿利•伯克”级驱逐舰红外目标特性研究工作。建立了“阿利•伯克”级驱逐舰三维物理模型，通过数值仿真分别计算了模型在巡航状态下的流场特性与中、长波红外辐射特性。研究结果表明，军事舰船目标烟囱排放羽流温度最高达688.5 K，舰体温度较低且均匀。舰船目标中波光谱辐射主要由高温羽流贡献，长波光谱辐射主要由常温舰体贡献。且中、长波光谱辐射特征受观测角度影响较大，中波光谱辐射强度最高值在正侧方探测角，达65000 W/(sr·μm)，长波光谱辐射强度最高值在正上方探测角，达18200 W/(sr·μm)。
Abstract:In naval warfare, infrared missiles are important weapons for destroying targets. To effectively demonstrate the ability of infrared missiles to destroy targets, it is necessary to first obtain the infrared target characteristics of typical destroyers. Therefore, this study investigated the infrared target characteristics of "Ali Burke" class destroyers. A three-dimensional physical model of the "Ali Burke" class destroyer was established. The flow field characteristics and medium- and long-wave infrared radiation characteristics of the model in the cruise state were calculated via numerical simulations. The results indicate that the plume temperature of the target chimney of the military ship is as high as 688.5 K, and the temperature of the ship body is low and uniform. The medium-wave spectral radiation of ship targets is primarily contributed by high-temperature plumes, whereas the long-wave spectral radiation is primarily contributed by the normal-temperature ship body. In addition, the characteristics of medium- and long-wave spectral radiation are significantly affected by the observation angle. The maximum value of the medium-wave spectral radiation intensity is 65000 W/(sr·μm), and the highest value of the long-wave spectral radiation intensity, 18200 W/(sr·μm), is reached at the detection angle directly above.
-
Keywords:
- ship /
- infrared target characteristics /
- flow field /
- spectral radiation intensity
-
0. 引言
红外成像技术通过接收目标的红外辐射，将肉眼不可见的红外辐射能量转化为电信号进行成像，可以探测到隐藏目标或伪装目标，具有很好的抗干扰能力，可以全天时全天候工作；可见光成像技术通过捕捉物体的反射光进行成像，可以很好地获取物体的细节信息、颜色信息和纹理特征[1, 2]，但易受天气、光照的影响。因此，可见光与红外图像具有较强的互补特性，只有将这两类图像信息进行融合，综合两类图像的互补信息，才能够有效提高成像系统的性能[3-4]。
目前，可见光与红外图像融合方法大致可以分为多尺度变换融合[5-6]、稀疏表示融合[7-11]、混合模型融合[12-14]、拟态融合[15]、深度学习融合[16-17]等。多尺度变换融合方法采用尺度变换数学模型描述图像特征，采用相应的融合规则进行合并，得到最终的融合图像。由于尺度变换数学模型往往对某一类图像特征敏感，易造成融合图像信息缺失。稀疏表示融合方法利用字典原子表征图像特征，通过合并稀疏系数重构得到最终的融合图像。Yang等[11]首次提出基于稀疏表示的图像融合算法，随后，Yu等[10]提出联合稀疏融合方法，Liu等[11]提出自适应稀疏融合方法，提高了稀疏融合的效果。混合模型融合方法采用多种算法进行嵌入、组合，如多尺度稀疏融合[12]、多尺度嵌入融合[13-14]等，混合模型融合方法在一定程度上克服了单一融合方法的局限性。拟态融合方法[15]是利用可见光与红外图像的差异特征，驱动选择相应的融合算法，为图像融合方法提供了新的研究思路。深度学习融合方法利用神经网络自主学习图像特征，获得了较好的融合效果[16-17]。
稀疏表示在过完备字典上对图像信号进行表示，可以有效表征图像特征，是目前图像融合的热点之一。但传统稀疏表示图像融合方法将图像分割成重叠的块，字典训练及稀疏分解过程均以图像块为单位进行，块与块之间是相互独立的，造成图像特征不能被准确提取，融合效果较差[18]。图像中的局部位置和非局部位置存在相似结构，这种特征被称为非局部相似性，通过非局部相似性可以建立图像相似块之间的联系。非局部相似性在图像去噪、图像去模糊等方面取得较好效果[12, 13]。
可见光与红外图像融合就是要将两类图像的显著特征进行有效合并，而实际上，图像中的显著特征不仅与局部结构有关，还与图像非局部结构有关。从两类图像的成像机理和成像特点可以看出，可见光图像具有丰富的结构信息，图像块反映了局部的几何结构，相似的结构在图像不同位置重复出现，具有非局部结构相似性。而红外图像中，绝大多数为背景区域，其灰度变化缓慢，背景区域的图像块之间存在较强相关性，也具有明显的非局部结构相似性。为了直观显示，如图 1所示，通过计算图像结构块的相似性，在可见光与红外图像分别自动选取了7个图像相似块，并进行局部放大，可以看出，图像相似块具有明显的非局部相似性，这也说明，可见光与红外图像的局部位置及非局部位置上均存在结构相似的图像块。因此，通过建立图像局部结构与非局部结构的内在联系，将相似的图像块构建成图像结构组，以结构组代替传统的图像块，利用结构组的重复结构和冗余信息，通过对结构组进行字典学习和稀疏分解，可以使字典原子和稀疏系数更有效地描述图像显著特征，从而有效解决传统稀疏表示融合方法存在的问题。
![]()
图 1 可见光与红外图像非局部相似性示意图Figure 1. Schematic diagram of non-local similarity for visible and infrared images本文针对传统稀疏表示融合方法字典表征能力不足，稀疏系数不准确的问题，通过非局部相似性建立相似图像块之间的联系，以欧氏距离为度量准则，依据图像的非局部相似性特征将相似图像块采用向量化、矩阵化的方式构建为结构组矩阵，在此基础上，构建组K-SVD训练模型，以结构组代替传统图像块，优化结构组字典学习和稀疏分解问题，将源图像的相似结构特征与融合过程充分结合。该方法建立了图像块之间的内在联系，以结构组矩阵为单位进行字典训练和稀疏分解，增强了字典原子的表征能力和稀疏系数的准确性，提高了图像的融合性能。
1. 模型构建
1.1 结构组矩阵构造
结构组矩阵构造基于图像的非局部相似性特征，首先采用滑动窗将源图分割成大小相同的图像块，然后根据相似性度量准则，从图像的确定邻域内提取相似块，对相似块采取向量化、矩阵化的方式构建成结构组矩阵，用于字典训练和融合过程，具体过程如图 2所示。
1）图像块提取
假设源图像I1、I2的大小均为M×N像素，采用滑动窗，滑动步长为1像素，按照图中所示的轨迹，自左向右、自上而下进行滑动。分别将源图像分解为大小为n×n的(M-n+1)·(N-n+1)个图像块。
2）结构组构建
在L×L的邻域内，以欧式距离作为度量准侧，计算每一个初始图像块xi最相似的s个图像块。将初始图像块xi与s个相似图像块构成一个结构组Gi，这样每个结构组中共有s+1个图像块。对于M×N的输入源图像，共有(M-n+1)·(N-n+1)个初始块，则可以构建(M-n+1)·(N-n+1)个结构组。
3）图像块向量化
对于初始图像块，按照列向量的顺序，对图像块进行逐像素排列，作向量化处理，可以得到图像块的向量表示${v_i} \in {R^{{n^2}}}$，i=1, 2, …, (M-n+1)·(N-n+1)。
4）结构组矩阵化：将结构组中的s+1个图像块作为矩阵的列进行组合，得到每个结构组对应的结构组矩阵${{\boldsymbol{M}}_i} \in {R^{{n^2} \times \left( {s + 1} \right)}}$，i=1, 2, …, (M-n+1)·(N-n+1)。
1.2 组K-SVD字典训练模型
假设信号y∈Rn，字典D∈Rn×m，当m大于n时，称字典D为冗余字典。信号y可以表示为y=Dα或者y≈Dα。α∈Rm为稀疏表示系数，α中的多数元素为0，只有少数元素不为0，称α是稀疏的，y可以在字典D上稀疏表示。因此，稀疏表示就是将信号表示为一个冗余字典中少数原子的线性组合。由于实际应用中存在噪声，当考虑噪声时，信号y在字典D上的稀疏表示求解可以表示为：
$$\begin{array}{*{20}{c}} {\arg \mathop {\min }\limits_{\mathit{\boldsymbol{D}}, \alpha } \left\| {y - \mathit{\boldsymbol{D}}\alpha } \right\|_2^2}&{{\rm{s}}{\rm{.t}}.}&{{{\left\| \alpha \right\|}_0} \leqslant \varepsilon } \end{array}$$ (1) 式中：$y \in {R^n}$为信号向量；D∈Rn×m（m＞n）为过完备字典；‖‖0是L0范数，可以计算向量中非零项的个数；ε≥0是允许的误差阈值，预先设定。
对于结构组矩阵M，在对应的字典D上的稀疏表示优化问题可以表示为如下的形式：
$$\begin{array}{*{20}{c}} {\arg \mathop {\min }\limits_{\mathit{\boldsymbol{D}}, \mathit{\boldsymbol{A}}} \left\| {\mathit{\boldsymbol{M}} - \mathit{\boldsymbol{DA}}} \right\|_{\rm{F}}^2}&{{\rm{s}}.{\rm{t}}.}&{{{\left\| {{\mathit{\boldsymbol{\alpha }}_i}} \right\|}_0} \le \mu , 1 \le i \le \left( {s + 1} \right)} \end{array}$$ (2) 式中：$\mathit{\boldsymbol{M}} \in {R^{{n^2} \times \left( {s + 1} \right)}}$为结构组矩阵；$\mathit{\boldsymbol{D}} \in {R^{{n^2} \times L}}$为字典；$\mathit{\boldsymbol{A}} \in {R^{L \times \left( {s + 1} \right)}}$为系数矩阵；αi是系数矩阵A的第i列；F表示Frobenius范数。
在字典学习中，字典更新阶段是在系数矩阵A不变的情况下，对字典D进行更新。可以表示为下面的优化问题：
$$\begin{array}{*{20}{c}} {\left\{ {\mathit{\boldsymbol{\hat D}}, \mathit{\boldsymbol{\hat A}}} \right\}{\rm{ = arg}}\mathop {{\rm{min}}}\limits_{\mathit{\boldsymbol{D}}, \mathit{\boldsymbol{A}}} \left\| {\mathit{\boldsymbol{M}}{\rm{ - }}\mathit{\boldsymbol{DA}}} \right\|_{{F}}^{\rm{2}}}&{{\rm{s}}.{\rm{t}}.}&{\mathit{\boldsymbol{A}} \odot \mathit{\boldsymbol{M}} = 0} \end{array}$$ (3) 式中：A☉M是两个大小相同矩阵之间的Schur乘积。掩膜矩阵M中只有0和1两种数值，通过M={|A|=0}得到，也就是，如果A(i, j)=0时，M(i, j)=1，其他的位置为0。因此，A☉M使得A中所有的零元素保持不变。虽然式(3)中的优化问题比求解整个字典简单，但是仍然是一个非凸的问题，求解困难。因此采用块坐标下降法，固定A并采用MOD算法最小化$\left\| {\mathit{\boldsymbol{M}} - \mathit{\boldsymbol{DA}}} \right\|_F^2$，对D进行求解。之后固定D对A进行更新，也就是通过解决下面的优化问题实现：
$$\begin{array}{*{20}{c}} {\mathit{\boldsymbol{\hat A}}{\rm{ = arg}}\mathop {{\rm{min}}}\limits_\mathit{\boldsymbol{A}} \left\| {\mathit{\boldsymbol{M}}{\rm{ - }}\mathit{\boldsymbol{DA}}} \right\|_\mathit{F}^{\rm{2}}}&{{\rm{s}}.{\rm{t}}.}&{\mathit{\boldsymbol{A}} \odot \mathit{\boldsymbol{M}} = 0} \end{array}$$ (4) 上式中的优化问题有一个封闭解，可以通过单独对A中的每一列求解，并且在求解过程中，只更新每一列的非零项，值为零的项保持不变，这样的优化问题可以表示为：
$${\hat \alpha _i}{\rm{ = }}\arg \mathop {\min }\limits_{{\alpha _i}} \left\| {{m_i} - {\mathit{\boldsymbol{D}}_i}{\alpha _i}} \right\|_2^2$$ (5) 式中：Di是D的子矩阵，包含在优化过程中参与表示的字典原子。$ {\hat \alpha _i}$是系数矩阵A第i列中的非零位置。
上述迭代算法先更新字典D，后更新系数矩阵A，可以得到式(3)优化问题的近似解。组K-SVD算法对迭代更新过程进行了改进，结构组矩阵M在字典D上的稀疏编码优化问题可以表示为：
$$\begin{array}{*{20}{c}} {\left\{ {\mathit{\boldsymbol{\hat D}}, \mathit{\boldsymbol{\hat A}}} \right\}{\rm{ = arg}}\mathop {{\rm{min}}}\limits_{\mathit{\boldsymbol{D}}, \mathit{\boldsymbol{A}}} \left\| {\mathit{\boldsymbol{M}}{\rm{ - }}\sum\limits_{{\rm{j = 1}}}^{\rm{n}} {{\mathit{\boldsymbol{d}}_\mathit{j}}} \mathit{\boldsymbol{a}}_\mathit{j}^{\rm{T}}} \right\|_\mathit{F}^{\rm{2}}}\\ {\begin{array}{*{20}{c}} {{\rm{s}}.{\rm{t}}.}&{\nabla 1 \le j \le L, {\mathit{\boldsymbol{m}}_j} \odot {\mathit{\boldsymbol{a}}_j} = 0} \end{array}} \end{array}$$ (6) 式中：$\mathit{\boldsymbol{M}} \in {R^{{n^2} \times \left( {s + 1} \right)}}$为结构组矩阵；dj是字典D的第j列；ajT是A的第j行；mjT是矩阵M的第j行，式中mj⊙aj是两个矩阵之间Schur乘积。采用块坐标下降法，通过对矩阵${\mathit{\boldsymbol{E}}_j} = \left( {\mathit{\boldsymbol{M}} - {\sum _{i \ne j}}{\mathit{\boldsymbol{d}}_i}\mathit{\boldsymbol{a}}_i^{\rm{T}}} \right) \odot \left( {{1_d} \cdot \mathit{\boldsymbol{m}}_j^{\rm{T}}} \right)$进行SVD分解来更新(dj, aj)，j=1, 2, …, L。掩膜矩阵${1_d} \cdot \mathit{\boldsymbol{m}}_j^T$是一个秩为1，大小为n2×(s+1)的矩阵，由d个mjT作为矩阵的行组成。掩膜矩阵有效地将$\mathit{\boldsymbol{X}} - \sum\nolimits_{i \ne j} {{\mathit{\boldsymbol{d}}_i}\mathit{\boldsymbol{a}}_i^{\rm{T}}} $相应的原子移除。
2. 融合方法
可见光与红外图像组K-SVD融合方法如图 3所示，具体步骤如下：
1）图像块提取：对于大小为M×N的输入源图像I1、I2，采用滑动窗技术，滑动步长为1像素，按照从左到右、自上而下的轨迹遍历整幅图像。设图像块大小为n×n，则可以分别提取源图像I1、I2的(M-n+1)·(N-n+1)个图像块。
2）结构组矩阵构建：对源图像I1，在L×L的邻域内，以欧式距离为度量准则，计算出初始图像块的最相似的s个图像块，并进行块向量化、结构组矩阵化处理，得到相应的结构组矩阵；对源图像I2，每个初始块对应的相似块与I1相同。
3）训练样本选取及字典学习：从源图像中随机选出x个初始块作为训练样本，采用组K-SVD字典训练方法求解源图像对应的字典。
4）稀疏系数求解：采用系数重用OMP（Orthogonal Matching Pursuit）算法，在得到的字典D上对源图像对应的结构组矩阵${\mathit{\boldsymbol{M}}_{1i}} \in {R^{{n^2} \times \left( {s + 1} \right)}}$，${\boldsymbol{M}_{2i}} \in {R^{{n^2} \times \left( {s + 1} \right)}}$，i=1, 2, …, (M-n+1)·(N-n+1)进行稀疏分解，得到对应的稀疏系数矩阵${\mathit{\boldsymbol{A}}_{1i}} \in {R^{L \times \left( {s + 1} \right)}}$，${\mathit{\boldsymbol{A}}_{2i}} \in {R^{L \times \left( {s + 1} \right) \cdot }}$。
5）稀疏系数融合：α1i(t)、α2i(t)，t=1, 2, …, s+1为系数向量，分别为稀疏系数矩阵A1、A2的第t列。对应的稀疏系数向量采用绝对值取大的规则进行融合：
$${\alpha _{{\rm{Fi}}}} = \left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{\alpha _{1i}}\left( t \right)}&{\left| {{\alpha _{1i}}\left( t \right)} \right| \geqslant \left| {{\alpha _{2i}}\left( t \right)} \right|} \\ {{\alpha _{2i}}\left( t \right)}&{\left| {{\alpha _{2i}}\left( t \right)} \right| < \left| {{\alpha _{1i}}\left( t \right)} \right|} \end{array}}&{} \end{array}} \right.t = 1, 2, \cdots , s + 1$$ (7) 最终得到融合系数矩阵AFi。
6）图像重构：将得到的系数矩阵AFi与字典D重构得到融合结构组矩阵，对结构组矩阵进行逆变换得到融合的图像块，之后将每一个图像块放回到对应的位置，相同位置像素进行均值化，得到最终的融合图像。
3. 实验结果及分析
3.1 设定
为了评价组K-SVD融合性能，实验采用4组可见光与红外图像，分别为“UNcape”、“Kaptein”、“Dune”和“Road”。对比的融合方法采用3种传统的稀疏融合方法，分别为稀疏表示融合方法（Sparse Representation, SR）[9]、联合稀疏融合方法（Joint Sparse Representation, JSR）[10]和自适应稀疏融合方法（Adaptive Sparse Representation, ASR）[11]。评价指标采用Q0[21]、Qw[22]、Qe[22]和Qab/f[23]四个评价指标。Q0从图像的亮度、对比度和结构失真度来衡量融合图像的退化程度，Qw描述源图像的显著信息转移到融合图像的合并程度，而Qe和Qab/f主要描述源图像到融合图像的边缘信息融合程度。评价指标数值越大，说明融合效果越好。
组K-SVD融合方法参数设定如下：①结构组构造模型参数，图像块大小为8×8，滑动窗大小为40×40，滑动步长为1像素。②结构组稀疏模型参数，训练样本图像块数量为10000，相似图像块个数设定为19，字典大小为64×256，字典原子稀疏度为5。
3.2 结果分析
4组可见光与红外图像的融合结果如图 4~图 7所示，其中，(a)和(b)分别表示可见光与红外图像，(c)、(d)、(e)分别是稀疏表示融合方法、自适应稀疏融合方法和联合稀疏融合方法得到的融合图像，(f)是本文方法得到的融合图像。
图 4为“Nato_camp”的融合实验结果。从图 4(a)和(b)可以看出，红外图像行人目标位置清晰，背景信息（如树木、道路和房屋等）细节不清楚，信息丢失较为严重；而可见光图像中，道路、树木和房屋等背景信息清晰，但行人目标缺失，可以看出这两类图像具有很大的互补性。图 4(c)是SR融合图像，可以看出，融合图像的边缘信息较为模糊，纹理信息缺失严重；图 4(d)是ASR融合图像，融合图像较好保留了两类图像的主要互补信息，但依然存在图像边缘信息模糊，存在一定的伪影现象；图 4(e)是JSR融合图像，行人目标清晰，对比度较高，但融合图像存在过增强现象，背景信息缺失严重。图 4(f)是本文方法的融合图像，从结果可以看出，融合图像行人目标位置、轮廓信息清楚，背景信息如树木、房屋、道路等保留较好，图像整体视觉效果好于传统稀疏融合方法。
图 5为“Kaptein”的融合实验结果。图 5(a)和(b)分别是可见光和红外图像，图 5(c)、(d)和(e)分别是SR融合图像、ASR融合图像和JSR融合图像。从结果可以看出，SR融合图像虽然保留了红外图像的目标信息，但图像背景细节信息丢失严重，例如地面、树木等；ASR融合图像视觉效果好于SR融合图像，行人目标、背景信息较为清楚，但边缘信息较为模糊；JSR融合图像行人目标的对比度较高，但融合图像过增强现象严重，如树木、门窗的边缘细节信息丢失。图 5(f)是本文方法的融合结果，可以看出，融合图像包含了两类图像的主要互补信息，行人、地面、树木等细节信息清楚，视觉效果较好。
图 6和图 7分别为“Duine”和“Road”的融合实验结果。从实验结果来看，SR和ASR的融合效果相对较差，JSR融合图像虽然对比度较高，但细节信息丢失严重，本文方法的融合图像视觉效果较好，图像中的细节信息及边缘轮廓信息保留完整。从主观视觉来看，4组实验结果评价结果一致。
4组可见光与红外图像融合结果的客观评价指标如表 1~表 4所示，表 5为客观评价指标的平均值。从结果来看，对于“Kaptein”图像，本文融合方法获得最好的融合效果；对于“Nato_camp”图像，本文的融合方法在评价指标Qw、Qe和Qab/f获得最优的结果，Q0指标略低于ASR的结果；对于“Duine”和“Road”图像，本文的融合方法在评价指标Qw和Qe获得最优的结果，而指标Q0和Qab/f低于ASR方法，这说明本文的融合方法在提取显著特征信息和边缘信息保留方面效果较好，总体来看，本文的融合方法优于传统的稀疏表示融合方法。
表 1 “Nato_camp”图像融合客观评价指标Table 1. Objective evaluation index of Nato_camp imagesMethod Q0 Qw Qe Qab/f SR 0.5857 0.7322 0.7151 0.3216 ASR 0.5966 0.7617 0.7439 0.4673 JSR 0.5649 0.7311 0.7140 0.3760 The proposed 0.5926 0.7737 0.7556 0.4871 表 2 “Kaptein”图像融合客观评价指标Table 2. Objective evaluation index of Kaptein imagesMethod Q0 Qw Qe Qab/f SR 0.5759 0.7365 0.7193 0.2489 ASR 0.5760 0.7659 0.7480 0.4206 JSR 0.5445 0.7231 0.7062 0.3400 The proposed 0.5794 0.7941 0.7755 0.4807 表 3 “Duine”图像融合客观评价指标Table 3. Objective evaluation index of Duine imagesMethod Q0 Qw Qe Qab/f SR 0.6426 0.8760 0.8555 0.2463 ASR 0.6781 0.9247 0.9031 0.5604 JSR 0.3215 0.7440 0.7266 0.2091 The proposed 0.6750 0.9312 0.9094 0.4783 表 4 “Road”图像融合客观评价指标Table 4. Objective evaluation index of road imagesMethod Q0 Qw Qe Qab/f SR 0.6909 0.7994 0.7807 0.5101 ASR 0.6941 0.8053 0.7865 0.6073 JSR 0.6346 0.7763 0.7581 0.5210 The proposed 0.6924 0.8099 0.7909 0.5815 表 5 不同融合方法评价指标的平均值Table 5. The average evaluation for different fused methodsMethod Q0 Qw Qe Qab/f SR 0.6238 0.7860 0.7677 0.3276 ASR 0.6358 0.8156 0.7965 0.5066 JSR 0.5264 0.7436 0.7262 0.3544 The proposed 0.6349 0.8272 0.8079 0.5053 为了进一步对融合效果进行评价，采用两组名为“Nato_camp”和“Duine”的图像序列进行实验。“Nato_camp”图像序列为同一场景中32幅在时间上连续的图片，“Duine”图像序列为同一场景中23幅在时间上连续的图片。采用Q0、Qw、Qe和Qab/f对融合效果进行评价。将得到的评价指标绘制为曲线图进行对比，具体的绘制方法为：分别将4种融合方法的Q0值绘制为曲线，并将这4条曲线整合到同一幅图像中；另外3种评价指标曲线图的绘制方法与Q0相同。
图 8和图 9分别为图像序列“Nato_camp”和“Duine”对应的曲线图。如图 8所示，(a)、(b)、(c)、(d)分别为评价指标Q0、Qw、Qe和Qab/f对应的曲线图。每幅图中，点状曲线对应SR方法的评价指标，圈型曲线对应ASR方法，块状曲线对应JSR方法，叉型曲线对应本文提出的方法，右下角为图像序列不同融合方法对应评价指标的平均值。通过两组分析曲线可知：在图 8中本文方法的Q0值低于ASR方法，Qw、Qe和Qab/f值均为最高；在图 9中本文方法的Q0、Qw和Qe值略低于ASR方法，Qab/f值为最高。
因此，通过对4组可见光与红外图像及两组图像序列进行融合实验，结合主观视觉效果和客观评价指标对实验结果判断，本文方法与另外3种方法相比，融合效果最好。
4. 结语
本文提出了可见光与红外图像组K-SVD融合方法，利用图像非局部相似性，以欧式距离为度量准则，构建图像结构组矩阵，以结构组代替传统的图像块，在此基础上，研究了组K-SVD的稀疏分解模型，以结构组矩阵作为训练样本，构建学习字典，并进行组稀疏分解，采用绝对值取大的合并规则进行合并，得到最终的融合图像。实验结果表明，该方法可以有效提取图像的显著特征信息，减少边缘特征信息丢失，融合图像视觉效果较好，在主观和客观评价上，都优于传统的稀疏表示融合方法。
4. 结语
本文提出了可见光与红外图像组K-SVD融合方法，利用图像非局部相似性，以欧式距离为度量准则，构建图像结构组矩阵，以结构组代替传统的图像块，在此基础上，研究了组K-SVD的稀疏分解模型，以结构组矩阵作为训练样本，构建学习字典，并进行组稀疏分解，采用绝对值取大的合并规则进行合并，得到最终的融合图像。实验结果表明，该方法可以有效提取图像的显著特征信息，减少边缘特征信息丢失，融合图像视觉效果较好，在主观和客观评价上，都优于传统的稀疏表示融合方法。
-
![]()
图 5 “阿利·伯克”级驱逐舰温度场计算结果
Figure 5. Calculation results of temperature field for the "Ali Burke" class destroyer
![]()
图 6 “阿利·伯克”级驱逐舰壁面温度计算结果
Figure 6. Calculation results of wall temperature for the Ali Burke class destroyer
![]()
图 7 “阿利·伯克”级驱逐舰速度场计算结果
Figure 7. Calculation results of velocity field for the "Ali Burke" class destroyer
![]()
图 8 “阿利·伯克”级驱逐舰CO2质量分数计算结果
Figure 8. Calculation results of CO2 mass fraction for the Ali Burke class destroyer
![]()
图 9 “阿利·伯克”级驱逐舰红外探测点分布示意图
Figure 9. Schematic diagram of the distribution of infrared detection points on the Ali Burke class destroyer
![]()
图 10 中波红外典型探测点示意图
Figure 10. Schematic diagram of typical detection points for medium-wave infrared
![]()
图 11 中波光谱辐射强度示意图
Figure 11. Schematic diagram of the radiation intensity of the Medium-wave spectrum
![]()
图 12 长波红外典型探测点示意图
Figure 12. Schematic diagram of typical long-wave infrared detection points
![]()
图 13 长波光谱辐射强度示意图
Figure 13. Schematic diagram of long-wave spectral radiation intensity
![]()
图 14 外场采集渔船红外成像图
Figure 14. Infrared imaging of fishing boats collected from external fields
![]()
图 15 渔船仿真结果与试验结果对比
Figure 15. Comparison between simulation and experimental results of fishing vessels
表 1 “阿利·伯克”级驱逐舰主要参数
Table 1 Main parameters of the "Ali Burke" class destroyer
Model "Ali Burke" class destroyer Length/m 155.29 Width/m 20.4 Draft/m 6 Tonnage/t 9238 Engine 4*LM-2500 gas turbine engine Highest Speed 31knots Endurance 4200 nautical miles（20 knots） 
下载: 导出CSV
表 2 红外仿真计算参数设置
Table 2 Infrared simulation calculation parameters setting
Infrared calculation setting parameters Value Detection range/km 5 Number of detection points 60 Gas components H2O, CO2 Gas spectral band model Narrow spectral band Solid emissivity 0.8 Calculation method for atmospheric
transmittanceModtran Spectral calculation range/μm Medium-wave: 3-5
Long-wave: 8-12
下载: 导出CSV
-
[1] 杨栋, 高德亮, 曹耀心, 等. 红外导引头抗诱饵干扰研究[J]. 飞控与探测, 2020, 3(3): 80-85. https://www.cnki.com.cn/Article/CJFDTOTAL-FKTC202003011.htm YANG Dong, GAO Deliang, CAO Yuexin, et al. Research on anti-decoy interference of infrared seeker[J]. Flight Control and Detection, 2020, 3(3): 80-85. https://www.cnki.com.cn/Article/CJFDTOTAL-FKTC202003011.htm
[2] 陈彬. 海面背景下舰船目标红外辐射特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2015. CHEN Bin. Research on Infrared Radiation Characteristics of Ship Targets against Sea Surface Background[D]. Harbin: Harbin Institute of Technology, 2015.
[3] 毛佳晨. 海洋背景下舰船目标的红外成像特征研究[D]. 西安: 西安电子科技大学, 2019. MAO Jiachen. Research on Infrared Imaging Characteristics of Ship Targets in Marine Background[D]. Xi'an: Xidian University, 2019.
[4] 张胜辉. 舰船目标红外中/长波特性分析与检测方法研究[D]. 长沙: 国防科学技术大学, 2015. ZHANG Shenghui. Analysis and Detection Methods of Infrared Medium/Long-wave Characteristics of Ship Targets[D]. Changsha: National University of Defense Technology, 2015.
[5] 庄亚平, 王建勋, 左超, 等. 舰船热尾流海面和大气传输红外辐射特性分析[J]. 船电技术, 2019, 39(7): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-CDJI201907019.htm ZHUANG Yaping, WANG Jianxun, ZUO Chao, et al. Analysis of infrared radiation characteristics of ship thermal wake transmitted by sea surface and atmosphere[J]. Ship Electric Technology, 2019, 39(7): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-CDJI201907019.htm
[6] 仇荣超, 娄树理, 李延军, 等. 多波段红外图像的海面舰船目标检测[J]. 光谱学与光谱分析, 2019, 39(3): 698-704. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201903007.htm QIU Rongchao, LOU Shuli, LI Yanjun, et al. Detection of ship targets on sea surface in multi-band infrared images[J]. Spectroscopy and Spectral Analysis, 2019, 39(3): 698-704. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201903007.htm
[7] 罗婷婷. 美国海军"阿利•伯克"Flight Ⅲ型驱逐舰解析[J]. 现代舰船, 2011, 4: 32-35. https://www.cnki.com.cn/Article/CJFDTOTAL-XDJC201104015.htm LUO Tingting. Analysis of the USN arleigh burke class flight Ⅲ destroyer[J]. Modern Warships, 2011, 4: 32-35. https://www.cnki.com.cn/Article/CJFDTOTAL-XDJC201104015.htm
[8] 周强. 美国海军"阿利•伯克"级"宙斯盾"驱逐舰[J]. 现代兵器, 2003, 8: 55. https://www.cnki.com.cn/Article/CJFDTOTAL-XDBQ200308016.htm ZHOU Qiang. The USN arleigh burke class aegis destroyer[J]. Modern Weapons, 2003, 8: 55. https://www.cnki.com.cn/Article/CJFDTOTAL-XDBQ200308016.htm
[9] 吴沿庆, 廖守亿, 张作宇, 等. 飞机发动机尾流流场数值模拟与红外特性计算[J]. 激光与红外, 2017, 47(11): 1380-1385. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201711011.htm WU Yanqing, LIAO Shouyi, ZHANG Zuoyu, et al. Numerical simulation of flow field and infrared characteristic calculation of aircraft engine wake vortex[J]. Laser & Infrared, 2017, 47(11): 1380-1385. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201711011.htm
[10] 阮鸿雁, 吕建军, 杨波. 激光切割超音速喷嘴的流场特性[J]. 中国激光, 2009, 36(5): 1233-1238. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ200905043.htm RUAN Hongyan, LYU Jianjun, YANG Bo. Flow field characteristics of supersonic nozzles for laser cutting[J]. Chinese Journal of Lasers, 2009, 36(5): 1233-1238. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ200905043.htm
[11] 刘滨, 洪文鹏, 王瑞芬, 等. 红外成像系统中带有距离修正的目标能量预测方法[J]. 飞控与探测, 2020, 3(2): 59-63. https://www.cnki.com.cn/Article/CJFDTOTAL-FKTC202002008.htm LIU Bin, HONG Wenpeng, WANG Ruifeng, et al. Target energy prediction method with distance correction in infrared imaging system[J]. Flight Control and Detection, 2020, 3(2): 59-63. https://www.cnki.com.cn/Article/CJFDTOTAL-FKTC202002008.htm
[12] 丁宣浩, 蔡如华. 光线方程的数值解与光线追迹[J]. 西南师范大学学报(自然科学版), 2005, 30(4): 667-671. https://www.cnki.com.cn/Article/CJFDTOTAL-XNZK200504017.htm DING Xuanhao, CAI Ruhua. Numerical solution of ray equation and ray tracing[J]. Journal of Southwest China Normal University (Natural Science Edition), 2005, 30(4): 667-671. https://www.cnki.com.cn/Article/CJFDTOTAL-XNZK200504017.htm
-
期刊类型引用(3)
1. 山清. 基于变跟踪点的红外场景仿真原理及关键技术研究. 激光与红外. 2020(07): 825-832 . 
百度学术
2. 姬爽，王艳奎，吴根水. 变跟踪点红外场景仿真建模时间误差分析. 激光与红外. 2019(01): 59-66 . 百度学术
3. 王艳奎，吴根水，姬爽. 导引头变跟踪点红外图像仿真建模研究. 红外技术. 2018(12): 1170-1175 . 本站查看
其他类型引用(3)
计量
- 文章访问数: 112
- HTML全文浏览量: 44
- PDF下载量: 56
- 被引次数: 6
