三维空间中高超声速弱目标RAE‑HT‑TBD积累检测方法与流程

本发明属于雷达目标检测领域，用于解决三维空间中高超声速弱目标的检测问题。

背景技术：

在对临近空间目标的预警探测研究中，一方面临近空间飞行器所具有的高超声速特性使其可在两小时内对全球任意目标实施打击；另一方面，临近空间飞行器高超声速运动所产生的激波等离子体会使得目标的雷达反射截面积发生严重衰减，导致目标回波十分微弱，从而使目标获得较强的隐身性能，进一步加大了雷达探测难度。临近空间飞行器的以上两点特性使得现有的预警系统很难对其进行有效拦截。因此，研究临近空间高超声速弱目标的检测跟踪问题对于增强国家空天安全具有重要意义。

基于hough的检测前跟踪(ht-tbd)是一种检测强杂波环境下微弱目标的有效技术，具有对随机噪声鲁棒性强、对目标位置不确定性及局部缺损不敏感等优点，在雷达探测领域得以广泛应用。tbd采用多帧数据批处理的思想，对多个采样周期数据进行非相参积累，提高了目标航迹信噪比，避免了传统单帧门限检测方法可能导致的目标潜在有用信息的丢失问题。目前，针对三维空间目标检测的ht-tbd积累检测方法主要采用直接变换法以及投影映射法两种方式：由于三维空间直线至少需要四个参数进行解析表示，因此直接进行三维hough变换对应的参数空间至少是四维空间，高维参数空间产生很大计算负担，不利于工程实际应用；投影映射方法将三维点迹投影至xyz坐标系下的多个直角坐标平面进行二维hough变换并通过融合各平面检测结果得到最终三维检测航迹，降维映射处理虽然减小了计算量，但是此种投影方式存在两个主要缺陷：一方面，这种投影变换机制极易在各直角坐标平面积累噪声，从而产生信噪比损失，并且信噪比损失随着观测区域的增大近似呈对数增加；另一方面，对于远距离临近空间目标，这种投影变换机制会产生距离量测与角度量测的耦合问题，即使角度量测误差很小，也会导致较大横向误差。例如，设雷达测角误差为0.2°，目标径向距离为600km，则目标横向误差标准差可达2.1km，根据高斯误差分布定律中的3σ准则，目标横向误差最大可达6.3km，使得利用直角坐标平面hough变换难以实现对临近空间远距离目标的有效检测。因此，本发明提出一种三维空间中高超声速弱目标rae-ht-tbd积累检测方法，在保证优良检测性能的同时具有较小计算量，可以适用于工程应用。

技术实现要素：

本发明的目的是利用三级降维分解的方式逐级完成点迹筛选，提出一种三维空间中高超声速弱目标rae-ht-tbd积累检测方法，通过依次在径向距离(range)-时间平面、方位角(azimuth)-时间平面、仰角(elevation)-时间平面内的二维hough变换可以解决临近空间内对于三维位置量测高超声速弱目标检测过程所具有的计算量大、检测概率低的问题。

本发明提出的三维空间中高超声速弱目标rae-ht-tbd积累检测方法，其特征在于，包括以下技术措施：

步骤一、在一个较高的虚警率pfa下设置初始门限η1，只保留功率超过门限的点迹，删除其余点迹，从而过滤掉大部分无关量测以消除部分杂波影响，得到过门限后的量测数据，包含量测点的距离、方位、时间和能量信息，设定其初始门限为：

η1＝-ln(pfa)

接下来，将进行基于hough变换的三重点迹筛选，由于各级点迹筛选的流程类似，下面将仅以第一级点迹筛选为例说明本发明提出的方法。

其中，较高的虚警率pfa是指10^-1量级。

步骤二、对于投影映射到距离-时间平面的点迹进行规格化处理：由于此时距离-时间平面中横纵坐标量级相差悬殊，会使得量级较小一维数据信息丢失，无法实现对目标的有效检测，因此需要求得下式所示的规格化系数μ，并得到规格化后的r-t数据(t,r/μ)：

其中，[lg(|rmax/tmax|)]表示大于lg(|rmax/tmax|)的最小整数，rmax,tmax分别表示r与t中的最大值。

步骤三、hough变换参数空间离散化：为了增大误差容限进而增强算法检测能力，使近似共线的量测点在参数空间内仍能交于同一单元，根据雷达和目标位置以及雷达的测距和测角误差对参量空间进行分割，确定分割的单元数；ρ-θ参数空间ω1将被离散成个参数单元，分辨率其中，l取雷达探测范围的两倍，离散化参数单元中心坐标：

步骤四、hough变换参数初始化：建立参数空间能量积累矩阵f1(i,j)以及投票点记录矩阵f2{i,j}，用于存储各参数单元能量积累数据以及各参数单元的投票点状态信息，同时预设部分内存空间用于存储过门限参数单元逆映射后的的数据空间点迹集合并对上述矩阵进行置零初始化。

步骤五、依次选取步骤二中的过初始门限数据，通过下式将距离-时间平面中的点变换到参数空间ω1，得到相应的参数曲线ξ：

ρ＝μtcosθ+rsinθ

其中，ρ表示数据点所在直线与原点的距离，θ表示数据点和原点连线与横轴正方向的夹角。

步骤六、在参数空间进行能量积累来提升参数空间积累性能，取参数空间ω1内功率积累最大值的χ⁽¹⁾倍作为非相参积累门限只提取出能量积累值超过的参数单元集合；

假设预处理后点迹集合a中第k帧第l个点迹被映射到参数空间对应的正弦曲线为则表征k时刻该曲线经过参数空间ω1中第i行第j列参数单元的数学判别式：

其中，l＝1,2,…,mk，k＝1,2,…,n；如果点迹满足上式，则利用该点对参数单元进行能量积累；否则，点迹对参数单元的能量积累无贡献。

步骤七、hough变换逆映射：经过步骤六，只有满足能量积累值超过的参数单元所对应的数据空间点迹集合内的点才被认为是可能的目标点迹，利用投票点记录矩阵f2(i,j)得到数据空间的逆映射点迹，同时删除数据空间内除了集合外的所有点迹，从而完成第一级基于径向距离-时间平面hough变换的点迹筛选，其中，k＝1,2,…,n，表示集合内点迹个数。

步骤八、依次进行方位角-时间以及仰角-时间平面内的后两级hough变换点迹筛选：由于接下来将要依次进行的第二级基于方位角-时间平面hough变换的点迹筛选以及第三级基于仰角-时间平面hough变换的点迹筛选的基本原理和流程与第一级基于径向距离-时间hough变换的点迹筛选类似，只需以前一级筛选出来的点迹集合为输入，将对应公式中的变量由径向距离r依次替换方成位角和仰角重复步骤三至七即可逐级完成点迹筛选，因此不再详细论述。

步骤九、将各帧筛选出的数据点按照时序进行关联得到可能航迹，并针对经过航迹检验后可能存在多条航迹对应同一个目标的情况对相似航迹进行航迹合并，具体措施为：

在对检测出的可能航迹进行同帧点迹两两比较的基础上，可以得到相同点迹的个数δ。如果δ超过某一门限δ0，则合并这两条航迹。门限δ0可以表示为：

其中，表示不大于的最大整数，τ为约束因子，取值为整数，满足

在进行合并航迹时，对于两条航迹中相同时刻的不同点迹，将只保存其中回波功率更大的量测点迹。

与背景技术相比，本发明提出的三维空间中高超声速弱目标rae-ht-tbd积累检测方法具有如下有益效果：

(1)由于三维空间直线含有四个自由度，因此直接进行三维hough变换对应的参数空间至少是四维空间，这将导致算法占据大量内存空间，巨大的计算量也会使算法失去实时性。本发明利用分级降维的思想将三维临近空间点迹映射至三个二维平面进行hough变换，本发明的计算复杂度介于o(qm)和3·o(qm)之间，远低于直接法o(qm³)的计算复杂度(其中，q为数据空间待处理单元数，并假设各参数均离散化为m种取值)。因此，当远距离条件下m取值较大时，本发明参数空间维数的减小显著降低了计算复杂度，易于工程实现；

(2)为最大限度地避免传统的将三维点迹投影到两个直角坐标平面进行hough变换方法存在的噪声积累问题以及距离量测、角度量测的耦合问题，考虑到目标径向距离和角度测量数据受距离影响很小，本发明依次将三维量测点投影至径向距离-时间平面、方位角-时间平面、仰角-时间平面进行三级hough变换检测，分别对径向距离、方位角以及仰角进行解耦处理，充分利用精度较高的距离量测信息进行点迹筛选，同时，利用角度信息有效避免了不同方向上的错误积累。例如，设雷达测角误差为0.2^°，目标径向距离为600km，则目标横向误差标准差可达2.1km，根据高斯误差分布定律中的3σ准则，目标横向误差最大可达6.3km，而解耦后的径向距离误差大约只有200m至400m。因此，本发明可以最大限度地减小噪声积累以及远距离条件下距离、角度量测耦合引起的较大横向误差所产生的不利影响，提升检测精度。

附图说明

附图1是本发明提出的三维空间中高超声速弱目标rae-ht-tbd积累检测方法的整体流程图；

附图2是本发明实施例中强杂波环境下近似直线运动的临近空间目标三维量测点迹图；

附图3是本发明实施例中各级点迹筛选流程图(以第一级点迹筛选为例)；

附图4是本发明实施例中各级点迹筛选后对应时间平面内的剩余点迹映射图；

附图5是本发明实施例中各级hough变换所对应的参数空间能量积累直方图；

附图6是本发明实施例中各级点迹筛选后的三维空间剩余点迹映射图；

附图7是本发明实施例中三重点迹筛选后按照时序关联的可能航迹图；

附图8是本发明实施例中经过航迹融合的最终检测航迹图。

具体实施方式

实施例条件：假设临近空间高超声速弱目标在三维空间内做近似匀速直线运动，初始位置为(350000m,350000m,40000m)，即量测向量为初始速度为(3000m/s,2500m/s,150m/s)，目标rcs＝2m²。相应地，传感器为位于坐标原点的三坐标雷达，雷达扫描周期t＝2s，发射功率pt＝10kw，天线增益g＝10⁴，工作波长λt＝0.15m，径向距离量测误差为200m，方位角量测误差为0.2°，仰角量测误差为0.2°，量测积累帧数n＝7。在每帧数据中，杂波密度λc＝800，杂波个数服从泊松分布，杂波位置在雷达视域内服从均匀分布，信噪比snr＝6db。参数空间ω1、ω2、ω3均离散化为180×300个分辨单元，能量积累门限均取各自参数空间ωi内功率积累最大值的0.85倍，即χ⁽¹⁾＝χ⁽²⁾＝χ⁽³⁾＝0.85，其中，i＝1,2,3。约束因子τ＝1，航迹合并门限δ0＝4。结合说明书附图1，本发明的具体处理流程如下所示：

步骤1：设置初始门限：

由于相参处理系统输出的原始量测的数据量非常大，如附图2所示，因此需要进行预处理；为了尽可能多的保留目标信息，通常在一个较高的虚警率下进行初始门限检测，过滤掉大部分无关量测以消除部分杂波影响，得到过门限后的量测数据；对于相参积累后经平方率检波得到的雷达数据，根据较高的虚警概率，设置一个第一门限，从而消除部分干扰点迹，设定其初始门限

η1＝-ln(pfa)＝-ln(0.01)≈-4.6

步骤2：规格化处理：

规格化系数

步骤3：hough变换参数空间离散化：

为了增大误差容限进而增强算法检测能力，使近似共线的量测点在参数空间内仍能交于同一单元，根据雷达和目标位置以及雷达的测距和测角误差将参数空间ω1将被离散成个参数单元，分辨率l≈600km，离散化参数单元中心坐标：

步骤4：hough变换参数初始化：

能量积累矩阵f1(i,j)＝zeros(180,300)；

投票点记录矩阵f2{i,j}＝zeros(180,300)。

步骤5：通过下式将距离-时间平面中的点变换到参数空间ω1，得到相应的参数曲线ξ：

ρ＝μtcosθ+rsinθ

步骤6：设立能量积累门限

能量积累门限为各自参数空间ωi内功率积累最大值的0.85倍，即χ⁽¹⁾＝χ⁽²⁾＝χ⁽³⁾＝0.85，其中，i＝1,2,3。

特别地，表征k时刻曲线经过参数空间ω1中第i行第j列参数单元的数学判别式：

其中，l＝1,2,…,mk，k＝1,2,…,n；如果点迹满足上式，则利用该点对参数单元进行能量积累；否则，点迹对参数单元的能量积累无贡献；hough变换后的各级点迹筛选流程图(以第一级点迹筛选为例)以及各级点迹筛选后对应时间平面内的剩余点迹映射图分别如附图3、附图4所示；经过参数空间能量积累之后的参数空间积累直方图如附图5所示。

步骤7：hough变换逆映射：

利用投票点记录矩阵f2{i,j}保存满足能量积累值超过的参数单元所对应的数据空间点迹集合完成第一级基于径向距离-时间平面hough变换的点迹筛选，各级点迹筛选后的三维空间剩余点迹映射图如附图6所示。

步骤8：按照上述流程依次进行方位角-时间以及仰角-时间平面内的后两级hough变换点迹筛选：由于接下来将要依次进行的第二级基于方位角-时间平面hough变换点迹筛选以及第三级基于仰角-时间平面hough变换点迹筛选的基本流程与第一级类似，将不再详细论述。

步骤9：将各帧筛选出的数据点按照时序进行关联得到可能航迹，三重点迹筛选后按照时序关联的可能航迹图如附图7所示；针对经过航迹检验后可能存在多条航迹对应同一个目标的情况对相似航迹进行航迹合并，相同点迹数合并门限δ0＝4，经过航迹融合的最终检测航迹图如附图8所示。

步骤10：重复步骤2～步骤9，直至所有n帧数据全部处理并存储完毕。

从附图8可以看出，在杂波密度λc＝800、信噪比snr＝6db的条件下，经过本发明的基于分级降维rae-ht-tbd的积累检测方法大部分目标点迹可以被准确检测出，经过时序关联可以准确起始航迹，经过1000次蒙特卡洛仿真实验，检测概率可以达到80％以上，适用于工程应用，易于进一步的目标跟踪等数据处理。