CN114488112B - 一种基于不同体制雷达的波形与点迹融合检测方法 - Google Patents
一种基于不同体制雷达的波形与点迹融合检测方法 Download PDFInfo
- Publication number
- CN114488112B CN114488112B CN202111592160.6A CN202111592160A CN114488112B CN 114488112 B CN114488112 B CN 114488112B CN 202111592160 A CN202111592160 A CN 202111592160A CN 114488112 B CN114488112 B CN 114488112B
- Authority
- CN
- China
- Prior art keywords
- target
- track
- ultra
- wideband radar
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional [2D] tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional [2D] tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/87—Combinations of radar systems, e.g. primary radar and secondary radar
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/50—Systems of measurement based on relative movement of target
- G01S17/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/66—Tracking systems using electromagnetic waves other than radio waves
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于不同体制雷达的波形与点迹融合检测方法,主要解决现有技术不能对波形数据和点迹数据融合检测的问题。其方案是:分别通过激光雷达和超宽带雷达获取目标的点迹数据和波形数据;对点迹数据信息进行轨迹生成,并按超宽带雷达检测帧时间对轨迹进行平滑滤波,得到基于点迹数据的目标轨迹;利用目标轨迹段反推此刻目标状态在波形数据中对应的目标回波,并在对应的波形数据中滤除该回波分量,再对其进行轨迹生成和区域平滑,得到基于波形数据的轨迹;将基于点迹数据的目标轨迹与基于波形数据的轨迹融合,得到融合后的目标轨迹。本发明能在检测时抑制虚警目标,改善精度稀释,减少多目标时的虚假目标,提高检测精度,可用于雷达探测。
Description
技术领域
本发明属于雷达技术领域,更进一步涉及一种目标融合检测方法,可用于雷达探测。
背景技术
在使用雷达传感器进行探测时,传感器根据测量数据类型的不同可以分为两类:测量值为波束的雷达与测量值为点迹的雷达,其中:
激光雷达Lidar是一种测量值为点迹的雷达,其能以较高的采样频率获取多个目标较为精确的二维位置,然而由于激光雷达发射和接收的激光信号波长短,其电磁波绕射能力弱,因而,在强光影响和障碍物遮挡的情况下,对某些方位不能得到有效测量值。如图1所示,阴影区域为雷达所不能探测到的区域,其中阴影区域I由干扰造成,阴影区域II由障碍物遮挡造成。
超宽带雷达UWB是一种测量值为波束的雷达,其通过多个雷达节点所得波形来进行二维定位,该雷达具有一维距离像分辨率高,作用距离远的特点。但在目标定位时有着两方面不足:
一是存在精度稀释问题,如图2a,图2b分别描述了一个由两组超宽带雷达传感器组成的一个目标检测系统,其中图2b中的目标相较图2a中的目标距离传感器更远,每个系统中的两个超宽带雷达传感器用不同的实框和空白框表示,阴影部分表示雷达系统对检测到的目标的预警区域。由于在超宽带雷达进行目标检测时,随着目标与传感器之间距离的增加,波形中属于目标部分的振幅会逐渐减小,使得靠近传感器的目标回波比远离传感器的目标回波强得多,造成超宽带雷达对远目标定位的精度稀释。从图2可看出,由于图2b中的目标相较图2a中的目标距离传感器更远,因而雷达系统对图2b中目标的检测区域相比于对图2a更加宽泛,表明其定位精度下降。
二是存在多目标检测时出现虚假目标问题,如图3所示。图3是由两组超宽带雷达传感器组成的一个检测系统在多目标检测情况下的示意图,两个超宽带雷达传感器的组成部分用不同的实框和空白框表示,阴影部分表示雷达系统对检测到的目标的预警区域,弧形的大小代表回波的强弱。当存在多个目标时,强目标的回波与弱目标的回波重叠区域会强于弱目标真实区域的回波,在图3中表示为虚警区域的强度大于目标二检测区域的强度。由此造成检测到两个强于弱目标的虚假目标。
由于使用单一种类传感器检测目标时存在上述局限性,为了实现更高的定位和测绘精度,利用不同种类传感器进行融合检测的方法近年来引起了越来越多的关注。融合来自超宽带雷达和激光雷达的数据具有两个方面的挑战性。第一,由于这两种传感器工作原理不同,它们采集的数据类型也不同,类型为波束的数据和类型为点迹的数据没有用于对应匹配和特征配准的共同特征。第二,对超宽带雷达世界坐标系与激光雷达世界坐标系之间的转换关系进行实时求解时缺乏测量环境的先验信息。基于上述两种原因,现有的目标融合检测技术均为基于点迹融合的处理方法,即将超宽带雷达网络所得点迹与激光雷达所得点迹进行联合处理。如申请号为CN106680806A,标题为一种多雷达点迹融合方法的专利,申请号为CN105093186A,标题为一种基于异构雷达传感网的多目标融合检测方法,其目标融合方法均采用点迹融合方式,这种方式在当雷达传感器测量值为波束时,无法改善精度稀释问题和多目标时的虚假目标问题。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种基于雷达波形与雷达点迹的目标融合检测方法,以减小精度稀释和多目标时的虚假目标,提高检测率。
为实现上述目的,本发明的技术方案包括如下:
(1)分别使用激光雷达传感器和超宽带雷达传感器获取目标的点迹数据信息和波形数据信息;
(2)使用基于三维图像处理的方法对目标的点迹数据信息进行轨迹生成,得到估计的目标区域,并按照超宽带雷达检测帧时间对该估计目标区域进行区域平滑,得到与超宽带雷达检测帧时间对齐的基于点迹数据的目标轨迹;
(3)利用该轨迹段反推此刻目标状态在超宽带雷达检测数据中对应的目标回波,并在对应的波形数据中滤除该回波分量;
(4)对回波分量滤除后的超宽带雷达波形数据进行轨迹生成和区域平滑处理,得到基于波形数据的目标轨迹;
(5)将基于点迹数据的目标轨迹和基于波形数据的目标轨迹融合,得到最终综合后的目标轨迹。
本发明于现有技术相比,具有如下优点:
1.本发明提出的技术涉及波形数据和点迹数据的融合,能够改善点迹融合技术所不能解决的精度稀释和多目标时的虚假目标问题。
2.本发明能通过波形抑制剔除在波形数据中测量值为点迹的雷达已检测到的目标所对应的回波,这种处理方法有利于在检测时抑制虚警目标。
3.本发明构建了一种联合处理波束数据传感器和点迹数据传感器数据的检测前跟踪框架,可以有效提高检测率。
附图说明
图1是现有激光雷达检测系统示意图;
图2是现有用两组超宽带雷达传感器组成的一个检测系统分别对相距传感器近距离目标和远距离目标进行探测的回波信号示意图;
图3是现有用两组超宽带雷达传感器组成的一个检测系统在多目标检测情况下的回波信号示意图;
图4是本发明的实现流程图;
图5是本发明仿真实验使用的场景说明图;
图6是本发明中激光雷达所得的点迹数据图;
图7是本发明中激光雷达所生成的目标轨迹图;
图8是本发明中超宽带雷达所得的波形数据图;
图9是本发明中超宽带雷达进行回波抑制后的波形数据图;
图10是本发明中超宽带雷达所得的目标轨迹图;
图11是用本发明融合处理后的综合目标轨迹图;
图12是用现有仅由超宽带雷达网络所得的目标轨迹图。
具体实施方式
以下结合附图,对本发明的实施例和效果进行详细描述:
参照图4,本实例的实现步骤如下:
步骤1.获取目标的点迹数据信息和波形数据信息。
目标的点迹数据信息和波形数据信息分别使用激光雷达和超宽带雷达获取,即:
使用激光雷达获得的点迹数据信息为:{xi,yi,ti|i=1,...,Ni,t1≤ti≤t2};
使用超宽带雷达获取的波束数据信息为:
其中,xi,yi表示第i个雷达点迹在搜索区域中的二维坐标,ti表示第i个雷达点迹在ti时刻采样,Ni表示激光雷达采样的点数;表示第s个雷达中,在第t帧中第j个距离分辨单元的幅度,每个距离分辨单位的长度为rz,Nc表示距离分辨单元的数目。
步骤2.获取动目标的测量值。
对激光雷达获取的点迹信息进行静态杂波抑制,将抑制后得到的点迹数据作为动目标的测量值:{xi,yi,ti|i=1,...,Nm,t1≤ti≤t2},其中,Nm为对激光雷达对动目标检测到的点迹数目。
步骤3.使用基于三维图像处理的方法对动目标的测量值进行轨迹生成。
(3.1)对空间时间进行三维划分并对点迹数据进行投票:
(3.1.1)将搜索区域按照x,y坐标轴分为Nx×Ny个宽度为Δxy的单元格,将采集数据的时间t1到t2分为宽度为Δt的Nt个单元格,即将空间和时间构成的三维区域划分为Nx×Ny×Nt个单元格,则该三维区域可用一个Nx×Ny×Nt的空矩阵表征;
(3.1.2)将激光雷达的每个点迹放入空矩阵中,完成对点迹数据的投票:
其中[·]表示四舍五入,表示一次投票后矩阵在位置的票数,表示一次投票前矩阵在位置的票数;
(3.2)用三维区域生长的方法检测目标区域:
(3.2.1)在三维矩阵中均匀选取种子单元格,以一个种子单元格作为中心向外生长,得到该单元格中完整的三维目标区域,并提取向该三维区域投票的点迹集合Rq:
Rq={xa,ya,ta|a=1,...,Nq},
其中,xa,ya表示ta时刻向三维区域投票的点迹在搜索区域的二维坐标,Nq表示该种子单元格中的点迹的数目;
(3.2.2)重复(3.2.1)直到获取全部的三维目标区域,将每一个包含点迹的三维目标区域看作一个目标轨迹,该区域包含的点迹即为该目标轨迹的测量值。
步骤4.按照超宽带雷达检测帧时间对估计目标区域进行区域平滑。
计算超宽带雷达检测帧的个数Nt:
其中,Δu为超宽带雷达网络帧间时间差,tb为雷达探测获取的点迹的时间戳;
按照超宽带雷达检测帧时间对估计目标区域进行区域平滑,得到时间平滑后的轨迹Tq:
Tq={xb,yb,tb|b=1,...,Nt}
其中,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标。
步骤5.利用点迹数据处理所得轨迹段反推目标状态在超宽带雷达检测数据中对应的目标回波。
(5.1)假设轨迹段Tq中包含NT个目标的轨迹,第s个超宽带雷达在各帧中到NT个目标的距离为
其中表示第s个超宽带雷达在第t帧中到第u个目标的距离,Nt表示超宽带雷达探测帧的个数,u为一个半径为du的圆形目标;
(5.2)提取超宽带雷达s检测到的波形数据中对应于目标的回波
其中,表示第s个超宽带雷达在第t帧中的第j个距离分辨单位的幅度,rz表示每个距离分辨单位的长度;
(5.3)根据提取出对应于目标的回波估计出目标u的回波分量强度
步骤6.在超宽带雷达检测到的波形数据中滤除对应于目标的回波分量。
在每个超宽带雷达波形对应于目标的回波中减去目标回波分量得到滤除已检测目标回波分量的波形测量值
其中,表示第s个超宽带雷达在第t帧中的第j个距离分辨单位的幅度,表示第s个超宽带雷达在第t帧中到第u个目标的距离,rz表示每个距离分辨单位的长度,du为圆形目标u的半径。
步骤7.使用滤除已检测目标回波分量的波形测量值进行轨迹生成和区域平滑处理。
(7.1)使用滤除已检测目标回波分量的波形测量值对搜索区域生成轨迹:
(7.1.1)将搜索区域按照x,y坐标轴分为Nx×Ny个宽度为Δxy的单元格,将采集数据的时间t1到t2分为宽度为Δu的Nt个单元格,每帧对应一个时间单元格,即将空间和时间构成的三维区域划分为Nx×Ny×Nt个单元格;
(7.1.2)对Nr个UWB雷达的第t帧数据的搜索区域成像:
设置检测距离检测门限Tr和波形数据幅度检测门限Tf,对搜索区域内距离j大与距离检测门限Tr的单元格进行投票,即当测量值大于幅度检测门限Tf时,对以雷达s为中心,jrz为半径的圆所经过的单元格中的票数加1;
遍历第t帧所有波形数据,得到一个Nx×Ny矩阵该矩阵表示雷达s在t帧时的搜索区域像,将Nr个雷达所得搜索区域像相乘,得到第t帧的综合区域像Φt:
(7.1.3)将超宽带雷达检测帧内的Nt个综合区域像{Φt|t=1,...,Nt}组成三维矩阵,使用三维区域生长方法对三维矩阵进行成像处理,得到超宽带雷达的目标区域;
(7.2)对每一个三维目标区域进行腐蚀和膨胀操作使其平滑,再对每帧目标区域进行凝聚操作,获取目标轨迹Tp:
Tp={xb,yb,tb|b=1,...,Nt}
其中,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标,Nt表示超宽带雷达检测帧数。
步骤8.将基于点迹数据的目标轨迹和基于波形数据的目标轨迹融合。
(8.1)对激光雷达目标轨迹与超宽带雷达目标轨迹从时间上进行外推,得到对应于{tb|b=Nt+1,...,Nt+w}的外推轨迹测量值,再将外推轨迹测量值与原轨迹关联后得到目标轨迹:
Tp={xb,yb,tb|b=1,...,Nt+w}
Tq={xb,yb,tb|b=1,...,Nt+w}
其中,Tp表示超宽带雷达目标轨迹,Tq表示激光雷达目标轨迹,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标,Nt表示超宽带雷达检测帧数,w表示外推的时间;
(8.2)枚举激光雷达和超宽带雷达轨迹所有可能的关联形式,并计算每一种可能的关联形式中激光雷达轨迹与超宽带雷达轨迹的最优次模式分配距离OSPA;
(8.3)将计算出的最优次模式分配距离最小的一对轨迹作为两种轨迹的最优关联形式,利用该最优关联形式对激光雷达和超宽带雷达轨迹进行融合,得到最终综合的目标轨迹。
本发明效果可通过以下实验效果作进一步描述:
一.实验条件:
构建实验场景如图5所示,其包括2个移动的行人目标,分别为目标1,目标2,三个超宽带雷达、一个激光雷达,使用三个超宽带雷达和一个激光雷达组成的雷达网络对一个与走廊相连接的大厅中两个行人目标进行监测,大厅的一边是玻璃,从玻璃中照入的阳光和大厅中的立柱均对激光雷达目标探测产生影响与遮挡。
二.实验的内容与结果分析:
仿真实验1,基于上述实验条件,用本发明场景中的激光雷达获取周围环境的点迹数据,并对点迹数据进行静态杂波抑制操作,结果如图6所示,其中虚框标出的部分分别对应为目标1和目标2的点迹数据,从图6可以看到,目标1的轨迹可以完整的探测到,但是由于遮挡作用,目标2只有部分轨迹可以被探测到;
仿真实验2,基于上述实验条件,将仿真实验1所得到的目标1和目标2这两个移动目标测量值分别经过基于三维区域生长的目标轨迹检测算法生成的目标轨迹,结果如图7所示。从图7可以看出目标2的轨迹有4部分的中断,如图7目标2虚线框标注,这是因为目标2所在部分区域由于干扰和遮挡无法被检测到。
仿真实验3,基于上述实验条件,用本发明场景中的三个超宽带雷达获取周围环境的波形数据,结果如图8所示,其中,图8a是第一个超宽带雷达1获取的数据,图8b是第二个超宽带雷达2获取的数据,图8c是第三个超宽带雷达3获取的数据。从图8可以看出,目标回波随着距离增大而减弱。
仿真实验4,基于上述实验条件,从仿真实验3中的波形数据中滤除已检测目标回波分量,得到回波分量抑制后的波形数据,如图9所示,其中,图9a是第一个超宽带雷达1经回波分量抑制后的波形数据,图9b是第二个超宽带雷达2经回波分量抑制后的波形数据,图9c是第三个超宽带雷达3经回波分量抑制后的波形数据。图9中各个雷达经回波分量抑制后的波形数据即对应激光雷达未能检测到的目标轨迹。
仿真实验5,基于上述实验条件,对仿真实验4得到的波形数据进行成像,生成轨迹,结果如图10所示,图10中虚框标出部分为利用回波抑制后的超宽带雷达波形数据所成像的轨迹,即目标2未被激光雷达检出的轨迹。
仿真实验6,基于上述实验条件,将仿真实验2中得到的激光雷达轨迹与仿真实验5中得到的超宽带雷达轨迹进行关联融合,得到目标1融合后的轨迹和目标2融合后的轨迹,结果如图11所示。其中图11中目标1和目标2的最终融合轨迹是由激光雷达所得轨迹与超宽带雷达所得轨迹融合后构成的,目标2轨迹中激光雷达漏检部分由超宽带雷达检测出并补偿,可以看出融合后轨迹与真实轨迹一致。
仿真实验7,单纯使用仿真实验3中超宽带雷达所得波形数据进行多帧联合检测生成目标1和目标2轨迹,结果如图12所示,从图12中可以看出,目标1和目标2的生成轨迹中存在一些虚假目标。
图11和图12的结果表明,本发明可以有效减少精度稀释和多目标时的虚假目标问题,并能在检测时抑制虚警,获得更精确的目标跟踪轨迹。
Claims (8)
1.一种基于不同体制雷达的波形与点迹融合检测方法,其特征在于,包括:
(1)分别使用激光雷达传感器和超宽带雷达传感器获取目标的点迹数据信息和波形数据信息;
(2)使用基于三维图像处理的方法对目标的点迹数据信息进行轨迹生成,得到估计的目标区域,并按照超宽带雷达检测帧时间对该估计目标区域进行区域平滑,得到与超宽带雷达检测帧时间对齐的基于点迹数据的目标轨迹;
(3)利用轨迹反推此刻目标状态在超宽带雷达检测数据中对应的目标回波,并在对应的波形数据中滤除该回波分量;
(4)对回波分量滤除后的超宽带雷达波形数据进行轨迹生成和区域平滑处理,得到基于波形数据的目标轨迹;
(5)将基于点迹数据的目标轨迹和基于波形数据的目标轨迹融合,得到最终综合后的目标轨迹。
2.如权利要求1所述的方法,其特征在于,(1)中获取到的点迹数据信息和波形数据信息,分别表示如下:
点迹数据信息为:{xi,yi,ti|i=1,...,Ni,t1≤ti≤t2};
波束数据信息为:
其中,xi,yi表示第i个激光雷达传感器点迹在搜索区域中的二维坐标,Ni表示激光雷达采样的点数;表示第s个超宽带雷达传感器中,在第t帧中第j个距离分辨单元的幅度,Nc表示距离分辨单元的数目,ti表示第i个激光雷达点迹在t时刻采样,t1、t2分别表示多帧联合处理中的起始时间和结束时间;Nr表示超宽带雷达传感器的数量。
3.如权利要求2所述的方法,其特征在于,(2)中使用基于三维图像处理的方法对目标的点迹数据信息进行轨迹生成,实现如下:
(2a)对空间时间进行三维划分并对点迹数据进行投票:
(2a1)将搜索区域按照x,y坐标轴分为Nx×Ny个宽度为Δxy的单元格,将采集数据的时间t1到t2分为宽度为Δt的Nt个单元格,即将空间和时间构成的三维区域划分为Nx×Ny×Nt个单元格,则该三维区域可用一个Nx×Ny×Nt的空矩阵表征;t1≤ti≤t2;
(2a2)将激光雷达的每个点迹放入空矩阵中,完成对点迹数据的投票:
其中[·]表示四舍五入,表示一次投票后矩阵在位置的票数,表示一次投票前矩阵在位置的票数;
(2b)用三维区域生长的方法检测目标区域:
(2b1)在三维矩阵中均匀选取种子单元格,以一个种子单元格作为中心向外生长,得到该单元格中完整的三维目标区域,并提取向该三维区域投票的点迹集合:
Rq={xa,ya,ta|a=1,...,Nq},
其中,xa,ya表示ta时刻向三维区域投票的点迹在搜索区域的二维坐标,Nq表示该种子单元格中的点迹的数目;
(2b2)重复(2b1)直到获取全部的三维目标区域,将每一个包含点迹的三维目标区域看作一个目标轨迹,该区域包含的点迹即为该目标轨迹的测量值。
4.如权利要求1所述的方法,其特征在于,(2)中按照超宽带雷达检测帧时间对估计目标区域进行区域平滑,通过如下公式进行:
Tq={xb,yb,tb|b=1,...,Nt}
tb=bΔu
其中,Tq为平滑后的轨迹,Δu为超宽带雷达网络帧间时间差,Nt表示超宽带雷达帧的个数,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标。
5.如权利要求4所述的方法,其特征在于,(3)中利用点迹数据处理所得轨迹段反推目标状态在超宽带雷达检测数据中对应的目标回波,实现如下:
(3a1)假设该轨迹段中包含NT个目标的轨迹,估计第s个超宽带雷达在各帧中到NT个目标的距离:
其中表示第s个超宽带雷达在第t帧中到第u个目标的距离,Nt表示超宽带雷达探测帧的个数,u为一个半径为du的圆形目标;
(3a2)提取超宽带雷达s检测到的波形数据中对应于目标的回波:
其中,表示第s个超宽带雷达在第t帧中的第j个距离分辨单位的幅度,rz表示每个距离分辨单位的长度;
(3a3)通过提取出的对应于目标的回波估计出目标u的回波分量强度
6.如权利要求1所述的方法,其特征在于,(3)中在超宽带雷达检测到的波形数据中滤除对应于目标的回波分量,是在每个超宽带雷达波形对应于目标的回波中减去目标回波分量得到滤除已检测目标回波分量的波形测量值
其中,表示第s个超宽带雷达在第t帧中的第j个距离分辨单位的幅度,表示第s个超宽带雷达在第t帧中到第u个目标的距离,rz表示每个距离分辨单位的长度,du为圆形目标u的半径。
7.如权利要求1所述的方法,其特征在于,(4)中对滤除已检测目标回波分量的波形测量值进行轨迹生成和区域平滑处理,实现如下:
(4a)使用回波抑制后的波形数据对搜索区域生成轨迹:
(4a1)将搜索区域按照x,y坐标轴分为Nx×Ny个宽度为Δxy的单元格,将采集数据的时间t1到t2分为宽度为Δu的Nt个单元格,每帧对应一个时间单元格,即将空间和时间构成的三维区域划分为Nx×Ny×Nt个单元格;
(4a2)对Nr个UWB雷达的第t帧数据的搜索区域成像:
设置检测距离检测门限Tr和波形数据幅度检测门限Tf,对搜索区域内距离j大与距离检测门限Tr的单元格进行投票,即当测量值大于幅度检测门限Tf时,对以雷达s为中心,jrz为半径的圆所经过的单元格中的票数加1;
遍历完第t帧所有波形数据后可以得到一个Nx×Ny矩阵该矩阵表示雷达s在t帧时的搜索区域像,将Nr个雷达所得搜索区域像相乘,得到第t帧的综合区域像Φt:
(4a3)将超宽带雷达检测帧内的Nt个综合区域像{Φt|t=1,...,Nt}组成三维矩阵,使用三维区域生长方法对三维矩阵进行处理,得到超宽带雷达的目标区域;
(4b)对每一个三维目标区域进行腐蚀和膨胀操作使其平滑,再对每帧目标区域进行凝聚操作,获取目标轨迹Tp:
Tp={xb,yb,tb|b=1,...,Nt}
其中,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标,Nt表示超宽带雷达检测帧数。
8.如权利要求1所述的方法,其特征在于,(5)中将基于点迹数据的目标轨迹和基于波形数据的目标轨迹融合,实现如下:
(5a)对激光雷达目标轨迹与超宽带雷达目标轨迹从时间上进行外推,得到对应于{tb|b=Nt+1,...,Nt+w}的外推轨迹测量值,得到目标轨迹:
Tp={xb,yb,tb|b=1,...,Nt+w}
Tq={xb,yb,tb|b=1,...,Nt+w}
其中,Tp表示超宽带雷达目标轨迹,Tq表示激光雷达目标轨迹,xb,yb表示对应于时间tb的目标轨迹段的测量值的二维坐标,Nt表示超宽带雷达检测帧数,w表示外推的时间;
(5b)枚举激光雷达和超宽带雷达轨迹所有可能的关联形式,并计算每一对激光雷达轨迹与超宽带雷达轨迹的最优次模式分配距离;
(5c)将最优次模式分配距离最小的一对轨迹作为两种轨迹的最优关联形式,使用最优关联形式对激光雷达和超宽带雷达轨迹进行融合,得到最终综合的目标轨迹。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111592160.6A CN114488112B (zh) | 2021-12-23 | 2021-12-23 | 一种基于不同体制雷达的波形与点迹融合检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111592160.6A CN114488112B (zh) | 2021-12-23 | 2021-12-23 | 一种基于不同体制雷达的波形与点迹融合检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114488112A CN114488112A (zh) | 2022-05-13 |
| CN114488112B true CN114488112B (zh) | 2024-10-22 |
Family
ID=81493218
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111592160.6A Active CN114488112B (zh) | 2021-12-23 | 2021-12-23 | 一种基于不同体制雷达的波形与点迹融合检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114488112B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115015905A (zh) * | 2022-05-25 | 2022-09-06 | 北京理工雷科雷达技术研究院有限公司 | 一种基于多帧轨迹关联的一维距离像船只检测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106680806A (zh) * | 2016-11-24 | 2017-05-17 | 清华大学 | 一种多雷达点迹融合方法 |
| CN109212514A (zh) * | 2018-09-29 | 2019-01-15 | 河北德冠隆电子科技有限公司 | 一种雷达检测设备对运动以及静止目标持续跟踪关联方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3664110B2 (ja) * | 2001-07-04 | 2005-06-22 | 日産自動車株式会社 | 物体種別判定装置及び物体種別判定方法 |
| US8717545B2 (en) * | 2009-02-20 | 2014-05-06 | Digital Signal Corporation | System and method for generating three dimensional images using lidar and video measurements |
| US8098193B2 (en) * | 2009-11-05 | 2012-01-17 | Honeywell International Inc. | Digitally controlled UWB millimeter wave radar |
| CN112285700B (zh) * | 2020-08-24 | 2023-12-15 | 江苏大学 | 基于激光雷达与毫米波雷达融合的机动目标跟踪方法 |
| CN112232248B (zh) * | 2020-10-22 | 2023-04-07 | 中国人民解放军战略支援部队信息工程大学 | 一种多线LiDAR点云数据平面特征的提取方法及装置 |
| CN112613532B (zh) * | 2020-11-26 | 2023-08-01 | 西安电子科技大学 | 基于雷达与循环神经网络补全红外融合的动目标跟踪方法 |
| CN113376648B (zh) * | 2021-04-26 | 2023-02-03 | 上海航天控制技术研究所 | 基于激光雷达探测的高速非合作目标轨迹提取方法 |
| CN113538410B (zh) * | 2021-08-06 | 2022-05-20 | 广东工业大学 | 一种基于3d激光雷达和uwb的室内slam建图方法 |
-
2021
- 2021-12-23 CN CN202111592160.6A patent/CN114488112B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106680806A (zh) * | 2016-11-24 | 2017-05-17 | 清华大学 | 一种多雷达点迹融合方法 |
| CN109212514A (zh) * | 2018-09-29 | 2019-01-15 | 河北德冠隆电子科技有限公司 | 一种雷达检测设备对运动以及静止目标持续跟踪关联方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114488112A (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106918807A (zh) | 一种雷达回波数据的目标点迹凝聚方法 | |
| CN112051568B (zh) | 一种两坐标雷达的俯仰测角方法 | |
| CN101963662B (zh) | 基于短时分数阶傅里叶域滤波的自聚焦预处理方法 | |
| CN101196562A (zh) | 一种基于改进的em算法的激光雷达波形数据分解的方法 | |
| CN116699602A (zh) | 一种基于毫米波雷达和相机融合的目标检测系统及方法 | |
| CN114814832A (zh) | 基于毫米波雷达的人体跌倒行为实时监测系统及监测方法 | |
| CN113552563B (zh) | 垂测信息和高频地波雷达杂波信息的对应性分析方法 | |
| CN118334736A (zh) | 基于毫米波雷达的多目标身份识别与行为监测方法 | |
| Li et al. | Pedestrian liveness detection based on mmwave radar and camera fusion | |
| CN109100697B (zh) | 一种基于地面监视雷达系统的目标凝聚方法 | |
| CN115061113A (zh) | 用于雷达的目标检测模型训练方法、装置及存储介质 | |
| CN108983194B (zh) | 一种基于地面监视雷达系统的目标提取及凝聚方法 | |
| CN120630199B (zh) | 一种井下环境激光雷达与毫米波雷达协同slam方法 | |
| Liu et al. | High-precision imaging and dense vehicle target detection method for airborne single-pass circular synthetic aperture radar | |
| CN115166764B (zh) | 一种毫米波雷达辅助的激光雷达运动目标实时去除方法 | |
| CN114488112A (zh) | 一种基于不同体制雷达的波形与点迹融合检测方法 | |
| CN111624605B (zh) | 一种基于角度维回波特征的航海雷达目标检测方法 | |
| CN117647788B (zh) | 一种基于人体3d点云的危险行为识别方法及设备 | |
| CN116879863B (zh) | 一种连续波4d毫米波雷达多目标测量方法及系统 | |
| CN117347995A (zh) | 一种基于车载传感器的目标融合方法及系统 | |
| CN114169402B (zh) | 一种基于门控循环单元的探地雷达地下空洞目标自动识别方法 | |
| CN118584474A (zh) | 一种基于毫米波雷达的无人机多目标检测系统及方法 | |
| CN113534096B (zh) | 基于样条函数的LiDAR信号几何特征提取方法及系统 | |
| CN117872341A (zh) | 一种基于距离-多普勒域的毫米波雷达人体存在检测方法 | |
| CN118112530A (zh) | 一种具有较强泛化能力的sar地面运动目标智能检测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |