基于PS-VF-TBD的地波雷达船只目标航迹起始方法
2022-04-29高昆鹏纪永刚王祎鸣孟俊敏
高昆鹏, 纪永刚, 王祎鸣, 孟俊敏
基于PS-VF-TBD的地波雷达船只目标航迹起始方法
高昆鹏1, 3, 纪永刚2, 王祎鸣1, 3, 孟俊敏1, 3
(1. 自然资源部第一海洋研究所, 山东 青岛 266061; 2. 中国石油大学(华东), 山东 青岛 266580; 3. 自然资源部海洋遥测技术创新中心, 山东 青岛 266061)
针对地波雷达低信噪比目标信号“闪烁”导致的航迹起始困难问题, 本文提出了一种基于PS-VF-TBD的航迹起始方法。该方法综合利用了地波雷达回波谱中目标信号的展宽特征和运动特性, 利用PS-VF-TBD方法对运动目标能量有效积累的特点, 提高对疑似目标多帧信号间的关联概率, 进而提高对真实目标的搜索概率和航迹起始成功率。通过地波雷达实测数据的航迹起始结果及长时间跟踪结果与同步AIS信息的对比分析, 表明本文提出的方法相对于基于DP-TBD的航迹起始方法可以给出更好的航迹起始结果, 提高了远端低信噪比下“闪烁”目标的航迹起始概率, 可提高目标跟踪性能, 并具有区分距离-多普勒谱中邻近目标的能力。
高频地波雷达; 航迹起始; PS-VF-TBD; 海上目标跟踪
高频地波雷达(high frequency surface wave radar, HFSWR)是海上舰船目标大范围连续探测的重要手段, 它利用高频电磁波(3～30 MHz)沿海面爬行的特性实现超视距的目标(海上船只、低空飞机等)探测, 并获取运动目标的位置及航速航向信息, 最远探测距离可达300 km, 因此又被称为地波超视距雷达[1-3]。在地波雷达船只目标跟踪过程中, 航迹起始是在目标跟踪过程中未进入稳定跟踪之前的航迹确认过程, 包括真实目标的确认和初始航迹段的确认, 是目标跟踪的首要环节, 因此航迹起始的性能对于海上目标航迹探测的最终结果有重要的影响。
目前常应用于地波雷达的航迹起始方法主要采用先检测后跟踪(detect before track, DBT)的处理流程, 如直观法、逻辑法[4]和Hough变换法等[5-6], 这些方法都是先在距离-多普勒(range-Doppler, RD)谱中进行恒虚警率检测(constant false alarm rate, CFAR)、目标测向和精确参数估计得到目标点迹信息, 实现目标的确认; 然后基于点迹数据开展航迹关联实现初始航迹段的确认[7-8]。这类航迹起始方法可以实现具有较高信噪比、能被持续稳定探测目标的航迹起始, 而对于尺寸较小或者远距离的船只目标, 目标回波强度较弱、信噪比较低, 难以保持持续稳定探测, 影响真实目标的确认和航迹关联过程, 导致航迹起始失败。检测前跟踪(track before detect, TBD)是近年来提出的针对低信噪比目标的检测方法, 其不对单帧回波谱数据进行目标有无的判断, 而是利用在时空序列上目标的关联性和杂波的随机性, 对多帧回波谱数据进行信号能量积累, 最大化积累目标信号能量从而提高目标的信噪比, 在目标检测的同时得到目标航迹。目前常应用于地波雷达的检测前跟踪方法主要有基于动态规划的TBD算法(DP-TBD)[9]和基于粒子滤波的TBD算法(PF- TBD)[10]。其中, PF-TBD算法需要大量的粒子保持鲁棒性, 计算复杂度较高, 实用性较弱, 因此通常较少用于工程应用; DP-TBD算法使用最优化策略在最大允许速度的可行区域内进行搜索, 可能会错误积累低信噪比条件下大于目标信号强度的噪声值, 导致弱目标信号能量积累错误等问题, 容易导致产生虚假航迹。
VF-TBD是近几年发展的一种新的TBD方法, 该方法假设目标在多帧处理期间具有恒定的速度, 当滤波器假设的速度与测量的速度相匹配时, 目标能量沿轨迹积分到最大, 具有提高信噪比和复杂度适中的优点。Zhou等[11]和Wang等[12]发展了基于VF-TBD的伪谱方法(PS-VF-TBD), 将该方法推广到距离-多普勒域, 通过构建点扩展函数的方式构建伪谱, 可以实现目标信号扩展情况下的信号能量精确积累, 仿真验证了在提高快速移动目标信噪比、进行弱目标提取方面具有优势。目前该方法主要用于常规多普勒雷达仿真研究, 还未应用于高频地波雷达。
相对于其他常规多普勒雷达而言, 高频地波雷达具有距离分辨率较低的特性, 而且船只目标具有慢速运动的特点, 因此, 有必要结合地波雷达实测数据特点发展适用于地波雷达的PS-VF-TBD航迹起始方法。本文首先结合实测数据给出了PS-VF- TBD的原理推导, 然后给出了基于PS-VF-TBD的航迹起始数据处理流程, 最后结合实测AIS数据对方法进行验证, 并与DP-TBD航迹起始方法进行了比较。
1 基于PS-VF-TBD的方法原理与处理流程
1.1 PS-VF-TBD原理
PS-VF-TBD依据多普勒雷达目标的展宽特性和运动特性, 通过构建伪谱和速度匹配滤波进行目标回波能量的有效积累, 伪谱构建示意图如图1所示。基于帧目标的距离-多普勒预测值使用如下公式构建伪谱:
其中,(n, k,n, k)为距离-多普勒单元格为(n, k,n, k)处的回波幅度值,、λ为目标的回波扩展参数,SS、SS为在距离维和多普勒维的单元格选取范围,为能量积累批次数。
若目标初始距离-多普勒为(n, 0,n, 0), 那么PS- VF-TBD的传递函数[11]可以表示为:
其中,,分别代表在距离和多普勒方向上从目标能量峰值点到量化单元格之间的距离。目标信号积累总能量可表示为:
其中,、代表能量积累区域在距离和多普勒方向上分别占据几个单元格, 如图2 所示所占单元格为n1—n8。最后将信号积累总能量与设定的阈值进行比较, 判断目标是否存在。
1.2 PS-VF-TBD航迹起始方法应用
PS-VF-TBD对于探测具有匀速直线运动的目标具有极佳的目标提取性能[13], 但是实际中的海上船只目标多为非匀速直线运动的慢速目标。再加上地波雷达具有较高的速度分辨率和较低的距离分辨率, 在实测数据中会存在目标在距离-多普勒谱中的速度抖动、距离变化缓慢的现象, 如果将该方法直接用于地波雷达实测数据, 则无法取得良好的目标提取效果。
图2 多帧能量积累示意图
针对以上问题, 本文采用低门限CFAR检测获取疑似目标点, 基于目标在距离-多普勒中的运动特征利用PS-VF-TBD对展宽目标信号准确积分的特性进行能量积累, 通过疑似目标点与信号积累点相结合的方式将目标的首、末帧特征信息进行关联, 实现真实目标的提取和初始航迹段的确认。一方面, 使用首帧疑似目标点的距离-多普勒坐标作为PS-VF-TBD的目标信号能量积累起始基准点以减少计算量; 另一方面, 通过末帧疑似目标点的距离-多普勒坐标与PS-VF- TBD的信号能量积累点的关联完成真实目标的确认。
其中, Δ、Δ分别为距离分辨率和多普勒速度分辨率,min为最小多普勒速度。对每个疑似目标点进行目标信号积累, 为排除其他目标信号的影响, 以疑似目标点的距离-多普勒坐标为基准限定距离-多普勒范围n_width、n_width基于假设的速度的平方V[11], 计算各单元格在第帧的距离-多普勒预测值:
进而得到对应的距离-多普勒单元格:
图3(a、b)展示了对仿真的疑似展宽目标3在限定区域内的PS-VF-TBD疑似目标信号能量积累。
图3 疑似目标信号积累与真实目标确认
对于远离雷达运动的目标(径向速度为负)在-1帧的搜索区域为:
其中,Z-1(,)代表在第-1帧RD谱中距离-多普勒坐标(,)处的幅度值。在航迹回溯中取A-1中回波幅度最大值对应的距离-多普勒坐标作为在-1中的目标搜索结果, 最终完成初始航迹段的确认。
1.3 航迹起始数据处理流程
基于本文提出的航迹起始方法的实测数据处理流程图如图4所示。
图4 数据处理流程图
步骤1: 利用多帧回波谱数据进行R-D-T三维数据构建。地波雷达数据往往有多个阵列从而形成多通道数据, 为了使航迹起始效果最优, 主要采用对强目标进行信噪比分析的方式进行不同通道的数据质量评估, 判别各通道信号接收质量的优劣, 通过对多通道的回波谱数据对比, 选取目标信噪比高, 背景更加纯净的通道来构建数据集。
步骤2: 数据预处理。该步包括对R-D-T数据进行杂波抑制和低门限CFAR检测。对于地波雷达目标探测而言, 地杂波和一阶海杂波都属于杂波干扰, 其较高的幅度对于步骤3中目标信号的积累会产生不良影响。本文采用文献[14]中基于信噪比的方法提取一阶海杂波并进行抑制, 基于地杂波固定在零频附近的特性进行提取并置为背景值, 对杂波抑制后的数据进行低门限CFAR检测。这样处理有利筛选出疑似目标点, 减少步骤3中PS-VF-TBD处理的计算量, 提高算法运行效率。
步骤3: 基于PS-VF-TBD的无向航迹起始。理论表明积累批次越多越有利于把握目标的运动规律, 疑似目标信号积累点的信噪比越高; 然而, 通过结合AIS目标运动特性分析, 发现积累批次越多目标速度发生变化的可能性越大, 会影响目标提取效果, 因此本文将航迹起始时间限定为10 min。结合回波谱数据中的目标展宽范围将展宽参数、设置为0.7; 综合考虑地波雷达分辨率和船只目标运动特点选取积累区域范围, 在10 min内按船速30节计算最远行驶9.3 km, 在回波谱中展宽目标的变化范围对应大约4～5个距离单元格, 因此积累区域范围设置为5; 综合考虑地波雷达较低的距离分辨率和航迹起始的准确性将关联参数设置为2。
步骤4: 目标航迹空间定位。对无向航迹进行测向、精确径向距离估计、精确径向速度估计, 实现起始航迹的精确参数估计, 然后结合雷达站的空间位置将目标航迹投影到地理坐标系, 完成目标航迹的空间定位。
步骤5: 航速航向估计。由于地波雷达测向误差较大, 在进行航迹空间定位后可能会出现航迹抖动现象, 为保持起始航迹的良好线性特征, 根据目标相对于雷达基站方位缓慢变化的规律性, 进行异常点筛除和方位角变化曲线拟合, 最后进行航速航向估计, 实现对起始航迹运动特征的把握。
2 基于实测数据的方法验证
采用2019年4月在威海海域获取的地波雷达目标探测实验数据, 验证本文提出的基于PS-VF-TBD的地波雷达航迹起始方法的性能。雷达工作中心频率4.7 MHz, 距离分辨率2.5 km, 接收阵元数为5个, 相干积累时间262.144 s, 选取其中10帧回波谱数据开展航迹起始, 相邻帧的时间间隔为1 min, 回波谱数据图5(a)所示。
图5 回波谱数据及AIS信息
通过回波谱数据可以看出船只目标多为面目标, 在多普勒维度上和距离维度上均存在一定的展宽。本文对船舶自动识别系统(automatic identification system, AIS)航迹数据进行了航迹分析, 典型AIS个例两个小时的航迹信息如图5(b、c)所示。可以看出目标在该时段航行期间(绿色部分)航速变化幅度达到0.5节, 在航迹起始阶段(红色部分)航速变化幅度为0.1节, 航迹起始阶段可以近似看作目标做匀速直线运动。
为比较本文提出方法的航迹起始性能, 与常规检测前跟踪方法DP-TBD算法进行了航迹起始结果对比并给出了长时间(60 min)跟踪结果。图6分别展现了目标在雷达坐标系下的航迹起始结果和长时间跟踪结果。通过图6可以看出, 两种方法起始的航迹数量较为接近, 航迹起始最远距离可以达到150 km左右, 在60～120 km距离处起始了相对较多的航迹。结合长时间跟踪结果可以看出从整体上看航迹结果大致相当, 由于两种航迹起始方法都是基于检测前跟踪思想, 因此从整体上无法直观看出方法的优劣及航迹起始质量。为了验证本文所提出方法的航迹起始效果, 采用同步AIS数据进行对比分析, 图7展示了具有同步AIS信息的部分长航迹目标跟踪个例结果。
本文选取典型个例进行具体分析, 图8展示了对个例1的分析结果:
通过目标航迹与AIS信息的对比可以看出, 本文方法实现了对目标56航迹的正确起始, 对该目标的最远探测距离可以达到114.5 km, 对地波雷达目标探测而言属于远端目标, 且起始航迹的径向速度误差小于0.5 km/h, 径向距离误差小于1.5 km, 均在地波雷达系统目标探测允许的误差范围之内; 而DP-TBD方法在航迹起始阶段的第6～8 min没有正确搜索到航迹点, 出现了航迹中断的问题, 导致航迹起始失败。下面结合目标的幅度信息和信噪比进行具体分析。
图6 不同方法下的结果对比
图 7 地波雷达航迹与AIS航迹结果对比
图8 个例1运动特性分析
通过图9可以看出, 在航迹起始阶段目标56的幅度强度较低, 整体上在189～194 dB范围内波动, 信噪比最低为3.5 dB最高为11.5 dB, 总体上在6 dB左右波动, 目标信号“闪烁”。该目标在前6时刻幅度相对较高, 在第6～8时刻幅度迅速降低, 由于DP-TBD算法是通过设定搜索范围通过搜索回波信号强度最大的位置确定目标, 另外该时段内目标信噪比较低(3～6 dB), 导致在连续批次的搜索中找不到目标点, 最终导致航迹起始失败; 而本文方法是依据目标的运动特性对航迹起始时段内的目标信号积累作为判断疑似目标点的依据, 虽然目标56在6～8时刻信号变弱, 但是多帧信号的有效积累使得目标信号突显, 如图9(c)所示, 实现了目标多帧信息的关联, 找到了起始航迹的始末帧点位并最终成功起始航迹。图10展示了该目标的航迹起始及长时间跟踪的V-T、R-T分布图。
图9 个例1幅度特性分析
图10 个例1 V-T、R-T分布
图11展示了对个例2的分析结果:
通过图11可以看出, 对于42号目标(AIS编号1881)和72号目标(AIS编号1890), 本文方法分别实现了两条航迹的起始, 并且通过长时间跟踪结果及AIS信息证明了航迹起始结果的正确性。基于DP-TBD的航迹起始方法则只完成其中一条航迹的起始, 并在航迹起始阶段存在航迹关联错误的现象, 在第1、2个时刻提取出了对应AIS编号1881的目标的径向距离、径向速度信息, 在第6～10个时刻提取出了对应AIS编号1890的目标的径向距离、径向速度信息。下面给出了回波谱数据及AIS投影结果进一步分析。
图11 个例2 R-V变化分析
图12 结合回波谱数据的个例2分析
通过回波谱数据及AIS投影结果可以看出, 这两个目标在距离-多普勒谱中为邻近目标, 相差2～3个多普勒单元格, 有相似的运动特征, 而且这两个目标在回波谱数据中能量较为分散, 目标中心点能量不突出, 在航迹起始阶段信噪比都是低于10 dB。通过这两个目标的幅度变化曲线可以看出, 目标42的幅度变化较为稳定, 幅度值在193 dB上下浮动; 而目标72的幅度变化较大, 总体上呈上升趋势, 幅度最低为180 dB左右, 幅度最高可达202 dB。在前5个时刻目标42的幅度值大于目标72, 在后5个时刻目标42的幅度值小于目标, 从而导致了基于DP-TBD起始的航迹在第3～6时刻航迹关联错误, 导致航迹起始失败, 最终影响了航迹起始质量; 而本文方法利用航迹起始阶段目标近似匀速直线运动的特性, 在设定的距离-多普勒范围内进行目标信号积累, 降低了目标信号积累错误的可能性, 从而实现了邻近弱目标的提取及区分。图13、图14通过V-T、R-T分布展示了上述两个目标起始航迹和长时间跟踪航迹结果。
图13 航迹起始阶段目标V-T、R-T分布图
图14 长时段目标航迹V-T、R-T分布图
图15 个例2在地理坐标系下的运动分析
图15展示了基于两种航迹起始方法提取的目标个例2在地理坐标系下的长时间跟踪结果及对应的AIS信息。
通过图15可以看出, 42号目标(AIS编号1881)和72号目标(AIS编号1890)在地理坐标系下也为邻近目标。通过目标的径向速度、径向距离、角度变化曲线可以看出, 这两个目标相对于雷达的径向距离、径向速度、方位角参数接近, 具有相似的动态信息特征, 使得在航迹起始阶段有相似的运动特征。由于地波雷达距离分辨率和方位分辨率较低, 再加上这两个目标在回波谱中呈现弱目标特性, 使得对此类目标的提取及区分存在一定的困难。对此类情况目标的区分进一步展现了本文方法对邻近目标具有一定的区分能力。
综上所述, 通过本文提出的航迹起始方法与基于DP-TBD的航迹起始方法对比, 证明了本文方法航迹起始的优越性能。实验结果表明, 相对于基于DP-TBD的航迹起始方法, 本文提出的航迹起始方法可以提高远端航迹的起始成功率, 并对于在距离-多普勒谱中具有相似运动特征的邻近目标有一定的区分能力。
3 结论
通过对高频地波雷达回波谱数据进行目标展宽特性分析, 结合AIS数据进行目标运动特性分析, 将PS-VF-TBD方法应用于地波雷达航迹起始, 通过对疑似目标首帧与末帧信息的关联实现虚假目标的剔除和真实航迹的确认, 最后与基于DP-TBD的航迹起始方法进行对比, 结果表明: 本文方法可实现较远端、非稳定观测的低信噪比目标的航迹起始, 进而提高航迹跟踪结果; 对于距离-多普勒谱中的邻近目标具有一定的区分能力。
[1] JI Y G, ZHANG J, WANG Y M, et al. Target monitoring using small-aperture compact high-frequency surface wave radar[J]. IEEE Aerospace & Electronic Systems Magazine, 2018, 33(3): 21-31.
[2] YANG Q S, JI Y G, LIU Y X, et al. Target tracking of surface wave over-the-horizon radars based on FJPDA algorithm[J]. Advances in Marine Science, 2019, 37(1): 8.
[3] 王瑞富, 孔祥超, 纪永刚, 等. 基于GIS的高频地波雷达回波信号的表达处理与显示[J]. 海洋科学, 2015, 39(7): 86-92.
WANG Ruifu, KONG Xiangchao, JI Yonggang, et al. Expression and display of echo signal of HF surface wave radar based on GIS[J]. Marine Sciences, 2015, 39(7): 86-92.
[4] 王健. 高频地波雷达舰船目标跟踪关键技术研究[D]. 哈尔滨: 哈尔滨工业大学, 2008.
WANG Jian. Research on Key Technologies of ship target tracking with high frequency ground wave radar[D]. Harbin: Harbin Institute of Technology, 2008.
[5] 王国宏, 苏峰, 毛士艺, 等. 杂波环境下基于Hough变换和逻辑的快速航迹起始算法[J]. 系统仿真学报, 2002, 14(7): 874-876.
WANG Guohong, SU Feng, MAO Shiyi, et al. Fast track initiation algorithm based on Hough transform and logic in clutter environment[J]. Journal of System Simulation, 2002, 14(7): 874-876.
[6] 金术玲, 梁彦, 潘泉, 等. 基于Hough变换和聚类的航迹起始算法[J]. 系统仿真学报, 2009, 21(8): 2362- 2364.
JIN Shuling, LIANG Yan, PAN Quan, et al. Track initiation algorithm based on Hough transform and clustering[J]. Journal of System Simulation, 2009, 21(8): 2362- 2364.
[7] PONSFORD A M, WANG J. A review of high frequency surface wave radar for detection and tracking of ships[J]. Turkish Journal of Electrical Engineering & Computer Sicences, 2010, 18(3): 409-428.
[8] VIVONE G, Braca P, Horstmann J. Knowledge- based multitarget ship tracking for HF surface wave radar systems[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(7): 3931-3949.
[9] YI W, MORELANDE M R, KONG L, et al. An efficient multi-frame track-before-detect algorithm for multi-target tracking[J]. IEEE Journal of Selected Topics in Signal Processing, 2013, 7(3): 421-434.
[10] 杨亚波, 夏永红, 匡华星, 等. 雷达微弱目标检测前跟踪技术研究[J]. 雷达与对抗, 2015, 35(2): 22-28.
YANG Yabo, XIA Yonghong, KUANG Huaxing, et al. Research on tracking technology before radar weak target detection[J]. RADAR & ECM, 2015, 35(2): 22-28.
[11] ZHOU Gongjian, WANG Liangliang. Pseudo-spectrum based speed square filter for track-before-detect in range- doppler domain[J]. IEEE Transactions on Signal Processing, 2019, 67(21): 5596-5610.
[12] WANG Liangliang, ZHOU Gongjian, LI Peiyuan. A complex pseudo-spectrum based velocity filtering method for track- before-detect[J]. Signal Processing, 2020, 174: 107651.
[13] ZHOU G J, WANG L L, Kirubarajan T. A pseudo- spectrum approach for weak target detection and tracking[J]. IEEE Transactions on Aerospace and Electronic Systems, 2019, 55(6): 3394-3412.
[14] 纪永刚, 张杰, 王彩玲, 等. 基于信噪比方法的天地波混合体制雷达一阶回波谱提取[J]. 电子与信息学报, 2015, 37(9): 2177-2182.
JI Yonggang, ZHANG Jie, WANG Cailing, et al. Extraction of first-order echo spectrum of heaven earth wave hybrid radar based on signal-to-noise ratio method[J]. Journal of Electronics and Information, 2015, 37(9): 2177-2182.
Ship track initiation method of a high-frequency surface wave radar based on PS-VF-TBD
GAO Kun-peng1, 3, JI Yong-gang2, WANG Yi-ming1, 3, MENG Jun-min1, 3
(1. First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China; 2. China University of Petroleum (East China), Qingdao 266580, China; 3. Technology Innovation Center for Ocean Telemetry, Ministry of Natural Resources, Qingdao 266061, China)
The high-frequency surface wave radar has a low signal-to-noise ratio, making the target signal “flicker” and tracking initiation difficult. To address this problem, this paper proposes a track initiation method based on PS-VF-TBD. This method makes comprehensive use of the broadening and motion characteristics of the target signal in the echo spectrum, as well as the PS-VF-TBD method that can effectively accumulate the energy of the moving target, to improve the correlation probability between multi-frame signals of suspected targets, thereby improving the search probability of the real target and the success rate of track initiation. The track initiation and long-time tracking performance based on the measured data and synchronous AIS information of the ground wave radar were compared and analyzed; the result reveals that the proposed method could give better track initiation results than the method based on DP-TBD, improve the track initiation probability of “flickering” targets with low signal-to-noise ratio at the far, and improve the target tracking performance. Moreover, it can distinguish adjacent targets in the range-Doppler spectrum.
high-frequency surface wave radar; track initiation; PS-VF-TBD; marine target tracking
Oct. 13, 2021
TN958.93
A
1000-3096(2022)04-0013-12
10.11759/hykx20211013002
2021-10-13;
2021-11-16
国家重点研发计划项目课题(2017YFC1405202); 国家自然科学基金(62031015)
[National Key R & D Program of China, No. 2017YFC1405202; the National Natural Science Foundation of China, No. 62031015]
高昆鹏(1996—), 男, 山东济宁人, 硕士, 主要从事高频地波雷达目标检测与跟踪方面研究, E-mail: gaokunpeng@fio.org.cn; 纪永刚(1977—),通信作者, 教授, 主要从事高频地波雷达海洋探测方面研究, E-mail: jiyonggang@upc.edu.cn
(本文编辑: 赵卫红)
