基于广义特征分解的全极化高分辨距离像目标检测方法
【技术领域】
[0001] 本发明属于雷达目标检测领域，涉及一种基于广义特征分解的全极化高分辨距离 像目标检测方法，适用于全极化雷达的目标检测。
【背景技术】
[0002] 传统的单极化宽带雷达只具有一个接收通道，在检测目标时，只能通过回波的幅 度信息来检测目标。随着雷达技术的不断发展，全极化雷达得到了越来越广泛的应用。由 于全极化雷达，具有四个接收通道，能提供丰富的极化信息，全极化雷达的目标检测性能优 于传统的单极化雷达。但现有的全极化雷达目标检测算法大多针对目标的SAR图像对目标 进行检测，如基于几何扰动的单目标检测器（STD)和部分目标检测器（PTD);而在全极化高 分辨距离像中，每个距离单元内都包含有杂波，因此现有的全极化雷达目标检测算法不适 用于全极化高分辨距离像的目标检测。
【发明内容】
[0003] 针对上述现有技术的不足，本发明的目的在于提出了一种基于广义特征分解的全 极化高分辨距离像目标检测方法，该方法利用全极化雷达的极化信息，通过广义特征分解 的方法抑制杂波，能够提高全极化高分辨距离像的目标检测率。
[0004] 本发明的实现思路是：在全极化雷达获得的全极化高分辨距离像（HRRP)中，目标 区的回波中包含了目标回波和杂波；本发明针对全极化高分辨距离像，先获得训练目标回 波和训练杂波作为训练数据，并通过广义特征分解的方法，找到训练数据中信杂比相对较 大的子空间；再获得全极化雷达的测试全极化高分辨距离像作为测试数据，根据训练数据 中信杂比相对较大的子空间，对测试数据进行检测。
[0005] 为达到上述技术目的，本发明采用以下技术方案予以实现。
[0006] -种基于广义特征分解的全极化高分辨距离像目标检测方法，其特征在于，包括 以下步骤：
[0007] 步骤1，根据已知的全极化雷达的回波，获取训练目标回波和训练杂波作为训练数 据；通过广义特征分解方法，从训练数据中找到信杂比相对较大的子空间，计算得到投影矩 阵P ;
[0008] 步骤2,获取全极化雷达的测试全极化高分辨距离像作为测试数据；根据投影矩 阵P对测试数据进行目标检测。
[0009] 本发明的特点和进一步改进在于：
[0010] (1)所述步骤1的具体子步骤为：
[0011] I. 1根据已知的全极化雷达的回波获取训练目标回波，将训练目标回波划分为N 片区域，并计算得到训练目标回波的第η片区域的相干向量，n = I... N ;
[0012] 1. 2根据已知的全极化雷达的回波获取训练杂波，将训练杂波划分为J片区域，并 计算得到训练杂波的第j片区域的相干向量，j = I... J ;
[0013] I. 3计算训练目标回波的相干向量的协方差矩阵Cw和训练杂波的相干向量的协 方差矩阵c(e);
[0014] 1. 4根据训练目标回波的相干向量的协方差矩阵Cw和训练杂波的相干向量的协 方差矩阵C?，计算得到投影矩阵P。
[0015] (2)所述步骤2的具体子步骤为：
[0016] 2. 1获取全极化雷达的测试全极化高分辨距离像作为测试数据；将测试数据划分 为L个距离单元，计算测试数据的第1个距离单元的相干矩阵，并提取测试数据的第1个距 离单元的相干向量，I = I... L ;
[0017] 2. 2将测试数据的第1个距离单元的相干向量kD⑴左乘投影矩阵P，得到测试数 据的第1个距离单元的重构相干向量k' D(l)为：
[0018] k' D(l) = PkD(l)
[0019] 2. 3求第1个距离单元的重构相干向量k'D(l)的2-范数I d(1) I I2;
[0020] 2. 4设定检测门限η，如果第1个距离单元的重构相干向量k'D(l)的2-范数 I k'D (I) I 12彡τι，贝IJ将测试数据判定为目标；否贝IJ将测试数据判定为杂波。
[0021] 本发明的有益效果是：本发明的基于广义特征分解的全极化高分辨距离像目标检 测方法，通过广义特征分解的方法对全极化雷达的全极化高分辨距离像进行目标检测，不 仅能够充分抑制杂波、增强信杂比，还能够检测到一些回波能量比较弱小的目标距离单元， 提高了全极化高分辨距离像的目标检测率。
【附图说明】
[0022] 下面结合【附图说明】和【具体实施方式】对本发明作进一步详细说明。
[0023] 图1是本发明的流程图。
[0024] 图2a是实验一所用的训练数据图；
[0025] 图2b是实验一所用的测试数据图。
[0026] 图3a是三种方法对图2b的测试数据在信杂比为OdB下检测时的接收机操作特性 (ROC)对比图，横坐标为虚警率，纵坐标为检测概率；
[0027] 图3b是三种方法对图2b的测试数据在信杂比为-IOdB到30dB间检测时的曲线 下面积（AUC)对比图，横坐标为信杂比，单位为（dB)，纵坐标为曲线下面积（AUC)。
[0028] 图4a是实验二所用的训练数据图；
[0029] 图4b是实验二所用的测试数据图。
[0030] 图5a是三种方法对图4b的测试数据在信杂比为OdB下检测时的接收机操作特性 (ROC)对比图，横坐标为虚警率，纵坐标为检测概率；
[0031] 图5b是本发明方法对图4b的测试数据进行杂波抑制后的波形图；
[0032] 图5c是用单极化非相干积累法对图4b的测试数据进行杂波抑制后的波形图；
[0033] 图5d是用全极化雷达四通道总能量的非相干积累方法对图4b的测试数据进行杂 波抑制后的波形图；
[0034] 图5b-图5d中，横坐标为距离单元，纵坐标为强度。
【具体实施方式】
[0035] 参照图I，本发明的基于广义特征分解的全极化高分辨距离像目标检测方法，包括 以下步骤：
[0036] 步骤1，根据已知的全极化雷达的回波，获取训练目标回波和训练杂波作为训练数 据；通过广义特征分解方法，从训练数据中找到信杂比相对较大的子空间，计算得到投影矩 阵P。
[0037] 步骤1的具体子步骤为：
[0038] I. 1根据已知的全极化雷达的回波获取训练目标回波，将训练目标回波划分为N 片区域，并计算得到训练目标回波的第η片区域的相干向量，具体步骤为：
[0039] I. I. 1将训练目标回波划分为N片区域（可以包含一定量的杂波），对训练目标回 波的每片区域连续探测R次，其中，对训练目标回波的第η片区域进行第r次探测时获得的 散射矩阵S w (r，η)为：
[0041] 其中，上标（0)表示训练目标回波，r代表第r次探测，r= 1...R，η代表第η片 区域，n = I. . . N ;
分别为已知的全极化雷达的 水平发射水平接收（HH)通道、垂直发射水平接收（HV)通道、水平发射垂直接收（VH)通道 以及垂直发射垂直接收（VV)通道接收的目标回波数据；
[0042] I. 1. 2根据对训练目标回波的第η片区域进行第r次探测时获得的散射矩阵Sw (r，η)，计算对训练目标回波的第η片区域进行第r次探测时的散射向量k? (r，η)为：
[0044] I. I. 3根据对训练目标回波的第η片区域进行第r次探测时的散射向量k? (r，η)， 计算训练目标回波的第η片区域的相干矩阵Dw (η)为：
[0046] 其中，上标H表示共辄转置；
[0047] 1. 1.4提取训练目标回波的第η片区域的相干矩阵Dw (η)的上三角元素，组成训 练目标回波的第η片区域的相干向量为：
[0049] 其中，下标D表示相干，上标T表示转置；
[0050] 1.2,根据已知的全极化雷达的回波获取训练杂波，将训练杂波划分为J片区域， 并计算得到训练杂波的第j片区域的相干向量，具体步骤为：
[0051] 1. 2. 1将训练杂波划分为