一种目标检测与跟踪一体化方法

摘要

本发明属于雷达信号处理技术领域，尤其涉及一种目标检测与跟踪一体化方法，包括：S1：选定满足强度要求的待检测单元范围，初步确定目标数目范围；S2：选定所述目标可能存在的速度通道，根据选定的所述速度通道，设计各目标运动模型；S3：初始化第1帧数据各目标对应单元的置信度值；S4：用k表示当前处理数据的编号，计算第k次各目标对应单元的置信度值；S5：根据S4中得到的第k次数据各目标对应单元的置信度值以及第k次数据的观测值，确定目标数目和目标可能存在的位置；S6：依次类推，在第k+1次数据到来时，重复S4和S5；S7：在判决为有目标的单元处用点迹进行标注，若有N个单元，则经过K次处理后，形成一幅K×N的检测后的点迹视图。

权利要求

1.一种目标检测与跟踪一体化方法，其特征在于，包括如下步骤：S1：选定满足强度要求的待检测单元范围，初步确定目标数目范围；

S2：选定所述目标可能存在的速度通道，根据选定的所述速度通道，设计各目标运动模型；

S3：初始化第1帧数据各目标对应单元的置信度值；

S4：用k表示当前处理数据的编号，在第k次数据到来时，若k＝1，则表示当前正处理第1次数据，则将步骤S3中的初始化置信度值作为第1次数据各目标对应单元的置信度值；若k≠1，则根据第k-1次各目标对应单元内的置信度值和步骤S2中所设定的运动模型，计算第k次各目标对应单元的置信度值；

S5：根据步骤S4中得到的第k次数据各目标对应单元的置信度值以及第k次数据的观测值，迭代计算第k次各目标对应单元的置信度值，然后选定阈值对置信度值进行门限处理，得到检测后的结果，确定目标数目和目标可能存在的位置；

S6：依次类推，在第k+1次数据到来时，重复步骤S4和S5；

S7：在判决为有目标的单元处用点迹进行标注，若有N个单元，则经过K次处理后，形成一幅K×N的检测后的点迹视图；

所述步骤S4和S5中的待检单元的置信度值正比于所述待检单元的原始置信度值与假定目标存在时该单元对应的观测似然函数的乘积。

2.根据权利要求1所述方法，其特征在于，将步骤S2中可能存在目标的运动模型建模为一阶马尔科夫模型，并建立目标的置信度值转移模型。

说明书

一种目标检测与跟踪一体化方法

技术领域

[0001] 本发明属于雷达信号处理中的目标检测与跟踪技术领域，尤其涉及一种目标检测与跟踪一体化方法。

背景技术

[0002] 复杂背景下的弱目标检测是一个具有共性的问题，在雷达、声呐、图像等领域中均具有重要的应用价值。其难点在于弱目标一般具有低的信噪比，背景存在很强、统计特性复杂的杂波和噪声，因此，传统的先检测后跟踪的处理方法往往处理性能不佳。

[0003] 为了提升在复杂背景下探测弱目标的信号处理能力，检测前跟踪(TBD-Track Before Detection)的处理方法被提出并得到了广泛的研究。TBD方法的主要思想是：如果能够首先跟踪目标的运动轨迹，并能够在多帧数据的基础上实现对目标能量的积累，那么在一定条件下，这种方法能够实现对弱目标信号能量的有效积累。到目前为止，TBD被认为是解决复杂背景下弱目标检测跟踪问题的重要技术手段。常用的TBD方法有Hough变换、动态规划、最大似然、粒子滤波以及基于离散网格(grid-based)的方法等。但是，Hough变换是沿着直线进行非相参积累，因此只能处理直线运动的目标，而动态规划、最大似然、粒子滤波等方法的运算复杂度都非常高，实时性难以得到保证，目前很难得到实用化。

发明内容

[0004] 针对上述问题，本发明提出了一种目标检测与跟踪一体化方法。包括如下步骤：

[0005] S1：选定满足强度要求的待检测单元范围，初步确定目标数目范围；

[0006] S2：选定所述目标可能存在的速度通道，根据选定的所述速度通道，设计各目标运动模型；

[0007] S3：初始化第1帧数据各目标对应单元的置信度值；

[0008] S4：用k表示当前处理数据的编号，在第k次数据到来时，若k＝1，则表示当前正处理第1次数据，则将步骤S3中的初始化置信度值作为第1次数据各目标对应单元的置信度值；若k≠1，则根据第k-1次各目标对应单元内的置信度值和步骤S2中所设定的运动模型，计算第k次各目标对应单元的置信度值；

[0009] S5：根据步骤S4中得到的第k次数据各目标对应单元的置信度值以及第k次数据的观测值，迭代计算第k次各目标对应单元的置信度值，然后选定阈值对置信度值进行门限处理，得到检测后的结果，确定目标数目和目标可能存在的位置；

[0010] S6：依次类推，在第k+1次数据到来时，重复步骤S4和S5；

[0011] S7：在判决为有目标的单元处用点迹进行标注，若有N个单元，则经过K次处理后，形成一幅K×N的检测后的点迹视图。

[0012] 优选的，将步骤S2中可能存在目标的运动模型建模为一阶马尔科夫模型，并建立目标的置信度值转移模型。

[0013] 优选的，步骤S4和S5中的待检单元的置信度值正比于所述待检单元的原始置信度值与假定目标存在时该单元对应的观测似然函数的乘积。

[0014] 本发明的有益效果在于：

[0015] 一种目标检测跟踪一体化方法，该方法不仅在处理高速机动目标和复杂的杂波背景时具有良好的检测跟踪性能和较好的鲁棒性，而且具有运算复杂度低的优点。

附图说明

[0016] 图1为本发明的算法流程框图；

[0017] 图2为本发明中所要处理的数据的到来和排列方式图；

[0018] 图3为仿真实验中三个目标的真实轨迹图；

[0019] 图4为在瑞利杂波仿真数据下本发明的检测与跟踪结果的点迹视图；

[0020] 图5为低海况实测海杂波下三个目标的检测与跟踪结果的点迹视图；

[0021] 图6为高海况实测海杂波下三个目标的检测与跟踪结果的点迹视图。

具体实施方式

[0022] 下面将结合附图对本发明加以详细说明，具体实施方式如下所述。

[0023] 本发明所述的基于递推贝叶斯原理的一种目标检测与跟踪一体化方法，如图1所示，包括以下步骤：

[0024] 第一步，选定满足强度要求的待检测单元范围，初步设定待检测单元中目标数目范围；

[0025] 考虑常规预警雷达系统，每一帧(即一次数据)有N个距离单元可能存在感兴趣目标，数据的排列与存储方式如图2所示。初步设定所选待检测单元范围内目标的数目范围，以便算法后续的顺利进行。

[0026] 第二步，选定感兴趣目标可能存在的速度通道，根据选定速度通道，设计各目标运动模型；

[0027] 目标运动模型中包括目标的速度、加速度等。

[0028] 第三步，初始化第1帧数据各目标对应单元的置信度值；

[0029] 设第1帧数据各目标对应单元的初始置信度值为Cinit。

[0030] 第四步，用k表示当前处理数据的编号，在第k次数据到来时，若k＝1，则表示当前正处理第1次数据，那么直接将第三步中的初始置信度值作为第1次数据各目标对应单元的置信度值；若k≠1，则根据第k-1次各目标对应单元的置信度值和第二步中所设定的运动模型，计算第k次各目标对应单元的原始置信度值；

[0031] 令符号 表示第k帧时第n个距离单元存在目标， 表示第k帧时第n个距离单元存在目标的原始置信度值。若k＝1，则有 若k≠1， 可以根据上一帧计算得到的 计算得到。

[0032] 第五步，根据第四步中得到的第k次数据各目标对应单元的置信度值以及第k次数据的观测值，迭代计算第k次各目标对应单元的置信度值，然后选定阈值对置信度值进行门限处理得到检测后的结果，确定目标数目和目标可能存在的位置；

[0033] 令 表示第k帧时第n个距离单元存在目标的更新置信度值，为目标存在时的观测似然函数，则根据原始置信度值和观测数据，第k帧各目标对应距离单元的更新置信度值可由式

[0034]

[0035] 来计算(其中，d为常数，可通过归一化等措施进行求解)。

[0036] 通过选取合适的阈值η，可对置信度值进行门限处理得到检测跟踪结果，即若则判为该距离单元有目标，若 则判为该距离单元无目标。

[0037] 第六步，依次类推处理，在第k+1次数据到来时，重复第四步和第五步；

[0038] 第七步，在判决为有目标的单元处用点迹进行标注，若有N个单元，则经过K次处理后，会形成一幅K×N的检测后的点迹视图。

[0039] 下面通过仿真实验进一步说明本发明所述方法的性能。仿真实验中预设了3个目标，两个目标进行直线运动v＝0.3个距离单元/帧，另一个目标作曲线运动，初始速度为v＝0.6个距离单元/帧，所有目标的信杂比均设置为5dB，三个目标的出现和消失的距离单元分别为：目标1起始于距离单元1，消失于距离单元300；目标2起始于距离单元20，消失于距离单元300；目标3起始于距离单元50，消失于距离单元250。仿真实验中三个目标的真实轨迹如图3所示。使用本发明所述方法对实验数据进行处理，选定速度通道v＝0.3个距离单元/帧，在瑞利杂波仿真数据下的检测与跟踪结果如图4所示。在低海况实测海杂波数据下的检测与跟踪结果如图5所示。在高海况实测海杂波数据下的检测与跟踪结果如图6所示。

[0040] 可以看出，本发明所述方法不仅能处理较低信杂比条件下的直线运动目标，还能准确跟踪较低信杂比条件下的曲线运动目标，且在高、低海况实测海杂波数据下也具有较好的检测跟踪效果。

[0041] 此实施例仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。