CN117192533A - Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search - Google Patents
Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search Download PDFInfo
- Publication number
- CN117192533A CN117192533A CN202311164769.2A CN202311164769A CN117192533A CN 117192533 A CN117192533 A CN 117192533A CN 202311164769 A CN202311164769 A CN 202311164769A CN 117192533 A CN117192533 A CN 117192533A
- Authority
- CN
- China
- Prior art keywords
- multipath
- point
- radar
- reflection
- ghost
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000003068 static effect Effects 0.000 claims abstract description 13
- 238000000926 separation method Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
技术领域Technical field
本发明属于车载毫米波多径信号处理技术领域,特别涉及一种多径鬼影识别技术。The invention belongs to the technical field of vehicle-mounted millimeter wave multipath signal processing, and particularly relates to a multipath ghost recognition technology.
背景技术Background technique
在利用车载毫米波雷达感知环境时,电磁波信号在不同物体之间的反射、衍射传播常常引起繁杂的多径信号。当多径信号与真实目标在距离多普勒谱中位于相同的分辨率单元时,会降低真实信号的回波能量,进而造成漏检;当多径信号与真实目标在距离多普勒谱中位于不相同的分辨率单元时,多径信号会形成鬼影,进而造成虚警。在自动驾驶场景下,多径信号造成的虚警(多径鬼影)现象更加常见,其带来的干扰也更加严重。When using vehicle-mounted millimeter wave radar to sense the environment, the reflection and diffraction propagation of electromagnetic wave signals between different objects often cause complex multipath signals. When the multipath signal and the real target are in the same resolution unit in the range Doppler spectrum, the echo energy of the real signal will be reduced, resulting in missed detection; when the multipath signal and the real target are in the range Doppler spectrum, When located in different resolution units, multipath signals will form ghosts and cause false alarms. In autonomous driving scenarios, false alarms (multipath ghosts) caused by multipath signals are more common, and the interference they cause is more serious.
近年来,针对车载毫米波雷达多径鬼影干扰问题,国内外研究机构开展了相关技术研究。2020年,C.Liu等人针对多径鬼影干扰问题,提出了一种基于微小距离差值的多径鬼影识别方法(C.Liu,S.Liu,C.Zhang,Y.Huang and H.Wang,"Multipath propagationanalysis and ghost target removal for FMCW automotive radars,"IETInternational Radar Conference(IET IRC 2020),Online Conference,2020,pp.330-334,doi:10.1049/icp.2021.0554.),然而,该方法应用场景十分受限,当目标距离雷达较近时,多径鬼影与真实目标相对于雷达距离的差值较大,上述方法无法有效识别多径鬼影;并且,针对相邻车道内距离相近、角度不同的两个真实目标也容易产生误识别。2023年,Yunda Li等人提出了一种基于DOA/DOD估计的多径鬼影识别方法(Y.Li and X.Shang,"Multipath Ghost Target Identification for Automotive MIMO Radar,"2022IEEE96th Vehicular Technology Conference(VTC2022-Fall),London,United Kingdom,2022,pp.1-5,doi:10.1109/VTC2022-Fall57202.2022.10012904.),通过区分天线的发射方向和接收方向的差异,进而识别由发射方向和入射方向不同的电磁传播路径引起的多径鬼影。然而,上述方法不仅无法解决接发射方向与接收方向相等情况下的多径鬼影,而且还受到天线布阵的严重影响。In recent years, domestic and foreign research institutions have carried out relevant technical research on the problem of multipath ghost interference in vehicle-mounted millimeter wave radar. In 2020, C. Liu et al. proposed a multipath ghost identification method based on small distance differences to solve the problem of multipath ghost interference (C. Liu, S. Liu, C. Zhang, Y. Huang and H .Wang, "Multipath propagation analysis and ghost target removal for FMCW automotive radars," IETInternational Radar Conference (IET IRC 2020), Online Conference, 2020, pp.330-334, doi:10.1049/icp.2021.0554.), however, this method The application scenario is very limited. When the target is close to the radar, the difference between the multipath ghost and the real target relative to the radar distance is large, and the above method cannot effectively identify the multipath ghost; and for adjacent lanes with similar distances , Two real targets with different angles are also prone to misrecognition. In 2023, Yunda Li and others proposed a multipath ghost identification method based on DOA/DOD estimation (Y.Li and X.Shang, "Multipath Ghost Target Identification for Automotive MIMO Radar," 2022IEEE96th Vehicular Technology Conference (VTC2022- Fall), London, United Kingdom, 2022, pp.1-5, doi:10.1109/VTC2022-Fall57202.2022.10012904.), by distinguishing the difference between the transmitting direction and the receiving direction of the antenna, and then identifying the differences between the transmitting direction and the incident direction. Multipath ghosting caused by electromagnetic propagation paths. However, the above method not only cannot solve the multipath ghosting when the transmitting and receiving directions are equal, but is also seriously affected by the antenna array.
总之,上述方法均存在一定缺陷,普适性不强。然而,多径鬼影若无法得到有效识别与抑制,多径鬼影引起的虚警将对场景中的真实目标探测带来严重干扰,不利于车载毫米波雷达的广泛应用。因此,深入研究普适性更强、效用更佳的车载毫米波雷达多径鬼影干扰识别方法极具重要价值。In short, the above methods all have certain shortcomings and are not universally applicable. However, if multipath ghosts cannot be effectively identified and suppressed, the false alarms caused by multipath ghosts will seriously interfere with the detection of real targets in the scene, which is not conducive to the widespread application of automotive millimeter wave radar. Therefore, it is of great value to conduct in-depth research on multipath ghost interference identification methods for vehicle-mounted millimeter wave radars that are more universal and more effective.
发明内容Contents of the invention
为解决上述技术问题,本发明提供一种基于反射点反演搜索的车载毫米波雷达多径鬼影识别方法。In order to solve the above technical problems, the present invention provides a vehicle-mounted millimeter wave radar multipath ghost identification method based on reflection point inversion search.
本发明采用的技术方案为:一种基于反射点反演搜索的车载毫米波雷达多径鬼影识别方法,包括:The technical solution adopted by the present invention is: a vehicle-mounted millimeter wave radar multipath ghost identification method based on reflection point inversion search, including:
S1、根据当前帧原始回波的ADC数据生成雷达点云;S1. Generate a radar point cloud based on the ADC data of the original echo of the current frame;
S2、对雷达点云进行动静点云分离,然后对分离得到的动态点云进行聚类,最后提取聚类得到的每一个簇的中心坐标与表征该中心坐标的多域信息;S2. Separate the dynamic and static point clouds of the radar point cloud, then cluster the separated dynamic point clouds, and finally extract the center coordinates of each cluster obtained by clustering and the multi-domain information characterizing the center coordinates;
S3、对步骤S2得到的簇中心坐标集合中任取两点进行关联,构成匹配集;S3. Correlate any two points from the cluster center coordinate set obtained in step S2 to form a matching set;
S4、从匹配集中取出一对关联点,基于反射点反演求解,得到搜索区域;S4. Take a pair of related points from the matching set, perform inversion and solve based on the reflection points, and obtain the search area;
S5、在步骤S4得到的搜索区域中搜索是否存在步骤S1得到的雷达点云,若存在,则执行步骤S6;否则返回步骤S4;S5. Search whether the radar point cloud obtained in step S1 exists in the search area obtained in step S4. If it exists, perform step S6; otherwise, return to step S4;
S6、对于经步骤S5处理后的这对关联点,基于多域信息进行多径鬼影判决;S6. For the pair of related points processed in step S5, perform multipath ghost judgment based on multi-domain information;
S7、遍历匹配集中的每一对关联点,重复步骤S4~S6,完成当前帧数据下的所有点云的多径鬼影识别。S7. Traverse each pair of associated points in the matching set, repeat steps S4 to S6, and complete multipath ghost recognition of all point clouds under the current frame data.
本发明的有益效果:本发明提出的一种基于反射点反演求解的车载毫米波雷达多径鬼影识别方法,可有效识别车载毫米波雷达回波中的多径鬼影。通过对MIMO毫米波雷达回波的预处理,生成雷达点云;通过雷达点云处理步骤,实现动静点云分离,动态点云的聚类,簇中心坐标以及多域信息获取;通过目标关联,将场景内的簇中心进行关联;计算三种假设下的反射点坐标;通过点云搜索步骤初步识别多径鬼影,最后通过多域特进一步判决多径鬼影。实测结果表明,本发明方法能有效识别车载毫米波雷达回波中的多径鬼影。Beneficial effects of the invention: The invention proposes a vehicle-mounted millimeter-wave radar multi-path ghost identification method based on reflection point inversion solution, which can effectively identify multi-path ghosts in vehicle-mounted millimeter-wave radar echoes. Through preprocessing of MIMO millimeter wave radar echo, radar point cloud is generated; through radar point cloud processing steps, dynamic and static point cloud separation, dynamic point cloud clustering, cluster center coordinates and multi-domain information acquisition are achieved; through target association, Correlate the cluster centers in the scene; calculate the reflection point coordinates under the three assumptions; initially identify multipath ghosts through the point cloud search step, and finally further determine the multipath ghosts through multi-domain characteristics. Actual measurement results show that the method of the present invention can effectively identify multipath ghosts in vehicle millimeter wave radar echoes.
附图说明Description of the drawings
图1为典型车载毫米波雷达多径信号传播示意图;Figure 1 is a schematic diagram of multipath signal propagation of a typical vehicle-mounted millimeter wave radar;
图2为点云动静分离的参考坐标系示例图;Figure 2 is an example diagram of the reference coordinate system for dynamic and static separation of point clouds;
图3为不同类型多径鬼影的几何特征示意图;Figure 3 is a schematic diagram of the geometric characteristics of different types of multipath ghosts;
其中,(a)为GS、T、P的几何特征示意图,(b)为GS、GM2、T'的几何特征示意图,(c)为O、P、T、GM1的几何特征示意图;Among them, (a) is a schematic diagram of the geometric characteristics of G S , T, and P, (b) is a schematic diagram of the geometric characteristics of G S , G M2 , and T', (c) is a schematic diagram of the geometric characteristics of O, P, T, and G M1 ;
图4为实验场景图;Figure 4 is the experimental scene diagram;
图5为实验数据处理结果;Figure 5 shows the experimental data processing results;
其中,(a)为实验场景下多帧累积的点云,(b)为运动目标点云及多径鬼影,(c)为多径鬼影识别结果,(d)为多径鬼影消除结果。Among them, (a) is the point cloud accumulated in multiple frames in the experimental scene, (b) is the moving target point cloud and multipath ghost, (c) is the multipath ghost recognition result, (d) is the multipath ghost elimination result.
具体实施方式Detailed ways
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。In order to facilitate those skilled in the art to understand the technical content of the present invention, the content of the present invention is further explained below with reference to the accompanying drawings.
图1为典型车载毫米波雷达多径信号传播示意图;雷达位置为点O,真实目标位置为点T,场景中存在一个反射点为点P;一般地,设电磁波反射为镜面反射;因此,电磁波的传播路径主要有以下几类:Figure 1 is a schematic diagram of multipath signal propagation of a typical vehicle-mounted millimeter wave radar; the radar position is point O, the real target position is point T, and there is a reflection point in the scene as point P; generally speaking, the electromagnetic wave reflection is assumed to be specular reflection; therefore, the electromagnetic wave The transmission paths mainly fall into the following categories:
直视路径:首先存在一条探测目标的双程直视路径,O→T→O;其次,还存在一条由反射点P后向散射引起的双程直视路径O→P→O;Direct-looking path: First, there is a two-way direct-looking path for detecting the target, O→T→O; secondly, there is also a two-way direct-looking path O→P→O caused by the backscattering of the reflection point P;
一阶多径:一阶路径指经过反射点P一次的电磁传播路径,在上述场景下存在两条一阶多径,一条是,O→T→P→O,另一条是O→P→T→O;First-order multipath: The first-order path refers to the electromagnetic propagation path that passes through the reflection point P once. In the above scenario, there are two first-order multipaths, one is O→T→P→O, and the other is O→P→T →O;
二阶多径:二阶多径指经过反射点P两次的电磁传播路径,具体的:O→P→T→P→O;Second-order multipath: Second-order multipath refers to the electromagnetic propagation path that passes through the reflection point P twice, specifically: O→P→T→P→O;
高阶多径:高阶多径指经过反射点P多次的电磁传播路径。一般的,由于电磁反射对信号的能量有较强的衰减,进多次反射的高阶路径也因信号能量微弱而可以忽略不计。High-order multipath: High-order multipath refers to the electromagnetic propagation path that passes through the reflection point P multiple times. Generally, since electromagnetic reflection has a strong attenuation of signal energy, the high-order path leading to multiple reflections can be ignored because the signal energy is weak.
电磁传播路径与检测目标的关系如表1所示:The relationship between electromagnetic propagation paths and detection targets is shown in Table 1:
表1电磁传播路径与检测目标的关系表Table 1 Relationship between electromagnetic propagation paths and detection targets
实验场景如图4所示,车载毫米波雷达向场景内持续地发射电磁波信号,接收天线接收回波信号,经混频、滤波、数字采样后得到原始的模数转换(Analog-to-DigitalConverter,ADC)数据。如图1所示,在典型车载毫米波雷达应用场景下,由反射点P所产生的多径鬼影主要有:The experimental scene is shown in Figure 4. The vehicle-mounted millimeter wave radar continuously emits electromagnetic wave signals into the scene, and the receiving antenna receives the echo signal. After mixing, filtering, and digital sampling, the original analog-to-digital converter (Analog-to-DigitalConverter, ADC) data. As shown in Figure 1, in a typical automotive millimeter-wave radar application scenario, the multipath ghosts generated by the reflection point P mainly include:
Typle-1型二次反射多径鬼影GM1,对应的电磁波传播路径为O→P→T→O;Type-1 secondary reflection multipath ghost G M1 , the corresponding electromagnetic wave propagation path is O→P→T→O;
Typle-2型二次反射多径鬼影GM2,对应的电磁波传播路径为O→T→P→O;Type-2 secondary reflection multipath ghost G M2 , the corresponding electromagnetic wave propagation path is O→T→P→O;
Typle-2型三次反射多径鬼影Gs,对应的电磁波传播路径为O→P→T→P→O;Type-2 triple reflection multipath ghost G s , the corresponding electromagnetic wave propagation path is O→P→T→P→O;
本领域技术人员应知本发明中的Typle-1型表示电磁波最终从真实目标处反射回雷达,Typle-2型表示电磁波最终从反射点处反射回雷达,图1中所示箭头为电磁波的传播方向。Those skilled in the art should know that Type-1 in the present invention means that the electromagnetic wave is finally reflected from the real target back to the radar, and Type-2 means that the electromagnetic wave is finally reflected from the reflection point back to the radar. The arrows shown in Figure 1 represent the propagation of electromagnetic waves. direction.
本发明的方法包括以下步骤:The method of the present invention includes the following steps:
S1:信号预处理;S1: signal preprocessing;
信号预处理步骤实现原始回波的ADC数据到雷达点云的生成过程。主要步骤包括快时间维快速傅里叶变换、慢时间维快速傅里叶变换、恒虚警检测、波达角度估计和坐标转换工作。假设处理一帧雷达ADC数据后,生成的雷达点云集合为P={p0,p1,p2,...,pNP},其中NP表示点云的个数,pi={xi,yi,vi,RCSi}表示第i个雷达点云及其多域信息。The signal preprocessing step realizes the generation process of the ADC data of the original echo to the radar point cloud. The main steps include fast time dimension fast Fourier transform, slow time dimension fast Fourier transform, constant false alarm detection, wave arrival angle estimation and coordinate conversion work. Assume that after processing one frame of radar ADC data, the generated radar point cloud set is P = {p 0 , p 1 , p 2 ,..., p NP }, where N P represents the number of point clouds, p i ={ x i , y i , vi , RCS i } represents the i-th radar point cloud and its multi-domain information.
进一步地,xi,yi,ri,vi,RCSi分别表示汽车坐标系下雷达点云的横坐标、汽车坐标系下雷达点云的纵坐标、雷达点云相对于雷达的距离、雷达点云相对于雷达的径向速度和雷达点云的雷达截面积计算值。预处理结果所得的多帧点云累积结果如图5(a)所示。Further, x i , y i , r i , vi , RCS i respectively represent the abscissa coordinate of the radar point cloud in the automobile coordinate system, the ordinate coordinate of the radar point cloud in the automobile coordinate system, the distance of the radar point cloud relative to the radar, The radial velocity of the radar point cloud relative to the radar and the calculated radar cross-sectional area of the radar point cloud. The multi-frame point cloud accumulation result obtained from the preprocessing result is shown in Figure 5(a).
S2:雷达点云处理;S2: Radar point cloud processing;
2.1动静点云分离2.1 Separation of dynamic and static point clouds
通过步骤S1所获得的雷达点云既包括了场景中运动目标的点云,又包括场景静物的点云。通过点云分离步骤,将原始雷达点云做动静分离处理,以便于后续处理。具体的,参考如图2所示的汽车坐标系,以汽车前进方向为y轴,垂直于y轴的水平方向为x轴,假设目标与雷达同高度,理论上地,车载毫米波雷达测得的静态物体(相对地面静止)的径向速度为:The radar point cloud obtained in step S1 includes both point clouds of moving targets in the scene and point clouds of still objects in the scene. Through the point cloud separation step, the original radar point cloud is separated from static and moving parts to facilitate subsequent processing. Specifically, refer to the car coordinate system shown in Figure 2, with the car's forward direction as the y-axis and the horizontal direction perpendicular to the y-axis as the x-axis. Assume that the target is at the same height as the radar. Theoretically, the vehicle-mounted millimeter wave radar measures The radial velocity of a static object (stationary relative to the ground) for:
其中,(xr,yr)为雷达在汽车坐标系下的位置,θr为雷达法线方向与汽车坐标系的夹角,vC,x,vC,y,w分别表示运动汽车的横纵向速度分量和转速。当雷达测得的点云速度与差异较大,则说明该点云是相对地面运动的。由此,点云的动静分离依据为:Among them, (x r , y r ) is the position of the radar in the car coordinate system, θ r is the angle between the normal direction of the radar and the car coordinate system, v C, x , v C, y and w respectively represent the movement of the car. Transverse and longitudinal velocity components and rotational speed. When the point cloud velocity measured by radar is the same as A large difference indicates that the point cloud is moving relative to the ground. Therefore, the basis for dynamic and static separation of point clouds is:
其中,|·|表示求绝对值,S,D分别表示静态点云集合和动态点云集合,△vsd为速度门限,本发明中,速度门限△vsd取为经验值0.8m/s。动静点云分离结果如图5(b)所示,从图5中可以看到,除了表征运动目标的真实点云外,还存在许多杂乱的多径鬼影。Among them, |·| means finding the absolute value, S and D respectively represent the static point cloud set and the dynamic point cloud set, and △v sd is the speed threshold. In the present invention, the speed threshold △v sd is taken as the empirical value of 0.8m/s. The result of separation of moving and static point clouds is shown in Figure 5(b). It can be seen from Figure 5 that in addition to the real point clouds representing moving targets, there are also many messy multipath ghosts.
2.2动态点云聚类2.2 Dynamic point cloud clustering
运动汽车等道路目标对于车载毫米波雷达是扩展的。因此,需要对雷达动态点云D进行聚类处理。考虑到汽车为刚体运动目标,本发明采用基于速度改进的DBSCAN聚类方法,对动态点云进行聚类处理。具体地,在传统的DBSCAN算法的核心点的邻域搜索过程中,对一个核心点邻域内的点做进一步的速度检查,将速度与核心点速度差异较大的点从该核心点的邻域内排除,在本实施例中取△vDBSCAN值为1m/s,即将速度与核心点速度差值大于1m/s的点从该核心点的邻域内排除,其余处理过程同传统的DBSCAN算法。本发明中,DBSCAN算法的两大重要参数取值为:邻域搜索半径Eps=2m,密度(邻域内最少的点数)Minpts=1。此外,用于核心点速度检查的速度门限取为经验值△vDBSCAN=1m/s。记,一帧雷达数据经过聚类处理后的结果为CLU={clu0,clu1,clu2,...,cluNcluster},其中,Ncluster表示聚类后的簇个数,clui={pj,...},pj∈D表示第i簇。Road targets such as sports cars are extended for on-board millimeter wave radar. Therefore, the radar dynamic point cloud D needs to be clustered. Considering that the car is a rigid body moving target, the present invention adopts the DBSCAN clustering method based on speed improvement to perform clustering processing on the dynamic point cloud. Specifically, in the neighborhood search process of the core point of the traditional DBSCAN algorithm, further speed checks are performed on the points in the neighborhood of a core point, and points with a large speed difference from the speed of the core point are removed from the neighborhood of the core point. To exclude, in this embodiment, the Δv DBSCAN value is set to 1m/s, that is, points with a speed difference greater than 1m/s from the core point are excluded from the neighborhood of the core point. The rest of the processing is the same as the traditional DBSCAN algorithm. In the present invention, the two important parameters of the DBSCAN algorithm are: neighborhood search radius Eps=2m, and density (the minimum number of points in the neighborhood) Minpts=1. In addition, the speed threshold used for core point speed check is taken as the empirical value Δv DBSCAN =1m/s. Note that the result of a frame of radar data after clustering processing is CLU={clu 0 ,clu 1 ,clu 2 ,...,clu Ncluster }, where N cluster represents the number of clusters after clustering, clu i = {p j ,...},p j ∈D represents the i-th cluster.
2.3簇中心坐标以及多域信息获取2.3 Cluster center coordinates and multi-domain information acquisition
聚类处理后,提取每一个簇clui={pj,...}的中心坐标,以及表征该中心点的多域信息(速度、距离、RCS)。After the clustering process, the center coordinates of each cluster clu i ={p j ,...} are extracted, as well as the multi-domain information (speed, distance, RCS) characterizing the center point.
首先,提取簇的中心坐标。对于一般目标,我们取簇中的点云的几何中心作为该簇的中心坐标(xc,yc):First, extract the center coordinates of the cluster. For general targets, we take the geometric center of the point cloud in the cluster as the center coordinates of the cluster (x c , y c ):
其中,(xi,yi)表示簇中的点云的坐标,N表示该簇内点云的个数。特别地,对于卡车等大型目标,由于表征该类目标的点云在空间上跨度很长,点云的几何中心并不能很好的表征目标造成多径鬼影的部位。考虑到RCS越大的部位造成多径鬼影的概率越大,本发明基于点云的RCS信息,对内部点云跨度超过10m的簇进行中心坐标提取。具体的,首先确定第i簇中的点云的最大RCS计算值而后,搜索RCS计算值与/>的差值在10dBsm之内的点云集,记为clus′i,/>最后,提取clus′i内点云的几何中心坐标作为簇clui的中心坐标。Among them, ( xi ,y i ) represents the coordinates of the point cloud in the cluster, and N represents the number of point clouds in the cluster. In particular, for large targets such as trucks, since the point cloud representing such targets spans a long space, the geometric center of the point cloud cannot well represent the parts of the target that cause multipath ghosting. Considering that parts with larger RCS have a greater probability of causing multipath ghosts, this invention extracts the center coordinates of clusters whose internal point cloud span exceeds 10m based on the RCS information of the point cloud. Specifically, first determine the maximum RCS calculated value of the point cloud in the i-th cluster. Then, search for the RCS calculated value and/> The point cloud set whose difference is within 10dBsm is recorded as clus′ i ,/> Finally, the geometric center coordinates of the point cloud in clus′ i are extracted as the center coordinates of cluster clu i .
进一步地,提取表征簇中心的多域特征。本发明取簇中心坐标到雷达的欧式距离作为簇中心到雷达的距离rc,取簇内点云的平均速度作为簇中心的速度vc,取簇内点云的最大RCS值作为簇中心的RCS,记为RCSc。Furthermore, multi-domain features characterizing the cluster centers are extracted. This invention takes the Euclidean distance from the cluster center coordinates to the radar as the distance r c from the cluster center to the radar, takes the average speed of the point cloud in the cluster as the speed v c of the cluster center, and takes the maximum RCS value of the point cloud in the cluster as the cluster center. RCS, denoted as RCS c .
进一步的,记录簇内点云的个数Nup。记2.3步骤处理的结果为其中/> Further, record the number Nup of point clouds in the cluster. Record the result of step 2.3 as Among them/>
S3:目标关联,取场景内任意两运动目标进行关联;S3: Target association, associate any two moving targets in the scene;
从S2生成的簇中心中任取不同的两点进行关联,构成一个匹配集。记关联运算为AS(·),则有匹配集合Γ:Cluster centers generated from S2 Take any two different points and associate them to form a matching set. Let the association operation be AS(·), then there is a matching set Γ:
Γ=AS(S)Γ=AS(S)
为降低计算量,上述步骤可进一步优化。由多径特征可知,真实目标与多径鬼影相对雷达距离差值较小,故,可选取一定的距离差内的目标点进行关联,以提高算法实时性。本领域技术人员应注意,这里的一定的距离差,在远距雷达、中距离雷达、近距雷达条件下可分别取值为5m、10m和15m。To reduce the amount of calculation, the above steps can be further optimized. It can be seen from the multipath characteristics that the relative radar distance difference between the real target and the multipath ghost is small. Therefore, target points within a certain distance difference can be selected for correlation to improve the real-time performance of the algorithm. Those skilled in the art should note that the certain distance difference here can take values of 5m, 10m and 15m respectively under the conditions of long-range radar, medium-range radar and short-range radar.
S4:反射点反演求解;S4: Reflection point inversion solution;
取匹配集Γ中一对关联点Γ0={cA,cB},Γ0∈Γ。其中,cA={xa,ya,ra,va,RCSa,Nupa},cB={xb,yb,rb,vb,RCSb,Nupb}。由于多径目标到雷达距离总是比真实目标到雷达的距离更长,故,可以做进一步假设ra<rb,即是说,cA为真实目标,cB为多径鬼影。设雷达所在位置为o(xo,yo),求解三种假设下的反射点。Take a pair of associated points Γ 0 = {c A ,c B }, Γ 0 ∈ Γ in the matching set Γ. Among them, c A = {x a , y a , r a , v a , RCS a , Nup a }, c B = {x b , y b , r b , v b , RCS b , Nup b }. Since the distance from the multipath target to the radar is always longer than the distance from the real target to the radar, a further assumption can be made that r a <r b , that is, c A is the real target and c B is the multipath ghost. Assume the location of the radar is o(x o , y o ), and solve the reflection points under the three assumptions.
4.1真实鬼影与Type-2型三次反射多径鬼影GS假设下的反射点求解4.1 Solution of reflection points under the assumption of real ghost and Type-2 three-reflection multipath ghost G S
如图3(a)所示,多径鬼影GS与真实目标关于反射点所在的反射面呈镜像对称。由此,真实目标T与多径鬼影GS连接线段的中垂线与多径鬼影GS与雷达O的连线的交点,恰好是反射点的位置。假设cA为真实目标T,cB为多径鬼影GS,反射点求解的具体过程为:As shown in Figure 3(a), the multipath ghost GS and the real target are mirror symmetrical with respect to the reflection surface where the reflection point is located. Therefore, the intersection point of the mid-perpendicular line connecting the real target T and the multipath ghost GS and the line connecting the multipath ghost GS and the radar O is exactly the position of the reflection point. Assume that c A is the real target T and c B is the multipath ghost G S. The specific process of solving the reflection point is:
首先计算直线 First calculate the straight line
进一步地,计算线段的中垂线lmid:Amidx+Bmidy+Cmid=0:Further, calculate the line segment The mid-perpendicular line l mid :A mid x+B mid y+C mid =0:
最终,中垂线lmid与直线的交点(xbm,ybm),也即是该假设下反射点PTS的理论坐标为:Finally, the mid-perpendicular line l mid is the same as the straight line The intersection point (x bm ,y bm ), that is, the theoretical coordinates of the reflection point P TS under this assumption are:
如图3所示,线段等效于图3(a)中T、GS所在的线段;As shown in Figure 3, the line segment Equivalent to the line segment where T and G S are located in Figure 3(a);
由于目标测量精度、信号噪声等因素的影响,目标的定位存在误差,进而影响反射点求解的精度,故,我们设置反射点搜索区域为以反射点PTS为圆心,半径为△R的区域,记为本实施例中△R的取值为1m.Due to the influence of target measurement accuracy, signal noise and other factors, there are errors in the positioning of the target, which in turn affects the accuracy of the reflection point solution. Therefore, we set the reflection point search area to be an area with the reflection point P TS as the center and a radius of △R, recorded as In this embodiment, the value of △R is 1m.
4.2真实鬼影与Type-2型二次反射多径鬼影GM2假设下的反射点求解4.2 Solution of reflection points under the assumption of real ghost and Type-2 secondary reflection multipath ghost G M2
如图3(b)所示,多径鬼影GM2为线段T'GS的中点,故,可以先求出T'坐标,再根据中点坐标公式确定二阶多径鬼影S,进而再利用步骤4.1中所述方法确定反射点。假设cA为真实目标T,cB为多径鬼影GM2,反射点求解的具体过程为:As shown in Figure 3(b), the multipath ghost G M2 is the midpoint of the line segment T'G S. Therefore, the coordinates of T' can be obtained first, and then the second-order multipath ghost S can be determined according to the midpoint coordinate formula. Then use the method described in step 4.1 to determine the reflection point. Assume that c A is the real target T and c B is the multipath ghost G M2 . The specific process of solving the reflection point is:
首先计算线段和线段/>的夹角为:First calculate the line segment and line segments/> The included angle is:
进一步的,将点cA绕点o旋转θ∠AOB可得到交点cA′,也即是图3(b)中的T',计算如下:Further, by rotating point c A around point o by θ ∠AOB , the intersection point c A′ can be obtained, which is T' in Figure 3(b). The calculation is as follows:
进一步地,根据中点公式可计算出多径鬼影GS坐标:Furthermore, the multipath ghost G S coordinates can be calculated according to the midpoint formula:
于此,再利用步骤4.1的方法,计算出反射点同理,本发明设置反射点搜索区域为以反射点/>为圆心,半径为△R的区域,记为/> Here, use the method in step 4.1 to calculate the reflection point In the same way, the present invention sets the reflection point search area to reflect points/> is the center of the circle and the area with radius △R, denoted as/>
4.3真实鬼影与Type-1型二次反射多径鬼影GM1假设下的反射点求解4.3 Solution of reflection points under the assumption of real ghost and Type-1 secondary reflection multipath ghost G M1
如图3(c)所示,多径鬼影GM1与真实目标两点存在两点特征,1)鬼影GM1、真实目标T以及雷达坐标三点共线;2)造成此多径鬼影GM1的反射点,落在以点O,点T为焦点,点到两焦距的距离和为|OT|+|PT|的椭圆轨迹上。假设cA为真实目标T,cB为多径鬼影GM2,反射点求解的具体过程为:As shown in Figure 3(c), the multipath ghost G M1 and the real target have two characteristics. 1) The three points of the ghost G M1 , the real target T and the radar coordinates are collinear; 2) The multipath ghost is caused by The reflection point of shadow G M1 falls on an elliptical trajectory with point O and point T as the focus, and the sum of the distances from the two focal lengths is |OT|+|PT|. Assume that c A is the real target T and c B is the multipath ghost G M2 . The specific process of solving the reflection point is:
首先获取计算线段和线段/>的夹角θ∠AOB。由于雷达测角精度以及噪声等因素的影响,多径鬼影GM1、真实目标T以及反射点P之间的在目标检测后并不一定共线。故我们设置一个合理的角度偏差门限,将两线段夹角θ∠AOB低于门限值△εθ时,即:First get the calculated line segment and line segments/> The angle θ ∠AOB . Due to the influence of radar angle measurement accuracy, noise and other factors, the multipath ghost G M1 , the real target T and the reflection point P are not necessarily collinear after target detection. Therefore, we set a reasonable angle deviation threshold, and when the angle θ ∠AOB between the two line segments is lower than the threshold value Δε θ , that is:
|θ∠AOB|≤△εθ |θ ∠AOB |≤△ε θ
认为o,cA,cB近似为三点共线。本实施例中△εθ取值为1.5°。It is considered that o, c A and c B are approximately three points collinear. In this embodiment, the value of Δε θ is 1.5°.
进一步地,求上述椭圆轨迹。椭圆的长半轴aE、短半轴bE、半焦距cE参数分别为:Further, find the above elliptical trajectory. The parameters of the ellipse’s semi-major axis a E , semi-minor axis b E , and half-focal length c E are respectively:
进一步的,椭圆与x轴正方向夹角为:Furthermore, the angle between the ellipse and the positive direction of the x-axis is:
最终,椭圆轨迹可表示为:Finally, the elliptical trajectory can be expressed as:
其中,(xcc,ycc)为线段的中点。Among them, (x cc ,y cc ) is a line segment the midpoint.
获取椭圆轨迹E后,同理,我们设置反射点搜索区域为以椭圆轨迹为轴心,缩放因子为scale1,scale2的椭圆环状区域,记为本实施例中缩放因子scale1,scale2取值分别为0.87、1.15。After obtaining the elliptical trajectory E, in the same way, we set the reflection point search area to be an elliptical annular area with the elliptical trajectory as the axis and the scaling factors as scale 1 and scale 2 , recorded as In this embodiment, the scaling factors scale 1 and scale 2 take values of 0.87 and 1.15 respectively.
当|θ∠AOB|>△εθ时,则排除该种假设,令 When |θ ∠AOB |>△ε θ , then this assumption is eliminated, let
S5:基于S4中求解的搜索域,在动静点云中搜索反射点,验证S4中的假设。S5: Based on the search domain solved in S4, search for reflection points in the dynamic and static point clouds to verify the hypothesis in S4.
通过S4的反射点区域的求解后,我们可以得到三个搜索区域于是,在上述区域内,在S1步骤中获取的动静点云集P搜索是否存在点云。若存在,则通过上述几何约束初步说明,Γ0={cA,cB}是一个由落在搜索域/>内的其他强RCS目标引起的真实目标与多径鬼影关联点对,且有,cA为真实目标,cB为多径鬼影。After solving the reflection point area of S4, we can get three search areas Therefore, in the above-mentioned area, the moving and static point cloud set P obtained in step S1 is searched for whether there is a point cloud. If it exists, it can be preliminarily explained through the above geometric constraints that Γ 0 = {c A , c B } is a result that falls in the search domain/> There are correlation point pairs between real targets and multipath ghosts caused by other strong RCS targets, and there are, c A is the real target, and c B is the multipath ghost.
S6:基于多域信息(距离、速度、RCS和点云个数)的多径鬼影判决。S6: Multipath ghost judgment based on multi-domain information (distance, speed, RCS and point cloud number).
为了提高多径鬼影识别的准确性,避免误识别,使用多域特征来进一步判决多径鬼影。由步骤S4已知,cA={xa,ya,ra,va,RCSa,Nupa}和cB={xb,yb,rb,vb,RCSb,Nupb}。判决器参考以下3个因子:In order to improve the accuracy of multipath ghost recognition and avoid misidentification, multi-domain features are used to further determine multipath ghosts. It is known from step S4 that c A ={x a ,ya ,ra ,va ,RCS a , Nup a } and c B ={x b ,y b ,r b , v b ,RCS b ,Nup b }. The decider refers to the following 3 factors:
因子δ1:速度,多径鬼影速度与真实目标速度存在特殊关系。Factor δ 1 : Speed. There is a special relationship between the multipath ghost speed and the real target speed.
其中,εv分别为速度门限,本实施例中取为经验值1m/s。Among them, ε v are speed thresholds respectively, which are taken as the empirical value 1m/s in this embodiment.
因子δ2:RCS,电磁波的反射将会使回波能量衰减,故真实目标的RCS计算值将大于多径鬼影的RCS计算值。Factor δ 2 : RCS. The reflection of electromagnetic waves will attenuate the echo energy, so the calculated RCS value of the real target will be greater than the calculated RCS value of the multipath ghost.
因子δ3:点云个数,表征真实目标的点云个数不低于由该真实目标对应的多径鬼影的点云个数。本实施例中,目标对应的点云个数即是表征该目标的簇中心对应簇内的点云个数。Factor δ 3 : The number of point clouds. The number of point clouds representing the real target is not less than the number of point clouds of multipath ghosts corresponding to the real target. In this embodiment, the number of point clouds corresponding to the target is the number of point clouds in the cluster corresponding to the cluster center representing the target.
综上,多径目标判决器可记为:In summary, the multipath target determiner can be recorded as:
当的计算结果等于3,则判定cB为多径鬼影;当/>的计算结果小于3,则判定cB不为多径鬼影;when The calculation result of is equal to 3, then it is determined that c B is a multipath ghost; when/> The calculation result is less than 3, then it is determined that c B is not a multipath ghost;
若某一簇中心被识别为多径鬼影,则其对应的簇内的点云均标记为多径鬼影。If a cluster center is identified as a multipath ghost, the point clouds in the corresponding cluster are all marked as multipath ghosts.
S7:遍历S3步骤获取的匹配集中的每一对关联对,重复步骤S4~S6,完成当前帧数据下的所有点云的多径识别任务。S7: Traverse each associated pair in the matching set obtained in step S3, repeat steps S4 to S6, and complete the multipath identification task of all point clouds under the current frame data.
最终的识别结果如图5(c)所示,标记为方框的雷达点云被识别为多径鬼影,从图5(c)中可以看到,通过所提发明技术处理后,多径鬼影能被有效识别,消除多径鬼影后的结果如图5(d)所示。实测实验验证了本发明的可行性与有效性。The final recognition result is shown in Figure 5(c). The radar point cloud marked as a box is identified as a multipath ghost. As can be seen from Figure 5(c), after being processed by the proposed inventive technology, the multipath Ghosts can be effectively identified, and the results after eliminating multipath ghosts are shown in Figure 5(d). Actual measurement experiments have verified the feasibility and effectiveness of the present invention.
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。Those of ordinary skill in the art will appreciate that the embodiments described here are provided to help readers understand the principles of the present invention, and it should be understood that the scope of the present invention is not limited to such specific statements and embodiments. Various modifications and variations of the present invention may occur to those skilled in the art. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention shall be included in the scope of the claims of the present invention.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311164769.2A CN117192533A (en) | 2023-09-11 | 2023-09-11 | Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311164769.2A CN117192533A (en) | 2023-09-11 | 2023-09-11 | Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117192533A true CN117192533A (en) | 2023-12-08 |
Family
ID=89001116
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311164769.2A Pending CN117192533A (en) | 2023-09-11 | 2023-09-11 | Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117192533A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN119337210A (en) * | 2024-10-16 | 2025-01-21 | 复睿智行智能科技(上海)有限公司 | A millimeter wave radar ghost suppression method and model based on machine learning |
WO2025025648A1 (en) * | 2023-08-03 | 2025-02-06 | 惠州市德赛西威智能交通技术研究院有限公司 | Multipath identification method and apparatus, vehicle, and medium |
-
2023
- 2023-09-11 CN CN202311164769.2A patent/CN117192533A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2025025648A1 (en) * | 2023-08-03 | 2025-02-06 | 惠州市德赛西威智能交通技术研究院有限公司 | Multipath identification method and apparatus, vehicle, and medium |
CN119337210A (en) * | 2024-10-16 | 2025-01-21 | 复睿智行智能科技(上海)有限公司 | A millimeter wave radar ghost suppression method and model based on machine learning |
CN119337210B (en) * | 2024-10-16 | 2025-05-27 | 复睿智行智能科技(上海)有限公司 | Millimeter wave radar ghost suppression method and model based on machine learning |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110427986B (en) | Target classification method of kernel support vector machine based on millimeter wave radar point cloud characteristics | |
CN117192533A (en) | Multipath ghost recognition method for vehicle-mounted millimeter wave radar based on reflection point inversion search | |
CN108761418B (en) | Radar multipath target eliminating method based on micro-Doppler characteristics | |
CN112240999B (en) | Method and system for processing radar reflections | |
CN110018461B (en) | Group target identification method based on high-resolution range profile and monopulse angle measurement | |
CN115079121B (en) | Multipath target identification method and device for vehicle-mounted radar and storage medium | |
CN112731307B (en) | RATM-CFAR detector based on distance-angle joint estimation and detection method | |
CN107340514A (en) | Hypersonic weak signal target RAE HT TBD integration detection methods in three dimensions | |
CN109031190B (en) | A passive time difference positioning method for high repetition frequency pulse signal | |
CN114966589B (en) | A multi-target detection method based on millimeter wave radar | |
CN115061113B (en) | Target detection model training method and device for radar and storage medium | |
WO2023108930A1 (en) | Point cloud speed-based millimeter wave radar big car identification method | |
CN111999726B (en) | Personnel positioning method based on millimeter wave radar | |
CN116990758A (en) | Multipath identification method, device, vehicle and medium | |
CN118011392A (en) | Point cloud high-resolution imaging method based on 4D millimeter wave radar | |
CN110967677B (en) | Time-frequency domain hierarchical multi-target distinguishing and measuring method | |
CN111707996B (en) | GEO satellite-borne SAR moving target detection method based on improved GRFT-STAP | |
CN113253251A (en) | FDA-MIMO radar detection method and system when target speed is unknown | |
CN102879774A (en) | Method and apparatus for synthesizing short flight paths | |
CN115144847A (en) | Parameter space multi-channel target search technology for weak target detection in sea clutter | |
CN113447901B (en) | Sea clutter identification and target detection method independent of probability model | |
CN113390406A (en) | Multi-target data association and positioning method based on passive multi-sensor system | |
CN115390028A (en) | Navigation radar multiple reflection false target identification technology based on arithmetic progression extraction | |
CN114839615A (en) | Target course angle fitting method for 4D millimeter wave radar and storage medium | |
CN113534066A (en) | Method and system for rejecting multi-reflection wild values of landing measurement radar in height direction |
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 |