本发明属于主动声纳水声目标探测和定位技术领域,具体地说,涉及一种主动声纳目标运动信息融合方法。
背景技术:
海洋声学目标探测技术对于维护国家主权,保障国家海洋环境安全,促进海洋探索与开发至关重要。特别是在民用领域,可用于海洋救险、鱼群探测和海上安全作业等;在军事和国防应用领域,可广泛应于军事上的各种防卫和攻击系统,例如海岸监视系统和反水雷反鱼雷系统等。
目前,水声目标一直向隐身化、长航时、大潜深、高攻击方向发展,给被动声纳探测技术带来了巨大挑战,同时促使主动声纳探测技术得到了足够的重视。当前,主动声纳探测技术的性能提升多依赖于低频大功率声源的研制和接收阵孔径的增大,其发展得益于电子元器件水平和海洋工程化能力的提高,而缺少理论上的突破,在融合目标回波信号多维度特性获取额外增益方面的研究较少,开展多维信息融合提高主动声纳目标定位精度和降低目标虚警率具有较高的实用价值。
现有的目标运动信息融合方法主要包括卡尔曼滤波、无迹卡尔曼滤波和粒子滤波,以及在现有的目标运动信息融合方法基础上的一些改进算法,现有的方法没有利用潜在目标的运动速度在各方位上投影分量对应的多普勒系数信息,存在目标跟踪速度慢、实时性差的问题,需要长时间多帧探测结果才能实现目标轨迹跟踪,而且虚警率高,常出现将混响干扰误判别为目标的问题。
技术实现要素:
为解决现有技术存在上述缺陷,本发明提出了一种主动声纳目标运动信息融合方法,在脉冲波主动声纳系统中,每发射一次脉冲信号,获得一帧探测结果,利用水声目标运动状态具有连续性和虚警目标空间位置分布具有随机性的特点,对近邻多帧探测结果进行数据融合,以提高主动声纳系统的水声目标的定位精度,并降低目标虚警率。
本发明提供了一种主动声纳目标运动信息融合方法,该方法包括:
步骤1)针对探测结果构成集合a={a1a2…ai…an},获得第i帧探测结果ai中的mi个潜在目标中的第j个潜在目标从时间ti到tk的、第i帧到第k帧的近邻帧之间的目标候选运动路径的路径数目集合pi,j,k′;
步骤2)对获得的近邻帧之间的目标候选运动路径的路径数目集合pi,j,k′,进行一次虚假路径剔除和目标运动信息融合,获得剩余候选运动路径的路径数目集合pi,j,k*和剩余候选运动路径的平均速度集合
步骤3)基于获得的剩余候选运动路径的路径数目集合
步骤4)对剩余运动路径的路径数目集合p1,进行二次虚假路径剔除和目标运动信息融合,获得第j个潜在目标的最终的运动路径和对应的运动轨迹信息,将其作为ai的目标运动信息;
步骤5)重复步骤1)-4),得到探测结果构成集合a中的每一帧探测结果的目标运动信息,并将每一帧探测结果的目标运动信息叠加在一起,形成主动声纳目标的运动路径。
作为上述技术方案的改进之一,所述步骤1)具体包括:
假设某一站点共连续发射n次脉冲信号,该站点连续发射的n次脉冲信号对应的发射时间分别为t={t1t2…ti…tn};其对应的各帧探测结果构成集合a={a1a2…ai…an};
其中,ai为第i帧探测结果;ai是一个mi行4列特征向量矩阵:
其中,mi为ai中包含的潜在目标点数;xi,j、yi,j、θi,j和
假设区域si,j,k为在ai中以第j个潜在目标的位置为中心,以ri,j,k为半径形成的区域,假设在区域si,j,k中,第j个潜在目标以第i帧为初始位置,查找ai中第j个潜在目标在第k帧中可能的位置;
其中,ri,j,k=(tk-ti)×vmax+dmax;
其中,ti为第i帧时刻;tk为第k帧时刻;tk大于ti;k>i,k和i为近邻帧;vmax为第j个潜在目标的最大运动速度;dmax为误差值;
根据上述查找结果和预先划定的区域si,j,k,获得在第k帧探测结果ak中,落在区域si,j,k中的可能的第j个潜在目标的目标数为pi,j,k,进而获得第j个潜在目标从时间ti到tk、第i帧到第k帧的候选运动路径的路径数目集合为pi,j,k′;集合pi,j,k′中包含的元素数目为pi,j,k。
作为上述技术方案的改进之一,所述步骤2)具体包括:
在ai中第j个潜在目标从时间ti到tk的、从第i帧运动到第k帧中的任一连续两帧探测结果形成的候选运动路径进行目标运动信息融合,得到ai和ak的平均运动速度矢量vi,j,k,h;
获得在候选运动路径的路径数目集合pi,j,k′中的各候选运动路径对应的平均运动速度向量构成平均速度集合vi,j,k:
vi,j,k=[vi,j,k,1vi,j,k,2…vi,j,k,h…vi,j,k,pi,j,k],
其中,ai和ak的平均运动速度矢量vi,j,k,h的大小为|vi,j,k,h|;
ai中第j个潜在目标的平均运动速度矢量vi,j,k,h的方向与方位角θi,j和θk,h的夹角分别为αi,j,k,h和βi,j,k,h;则ai中第j个潜在目标的平均运动速度矢量vi,j,k,h在方位角θi,j和θk,h上投影分量对应的多普勒系数分别为
其中,c为水中声波传播速度;
假设
其中,
如果εi,j,k,h≥εmax;则vi,j,k,h对应的候选运动路径为虚假路径,将其从候选运动路径的路径数目集合pi,j,k′中剔除;
对ai中第j个潜在目标从时间ti到tk的、从第i帧运动到第k帧中的剩余的每一连续两帧探测结果,重复上述过程,进行一次虚假路径剔除,将在候选运动路径的路径数目集合pi,j,k′中剔除虚假路径后,获得剩余候选运动路径的路径数目集合
同时,在剩余候选运动路径的路径数目集合
其中,
作为上述技术方案的改进之一,所述步骤3)具体包括:
根据构成的剩余运动路径的平均速度集合
将
其中,
该第j个潜在目标沿着
其中,l>k;k和l为近邻帧;
假设剔除区域si,j,k,h,l为在al中以坐标
根据上述查找结果和预先划定的剔除区域si,j,k,h,l,获得在第l帧探测结果al中,落在剔除区域si,j,k,h,l中的可能的第j个潜在目标的目标数为
同时,在ai中第j个潜在目标从时间tk到tl的、从第k帧运动到第l帧中的任一连续三帧探测结果形成的剩余运动路径进行目标运动信息融合,
得到第j个潜在目标从第k帧运动到第l帧的、ak和al的剩余运动路径的平均运动速度矢量
进而构成剩余运动路径的平均速度集合
其中,
作为上述技术方案的改进之一,所述步骤4)具体包括:
根据构成的剩余运动路径的平均速度集合
将
其中,剩余运动路径的平均运动速度矢量
和
假设
如果
对ai中第j个潜在目标从时间tk到tl的、从第k帧运动到第l帧中的剩余的每一连续三帧探测结果,重复上述过程,进行二次虚假路径剔除,将在剩余运动路径的路径数目集合p1中剔除虚假路径后,
得到第j个潜在目标的最终的运动路径,则该运动路径为当
进而获得第j个潜在目标的最终的运动路径的运动轨迹信息;
其中,第j个潜在目标的最终的运动路径为当
所述最终的运动路径的运动轨迹信息包括:在时刻ti、tk和tl对应的ai、ak和al的坐标依次为(xi,yi)、
本发明与现有技术相比的有益效果是:
本发明的方法充分利用水声目标运动状态具有连续性,虚警目标空间位置分布具有随机性,提出目标运动信息融合方法,提高了主动声纳目标定位精度,同时,降低了目标虚警率。
附图说明
图1是本发明的一种主动声纳目标运动信息融合方法的流程图;
图2是本发明的一种主动声纳目标运动信息融合方法的一个具体实施例的合作目标的运动轨迹图;
图3是未采用本发明的一种主动声纳目标运动信息融合方法的虚警目标和合作目标的运动轨迹示意图。
具体实施方式
现结合附图对本发明作进一步的描述。
在本发明提供的方法中,对如下表述进行了具体限定:
近邻帧是指从第i帧到第k帧,且k>i,第i帧到第k帧之间存在多帧。
近邻帧之间目标候选运动路径是指第j个潜在目标从时间ti到tk的、从第i帧到第k帧、记录第j个潜在目标的运动信息的运动路径;
其中,所述潜在目标的运动信息包括:第i帧探测结果ai中第j个潜在目标的x轴坐标xi,j、y轴坐标yi,j、波束形成对应的方位角θi,j、第j个潜在目标的运动速度在θi,j方位上投影分量对应的多普勒系数
候选运动路径是指潜在目标的运动路径;
以第i帧探测结果ai中的mi个潜在目标中的每一个潜在目标的位置为圆心,以ri,j,k为半径形成区域si,j,k,在该区域内获得多条针对该潜在目标的候选运动路径,再利用本发明的方法进行一次虚假路径剔除和二次虚假路径剔除,最后得到一个针对该潜在目标的运动路径,并获取该潜在目标的运动轨迹信息;对每一帧探测结果均重复上述过程,获得对应的潜在目标的运动轨迹,及其对应的潜在目标的运动轨迹信息。
连续两帧潜在目标的运动信息融合是指针对第j个潜在目标,从时间ti到tk内、从第i帧到第k帧中任选连续两帧,获取该连续两帧之间形成的运动路径对应的平均运动速度向量,根据预先设定的多普勒判别阈值,判断该运动路径是否为虚假路径,进行一次虚假路径剔除,对非虚假路径的运动路径进行目标运动信息融合。
每连续多帧进行潜在目标的运动信息融合是指针对第j个潜在目标,从时间ti到tk内、从第i帧到第k帧中任选至少连续三帧,获取该至少连续三帧之间形成的运动路径对应的平均运动速度向量,根据预先设定的多普勒判别阈值,判断该运动路径是否为虚假路径,对非虚假路径的运动路径进行目标运动信息融合。
实施例1.
如图1所示,本发明提供了一种主动声纳目标运动信息融合方法,该方法包括:
步骤1)针对探测结果构成集合a={a1a2…ai…an},获得第i帧探测结果ai中的mi个潜在目标中的第j个潜在目标从时间ti到tk的、第i帧到第k帧的近邻帧之间的目标候选运动路径的路径数目集合pi,j,k′;
具体地,假设某一站点共连续发射n次脉冲信号,该站点连续发射的n次脉冲信号对应的发射时间分别为t={t1t2…ti…tn};其对应的各帧探测结果构成集合a={a1a2…ai…an};
其中,ai为第i帧探测结果;ai是一个mi行4列特征向量矩阵:
其中,mi为ai中包含的潜在目标点数;xi,j、yi,j、θi,j和
假设区域si,j,k为在ai中以第j个潜在目标的位置为中心,以ri,j,k为半径形成的区域,假设在区域si,j,k中,第j个潜在目标以第i帧为初始位置,查找ai中第j个潜在目标在第k帧中可能的位置,即第j个潜在目标从第i帧运动到第k帧的可能位置;
其中,ri,j,k=(tk-ti)×vmax+dmax;
其中,ti为第i帧时刻;tk为第k帧时刻;tk大于ti;k>i,k和i为近邻帧;vmax为第j个潜在目标的最大运动速度;dmax为误差值;其中,该误差值为声速估计误差、测时误差和测向误差导致的定位误差最大值;
根据上述查找结果和预先划定的区域si,j,k,获得在第k帧探测结果ak中,落在区域si,j,k中的可能的第j个潜在目标的目标数为pi,j,k,进而获得第j个潜在目标从时间ti到tk、第i帧到第k帧的候选运动路径的路径数目集合为pi,j,k′;集合pi,j,k′中包含的元素数目为pi,j,k。
其中,pi,j,k′中包含多条候选运动路径;上述候选运动路径为ai中第j个潜在目标从ai运动到ak可能的第j个潜在目标的运动路径。
步骤2)对获得的近邻帧之间的目标候选运动路径的路径数目集合pi,j,k′,进行一次虚假路径剔除和目标运动信息融合,获得剩余候选运动路径的路径数目集合pi,j,k*和剩余候选运动路径的平均速度集合
具体地,在ai中第j个潜在目标从时间ti到tk的、从第i帧运动到第k帧中的任一连续两帧探测结果形成的候选运动路径进行目标运动信息融合,得到ai和ak的平均运动速度矢量vi,j,k,h;
获得在候选运动路径的路径数目集合pi,j,k′中的各候选运动路径对应的平均运动速度向量构成平均速度集合vi,j,k:
vi,j,k=[vi,j,k,1vi,j,k,2…vi,j,k,h…vi,j,k,pi,j,k],
其中,ai和ak的平均运动速度矢量vi,j,k,h的大小为|vi,j,k,h|;
ai中第j个潜在目标的平均运动速度矢量vi,j,k,h的方向与方位角θi,j和θk,h的夹角分别为αi,j,k,h和βi,j,k,h;则ai中第j个潜在目标的平均运动速度矢量vi,j,k,h在方位角θi,j和θk,h上投影分量对应的多普勒系数分别为
其中,c为水中声波传播速度;
假设
其中,
如果εi,j,k,h≥εmax;则vi,j,k,h对应的候选运动路径为虚假路径,将其从候选运动路径的路径数目集合pi,j,k′中剔除;
对ai中第j个潜在目标从时间ti到tk的、从第i帧运动到第k帧中的剩余的每一连续两帧探测结果,重复上述过程,讲行一次虚假路径剔除,将在候选运动路径的路径数目集合pi,j,k′中剔除虚假路径后,获得剩余候选运动路径的路径数目集合
同时,在剩余候选运动路径的路径数目集合
其中,
所述目标信息融合是指ai中mi个潜在目标中的每一个潜在目标均进行一次虚假路径剔除,得到满足判决条件:εi,j,k,h<εmax的每一个潜在目标对应的x轴坐标、y轴坐标、波束形成对应的方位角、该潜在目标的运动速度在方位角方位上投影分量对应的多普勒系数的运动信息,形成一个新的
步骤3)基于获得的剩余候选运动路径的路径数目集合
具体地,根据构成的剩余运动路径的平均速度集合
将
其中,
该第j个潜在目标沿着
其中,l>k;k和l为近邻帧;
假设剔除区域si,j,k,h,l为在al中以坐标
根据上述查找结果和预先划定的剔除区域si,j,k,h,l,获得在第l帧探测结果al中,落在剔除区域si,j,k,h,l中的可能的第j个潜在目标的目标数为
同时,在ai中第j个潜在目标从时间tk到tl的、从第k帧运动到第l帧中的任一连续三帧探测结果形成的剩余运动路径进行目标运动信息融合,
得到第j个潜在目标从第k帧运动到第l帧的、ak和al的剩余运动路径的平均运动速度矢量
进而构成剩余运动路径的平均速度集合
其中,
步骤4)对剩余运动路径的路径数目集合p1,进行二次虚假路径剔除和目标运动信息融合,获得第j个潜在目标的最终的运动路径和对应的运动轨迹信息,将其作为ai的目标运动信息;
具体地,根据构成的剩余运动路径的平均速度集合
将
其中,剩余运动路径的平均运动速度矢量
和
假设
如果
对ai中第j个潜在目标从时间tk到tl的、从第k帧运动到第l帧中的剩余的每一连续三帧探测结果,重复上述过程,进行二次虚假路径剔除,将在剩余运动路径的路径数目集合p1中剔除虚假路径后,
得到第j个潜在目标的最终的运动路径,则该运动路径为当
进而获得第j个潜在目标的最终的运动路径的运动轨迹信息;
其中,第j个潜在目标的最终的运动路径为当
所述最终的运动路径的运动轨迹信息包括:在时刻ti、tk和tl对应的ai、ak和al的坐标依次为(xi,yi)、
步骤5)重复步骤1)-4),得到探测结果构成集合a中的每一帧探测结果的目标运动信息,并将每一帧探测结果的目标运动信息叠加在一起,形成主动声纳目标的运动路径,为实现低误判率主动声纳目标探测系统做准备。
实施例2.
采用2018年4月在某一海域开展的一次浅海双基地主动声纳探测试验进行说明。海试中,采用吊放换能器作为声源,在某站点共发射8次脉冲信号,各脉冲信号相对于第1次脉冲信号发射时刻的时间t分别为0秒、120秒、341秒、541秒、891秒、1231秒、1731秒和2301秒;采用“实验1号”双体船的水下潜体作为合作目标;接收信号的基阵位于海底。经过预处理,各帧探测结果构成集合为:
a={a1a2…ai…a8};
在各帧探测结果中包含的潜在目标数分别为51、40、69、81、294、271、84、103;
假设潜在目标最大运动速度为vmax=8节;由于声速估计误差、测时误差和测向误差等导致的定位误差最大值为dmax=150米;
本发明的具体方法如下:
步骤1)针对探测结果构成集合a={a1a2…ai…a8},获得第1帧探测结果a1中的51个潜在目标中的第j个潜在目标从时间t1到t2的、第1帧到第2帧的近邻帧之间的目标候选运动路径的路径数目集合;
具体地,假设区域s1,j,2为在a1中以第j个潜在目标的位置为中心,以r1,j,2=120×8+150=1110米为半径形成的区域,假设在区域s1,j,k中,第j个潜在目标以第1帧为初始位置,查找a1中第j个潜在目标在第2帧中可能的位置;
根据上述查找结果和预先划定的区域s1,j,2,获得在第2帧探测结果a2中,落在区域s1,j,2中的可能的第j个潜在目标的目标数为p1,j,2,其中,j取值为37、39、40、41时,p1,j,2>0;j取其他值时,p1,j,2=0。其中p1,37,2=1、p1,39,2=2、p1,40,2=2、p1,41,2=2。进而获得第j个潜在目标从时间t1到t2的、第1帧到第2帧的候选运动路径的路径数目集合为p1,j,2′,p1,j,2′中包含7个元素,即有7条符合要求的运动路径。
步骤2)从获得的近邻帧之间的目标候选运动路径的路径数目集合中,进行一次虚假路径剔除和目标运动信息融合,获得剩余候选运动路径的路径数目集合p1,j,2*和剩余运动路径的平均速度集合
具体地,多普勒判别阈值εmax=2e-4;如果ε1,j,2,h≥εmax;则v1,j,2,h对应的候选运动路径为虚假路径,将其从候选运动路径的路径数目集合p1,j,2′中剔除;
对a1中第j个潜在目标从时间t1到t2的、从第1帧运动到第2帧中的剩余的每一连续两帧探测结果,重复上述过程,进行一次虚假路径剔除,将在候选运动路径数目集合p1,j,2′中剔除虚假路径后,获得剩余候选运动路径的路径数目集合
同时,在a1中第j个潜在目标从时间t1到t2的剩余候选运动路径的路径数目集合
根据构成的剩余运动路径的平均速度集合
步骤3)基于获得的剩余候选运动路径的路径数目集合
该第j个潜在目标沿着
假设区域s1,j,2,1,3为在a3中以坐标
根据上述查找结果和预先划定的区域s1,j,2,1,3,获得在第3帧探测结果a3中,落在区域s1,j,2,1,3中的可能的第j个潜在目标的目标数为
则第j个潜在目标从时间t2到t3、第2帧到第3帧的各剩余运动路径对应的平均速度集合
步骤4)对剩余运动路径的路径数目集合p1,进行二次虚假路径剔除和目标运动信息融合,获得第j个潜在目标的最终的运动路径和对应的运动轨迹信息,将其作为ai的目标运动信息;
具体地,根据构成的剩余路径的平均速度集合
第j个潜在目标的最终的运动路径为在时刻t1、t2和t3对应的a1、a2和a3的坐标依次顺序连接(-138.09米,-1582.3米)、(-276.00米,-1581.2米)和(-547.65米,-1568.3米)形成的运动路径;
第j个潜在目标的最终的运动路径的运动轨迹信息包括:在时刻t1、t2和t3对应的a1、a2和a3的坐标依次为(-138.09米,-1582.3米)、(-276.00米,-1581.2米)和(-547.65米,-1568.3米),从时刻t1到t2的平均运动速度矢量为
将第j个潜在目标的最终的运动路径及其对应的运动轨迹信息作为第1帧探测结果a1的目标运动信息;
步骤5)重复步骤1)-4),得到该站点连续发射8次脉冲信号的0-2301s时间段内,探测结果构成集合a中剩下的a2、a3、a4、a5、a6、a7、a8中每一帧探测结果的目标运动信息,其中,由于a7、a8中的连续帧数少于2帧,不满足本发明的处理条件。因此,本发明的方法不对a7、a8探测结果进行处理,探测结果构成集合a中剩下的a2、a3、a4、a5、a6处理结果如下:
a2的目标运动信息包括:在时刻t2、t3和t4对应的a2、a3和a4的坐标依次为(-138.08米,-1582.3米)、(-276.00米,-1581.2米)和(-547.51米,-1567.9米),从时刻t2到t3的平均运动速度矢量为
a3的目标运动信息包括:在时刻t3、t4和t5对应的a3、a4和a5的坐标依次为(-547.65米,-1568.3米)、(-272.71米,-1562.3米)和(136.17米,-1560.4米),从时刻t3到t4的平均运动速度矢量为
a4的目标运动信息包括:目标1)在时刻t4、t5和t6对应的a4、a5和a6的坐标依次为(-136.37米,-1562.6米)、(-544.89米,-1560.4米)和(-812.16米,-1549.7米),从时刻t4到t5的平均运动速度矢量为
a5的目标运动信息包括:目标1)在时刻t5、t6和t7对应的a5、a6和a7的坐标依次为(-544.89米,-1560.4米)、(-812.16米,-1549.7米)和(-1214.9米,-1543.7米),从时刻t5到t6的平均运动速度矢量为
a6的目标运动信息包括:目标1)在时刻t6、t7和t8对应的a6、a7和a8的坐标依次为(-677.01米,-1550.6米)、(-1214.9米,-1543.7米)和(-1750.0米,-1536.0米),从时刻t6到t7的平均运动速度矢量为
经过上述结果,可以确定,目标1)、目标2)、目标3)、目标4)各自对应的运动路径为每个潜在目标的最终的运动路径,
如图2所示,将上述每一帧探测结果的4个潜在目标的最终的运动路径进行叠加,形成图2所示的合作目标运行轨迹,完成主动声纳目标运动信息融合,为后续处理做好准备。其中,图2所示的合作目标运行轨迹是将每一帧探测结果的运动路径进行叠加的结果。
图2为采用本发明方法后,在某站点进行探测的结果图,图2以接收阵所在位置,即接收阵位置为坐标原点(0,0),横坐标表示该站点探测区域的x轴坐标,纵坐标表示该站点探测区域的y轴坐标,发射船位置坐标为(1071.9米,-4257.9米),图2中的合作目标运行轨迹每一帧探测结果ai的最终的运动路径叠加后形成的运动路径,图2中发现1个合作目标(实验1号)和1个非合作目标,其中,所发现的合作目标的真实运动轨迹采用gps进行记录,本发明方法所发现的合作目标的运动轨迹,即定位结果,与gps记录的真实运动轨迹相比,平均定位误差为135米,误差率为0.87%;1个非合作目标为虚警目标,非合作目标运动轨迹为经过本发明的方法剔除后剩余的虚假运动路径。
图3为该站点未采用本发明的方法进行运动信息融合的探测结果图,图3以接收阵所在位置为坐标原点(0,0),横坐标表示该站点探测区域的x轴坐标,纵坐标表示该站点探测区域的y轴坐标,发射船位置坐标为(1071.9米,-4257.9米),图3中包含大量虚警目标,不能有效的获取合作目标运动轨迹信息。
图2与图3相比,通过对近邻多帧探测结果进行数据融合,将潜在目标数目从993个降低到13个,其中保留下来的12个潜在目标是合作目标(实验1号)及其多途产生的,虚假目标得到较大程度的剔除。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。