SAR ATR method based on canonical correlations analysis of features extracted by 2D random projection

合成孔径雷达（Synthetic aperture radar，SAR）具有很强的数据测量和成像能力，因而在不同领域发挥着重要作用。以SAR自动目标识别（Automatic target recognition，ATR）为代表的关键技术成为情报侦察、情报解译以及导弹制导等军事应用的重要支撑^[1]。伴随高分辨SAR成像技术的发展和成熟，基于SAR图像的战术目标识别方法应运而生。这些方法主要采用“特征提取+分类器”的两阶段思路。特征提取从原始SAR图像出发，秉承去粗存精的思维抽取有价值的目标特征，如几何形状、灰度分布以及散射特性等。参考文献[2-3]通过目标区域等几何特性进行分析和匹配识别，体现了目标物理外观信息。然而，受SAR图像中的噪声、杂波等影响以及可能出现的人为干扰，几何形状等信息的抽取和分析难度较大。参考文献[4-10]以主成分分析（Principal component analysis，PCA）、单演信号、模态分解等数学投影或信号分解算法为基本手段获取SAR图像矢量化特征。这些特征具有良好的一致性，并且提取效率很高。不足的是，它们往往不能反映目标的物理层信息。描述SAR目标散射特性的特征包括极化信息^[11]、局部散射中心^[12-13]等。其中，散射中心特征应用更为普遍，可通过属性参数（位置、长度、结构等）的匹配进行类别决策。然而，由于散射中心模型一般较为复杂，散射中心参数估计难度较大，在一定程度影响了识别精度。SAR目标识别中使用的分类器大多继承于传统模式识别领域，或是根据SAR图像的特点进行了优化改进，如参考文献[14]采用的K近邻，参考文献[15-16]使用的支持向量机（Support vector machine，SVM），参考文献[17-19]采用的稀疏表示分类（Sparse representation-based classification，SRC）以及近年来常见的卷积神经网络（Convolutional neural network，CNN）^[20-22]等。

根据上述介绍，设计或选取有效的特征是影响SAR目标识别性能的关键因素之一。文中提出基于二维随机投影特征多重集典型相关分析（Multiset canonical correlations analy-sis，MCCA）^[23-24]融合的识别方法。参考文献[25]根据压缩感知的原理提出了基于随机投影的特征提取算法，并在SAR目标识别中进行了应用。其结果表明随机投影作为一种非参数化的降维算法，具有很高的效率并且可以保持特征的鉴别力。考虑到一维随机投影对于SAR图像二维结构信息的潜在破坏性，参考文献[26-27]对其进行了扩展，提出二维随机投影算法。经验证，二维随机投影针对SAR目标识别问题具有性能优势。为此，文中以二维随机投影作为特征提取的基础算法，通过若干个不同的投影矩阵获取SAR图像的多层次特征描述。不同的投影结果之间具有良好的互补性，同时也存在一定的冗余或干扰。为此，文中进一步采用MCCA对不同投影矩阵下获得的特征进行融合分析。融合的特征矢量既可以保持不同层次特征之间的关联性，又可以有效去除冗余和干扰成分，从而为后续分类提供高鉴别力的特征输入。在分类阶段，选用稳健性较强的SRC作为基础分类器，通过融合特征矢量的训练和处理获取最终分类结果。

根据压缩感知原理，稀疏信号可基于少量的观测进行重构。因此，存在冗余的信号可按此原理进行有效降维。参考文献[25]根据这一理论设计了随机投影算法并在特征提取和目标识别中进行了应用。为进一步处理图像等二维数据，参考文献[26-27]将一维的投影矩阵延伸到二维。以图像分析为例，输入图像$X \in {{{R}}^{{n_1} \times {n_2}}}$经过两个随机投影矩阵${A} \in {{{R}}^{{m_1} \times {n_1}}}$和${B} \in {{{R}}^{{m_2} \times {n_2}}}$，按照$ {Y} = {AX}{{B}^{\text{T}}} $的方式变换得到特征矩阵${Y} \in {{{R}}^{{m_1} \times {m_2}}}$，其中$ {m_1} $远小于$ {n_1} $，$ {m_2} $远小于$ {n_2} $，实现高效降维。

根据参考文献[26]相关论述，为了使得二维投影后的特征矩阵能够有效保持图像特性，要求通过低维度的投影矩阵$ Y $能够对原始图像进行高精度重构。此时，需要对两个投影矩阵进行如下约束：

式中：$ {\left\| {X} \right\|_0} $计算$ {X} $的$ {\ell _0} $范数，即为矩阵中非零元素的个数；函数$ {\text{spark}} $获取当前矩阵中相互独立的列数。对于矩阵${A} \in {{{R}}^{{m_1} \times {n_1}}}$和${B} \in {{{R}}^{{m_2} \times {n_2}}}$，满足$ {\text{spark}}({A}) = {m_1} + 1 $，$ {\text{spark}}({B}) = {m_2} + 1 $。

基于上述准则构造的二维随机投影矩阵可有效获取SAR图像的鉴别力特征，并已成功应用于目标识别。然而，这种特征提取方法由于投影矩阵的不唯一性而具有一定的随机性，导致识别结果存在不稳定的问题。为此，文中采用若干个二维随机投影矩阵对原始图像进行处理。不同投影矩阵下的结果具有互补性，可有效弥补单一投影矩阵存在的不稳定。因此，联合不同投影矩阵下的结果有利于提供最终的识别性能。

传统典型相关分析（Canonical correlations analysis，CCA）主要用于比较两组随机变量，分析它们之间的关联性和差异。MCCA则是对传统CCA进行延伸扩展，主要用于分析多组不同随机变量之间的关联性和差异^[23-24]。记一组随机变量为$ {X_1},{X_2}, \cdots ,{X_n} $，它们的维度分别为$ {m_1},{m_2}, \cdots ,{m_n} $（$ {m_1} $为最小值）。MCCA算法首先对各个随机变量进行去均值和中心化。然后，设计如公式（2）所示的准则函数：

式中：$ {S_{ij}} = E({X_i}X_j^{\text{T}}) $计算$ {X_i} $与$ {X_j} $两个随机变量的互协方差矩阵。为了对这一准则函数进行最大化，开展以下优化过程：

对于公式（3）的优化过程，可利用Lagrange乘子法进行求解：

根据对公式（4）的求解结果进一步获取变换矩阵$ A = {[{\alpha _1},{\alpha _2}, \cdots ,{\alpha _n}]^{\text{T}}} $，其中：

对应地，从公式（5）中可以得到不同变量$ {X_i} $对应的投影方向，采用公式（6）对不同随机变量加权融合，得到结果如下：

文中采用MCCA对同一幅SAR图像的多个不同特征进行融合。最终获取的特征矢量不仅能保持不同投影矩阵的独立性，也可以去除可能存在的冗余或干扰成分。因此，融合后的特征矢量具有更强的鉴别力和简洁性，有利于提高识别算法的整体性能。

自稀疏表示理论提出和应用以来，其衍生工具在信号、图像处理等领域得到十分广泛的运用。SRC就是一个典型代表，主要用于模式识别中监督分类问题^[17-19]。对于多类别的训练样本，可根据它们建立全局字典$A = [{A_1},{A_2}, \cdots ,{A_C}] \in {{{R}}^{d \times N}}$，其中$ {A_i} $对应来自第$ i $类中$ {N_i} $个训练样本。测试样本$ y $在全局字典上进行线性表示，求解线性表示系数：

式中：$ x $为需要求解的系数矢量；$ \varepsilon $为误差门限。

根据现有文献，上述非凸优化问题可通过$ {\ell _{\text{1}}} $范数替代$ {\ell _0} $范数进行转换，从而简化优化问题，获取近似解。此外，也可采用正交匹配追踪等贪婪算法进行寻优。根据系数求解结果，可以按照公式（8）分别阶段单个类别的重构误差大小：

式中：$ {x_i} $为第$ i $类上的系数矢量；$ r(i){\text{ }} $即为$ i $类别的最终误差。最后，可以根据误差判定测试样本与哪一类训练样本最为相似，从而获得类别决策结果。

根据前文思想，文中设计如图1所示的目标识别流程，实线部分代表训练流程，虚线部分代表测试流程，中间部分为相同流程，具体实施描述如下。首先，按照二维随机投影的构造原则获取若干对随机投影矩阵，用于对SAR图像的特征提取。对于训练样本，分别采用不同的投影矩阵进行特征提取，进一步采用MCCA对它们进行融合，获得单一的特征矢量，存入全局字典。据此可构建与训练样本规模一致的全局字典。对于测试样本，按照上述步骤经随机投影特征提取和MCCA特征融合获得对应的特征矢量。最终，采用SRC对测试样本的特征矢量进行分类，判决目标所属类别。

Figure 1.  Procedure of ATR based on MCCA of 2D random projection features

二维随机投影矩阵在无需先验知识的情况下，保持了SAR图像的二维结构信息。同时，结合多个不同二维随机投影矩阵的结果，能够实现对目标特性互补的描述，进一步增强特征层的鉴别力。基于MCCA的特征融合有效继承了多层次二维随机投影特征的优势，同时去除了存在的冗余分量。SRC中的稀疏重构分类机制具有较强的噪声、遮挡稳健性，采用其进行最终决策有利于提高识别结果的稳健性。因此，文中方法通过结合特征提取和分类决策的优势提升整体识别性能。

在SAR目标识别领域，MSTAR数据集具有很强的代表性，是当前绝大多数相关方法开展验证的基础数据集。MSTAR计划中，研究人员在多种场景和操作条件下采集了图2所示的10类目标的SAR图像，分辨率达到0.3 m。据此，可分别构建训练和测试样本，验证目标识别算法的有效性。根据现有文献报道，在MSTAR数据集的样本基础上可设置诸如标准操作条件、目标型号差异、俯仰角差异等测试条件。此外，通过对原始MSTAR样本进行简单处理，可进一步构建噪声干扰、部分遮挡等测试条件。文中在上述条件中重点对标准操作条件、俯仰角差异以及噪声干扰三种场景进行针对性实验，检验提出方法的有效性和稳健性。

Figure 2.  Images of the vehicle targets in the MSTAR dataset

实验中还设置了一些对比方法，它们包括参考文献[25]基于随机投影的方法（Random projection），该方法采用一维随机投影算法进行SAR图像降维；参考文献[26]中二维随机投影的方法（记为2D random projection），该方法将一维投影扩展到二维，从而更好地保持图像结构信息；参考文献[8]采用单演信号特征的方法（记为Monogenic），结合不同成分特征进行决策；参考文献[20]中基于CNN的方法（记为CNN），该方法针对SAR目标识别设计了特定结构的CNN。

标准操作条件是SAR目标识别问题中一项基础且具有代表性意义的实验场景。根据MSTAR数据集包含的目标类别和不同条件下的SAR图像数据，设置表1所示的测试条件，且每个目标只有一种型号。测试集为各类目标15°俯仰角下的样本，其中BMP2和T72两类目标各包含3个型号。对比训练和测试集，两种样本虽然存在一定差异但整体相似度较高，可认为是标准操作条件。对于10类目标的识别结果，采用平均识别率进行定量评价，即正确识别样本数所占的比例。据此，表2对比了不同方法的识别性能。文中方法以99.08%的平均识别率优于四类对比方法。CNN方法性能排名第二，显示了深度学习网络优越性能。二维随机投影方法的性能优于传统的一维随机投影方法，体现了二维特征提取的必要性。相比单演信号方法，二维随机投影方法略有劣势，主要是特征鉴别力的差距。然而，单演信号特征分解的复杂度要远远高于二维随机投影。文中方法在多层次二维随机投影特征的基础上进行融合处理，经过MCCA过后的特征矢量显著提升了识别性能，反映了方法设计思路的科学性。

在上述标准操作条件中，一个重要的因素是测试样本与训练样本的俯仰角差异很小。实际上，测试样本很可能来自于测试样本差距较大俯仰角，此时两者的图像差异也会较为明显。根据MSTAR数据集中的SAR图像数据，设置如表3所示的训练和测试集。训练集包含了2S1、BRDM2以及ZSU23/4共3类目标在17°俯仰角下的SAR图像。测试样本来自差异较大的俯仰角，可进一步区分为两个子集，分别对应30°俯仰角和45°俯仰角。两个测试子集反映了不同程度的俯仰角差异，能够更为全面地反映方法在不同俯仰角差异程度条件下的识别性能。基于上述设置，获取不同方法的识别性能，统计如表4所示。从两个测试子集的角度来看，30°俯仰角下的整体性能要明显优于45°，反映了大俯仰角差异带来的识别难度急剧增加。根据上述结果，文中方法对于不同程度俯仰角差异的适应性和稳健性更强，能够在俯仰角差异的情形下保持良好性能。相比标准操作条件，CNN方法在此条件下的性能下降最为剧烈，表明在训练样本不足或与测试样本差异较大时，深度学习算法也存在很大的局限性。文中方法通过MCCA对多层次的二维随机投影特征进行融合，可获取稳健的识别性能，确保了在此场景下的性能。

图像噪声水平直接影响图像质量，低信噪比（Signal-to-noise ratio，SNR）的图像会给识别问题带来很大难度。原始MSTAR数据集的SAR图像整体信噪比较高，据此建立的训练和测试集具有相接近的噪声惠普，因此难以评价方法在测试样本不同信噪比条件下的识别性能。为此，文中按照参考文献[12]的思路和方法对表1中的测试样本进行噪声条件，以原始图像为参照分别构建多个不同信噪比（含−10 dB、−5 dB、0 dB、5 dB和10 dB）的样本测试集。然后，就不同信噪比下的测试样本对各类方法进行测试，获得如图3所示的平均识别率曲线。伴随信噪比的降低，各类方法的平均识别率均出现了较为明显的下降，充分反映了噪声干扰的影响。对比几类方法，文中方法可在不同噪声水平保持最佳性能，具有最强的噪声稳健性。与俯仰角差异的情况类似，CNN方法由于样本不足，导致性能下降十分剧烈。与二维随机投影方法相比，文中通过多个随机投影矩阵的联合使用和有效的特征融合显著提升了对于噪声干扰的稳健性。

Figure 3.  Comparison results under noise corruption

论文针对SAR图像目标识别问题提出结合二维随机投影和MCCA的特征提取方法，进而采用SRC进行分类决策。利用多个不同的二维随机投影矩阵获取原始SAR图像多层次的特征描述。进一步，采用MCCA对多层次的二维随机投影特征进行融合处理，获取一致性更强、维度更低、鉴别力更优的特征矢量。最后，通过SRC对融合的特征矢量进行分类，获取识别结果。在MSTAR数据集上开展了3类典型实验，覆盖了标准操作条件和代表性扩展操作条件。经实验对比分析，提出方法对于SAR ATR可以取得良好性能。