基于低秩矩阵恢复和稀疏表示的sar目标识别方法

基于低秩矩阵恢复和稀疏表示的sar目标识别方法
【专利摘要】本发明提供了一种基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，该方法特别适用于针对低图像质量的SAR图像进行目标识别，其先通过对SAR图像进行低秩矩阵恢复，得到对应的SAR低秩图像，使得绝大部分的SAR目标信息保留在SAR低秩图像中而去除噪声，并结合SAR图像中目标本身的散射特征信息仍处于低维空间、能够通过SAR低秩图像得以保留的特点，结合稀疏表示方法进行图像的目标识别，通过理论分析和实验数据均证明，能够有效的帮助提升识别正确率，解决现有技术中针对低图像质量SAR图像进行目标识别的准确率较低的问题。
【专利说明】
基于低秩矩阵恢复和稀疏表示的SAR目标识别方法
技术领域
[0001 ]本发明涉及雷达目标识别技术领域，尤其涉及一种基于低秩矩阵恢复和稀疏表示 的SAR目标识别方法。
【背景技术】
[0002] 合成孔径雷达（Synthetic Aperture Radar，简称SAR)技术，是采用搭载在卫星或 飞机上的移动雷达，获得高精度地理区域雷达目标图像的一种脉冲雷达技术。合成孔径雷 达是一种主动微波成像系统，通过对目标区域进行电磁波照射，并对回波信号进行信号解 析，它可以提供目标区域的高分辨率图像。它具有全天候、全天时的工作能力和一定的穿透 能力。鉴于它的这些优点，它被广泛应用于矿产探测、海洋环境监测和军事防御等领域。在 军事防御领域中对目标的识别的研究又是最广泛的，因此SAR目标自动识别（Automatic Target Recognition，ATR)的研究受到了国内外学者的广泛关注。
[0003] 然而，由于SAR的主动成像特点和成像过程中的复杂散射机制，SAR图像中的目标 特性和光学图像差别很大，这给目标特征提取和识别带来了很多困难。例如，在实际应用 中，由于SAR传感器探测的大都是机动的非合作目标，而目标往往采用了伪装、隐身等措施 来降低可探测性，因此获取的SAR目标图像往往存在着低信噪比情况，这种图像通常称为低 图像质量SAR图像。由于低信噪比的原因，若直接拿这些低图像质量的SAR图像来进行目标 识别，识别的正确率会大大降低。

【发明内容】

[0004] 针对现有技术中存在的上述问题，本发明提供了一种基于低秩矩阵恢复和稀疏表 示的SAR目标识别方法，特别适用于针对低图像质量的SAR图像进行目标识别，帮助提升其 识别正确率，用以解决现有技术中针对低图像质量SAR图像进行目标识别的准确率较低的 问题。
[0005] 为实现上述目的，本发明采用了如下技术手段：
[0006] 基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，包括如下步骤：
[0007] 1)对于多类不同的已知雷达目标，分别针对每一类采集多个已知雷达目标的SAR 图像作为训练样本，并分别获取各个类别中各个训练样本的SAR图像像素值向量;每个训练 样本的SAR图像像素值向量，是由作为训练样本的SAR图像的各个像素值排列构成；
[0008] 2)分别将每一类别中各个训练样本的SAR图像像素值向量组合构成相应类别的训 练样本图像数据矩阵，并采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩 阵所对应的训练样本图像低秩矩阵，进而根据每个训练样本图像低秩矩阵中对应于各个训 练样本的低秩图像像素值向量，分别恢复得到各个类别中各个训练样本对应的SAR低秩图 像；
[0009] 3)分别提取各个类别中各个训练样本对应的SAR低秩图像的稀疏特征矩阵，将针 对每个训练样本所提取的稀疏特征矩阵作为一个稀疏特征训练样本，从而由提取的各个类 别各个训练样本的稀疏特征矩阵的集合构成稀疏特征训练样本集；
[0010] 4)针对各个待测雷达目标，采集各个待测雷达目标的SAR图像作为测试样本，并分 别获取各个测试样本的SAR图像像素值向量;每个测试样本的SAR图像像素值向量，是由作 为测试样本的SAR图像的各个像素值排列构成；
[0011 ] 5)分别将各个测试样本的SAR图像像素值向量组合构成测试样本图像数据矩阵， 并采用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本图像低秩矩阵，进 而根据测试样本图像低秩矩阵中对应于各个训测试样本的低秩图像像素值向量，分别恢复 得到各个测试样本对应的SAR低秩图像；
[0012] 6)分别提取各个测试样本对应的SAR低秩图像的稀疏特征矩阵；
[0013] 7)针对每个测试样本的稀疏特征矩阵，利用稀疏特征训练样本集中的各个稀疏特 征训练样本建立稀疏线性方程，对相应测试样本的稀疏特征矩阵进行稀疏线性表示，并求 解得到该稀疏线性方程的系数向量，作为相应测试样本的稀疏系数向量;从而得到各个测 试样本的稀疏系数向量；
[0014] 8)针对每个测试样本的稀疏系数向量，分别提取其稀疏系数向量中对应于稀疏特 征训练样本集中每一类已知雷达目标的类别稀疏系数向量，然后分别计算利用每一类已知 雷达目标对应的类别稀疏系数向量通过稀疏线性方程对相应测试样本的稀疏特征矩阵进 行稀疏重构的重构误差，将重构误差最小的类别稀疏系数向量所对应的一个已知雷达目标 类别判定为相应测试样本对应的待测雷达目标所属的雷达目标类别；由此，实现对各个待 测雷达目标的类别识别。
[0015] 上述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法中，具体而言，所述步骤2) 中，采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩阵所对应的训练样本 图像低秩矩阵的具体方式为：
[0016] 针对于第i个类别的训练样本图像数据矩阵表示所采集的已知雷 达目标的类别总数，建立如下的低秩矩阵恢复目标函数：
[0017]
[0018]其中，Ax,4ΡΕΧ>1分别表示对第i个类别的训练样本图像数据矩阵D x,i进行低秩矩阵 恢复所得到的训练样本图像低秩矩阵和训练样本图像误差矩阵为针对第i个类别的训 练样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数，且(XA 1SUl I · I I*为核范数 运算符，M · I |ι为LI范数运算符；
[0019] 通过对所述低秩矩阵恢复目标函数求解，即可得到第i个类别的训练样本图像数 据矩阵Dx,i所对应的训练样本图像低秩矩阵A x,1;由此，分别求解得到各个类别的训练样本 图像数据矩阵所对应的训练样本图像低秩矩阵。
[0020] 上述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法中，作为优选方案，针对第 i个类别的训练样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数取值为
实中，M表示训练样本的SAR图像像素值向量的维数，Ni表示针对第i个 类别已知雷达目标所采集的训练样本的个数。
[0021] 上述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法中，具体而言，所述步骤5) 中，采用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本图像低秩矩阵具 体方式为：
[0022] 针对于测试样本图像数据矩阵Dy，建立如下的低秩矩阵恢复目标函数：
[0023]
[0024]其中，Ay和Ey分别表示对测试样本图像数据矩阵Dy进行低秩矩阵恢复所得到的训 练样本图像低秩矩阵和训练样本图像误差矩阵;λγ为针对测试样本图像数据矩阵的低秩矩 阵恢复所设定的正则化参数，且〇<x y<i; I I · I I*为核范数运算符，I I · I |ι为Li范数运算 符；
[0025] 通过对所述低秩矩阵恢复目标函数求解，即可得到测试样本图像数据矩阵Dy所对 应的测试样本图像低秩矩阵Ay。
[0026] 上述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法中，作为优选方案，针 对测试样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数~7的取值为 其中，M表示测试样本的SAR图像像素值向量的维数，Ny表示测试样本 的个数。
[0027] 上述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法中，具体而言，所述步骤8) 具体为：
[0028] 81)对于第j个测试样本的稀疏系数向量a」，分别提取其中对应于稀疏特征训练样 本集中每一类已知雷达目标的类别稀疏系数向量;其中，稀疏系数向量A对应于第i个类别 已知雷达目标的类别稀疏系数向量S 1(A)为：
[0029]
[0030]其中，类别稀疏系数向量S1(A)的维数与第j个测试样本的稀疏系数向量A的维数 相同，表示第j个测试样本的稀疏系数向量A中对应于第i个类别第m个训练样本的稀 疏特征矩阵的稀疏系数，且类别稀疏系数向量&(A)中对应于非第i个类别已知雷达目标的 稀疏系数的值均为零；1彡i<K，K表示所采集的已知雷达目标的类别总数；1彡m彡N i ,Ni表 示第i个类别包含的训练样本的个数；
[0031] 82)分别计算利用每一类已知雷达目标对应的类别稀疏系数向量通过稀疏线性方 程对相应测试样本的稀疏特征矩阵进行稀疏重构的重构误差，将重构误差最小的类别稀疏 系数向量所对应的一个已知雷达目标类别判定为相应测试样本对应的待测雷达目标所属 的雷达目标类别;其中，确定第j个测试样本对应的待测雷达目标所属的雷达目标类别的具 体计算方式为：
[0032]
[0033] 其中，Ij表示第j个测试样本对应的待测雷达目标所属的雷达目标类别;^表示第j 个测试样本的稀疏特征矩阵;X表示稀疏特征训练样本集，X= [Xi，X2，…，Xi，…，Χκ]，Xi表示 第i个类别的训练样本的稀疏特征矩阵所组合构成的稀疏特征训练样本子集，l<i<K，K表 示所采集的已知雷达目标的类别总数，且^[.^^???^,.,,???，.\:\]，5,.表示第1个类别的 第m个训练样本的稀疏特征矩阵，KmSN 1, N1表示第i个类别包含的训练样本的个数;Χδ, (Ctj)则表示利用稀疏系数向量对应于第i个类别已知雷达目标的类别稀疏系数向量S 1 (A)通过稀疏线性方程对第j个测试样本的稀疏特征矩阵进行稀疏重构的重构式；I I · I |2 为L2范数运算符；
[0034] 由此，实现对各个待测雷达目标的类别识别。
[0035] 相比于现有技术，本发明具有如下有益效果：
[0036] 本发明基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，特别适用于针对低图 像质量的SAR图像进行目标识别，其先通过对SAR图像进行低秩矩阵恢复，得到对应的SAR低 秩图像，使得绝大部分的SAR目标信息保留在SAR低秩图像中而去除噪声，并结合SAR图像中 目标本身的散射特征信息仍处于低维空间、能够通过SAR低秩图像得以保留的特点，结合稀 疏表示方法进行图像的目标识别，通过理论分析和实验数据均证明，能够有效的帮助提升 识别正确率，解决现有技术中针对低图像质量SAR图像进行目标识别的准确率较低的问题。
【附图说明】
[0037] 图1为本发明基于低秩矩阵恢复和稀疏表示的SAR目标识别方法的流程图。
[0038] 图2为BMP2、BTR70、T72三个不同类别雷达目标的可见光图像。
[0039] 图3为BMP2、BTR70、T72三个不同类别雷达目标的SAR图像。
[0040] 图4为Odb噪声SAR图像与其进行低秩矩阵恢复所得低秩图像的对比图。
[0041] 图5为IOdb噪声SAR图像与其进行低秩矩阵恢复所得低秩图像的对比图。
[0042]图6为20db噪声SAR图像与其进行低秩矩阵恢复所得低秩图像的对比图。
【具体实施方式】
[0043]下面结合附图和实施例对本发明的技术方案作进一步的说明。
[0044]为了解决现有技术中针对低图像质量SAR图像进行目标识别的准确率较低的问 题，需要找到适合于SAR目标识别的SAR图像降噪处理方案。近年来，E.Candes等人提出了一 种矩阵低秩恢复算法(Low Rank Matrix Recovery，LRMR)，其能够将训练样本矩阵分解为 一个低秩逼近矩阵和一个稀疏误差矩阵。受这一思路的启发，通过分析得知，虽然低图像质 量的SAR图像具有较强噪声，但是SAR图像中目标本身的散射特征信息仍处于低维空间，因 此可以采用低秩矩阵恢复的办法将低图像质量的SAR图像分解为低秩图像和误差图像两部 分;这样低秩图像中保留了绝大部分的SAR目标信息，而误差图像中主要是噪声。然后，用恢 复出来的测试样本低秩图像进行SAR目标识别，就可以帮助提高识别正确率。然而，通过恢 复提取出低图像质量SAR图像的低秩图像后，结合什么识别方式来进行SAR目标识别能够体 现低秩图像中所保留的SAR目标信息呢?近年来，随着压缩感知的理论的发展，基于压缩感 知的稀疏表示（sparse representation，SR)引起了信号处理和模式识别领域的很多研究 人员的关注。稀疏表示理论表明信号可以通过字典中的原子的线性组合来表示，并且这些 原子的分布是稀疏的，即大部分的系数是零或者接近零，只有与输入信号有较大相关性的 原子的对应的系数才不为零，稀疏系数的这种稀疏性蕴含着鉴别信息，并且能够在非零数 据元素较少的前提下描述目标的最主要特征信息。因此稀疏表示和稀疏重构理论在人脸识 另IJ、医学肿瘤识别和SAR图像目标识别等领域中都得到了较为广泛的应用。基于SAR图像中 目标本身的散射特征信息仍处于低维空间、能够通过SAR图像的低秩图像的以保留的特点， 适合采用基于稀疏表示的SAR目标识别方法。
[0045] 一、低秩矩阵恢复理论。
[0046] 对于某个目标类别，将其所有SAR图像测试样本表示为列矢量V1，构成数据矩阵D =[VI，V2，…，Vn] eRmXn，i = 1，2，…，n，其中m表示测试样本特征维数，η表示该类的测试样本 数目。对数据矩阵D进行分解，有效地恢复低秩矩阵AeR mxn和误差矩阵EeRmxn，因此该数据 矩阵可以表示为：
[0047] D=A+E； (1)
[0048] 其中，A代表没有散焦或者没有噪声的目标SAR图像，它具有目标特性决定的低秩 子空间结构，E代表误差项，表示测试图像中的噪声、散焦部分，这些误差项破坏了 SAR图像 的低秩结构。可以假设每一幅测试图像中某一小部分像素代表上述误差，也就是说误差项 可以比较大，但是约定在空间域里是稀疏的。因此可以将数据矩阵D分解为低秩矩阵和稀疏 误差矩阵。利用拉格朗日公式，可以将上述问题转化为下面的优化问题：
[0049]
⑵.
[0050] 其中，I I · I Io代表LO范数，表示矩阵中的非零项数目。λ>0是正则化参数，用来权 衡矩阵A的秩和误差矩阵E的稀疏性之间的关系。给定一个合适的λ，希望能够有效地恢复一 对(A，E)。然而，公式(5)是一个NP-hard问题，在满足特定的几何条件下LO范数可以松弛为 Ll范数的优化问题，同时秩可以转化为核范数，因此公式（5)可以转化为下面的凸优化问 题：
[0051 ]
.⑴.
[0052]其中，I I · I ^代表以范数，表示矩阵中所有元素绝对值之和。U · I I*表示核范数， 表示矩阵的奇异值之和。λ>0是正则化参数，用来权衡矩阵A的秩和误差矩阵E的稀疏性之 间的关系。
[0053]合理地选取正则化参数λ是有效解决式(6)的前提保证。对于实际的问题，可以根 据先验知识来选择λ以改善算法的性能。例如已知E是稀疏的，增大的λ可以允许恢复出一个 具有更大秩的矩阵Α。文献"ChenY,Zhang S，Zhao X.FacialExpression Recognition via Non-Negative Least-Squares Sparse Coding[J]. Information,2014,5(2):305-318?ΙΕ 明并确定了一个合理的选择，S
，允许通过轻微调整来获得最佳的结果； 也就是说，公式(3)能够有效地恢复低秩矩阵A和稀疏误差矩阵E，可以利用增广拉格朗日乘 子(Augmented Lagrange Multiplier，ALM)解决式（6)中的优化问题。
[0054] 二、稀疏表示理论。
[0055] 稀疏表示的理论表明输入信号可以通过一组基向量的线性组合来表示。给定k类 目标的足够的样本，构成
，任何一个新的测试样本y e Rm都可 以很好的通过属于它那一类的训练样本来很好的表示：
[0056] ； ⑷
[0057] >由于一个新的测试样本属于哪一类是不知道的， 所以k类不同类别的η个训练样本被用来作为基向量。
[0058] 衣后y就可以用所有训练样本来 线性表示：
[0059]
[0060] 其中(61 = 1;?,?, . . .aK]eRn。由于m<n，因此(5)式的解就不是唯一的。通行的做法 是找到最稀疏的解：
[0061 ]
[0062] 其中，I I · I Io和I I · I |2分别是LO范数和L2范数，ε是误差阈值，方程(6)的求解是 NP-hard问题。压缩感知理论的发展表明，如果解沒足够的稀疏，那么（6)式就可以被看成是 一个Ll范数最小的优化问题：
[0063]
(7>
[0064]其中I I · I |ι是Ll范数。求解范数的这一最优化问题可以采用正交匹配追踪 (Orthogonal Matching Pursuit，0ΜΡ)等算法。
[0065] 而基于稀疏表示的SAR图像目标识别方法的思路是：用所有类别已知雷达目标的 全部训练样本的稀疏特征矩阵的集合构成字典，对测试目标SAR图像稀疏特征矩阵在该字 典下进行稀疏表示，求得稀疏编码系数向量;根据类内系数对测试样本进行重构后与测试 样本的误差最小来进行识别。也就是说，得到测试样本的稀疏编码向量&后，将其判决为各 类子空间系数重构误差最小的类。即：
[0066] min Vi(j) ^|| y -XSi(ft) [],：; :(.:8) i
[0067] 第i类子空间系数AM)是和稀疏编码向量?:维数相同的矢量，且其中除第i类训练 样本对应的系数与&中对应的系数相同外，其余系数为零。
[0068] 根据上述的技术思路，本发明提出了一种基于低秩矩阵恢复和稀疏表示的SAR目 标识别方法。该方法流程如图1所示，包括如下步骤：
[0069] 1)对于多类不同的已知雷达目标，分别针对每一类采集多个已知雷达目标的SAR 图像作为训练样本，并分别获取各个类别中各个训练样本的SAR图像像素值向量;每个训练 样本的SAR图像像素值向量，是由作为训练样本的SAR图像的各个像素值排列构成。
[0070] 2)分别将每一类别中各个训练样本的SAR图像像素值向量组合构成相应类别的训 练样本图像数据矩阵，并采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩 阵所对应的训练样本图像低秩矩阵，进而根据每个训练样本图像低秩矩阵中对应于各个训 练样本的低秩图像像素值向量，分别恢复得到各个类别中各个训练样本对应的SAR低秩图 像。
[0071] 该步骤中，采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩阵所 对应的训练样本图像低秩矩阵的具体方式为：
[0072]针对于第i个类别的训练样本图像数据矩阵表示所采集的已知雷 达目标的类别总数，建立如下的低秩矩阵恢复目标函数：
[0073]
[0074]其中，Ax,4ΡΕΧ>1分别表示对第i个类别的训练样本图像数据矩阵D x,i进行低秩矩阵 恢复所得到的训练样本图像低秩矩阵和训练样本图像误差矩阵为针对第i个类别的训 练样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数，且(XA 1SUl I · I I*为核范数 运算符，M · I Ii为Li范数运算符；
[0075] 通过对所述低秩矩阵恢复目标函数求解，即可得到第i个类别的训练样本图像数 据矩阵Dx,i所对应的训练样本图像低秩矩阵A x,1;由此，分别求解得到各个类别的训练样本 图像数据矩阵所对应的训练样本图像低秩矩阵。
[0076] 而作为优选方式，针对第i个类别的训练样本图像数据矩阵的低秩矩阵恢复所设 定的正则化参数Ai的取值为i,. = 1/>腦(从，％);其中1表示训练样本的SAR图像像素值向 量的维数，N1表示针对第i个类别已知雷达目标所采集的训练样本的个数。
[0077] 3)分别提取各个类别中各个训练样本对应的SAR低秩图像的稀疏特征矩阵，将针 对每个训练样本所提取的稀疏特征矩阵作为一个稀疏特征训练样本，从而由提取的各个类 别各个训练样本的稀疏特征矩阵的集合构成稀疏特征训练样本集。
[0078] 该步骤中基于雷达目标图像的何种特征数据来建立其稀疏特征，是需要根据实际 应用个情况进行选择问题。对于一幅雷达目标SAR图像，根据物理光学近似的原理，目标的 散射回波近似成多个散射中心的回波响应总和(参见现有技术文献"Potter，L.C.;Ertin， E · ; Parker，J·T· ; Cetin，M·Sparsity and compressed sensing in radar imaging.Proceedings of the IEEE 2010,98，1006_1020."），这些稀疏散射中心提供了简 明的、物理相关的目标特性描述（参见现有技术文献" Petin, M.;Karl，W.C.;Castanon， D.A.Evaluation of a regularized SAR imaging technique based on recognition-oriented features. In Proc.SPIE 4053,Algorithms for Synthetic Aperture Radar Imagery VII，0rlando，FL，USA，24 April 2000;pp.40-51"）。经典的散射中心参数建模方 法是基于近似的物理光学模型，其参数包括像素值特征、像素频谱、入射角、接收角和极化 等;这些散射中心的参数都可以选择作为SAR目标识别的有效稀疏特征，用以建立稀疏特征 矩阵。
[0079] 4)针对各个待测雷达目标，采集各个待测雷达目标的SAR图像作为测试样本，并分 别获取各个测试样本的SAR图像像素值向量;每个测试样本的SAR图像像素值向量，是由作 为测试样本的SAR图像的各个像素值排列构成。
[0080] 5)分别将各个测试样本的SAR图像像素值向量组合构成测试样本图像数据矩阵， 并采用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本图像低秩矩阵，进 而根据测试样本图像低秩矩阵中对应于各个训测试样本的低秩图像像素值向量，分别恢复 得到各个测试样本对应的SAR低秩图像。
[0081 ]该步骤中，采用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本 图像低秩矩阵具体方式为：
[0082] 针对于测试样本图像数据矩阵Dy，建立如下的低秩矩阵恢复目标函数：
[0083]
[0084] 其中，Ay和Ey分别表示对测试样本图像数据矩阵Dy进行低秩矩阵恢复所得到的训 练样本图像低秩矩阵和训练样本图像误差矩阵;λ γ为针对测试样本图像数据矩阵的低秩矩 阵恢复所设定的正则化参数，且〇<xy<i; I I · I I*为核范数运算符，I I · I |ι为Li范数运算 符；
[0085] 通过对所述低秩矩阵恢复目标函数求解，即可得到测试样本图像数据矩阵Dy所对 应的测试样本图像低秩矩阵Ay。
[0086] 而作为优选方式，针对测试样本图像数据矩阵的低秩矩阵恢复所设定的正则化参 数取值为
其中，M表示测试样本的SAR图像像素值向量的维数，Ny 表示测试样本的个数。
[0087] 需要说明的是，该步骤中可以选择不同类别的待测雷达目标的测试样本组合来构 建测试样本图像数据矩阵;但相对而言，如果构建测试样本图像数据矩阵所采用的各个待 测雷达目标的测试样本均是来自同一个雷达目标类别(虽然具体所属类别尚未知），则更有 助于提升对测试样本图像数据矩阵进行地址矩阵恢复后的去噪效果，进而更有助于提高对 测试样本的目标识别准确率。
[0088] 6)分别提取各个测试样本对应的SAR低秩图像的稀疏特征矩阵。同样，该步骤中， 测试样本图像的像素值特征、像素频谱、入射角、接收角和极化等散射中心参数都可以选择 作为SAR目标识别的有效稀疏特征，用以建立稀疏特征矩阵。只是，针对训练样本和测试样 本所选择建立稀疏特征矩阵的参数应当相同。
[0089] 7)针对每个测试样本的稀疏特征矩阵，利用稀疏特征训练样本集中的各个稀疏特 征训练样本建立稀疏线性方程，对相应测试样本的稀疏特征矩阵进行稀疏线性表示，并求 解得到该稀疏线性方程的系数向量，作为相应测试样本的稀疏系数向量;从而得到各个测 试样本的稀疏系数向量。
[0090] 该步骤中建立的稀疏线性方程可以表示为：
[0091 ] yj = Xaj = Xiaji+X2aj2+· · -+Xiaji+- --+XkOjk ；
[0092]其中，^表示第j个测试样本的稀疏特征矩阵;X表示稀疏特征训练样本集，X= [X1, X2，…，X1，…，Xk] ,X1表示第i个类别的训练样本的稀疏特征矩阵所组合构成的稀疏特征训 练样本子集，I < i < K，K表示所采集的已知雷达目标的类别总数，且 木=[χ?ρχ,,ρ…，，气表示第i个类别的第I ll个训练样本的稀疏特征矩阵，1彡Ill SN1, N1表示第i个类别包含的训练样本的个数;表示第j个测试样本的稀疏系数向量，(^ =[aji，a j2，…，CJji，…，ajK]，α#表示第j个测试样本的稀疏系数向量CX j对应于第i个类别稀疏 特征训练样本子集的稀疏系数子向量，…，α??,?, ι'??,Λ, 表不弟j个测 试样本的稀疏系数向量A中对应于第i个类别第m个训练样本的稀疏特征矩阵的稀疏系数。 [0093] 根据压缩感知理论，该稀疏线性方程可以采用正交匹配追踪（Orthogonal Matching Pursuit，OMP)等算法求解，从而得到测试样本的稀疏系数向量。
[0094] 8)针对每个测试样本的稀疏系数向量，分别提取其稀疏系数向量中对应于稀疏特 征训练样本集中每一类已知雷达目标的类别稀疏系数向量，然后分别计算利用每一类已知 雷达目标对应的类别稀疏系数向量通过稀疏线性方程对相应测试样本的稀疏特征矩阵进 行稀疏重构的重构误差，将重构误差最小的类别稀疏系数向量所对应的一个已知雷达目标 类别判定为相应测试样本对应的待测雷达目标所属的雷达目标类别；由此，实现对各个待 测雷达目标的类别识别。
[0095] 该步骤具体为：
[0096] 81)对于第j个测试样本的稀疏系数向量a」，分别提取其中对应于稀疏特征训练样 本集中每一类已知雷达目标的类别稀疏系数向量;其中，稀疏系数向量A对应于第i个类别 已知雷达目标的类别稀疏系数向量S1(A)为：
[0097]
[0098] 兵〒，关別怖龅系双|nj重OUQj爛维双与弟J个测讯秤本的稀疏系数向量Qj的维数 相同，0V,表示第j个测试样本的稀疏系数向量^中对应于第i个类别第m个训练样本的稀 疏特征矩阵的稀疏系数，且类别稀疏系数向量&(A)中对应于非第i个类别已知雷达目标的 稀疏系数的值均为零；1彡i<K，K表示所采集的已知雷达目标的类别总数；1彡m彡N i ,Ni表 示第i个类别包含的训练样本的个数；
[0099] 82)分别计算利用每一类已知雷达目标对应的类别稀疏系数向量通过稀疏线性方 程对相应测试样本的稀疏特征矩阵进行稀疏重构的重构误差，将重构误差最小的类别稀疏 系数向量所对应的一个已知雷达目标类别判定为相应测试样本对应的待测雷达目标所属 的雷达目标类别;其中，确定第j个测试样本对应的待测雷达目标所属的雷达目标类别的具 体计算方式为：
[0100]
[0101 ]其中，Ij表示第j个测试样本对应的待测雷达目标所属的雷达目标类别;^表示第j 个测试样本的稀疏特征矩阵;X表示稀疏特征训练样本集，X= [Xi，X2，…，Xi，…，Χκ]，Xi表示 第i个类别的训练样本的稀疏特征矩阵所组合构成的稀疏特征训练样本子集，l<i<K，K表 示所采集的已知雷达目标的类别总数，且义=[' 141^_，',,，~4.\]，1.表示第1个类别的 第ru个训练样本的稀疏特征矩阵，Kn 1SN1,N1表示第i个类别包含的训练样本的个数； ? I |2为L2范数运算符;XS1(Cij)则表示利用稀疏系数向量a」对应于第i个类别已知雷达目 标的类别稀疏系数向量^"七)通过稀疏线性方程对第j个测试样本的稀疏特征矩阵进行稀 疏重构的重构式，即：
[0102]
[0103] 由此，实现对各个待测雷达目标的类别识别。
[0104] 本发明基于低秩矩阵恢复和稀疏表示的SAR目标识别方法可以应用于基于计算机 编程自运行的雷达目标识别系统，实现自动化的雷达目标识别。
[0105] 下面通过实施例对本发明的技术方案作进一步的说明。
[0106] 实施例：
[0107]本实施例使用MSTAR数据库来做实验，它是由美国圣地亚哥国家实验室X波段的 SAR系统获取的实测数据，它具有0.3m X 0.3m的分辨率，分别在15度和17度的俯仰角下获得 的。本实施例中使用MSTAR数据库中BMP2(步兵坦克）、BTR70(装甲运兵车）、T72(T-72型主站 坦克)这三类目标来进行实验。BMP2、BTR70、T72三个不同类别雷达目标的可见光图像分别 如图2中（2a)、（ 2b)、（ 2c)所示，而BMP2、BTR70、T72三个不同类别雷达目标的SAR图像则分别 如图3中（3a)、（3b)、（3c)所示，可以看到，可见光图像与SAR雷达图像有明显的差别，并且 SAR图像普遍具有较低的信噪比，对雷达目标的识别造成了不利影响。MSTAR数据库中每一 幅SAR图像的像素密度都是128行X 128列，在本实施例中用17度俯仰角的SAR图像数据来用 作训练样本数据，15度俯仰角的SAR图像数据来做测试样本数据，训练样本和测试样本的数 量如表1所示：
[0108]表1
[0110] 同时，分别针对本实施例中的训练样本和测试样本分别添加噪声形成信噪比为 0db、10db、20db、30db的噪声训练样本和噪声测试样本，然后采用本发明基于低秩矩阵恢复 和稀疏表示的SAR目标识别方法，利用各不同信噪比的噪声训练样本对不同信噪比的噪声 测试样本进行目标类别的识别;在采用本发明方法进行目标识别处理的过程中，图4、图5、 图6分别示出了对Odb、IOdb、20db噪声图像进行低秩矩阵恢复的效果对比情况，其中图4的 (4&)、图5的（5 &)、图6的(6&)分别为0(113、10(113、20(113噪声3厶1?图像，而图4的(413)、图5的（513)、 图6的（6b)分别为0db、10db、20db噪声SAR图像通过低秩矩阵恢复后得到的低秩图像，可以 看到通过低秩矩阵恢复后得到的低秩图像的信噪比均得以明显增强。同时，为了作为对比， 还采用了不使用低秩矩阵恢复的稀疏表示分类识别方法，利用各不同信噪比的噪声训练样 本对不用信噪比的噪声测试样本进行目标类别的识别，与本发明方法进行识别率对比，统 计得到两种方法的识别率如表2所示。
[0111] 表2
[0113]可以看到，针对于不同信噪比的噪声SAR图像目标识别，采用本发明基于低秩矩阵 恢复和稀疏表示的SAR目标识别方法的识别正确率均有明显的提高。
[0114]综上所述，本发明基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，特别适用于 针对低图像质量的SAR图像进行目标识别，其先通过对SAR图像进行低秩矩阵恢复，得到对 应的SAR低秩图像，使得绝大部分的SAR目标信息保留在SAR低秩图像中而去除噪声，并结合 SAR图像中目标本身的散射特征信息仍处于低维空间、能够通过SAR低秩图像得以保留的特 点，结合稀疏表示方法进行图像的目标识别，通过理论分析和实验数据均证明，能够有效的 帮助提升识别正确率，解决现有技术中针对低图像质量SAR图像进行目标识别的准确率较 低的问题。
[0115]最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较 佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技 术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本 发明的权利要求范围当中。
【主权项】
1. 基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在于，包括如下步骤： 1) 对于多类不同的已知雷达目标，分别针对每一类采集多个已知雷达目标的SAR图像 作为训练样本，并分别获取各个类别中各个训练样本的SAR图像像素值向量;每个训练样本 的SAR图像像素值向量，是由作为训练样本的SAR图像的各个像素值排列构成； 2) 分别将每一类别中各个训练样本的SAR图像像素值向量组合构成相应类别的训练样 本图像数据矩阵，并采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩阵所 对应的训练样本图像低秩矩阵，进而根据每个训练样本图像低秩矩阵中对应于各个训练样 本的低秩图像像素值向量，分别恢复得到各个类别中各个训练样本对应的SAR低秩图像； 3) 分别提取各个类别中各个训练样本对应的SAR低秩图像的稀疏特征矩阵，将针对每 个训练样本所提取的稀疏特征矩阵作为一个稀疏特征训练样本，从而由提取的各个类别各 个训练样本的稀疏特征矩阵的集合构成稀疏特征训练样本集； 4) 针对各个待测雷达目标，采集各个待测雷达目标的SAR图像作为测试样本，并分别获 取各个测试样本的SAR图像像素值向量;每个测试样本的SAR图像像素值向量，是由作为测 试样本的SAR图像的各个像素值排列构成； 5) 分别将各个测试样本的SAR图像像素值向量组合构成测试样本图像数据矩阵，并采 用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本图像低秩矩阵，进而根 据测试样本图像低秩矩阵中对应于各个训测试样本的低秩图像像素值向量，分别恢复得到 各个测试样本对应的SAR低秩图像； 6) 分别提取各个测试样本对应的SAR低秩图像的稀疏特征矩阵； 7) 针对每个测试样本的稀疏特征矩阵，利用稀疏特征训练样本集中的各个稀疏特征训 练样本建立稀疏线性方程，对相应测试样本的稀疏特征矩阵进行稀疏线性表示，并求解得 到该稀疏线性方程的系数向量，作为相应测试样本的稀疏系数向量;从而得到各个测试样 本的稀疏系数向量； 8) 针对每个测试样本的稀疏系数向量，分别提取其稀疏系数向量中对应于稀疏特征训 练样本集中每一类已知雷达目标的类别稀疏系数向量，然后分别计算利用每一类已知雷达 目标对应的类别稀疏系数向量通过稀疏线性方程对相应测试样本的稀疏特征矩阵进行稀 疏重构的重构误差，将重构误差最小的类别稀疏系数向量所对应的一个已知雷达目标类别 判定为相应测试样本对应的待测雷达目标所属的雷达目标类别；由此，实现对各个待测雷 达目标的类别识别。2. 根据权利要求1所述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在 于，所述步骤2)中，采用低秩矩阵恢复方法分别求得每一类别的训练样本图像数据矩阵所 对应的训练样本图像低秩矩阵的具体方式为： 针对于第i个类别的训练样本图像数据矩阵Dx,n l<i<K，K表示所采集的已知雷达目 标的类别总数，建立如下的低秩矩阵恢复目标函数：其中，Αχ,4ΡΕΧ>1分别表示对第i个类别的训练样本图像数据矩阵Dx>1进行低秩矩阵恢复 所得到的训练样本图像低秩矩阵和训练样本图像误差矩阵;针对第i个类别的训练样 本图像数据矩阵的低秩矩阵恢复所设定的正则化参数，且0<M<1; | | · | I*为核范数运算 符，11 · I li为L1范数运算符； 通过对所述低秩矩阵恢复目标函数求解，即可得到第i个类别的训练样本图像数据矩 阵Dx>1所对应的训练样本图像低秩矩阵Ax,1;由此，分别求解得到各个类别的训练样本图像 数据矩阵所对应的训练样本图像低秩矩阵。3. 根据权利要求2所述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在 于，针对第i个类别的训练样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数的取;其中，Μ表示训练样本的SAR图像像素值向量的维数，Ni表示针对第 i个类别已知雷达目标所采集的训练样本的个数。4. 根据权利要求1所述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在 于，所述步骤5)中，采用低秩矩阵恢复方法求得测试样本图像数据矩阵所对应的测试样本 图像低秩矩阵具体方式为： 针对于测试样本图像数据矩阵Dy，建立如下的低秩矩阵恢复目标函数：其中，Ay和Ey分别表示对测试样本图像数据矩阵Dy进行低秩矩阵恢复所得到的训练样 本图像低秩矩阵和训练样本图像误差矩阵;λγ为针对测试样本图像数据矩阵的低秩矩阵恢 复所设定的正则化参数，且〇<\<1; 11 · 11*为核范数运算符，11 · | |ι为L1范数运算符； 通过对所述低秩矩阵恢复目标函数求解，即可得到测试样本图像数据矩阵Dy所对应的 测试样本图像低秩矩阵A y。5. 根据权利要求4所述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在 于，针对测试样本图像数据矩阵的低秩矩阵恢复所设定的正则化参数~7的取值为;其中，Μ表示测试样本的SAR图像像素值向量的维数，Ny表示测试样本 的个数。6. 根据权利要求1所述基于低秩矩阵恢复和稀疏表示的SAR目标识别方法，其特征在 于，所述步骤8)具体为： 81) 对于第j个测试样本的稀疏系数向量α」，分别提取其中对应于稀疏特征训练样本集 中每一类已知雷达目标的类别稀疏系数向量;其中，稀疏系数向量^对应于第i个类别已知 雷达目标的类别稀疏系数向量"㈧为：其中，类别稀疏系数向量的维数与第j个测试样本的稀疏系数向量^的维数相 同表示第j个测试样本的稀疏系数向量^中对应于第i个类别第m个训练样本的稀疏 特征矩阵的稀疏系数，且类别稀疏系数向量&(A)中对应于非第i个类别已知雷达目标的稀 疏系数的值均为零；1彡i<K，K表示所采集的已知雷达目标的类别总数；1彡mSNnNi表示 第i个类别包含的训练样本的个数； 82) 分别计算利用每一类已知雷达目标对应的类别稀疏系数向量通过稀疏线性方程对 相应测试样本的稀疏特征矩阵进行稀疏重构的重构误差，将重构误差最小的类别稀疏系数 向量所对应的一个已知雷达目标类别判定为相应测试样本对应的待测雷达目标所属的雷 达目标类别;其中，确定第j个测试样本对应的待测雷达目标所属的雷达目标类别的具体计 算方式为：其中，1』表示第j个测试样本对应的待测雷达目标所属的雷达目标类别;η表示第j个测 试样本的稀疏特征矩阵;X表示稀疏特征训练样本集，x= [Χι，X2，…，Xi，…，Xk]，Xi表示第i个 类别的训练样本的稀疏特征矩阵所组合构成的稀疏特征训练样本子集，l<i<K，K表示所 采集的已知雷达目标的类别总数，且不=[~^.^"，\",，"％.^,]，\表示第丨个类别的第1 11 个训练样本的稀疏特征矩阵，表示第i个类别包含的训练样本的个数;Χδ^α」） 则表示利用稀疏系数向量Α对应于第i个类别已知雷达目标的类别稀疏系数向量通 过稀疏线性方程对第j个测试样本的稀疏特征矩阵进行稀疏重构的重构式；I I · I |2为L2范 数运算符； 由此，实现对各个待测雷达目标的类别识别。
【文档编号】G06K9/00GK106056070SQ201610368234
【公开日】2016年10月26日
【申请日】2016年5月26日
【发明人】张新征, 王亦坚, 常云鹤
【申请人】重庆大学