[go: up one dir, main page]

CN114779235A - Road boundary detection method and system based on vehicle-mounted millimeter-wave radar - Google Patents

Road boundary detection method and system based on vehicle-mounted millimeter-wave radar Download PDF

Info

Publication number
CN114779235A
CN114779235A CN202210072589.0A CN202210072589A CN114779235A CN 114779235 A CN114779235 A CN 114779235A CN 202210072589 A CN202210072589 A CN 202210072589A CN 114779235 A CN114779235 A CN 114779235A
Authority
CN
China
Prior art keywords
boundary
point
starting point
points
vehicle
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
Application number
CN202210072589.0A
Other languages
Chinese (zh)
Inventor
杜红灯
石志轩
刘鹏耀
顾翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Hanrun Automobile Electronics Co ltd
Original Assignee
Jiangsu Hanrun Automobile Electronics Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu Hanrun Automobile Electronics Co ltd filed Critical Jiangsu Hanrun Automobile Electronics Co ltd
Priority to CN202210072589.0A priority Critical patent/CN114779235A/en
Publication of CN114779235A publication Critical patent/CN114779235A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

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

The invention discloses a road boundary detection method and a system based on a vehicle-mounted millimeter wave radar.

Description

基于车载毫米波雷达的道路边界检测方法和系统Road boundary detection method and system based on vehicle millimeter wave radar

技术领域technical field

本发明涉及道路检测技术领域,尤其涉及一种基于车载毫米波雷达的道路边界检测方法和系统。The invention relates to the technical field of road detection, in particular to a road boundary detection method and system based on a vehicle-mounted millimeter-wave radar.

背景技术Background technique

传统道路边界检测一般是利用激光雷达来实现,但是激光雷达价格昂贵,无法得到普及,而车载毫米波雷达价格低廉,并且对目标有较好的测距测速能力,对雨雾有较好的穿透能力,以及不被光照强度影响,成为了智能驾驶方案中不可替代的传感器选择。Traditional road boundary detection is generally realized by using lidar, but lidar is expensive and cannot be popularized, while vehicle-mounted millimeter-wave radar is cheap, has good ranging and speed measurement capabilities for targets, and has good penetration of rain and fog. Capability, and not being affected by light intensity, has become an irreplaceable sensor choice in intelligent driving solutions.

由于高速道路与城市道路上存在车辆等干扰,利用现有检测方法根据车载毫米波雷达采集的雷达数据获取的道路边界线输出往往在时域上不稳定,从而导致检测精确度较低。Due to the interference of vehicles on highways and urban roads, the output of road boundary lines obtained by existing detection methods based on radar data collected by vehicle-mounted millimeter-wave radars is often unstable in the time domain, resulting in low detection accuracy.

发明内容SUMMARY OF THE INVENTION

为了解决上述技术问题,本发明提供了一种基于车载毫米波雷达的道路边界检测方法和系统,其解决了现有技术中根据车载毫米波雷达采集的雷达数据获取的道路边界线输出在时域上不稳定,检测精确度较低的问题。In order to solve the above technical problems, the present invention provides a road boundary detection method and system based on a vehicle-mounted millimeter-wave radar, which solves the problem in the prior art that the output of the road boundary line obtained according to the radar data collected by the vehicle-mounted millimeter-wave radar is in the time domain. It is unstable and the detection accuracy is low.

一方面,本发明提供了一种基于车载毫米波雷达的道路边界检测方法,包括:In one aspect, the present invention provides a road boundary detection method based on a vehicle-mounted millimeter-wave radar, comprising:

获取车载毫米波雷达输出的当前帧所有的静态目标点;Obtain all static target points of the current frame output by the vehicle-mounted millimeter-wave radar;

根据初始筛选范围从所有的所述静态坐标点中筛选出多个稳定静态目标点,所述初始筛选范围根据雷达探测能力确定;Screen a plurality of stable static target points from all the static coordinate points according to the initial screening range, and the initial screening range is determined according to the radar detection capability;

将位于预设筛选范围内的所述稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为所述预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行;The stable static target point located in the preset screening range is used as the starting point of the alternative boundary; wherein, taking the central axis of the radar as the symmetry axis, two search ranges are set within the field of view of the vehicle-mounted millimeter-wave radar, and the The two set search ranges are used as the preset screening ranges, and each search range takes the sector radius boundary corresponding to the field of view range of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other boundary is parallel to the sector radius boundary;

对所述备选边界起始点按照与雷达原点之间距离的由小到大进行排序,将距离最小的所述备选边界起始点作为边界起始点;Sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point;

以所述边界起始点为起始点,在预设的查找范围内从所述稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止;Take the boundary starting point as the starting point, search for the boundary point from the stable static target point within the preset search range, and use the found boundary point as the new starting point to continue searching for the next boundary point until the search fails. until the next boundary point;

利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数;Use the least squares method to perform quadratic curve fitting on all the obtained boundary points to obtain the quadratic curve coefficient of the current frame;

利用卡尔曼滤波算法对所述当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。The quadratic curve coefficients of the current frame and the quadratic curve coefficients of the previous frame are filtered by using the Kalman filtering algorithm, so as to obtain the target quadratic curve coefficients and the target road boundary curve.

可选地,还包括:如果在预设的查找范围内查找到多个边界点,则求取每个边界点与雷达原点之间的距离;Optionally, it also includes: if a plurality of boundary points are found within the preset search range, calculating the distance between each boundary point and the origin of the radar;

将所有边界点按照求得的距离由小到大的顺序进行排序,并将距离最小值对应的边界点作为新的起始点。Sort all the boundary points in ascending order of the obtained distance, and use the boundary point corresponding to the minimum distance as the new starting point.

可选地,边界点与雷达原点之间的距离计算公式表示为:Optionally, the distance calculation formula between the boundary point and the radar origin is expressed as:

d2=αPx 2+βPy 2---(1)d 2 =αP x 2 +βP y 2 ---(1)

式中:d为计算得到的边界点与雷达原点之间的距离,α为横向的比例因子,β为纵向的比例因子,Px为边界点的横坐标,Py为边界点的纵坐标。In the formula: d is the distance between the calculated boundary point and the origin of the radar, α is the horizontal scale factor, β is the vertical scale factor, P x is the abscissa of the boundary point, and P y is the ordinate of the boundary point.

可选地,还包括:Optionally, also include:

如果查找得到的所有边界点满足设定的数量阈值和长度阈值,则转入利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数的步骤;If all the boundary points obtained by the search meet the set quantity threshold and length threshold, then transfer to the step of using the least squares method to perform quadratic curve fitting on all the boundary points obtained by the search to obtain the quadratic curve coefficient of the current frame;

否则,去除当前的所述边界起始点,并将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点。Otherwise, the current boundary starting point is removed, and the candidate boundary starting point located next to the current boundary starting point in the sorting is used as a new boundary starting point.

可选地,所述去除当前的所述边界起始点,并将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点的步骤包括:Optionally, the step of removing the current boundary starting point and using the candidate boundary starting point that is located next to the current boundary starting point in the sorting as a new boundary starting point includes:

统计剩余的所述备选边界起始点的数量;count the number of the remaining candidate boundary starting points;

如果所述剩余的所述备选边界起始点的数量大于零,则将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点。If the number of the remaining candidate boundary starting points is greater than zero, the candidate boundary starting point located next to the current boundary starting point in the sorting is used as a new boundary starting point.

可选地,在雷达坐标系下,通过设置横向阈值与纵向阈值得到所述初始筛选范围;其中,所述横向阈值根据真实场景下道路宽度、雷达探测能力得到;所述纵向阈值根据雷达探测能力确定。Optionally, in the radar coordinate system, the initial screening range is obtained by setting a horizontal threshold and a vertical threshold; wherein, the horizontal threshold is obtained according to the road width and radar detection capability in a real scene; the vertical threshold is obtained according to the radar detection capability. Sure.

可选地,利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数,包括:Optionally, use the least squares method to perform quadratic curve fitting on all the obtained boundary points to obtain the quadratic curve coefficients of the current frame, including:

获取上一帧查找得到的所有边界点;Get all the boundary points found in the previous frame;

将当前帧查找得到的所有边界点和上一帧查找得到的所有边界点进行融合;Fusion of all boundary points found in the current frame and all boundary points found in the previous frame;

利用最小二乘法对融合得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数。The quadratic curve fitting is performed on all the fused boundary points using the least squares method to obtain the quadratic curve coefficients of the current frame.

可选地,起始点的纵坐标越大,对应的预设的查找范围的横向宽度和纵向宽度越大。Optionally, the larger the vertical coordinate of the starting point, the larger the horizontal width and vertical width of the corresponding preset search range.

另一方面,本发明还提供了一种基于车载毫米波雷达的道路边界检测系统,包括:On the other hand, the present invention also provides a road boundary detection system based on a vehicle-mounted millimeter-wave radar, comprising:

目标点获取模块,被配置为获取车载毫米波雷达输出的当前帧所有的静态目标点;The target point acquisition module is configured to acquire all static target points of the current frame output by the vehicle-mounted millimeter wave radar;

初步筛选模块,被配置为根据初始筛选范围从所有的所述静态坐标点中筛选出多个稳定静态目标点,所述初始筛选范围根据雷达探测能力确定;a preliminary screening module, configured to screen out a plurality of stable static target points from all the static coordinate points according to an initial screening range, the initial screening range being determined according to the radar detection capability;

备选起始点获取模块,被配置为将位于预设筛选范围内的所述稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为所述预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行;The alternative starting point acquisition module is configured to use the stable static target point located within the preset screening range as the alternative boundary starting point; wherein, taking the radar center axis as the axis of symmetry, at the field of view of the vehicle-mounted millimeter-wave radar Two search ranges are set within the range, and the two set search ranges are used as the preset screening ranges, and each search range takes the sector radius boundary corresponding to the field of view angle range of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other a boundary parallel to the sector radius boundary;

初始起始点确定模块,被配置为对所述备选边界起始点按照与雷达原点的之间距离由小到大进行排序,将距离最小的所述备选边界起始点作为边界起始点;The initial starting point determination module is configured to sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point;

查找模块,被配置为以所述边界起始点为起始点,在预设的查找范围内从所述稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止;The search module is configured to take the boundary starting point as a starting point, search for a boundary point from the stable static target point within a preset search range, and use the found boundary point as a new starting point to continue searching for the next Boundary point until the next boundary point cannot be found;

曲线拟合模块,被配置为利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数;a curve fitting module, configured to perform quadratic curve fitting on all boundary points obtained by the least squares method to obtain quadratic curve coefficients of the current frame;

以及,滤波处理模块,被配置为利用卡尔曼滤波算法对所述当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。And, the filtering processing module is configured to perform filtering processing on the quadratic curve coefficient of the current frame and the quadratic curve coefficient of the previous frame by using the Kalman filtering algorithm, so as to obtain the target quadratic curve coefficient and the target road boundary curve.

再一方面,本发明还提供了一种电子设备,包括:处理器和存储器,所述存储器上存储有计算机可读指令,所述计算机可读指令被所述处理器执行时实现如上述所述的基于车载毫米波雷达的道路边界检测方法。In another aspect, the present invention also provides an electronic device, comprising: a processor and a memory, where computer-readable instructions are stored in the memory, and the computer-readable instructions are implemented as described above when executed by the processor A road boundary detection method based on vehicle-mounted millimeter-wave radar.

又一方面,本发明还提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如上述所述的基于车载毫米波雷达的道路边界检测方法。In another aspect, the present invention also provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, implements the above-mentioned road boundary detection method based on a vehicle-mounted millimeter-wave radar.

本发明的基于车载毫米波雷达的道路边界检测方法首先将输入的车载毫米波雷达静态目标点进行初步筛选,获得比较稳定的静态目标点,通过设置查找范围利用首个起始点进行聚类找出所有边界点的策略在稳定的静态目标点中筛选所有可能的道路边界点,然后,结合雷达上一帧筛选得到的边界点利用最小二乘法进行二次曲线拟合,最后,考虑道路边界不会短时间发生跳变,利用卡尔曼滤波算法对得到的二次曲线中三个系数进行滤波处理,得到时域中较为稳定的曲线输出,进而提高了道路边界检测精确度。The road boundary detection method based on the vehicle-mounted millimeter-wave radar of the present invention first performs preliminary screening of the input vehicle-mounted millimeter-wave radar static target points to obtain relatively stable static target points, and uses the first starting point to cluster to find out by setting the search range. The strategy of all boundary points selects all possible road boundary points in the stable static target points, and then uses the least squares method to fit the quadratic curve combined with the boundary points screened in the previous frame of the radar. Finally, considering that the road boundary will not If a jump occurs in a short time, the Kalman filter algorithm is used to filter the three coefficients in the obtained quadratic curve, and a relatively stable curve output in the time domain is obtained, thereby improving the detection accuracy of road boundaries.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings required in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some of the present invention. In the embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without any creative effort.

图1为本发明的基于车载毫米波雷达的道路边界检测方法的一些实施例的流程示意图;1 is a schematic flowchart of some embodiments of a road boundary detection method based on a vehicle-mounted millimeter-wave radar of the present invention;

图2是本发明的基于车载毫米波雷达的道路边界检测方法的另一些实施例的流程示意图;2 is a schematic flowchart of other embodiments of a road boundary detection method based on a vehicle-mounted millimeter-wave radar of the present invention;

图3是本发明的基于车载毫米波雷达的道路边界检测方法的再一些实施例的流程示意图;FIG. 3 is a schematic flowchart of some further embodiments of a road boundary detection method based on a vehicle-mounted millimeter-wave radar of the present invention;

图4是本发明的基于车载毫米波雷达的道路边界检测方法的道路边界点跟踪示意图;Fig. 4 is the road boundary point tracking schematic diagram of the road boundary detection method based on the vehicle-mounted millimeter wave radar of the present invention;

图5是本发明的基于车载毫米波雷达的道路边界检测方法的卡尔曼滤波算法流程图;Fig. 5 is the Kalman filter algorithm flow chart of the road boundary detection method based on vehicle-mounted millimeter wave radar of the present invention;

图6为本发明的基于车载毫米波雷达的道路边界检测系统的一些实施例的结构框图。FIG. 6 is a structural block diagram of some embodiments of the road boundary detection system based on the vehicle-mounted millimeter wave radar of the present invention.

具体实施方式Detailed ways

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art fall within the protection scope of the present invention.

一方面,参见图1所示,本发明实施例提供了一种基于车载毫米波雷达的道路边界检测方法,包括:On the one hand, as shown in FIG. 1 , an embodiment of the present invention provides a road boundary detection method based on a vehicle-mounted millimeter-wave radar, including:

步骤100:获取车载毫米波雷达输出的当前帧所有的静态目标点。Step 100: Acquire all static target points of the current frame output by the vehicle-mounted millimeter-wave radar.

本步骤中车载毫米波雷达可以为具有直接输出静态坐标点功能的毫米波雷达;当然也可以为不具有该功能的车载毫米波雷达,此时需要首先获取车载毫米波雷达输出的当前帧雷达回波信号,再从当前帧雷达回波信号中分析得到当前帧所有的静态目标点,该静态目标点即静态目标点的位置数据/坐标数据。In this step, the vehicle-mounted millimeter-wave radar can be a millimeter-wave radar with the function of directly outputting static coordinate points; of course, it can also be a vehicle-mounted millimeter-wave radar without this function. In this case, it is necessary to first obtain the current frame radar response output by the vehicle-mounted millimeter-wave radar. wave signal, and then analyze all the static target points of the current frame from the radar echo signal of the current frame, and the static target point is the position data/coordinate data of the static target point.

步骤200:根据初始筛选范围从所有的静态坐标点中筛选出多个稳定静态目标点,初始筛选范围根据雷达探测能力确定。Step 200: Screen out a plurality of stable static target points from all the static coordinate points according to the initial screening range, and the initial screening range is determined according to the radar detection capability.

本步骤中直接将位于初始筛选范围的静态目标点筛选出来就得到多个稳定静态目标点;其中初始筛选范围根据雷达探测能力来确定。例如:车载毫米波雷达在L米范围内探测的静态目标点的数据可靠度高,则将位于L米范围内探测的静态目标点筛选出来作为稳定静态目标点。这里初始筛选范围可以是矩形,也可以是三角形,当为矩形时,将初始筛选范围的纵向距离定为L米,将初始筛选范围的横向距离定为预设宽度,将预设宽度范围内的静态目标点筛选出来即可,无需过宽;也可以为三角形,只设定纵向距离阈值。一般地,预设宽度能够覆盖真实场景下全部或者绝大部分带护栏的道路即可。In this step, the static target points located in the initial screening range are directly screened to obtain a plurality of stable static target points; wherein the initial screening range is determined according to the radar detection capability. For example, if the data reliability of the static target points detected by the vehicle-mounted millimeter-wave radar within the range of L meters is high, the static target points detected within the range of L meters are screened out as stable static target points. Here, the initial screening range can be a rectangle or a triangle. When it is a rectangle, the vertical distance of the initial screening range is set as L meters, the horizontal distance of the initial screening range is set as the preset width, and the The static target point can be screened out, and it does not need to be too wide; it can also be a triangle, and only the vertical distance threshold is set. Generally, the preset width can cover all or most of the roads with guardrails in the real scene.

步骤300:将位于预设筛选范围内的稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行。Step 300: Use the stable static target point within the preset screening range as the starting point of the alternative boundary; wherein, taking the radar center axis as the axis of symmetry, two search ranges are set within the field of view of the vehicle-mounted millimeter-wave radar, The two set search ranges are used as the preset screening ranges, and each search range takes the sector radius boundary corresponding to the field of view range of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other boundary is parallel to the sector radius boundary.

本步骤中,每个查找范围两个边界间的距离应能满足护栏等边界的检测需求,查找范围在扇形半径方向的延伸长度与各种不同道路的边界特征相关,一般为经验值,通常地,查找范围在扇形半径方向的延伸长度在雷达中心轴线方向的投影长度小于初始筛选范围在雷达中心轴线方向对应的长度。在实际应用中,在获得稳定静态目标点之后,可按照纵坐标由小到大排序,方便后续查找边界点。In this step, the distance between the two boundaries of each search range should meet the detection requirements of borders such as guardrails, and the extension length of the search range in the direction of the sector radius is related to the boundary features of various roads, and is generally an empirical value. , the projection length of the extension length of the search range in the direction of the sector radius in the direction of the radar center axis is smaller than the corresponding length of the initial screening range in the direction of the radar center axis. In practical applications, after obtaining stable static target points, it can be sorted according to the ordinate from small to large, which is convenient for subsequent search of boundary points.

步骤400:对备选边界起始点按照与雷达原点的之间距离由小到大进行排序,将距离最小的备选边界起始点作为边界起始点。Step 400: Sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point.

在进行边界点查找前,计算各备选边界起始点与雷达原点之间的距离,并将计算出的距离有小到大进行排序,将距离最小的备选边界起始点作为边界起始点。Before searching for boundary points, calculate the distance between each candidate boundary starting point and the radar origin, sort the calculated distances from small to large, and take the candidate boundary starting point with the smallest distance as the boundary starting point.

步骤500:以边界起始点为起始点,在预设的查找范围内从稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止。Step 500: Take the boundary starting point as the starting point, search for the boundary point from the stable static target points within the preset search range, and use the found boundary point as the new starting point to continue to search for the next boundary point until no more can be found. until the next boundary point.

本步骤中首先确定一个边界起始点为开始点,在预设的查找范围内查找下一个边界起始点,当找到下一个边界起始点时,以此边界起始点作为新的起始点,查找下一个边界点,也即将预设的查找范围内查找到的稳定静态目标点作为边界点,实现道路边界点的跟踪。需要说明的是,每个起始点对应一个预设的查找范围。In this step, first determine a boundary starting point as the starting point, search for the next boundary starting point within the preset search range, when the next boundary starting point is found, use this boundary starting point as the new starting point, and search for the next The boundary point, that is, the stable static target point found within the preset search range is used as the boundary point to realize the tracking of the road boundary point. It should be noted that each starting point corresponds to a preset search range.

步骤600:利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数。Step 600: Use the least squares method to perform quadratic curve fitting on all the obtained boundary points to obtain quadratic curve coefficients of the current frame.

本步骤中利用最小二乘法对步骤400得到的点进行二次曲线拟合,拟合式子如下式所示:In this step, the least squares method is used to perform quadratic curve fitting on the points obtained in step 400, and the fitting formula is as follows:

Px=aPy 2+bPy+c---(2)P x =aP y 2 +bP y +c---(2)

式中:a为拟合的函数的二次项系数,b为拟合的函数的一次项系数,c为拟合的函数的常数项。Px为边界点的横坐标,Py为边界点的纵坐标。In the formula: a is the quadratic term coefficient of the fitted function, b is the linear term coefficient of the fitted function, and c is the constant term of the fitted function. P x is the abscissa of the boundary point, and P y is the ordinate of the boundary point.

步骤700:利用卡尔曼滤波算法对当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。Step 700 : filter the quadratic curve coefficients of the current frame and the quadratic curve coefficients of the previous frame by using the Kalman filtering algorithm to obtain the target quadratic curve coefficients and the target road boundary curve.

本步骤中利用当前帧与上一帧得到的a,b,c三个系数利用卡尔曼滤波算法对其进行滤波处理,得到时域上比较稳定的系数输出。In this step, the three coefficients a, b, and c obtained from the current frame and the previous frame are used to filter and process the three coefficients of the current frame and the previous frame, and the Kalman filtering algorithm is used to obtain a relatively stable coefficient output in the time domain.

本发明实施例的基于车载毫米波雷达的道路边界检测方法首先将输入的车载毫米波雷达静态目标点进行初步筛选,获得比较稳定的静态目标点,通过设置查找范围利用首个起始点进行聚类找出所有边界点的策略在稳定的静态目标点中筛选所有可能的道路边界点,然后,结合雷达上一帧筛选得到的边界点利用最小二乘法进行二次曲线拟合,最后,考虑道路边界不会短时间发生跳变,利用卡尔曼滤波算法对得到的二次曲线中三个系数进行滤波处理,得到时域中较为稳定的曲线输出,进而提高了道路边界检测精确度。本发明实施例的基于车载毫米波雷达的道路边界检测方法可在无人驾驶、辅助驾驶中具有十分重要的工程应用。The road boundary detection method based on the vehicle-mounted millimeter-wave radar according to the embodiment of the present invention first performs preliminary screening of the input vehicle-mounted millimeter-wave radar static target points to obtain relatively stable static target points, and uses the first starting point to perform clustering by setting the search range The strategy of finding all boundary points is to screen all possible road boundary points in the stable static target points, and then, combined with the boundary points screened in the previous frame of the radar, use the least squares method to perform quadratic curve fitting, and finally, consider the road boundary There will be no jump in a short time, and the Kalman filter algorithm is used to filter the three coefficients in the obtained quadratic curve to obtain a relatively stable curve output in the time domain, thereby improving the accuracy of road boundary detection. The road boundary detection method based on the vehicle-mounted millimeter-wave radar according to the embodiment of the present invention can have very important engineering applications in unmanned driving and assisted driving.

在一些实施例中,本发明的基于车载毫米波雷达的道路边界检测方法还包括:In some embodiments, the road boundary detection method based on the vehicle-mounted millimeter-wave radar of the present invention further comprises:

如果在预设的查找范围内查找到多个边界点,则求取每个边界点与雷达原点之间的距离;If multiple boundary points are found within the preset search range, the distance between each boundary point and the origin of the radar is calculated;

将所有边界点按照求得的距离由小到大的顺序进行排序,并将距离最小值对应的边界点作为新的起始点。Sort all the boundary points in ascending order of the obtained distance, and use the boundary point corresponding to the minimum distance as the new starting point.

本实施例中通过预设的查找范围一般情况下会找到一个边界点,当在同一个预设的查找范围查找到多个边界点时通过计算每个边界点与雷达原点之间的距离d,将距离雷达原点最近的边界点作为新的起始点。另外,在实际梳理过程中也可以通过程序进行设定即找到一个就会停止搜索,然后再利用此点进行下一步的筛选。In this embodiment, a boundary point is generally found through the preset search range. When multiple boundary points are found in the same preset search range, the distance d between each boundary point and the origin of the radar is calculated, Take the boundary point closest to the radar origin as the new starting point. In addition, in the actual carding process, it can also be set through the program, that is, if one is found, the search will be stopped, and then this point will be used for the next step of screening.

在一些实施例中,本发明的基于车载毫米波雷达的道路边界检测方法中边界点与雷达原点之间的距离d计算公式表示为:In some embodiments, the calculation formula of the distance d between the boundary point and the origin of the radar in the road boundary detection method based on the vehicle-mounted millimeter wave radar of the present invention is expressed as:

d2=αPx 2+βPy 2---(1)d 2 =αP x 2 +βP y 2 ---(1)

式中:d为计算得到的边界点与雷达原点之间的距离,α为横向的比例因子,β为纵向的比例因子,α,β取值范围均为[0,1],并且α=1-β,α,β可以为经验值或标定值,Px为边界点的横坐标,Py为边界点的纵坐标。In the formula: d is the distance between the calculated boundary point and the origin of the radar, α is the horizontal scale factor, β is the vertical scale factor, the value ranges of α and β are both [0, 1], and α=1 -β, α, β can be empirical values or calibration values, P x is the abscissa of the boundary point, P y is the ordinate of the boundary point.

本实施例中考虑边界点的特点,在求取距离时计算公式加上了一个影响因子,可以利用α与β的值来调节Px与Py对d的影响。In this embodiment, considering the characteristics of boundary points, an influence factor is added to the calculation formula when calculating the distance, and the values of α and β can be used to adjust the influence of Px and Py on d.

在一些实施例中,参见图2所示,本发明的基于车载毫米波雷达的道路边界检测方法还包括:In some embodiments, as shown in FIG. 2 , the road boundary detection method based on the vehicle-mounted millimeter-wave radar of the present invention further includes:

步骤800:如果查找得到的所有边界点满足设定的数量阈值和长度阈值,则转入利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数的步骤;Step 800: If all the boundary points obtained by the search satisfy the set quantity threshold and length threshold, then proceed to the step of using the least squares method to perform quadratic curve fitting on all the boundary points obtained by the search to obtain the quadratic curve coefficient of the current frame. ;

步骤900:否则,去除当前的边界起始点,并将排序中位于当前的边界起始点的下一位的备选边界起始点作为新的边界起始点。Step 900: Otherwise, remove the current boundary starting point, and use the candidate boundary starting point that is located next to the current boundary starting point in the sorting as a new boundary starting point.

本实施例中对于步骤500中查找得到的所有边界点,需要进行有效性判断,即判断所有边界点是否满足设定的数量阈值和长度阈值,若不满足,则去除当前的边界起始点,若找到的边界点满足要求,则进行步骤600。需要说明的是,这里设定的数量阈值和长度阈值为经验值或标定值。在其他一些实施例中,因为采用的是二次曲线拟合,因此边界点数据最少为3个,也可结合毫米波雷达输出点云个数综合确定数量阈值,例如雷达分辨率较高,输出的点云个数较多,此阈值可以适量放大,长度阈值的设定与雷达分辨率以及实际场景相关,此处设置阈值为8m。在其他一些实施例中,也可以在查找得到的所有边界点满足设定的数量阈值时,转入步骤600。In this embodiment, for all the boundary points found in step 500, a validity judgment needs to be performed, that is, it is judged whether all the boundary points satisfy the set quantity threshold and length threshold. If not, the current boundary starting point is removed. If the found boundary points meet the requirements, step 600 is performed. It should be noted that the quantity threshold and length threshold set here are empirical values or calibration values. In some other embodiments, because quadratic curve fitting is used, the boundary point data is at least 3, and the number threshold can also be comprehensively determined in combination with the number of output point clouds of the millimeter wave radar. The number of point clouds is large, and this threshold can be appropriately enlarged. The setting of the length threshold is related to the radar resolution and the actual scene, and the threshold is set to 8m here. In some other embodiments, it is also possible to go to step 600 when all the boundary points obtained by searching meet the set quantity threshold.

在一些实施例中,参见图2所示,本发明的基于车载毫米波雷达的道路边界检测方法中步骤900包括:In some embodiments, as shown in FIG. 2 , step 900 in the road boundary detection method based on the vehicle-mounted millimeter-wave radar of the present invention includes:

步骤901:统计剩余的备选边界起始点的数量;Step 901: Count the number of remaining candidate boundary starting points;

步骤902:如果剩余的备选边界起始点的数量大于零,则将排序中位于当前的边界起始点的下一位的备选边界起始点作为新的边界起始点。Step 902: If the number of remaining candidate boundary starting points is greater than zero, use the candidate boundary starting point located next to the current boundary starting point in the sorting as a new boundary starting point.

本实施例中若查找得到的所有边界点不满足设定的数量阈值和长度阈值,则去除当前的边界起始点,判断步骤400中获得的剩余可能的备选边界起始点是否为空,如果不为空,则选择下一个最近的点作为边界起始点,回到步骤500,如果为空,则输出为空,之后直接转入步骤100,处理下一帧车载毫米波雷达输出的静态目标点。In this embodiment, if all the boundary points obtained by the search do not meet the set quantity threshold and length threshold, the current boundary starting point is removed, and it is determined whether the remaining possible candidate boundary starting points obtained in step 400 are empty. If it is empty, select the next closest point as the starting point of the boundary, and go back to step 500. If it is empty, the output is empty, and then directly go to step 100 to process the static target point output by the vehicle-mounted millimeter-wave radar in the next frame.

可选地,本发明实施例的基于车载毫米波雷达的道路边界检测方法中在雷达坐标系下,通过设置横向阈值与纵向阈值得到初始筛选范围。其中,横向阈值根据真实场景下道路宽度、雷达探测能力得到;纵向阈值根据雷达探测能力确定。也即筛选出雷达探测能力对应的静态目标点,以得到可靠的静态目标点作为稳定静态目标点。需要说明的是,本实施例中横向阈值应能覆盖真实场景下全部或者大部分带护栏等边界的道路,纵向阈值应能覆盖雷达能够准确探测的距离范围。Optionally, in the road boundary detection method based on the vehicle-mounted millimeter-wave radar in the embodiment of the present invention, in the radar coordinate system, the initial screening range is obtained by setting a horizontal threshold and a vertical threshold. Among them, the horizontal threshold is obtained according to the road width and radar detection ability in the real scene; the vertical threshold is determined according to the radar detection ability. That is, the static target points corresponding to the radar detection capability are screened out, so as to obtain reliable static target points as stable static target points. It should be noted that in this embodiment, the horizontal threshold should cover all or most of the roads with borders such as guardrails in the real scene, and the vertical threshold should cover the distance range that the radar can accurately detect.

在一些实施例中,参见图3所示,本发明的基于车载毫米波雷达的道路边界检测方法中步骤600包括:In some embodiments, referring to FIG. 3 , step 600 of the road boundary detection method based on the vehicle-mounted millimeter-wave radar of the present invention includes:

步骤601:获取上一帧查找得到的所有边界点;Step 601: Acquire all boundary points obtained by searching in the previous frame;

步骤602:将当前帧查找得到的所有边界点和上一帧查找得到的所有边界点进行融合;Step 602: fuse all the boundary points obtained by the current frame search with all the boundary points obtained by the previous frame search;

步骤603:利用最小二乘法对融合得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数。Step 603: Use the least squares method to perform quadratic curve fitting on all the fused boundary points to obtain quadratic curve coefficients of the current frame.

本实施例中每个静态目标点都有一个ID或标记,通过融合来去除ID或标记重复的静态目标点。In this embodiment, each static target point has an ID or a mark, and static target points with repeated IDs or marks are removed by fusion.

可选地,本发明实施例的基于车载毫米波雷达的道路边界检测方法中起始点的纵坐标越大,对应的预设的查找范围的横向宽度和纵向宽度越大。需要说明的是纵坐标与预设的查找范围的横向宽度和纵向宽度之间的对应关系可根据经验或标定获得。本实施例中预设的查找范围是一个阈值框,起始就是边界起始点,在其左右和向上设置一个阈值框,在框内查找下一点。Optionally, the larger the ordinate of the starting point in the road boundary detection method based on the vehicle-mounted millimeter-wave radar in the embodiment of the present invention, the larger the horizontal and vertical widths of the corresponding preset search range. It should be noted that the correspondence between the ordinate and the preset horizontal width and vertical width of the search range can be obtained based on experience or calibration. The preset search range in this embodiment is a threshold value box, and the starting point is the starting point of the boundary, and a threshold value box is set on the left and right and upward of the threshold value box, and the next point is searched in the box.

下面具体说明一下基于车载毫米波雷达的道路边界检测方法的实现流程,参见图4所示:The following is a detailed description of the implementation process of the road boundary detection method based on the vehicle-mounted millimeter-wave radar, as shown in Figure 4:

1)获取雷达当前帧所有的静态目标点;1) Obtain all static target points in the current frame of the radar;

2)在雷达坐标系下,通过设置横向与纵向的阈值,得到一个初始筛选范围,初步筛选所有输入的静态目标点,获得稳定的静态目标点;其中,横向与纵向的阈值根据真实场景下道路宽度、雷达探测能力等指标确定;如横向阈值应能覆盖真实场景下全部或者大部分带护栏等边界的道路,纵向阈值应能覆盖雷达能够准确探测的距离范围,获得稳定静态目标点之后对其按照纵坐标的大小进行排序。2) In the radar coordinate system, by setting the horizontal and vertical thresholds, an initial screening range is obtained, and all input static target points are preliminarily screened to obtain stable static target points; among them, the horizontal and vertical thresholds are based on the road in the real scene. Width, radar detection ability and other indicators are determined; for example, the horizontal threshold should cover all or most of the roads with borders such as guardrails in the real scene, and the vertical threshold should cover the distance range that the radar can accurately detect. Sort by the size of the ordinate.

3)在初始筛选范围内确定出位于雷达的视场角(FOV,Field of view)范围内的区域,以雷达中心轴线为对称轴在该区域的两侧分别设定一个查找范围,将查找范围内的静态目标点定义为道路的备选边界起始点,在备选边界起始点为非空的情形下,执行步骤4);其中,每个查找范围以视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行,两个边界间的距离应能满足护栏等边界的检测需求;3) Determine the area within the range of the radar's field of view (FOV, Field of view) within the initial screening range, and set a search range on both sides of the area with the central axis of the radar as the symmetry axis. The static target point within is defined as the alternative boundary starting point of the road, and in the case that the alternative boundary starting point is non-empty, step 4) is performed; wherein, each search range takes the sector radius boundary corresponding to the field of view angle range as One of the borders and the other border are parallel to the sector radius border, and the distance between the two borders should meet the detection requirements of borders such as guardrails;

4)求取所有获得的备选边界起始点距离雷达原点的距离,由于边界起始点的特点,因此在求取距离时本发明加上了一个影响因子,其计算公式如式(1)所示:4) Calculate the distances from all the obtained alternative boundary starting points to the radar origin. Due to the characteristics of the boundary starting points, the present invention adds an influence factor when calculating the distance, and its calculation formula is shown in formula (1). :

d2=αPx 2+βPy 2---(1)d 2 =αP x 2 +βP y 2 ---(1)

式中:d为计算得到的备选边界起始点距离雷达原点的距离,α为横向的比例因子,β为纵向的比例因子,Px为备选边界起始点的横坐标,Py为备选边界起始点的纵坐标,因此可以利用α与β的值来调节Px与Py对d的影响。将所有备选边界起始点的距离d的最小值作为边界起始点;In the formula: d is the calculated distance from the starting point of the alternative boundary to the origin of the radar, α is the horizontal scale factor, β is the vertical scale factor, P x is the abscissa of the starting point of the alternative boundary, and P y is the alternative boundary The ordinate of the starting point of the boundary, so the value of α and β can be used to adjust the influence of P x and P y on d. Take the minimum value of the distance d from all alternative boundary starting points as the boundary starting point;

5)利用获得边界起始点,在其横向与纵向某个范围内,在稳定静态目标点中查找下一个边界点,然后以找到的边界点作为新的起始点,进而再查找下一个边界点,以此类推,在查找不到下一个边界点时,进入步骤6);5) Use the obtained boundary starting point to find the next boundary point in the stable static target point within a certain horizontal and vertical range, and then use the found boundary point as a new starting point, and then search for the next boundary point, By analogy, when the next boundary point cannot be found, go to step 6);

6)判断步骤5查找得到的边界点的有效性,判断步骤5)查找到的边界点的数量与长度是否满足设定的阈值;若不满足,则去除当前的备选边界起始点,判断步骤(3)中获得的剩余可能的备选边界起始点是否为空,如果不为空,则选择下一个最近的备选边界起始点作为边界起始点,回到步骤5),如果为空,输出为空。若找到的边界点满足要求,则进行步骤7);6) Judging the validity of the boundary points found in step 5, determine whether the number and length of the boundary points found in step 5) meet the set threshold; if not, remove the current alternative boundary starting point, and determine the step Whether the remaining possible candidate boundary starting points obtained in (3) are empty, if not, select the next closest candidate boundary starting point as the boundary starting point, go back to step 5), if it is empty, output Is empty. If the found boundary points meet the requirements, then proceed to step 7);

7)利用最小二乘法对步骤6)得到的点进行二次曲线拟合,拟合式子如式(2)所示:7) Use the least squares method to perform quadratic curve fitting on the points obtained in step 6), and the fitting formula is shown in formula (2):

Px=aPy 2+bPy+c---(2)P x =aP y 2 +bP y +c---(2)

式中a为拟合的函数的二次项系数,b为拟合的函数的一次项系数,c为拟合的函数的常数项;Px为边界点的横坐标,Py为边界点的纵坐标。where a is the quadratic term coefficient of the fitted function, b is the linear term coefficient of the fitted function, and c is the constant term of the fitted function; P x is the abscissa of the boundary point, and P y is the Y-axis.

8)然后利用当前帧与上一帧得到的a,b,c三个系数利用卡尔曼滤波的思想对其进行滤波处理,得到时域上比较稳定的系数输出。8) Then use the three coefficients a, b, and c obtained from the current frame and the previous frame to filter them with the idea of Kalman filtering to obtain a relatively stable coefficient output in the time domain.

参见图5所示,以道路中某一侧边界为例子,讲解卡尔曼滤波的原理,设k时刻与k-1时刻道路边界拟合的二次曲线如式(3)所示:Referring to Figure 5, taking a certain side boundary of the road as an example to explain the principle of Kalman filtering, let the quadratic curve fitting the road boundary at time k and time k-1 is shown in formula (3):

Figure BDA0003482651200000111
Figure BDA0003482651200000111

式中ak,bk,ck与ak-1,bk-1,ck-1分别为k时刻与k-1时刻拟合得到多项式系数,Pk x,Pk -1 x分别为k时刻与k-1时刻边界点的横坐标,Pk y,Pk-1 y分别为k时刻与k-1时刻边界点的纵坐标。因此可以设状态向量x=[a b c]T,由于车载毫米波雷达可以测量得到道路边界点的坐标信息,因此可以设测量向量

Figure BDA0003482651200000112
为每一帧聚类得到的所有边界点的横坐标值,n为边界点的数量。得到两个状态向量与测量向量之后,其状态方程与观测方程分别如式(4)与式(5)所示:In the formula, a k , b k , c k and a k-1 , b k-1 , c k-1 are the polynomial coefficients obtained by fitting at time k and time k-1 respectively, P k x , P k -1 x respectively is the abscissa of the boundary point at time k and time k-1, and P k y and P k-1 y are the ordinate of the boundary point at time k and time k-1, respectively. Therefore, the state vector x=[abc] T can be set. Since the vehicle-mounted millimeter-wave radar can measure the coordinate information of the road boundary point, the measurement vector can be set.
Figure BDA0003482651200000112
The abscissa values of all boundary points obtained by clustering for each frame, n is the number of boundary points. After two state vectors and measurement vectors are obtained, the state equation and observation equation are shown in equations (4) and (5), respectively:

xk=Axk-1+wk-1---(4)x k =Ax k-1 +w k-1 ---(4)

zk=Hxk+vk---(5)z k = Hx k +v k ---(5)

式中,A为状态转移矩阵,Wk-1为预测噪声的协方差矩阵,H为观测矩阵,Vk为测量噪声的协方差矩阵,Xk与Xk-1分别为k时刻与k-1时刻的状态向量,Zk与Zk-1分别为k时刻与k-1时刻的测量向量。在汽车行驶的过程中,由于k时刻与k-1雷达坐标系会存在一定转角θ(逆时针为正,顺时针为负),需要进行坐标系旋转变换,得到k-1时刻在k时刻坐标系下的曲线方程。这里xk-1我们取经过坐标变换之后的状态向量,因此,状态转移矩阵A为单位矩阵。由观测向量可得到观测矩阵如式(6)所示:In the formula, A is the state transition matrix, W k-1 is the covariance matrix of the predicted noise, H is the observation matrix, V k is the covariance matrix of the measurement noise, X k and X k-1 are the k time and k- The state vector at time 1, Z k and Z k-1 are the measurement vectors at time k and time k-1, respectively. In the process of driving the car, since there will be a certain rotation angle θ between the k-1 radar coordinate system and the k-1 radar coordinate system (counterclockwise is positive, clockwise is negative), the coordinate system needs to be rotated and transformed to obtain the coordinates of k-1 at time k. The curve equation under the system. Here x k-1 we take the state vector after coordinate transformation, therefore, the state transition matrix A is the identity matrix. The observation matrix can be obtained from the observation vector as shown in formula (6):

Figure BDA0003482651200000113
Figure BDA0003482651200000113

式中,Py为每帧边界点聚类得到的所有边界点的纵向坐标,由此得到了卡尔曼滤波所必需的参数,其具体流程如图5所示。图中

Figure BDA0003482651200000121
为k时刻状态向量的先验估计,Pk|k-1为k时刻先验估计的协方差矩阵,Pk-1为k-1时刻后验估计误差的协方差矩阵,Q为状态估计噪声的协方差矩阵,R为观测噪声的协方差矩阵。卡尔曼滤波的详细步骤为:In the formula, P y is the longitudinal coordinates of all boundary points obtained by clustering the boundary points of each frame, thus obtaining the necessary parameters of the Kalman filter. The specific process is shown in Figure 5. pictured
Figure BDA0003482651200000121
is the prior estimation of the state vector at time k, P k|k-1 is the covariance matrix of the prior estimation at time k, P k-1 is the covariance matrix of the posterior estimation error at time k-1, and Q is the state estimation noise The covariance matrix of , R is the covariance matrix of the observation noise. The detailed steps of the Kalman filter are:

判断上一时刻(k-1)状态向量是否为空,若为空,则判断当前时刻(k)选取的边界点是否满足阈值,若不满足,直接结束,否则,进行二次函数拟合,并输出特征;Determine whether the state vector at the last moment (k-1) is empty, if it is empty, then determine whether the boundary point selected at the current moment (k) meets the threshold, if not, end directly, otherwise, perform quadratic function fitting, and output features;

如果(k-1)时刻的状态向量不为空,则判断(k)时刻选取的边界点是否满足阈值,若不满足,直接结束,否则,进行步骤4);If the state vector at time (k-1) is not empty, then judge whether the boundary point selected at time (k) satisfies the threshold, if not, end directly, otherwise, go to step 4);

进行卡尔曼滤波预测过程,如式(7)和(8)所示:The Kalman filter prediction process is performed, as shown in equations (7) and (8):

Figure BDA0003482651200000122
Figure BDA0003482651200000122

Pk|k-1=APk-1AT+Q---(8)P k|k-1 =AP k-1 A T +Q---(8)

利用当前时刻筛选得到的边界点,生成当前时刻的观测矩阵H和状态向量zkUse the boundary points screened at the current moment to generate the observation matrix H and the state vector z k at the current moment;

卡尔曼滤波的校正阶段,公式如式(9),(10)和(11):The correction stage of Kalman filter, the formula is as formula (9), (10) and (11):

Kk=Pk|k-1HT(HPk|k-1HT+R)-1---(9)K k =P k|k-1 H T (HP k|k-1 H T +R) -1 ---(9)

Figure BDA0003482651200000123
Figure BDA0003482651200000123

Pk=(I-KkH)Pk|k-1---(11)P k =(IK k H)P k|k-1 ---(11)

式中,Kk为卡尔曼增益,

Figure BDA0003482651200000124
为卡尔曼滤波得到的最优估计值,Pk为k时刻后验估计误差的协方差矩阵,I为1的矩阵。where K k is the Kalman gain,
Figure BDA0003482651200000124
is the optimal estimated value obtained by Kalman filtering, P k is the covariance matrix of the posterior estimation error at time k, and I is a matrix of 1.

最后,将卡尔曼滤波得到的最优估计值输出。Finally, the optimal estimated value obtained by Kalman filtering is output.

另一方面,参见图6所示,本发明实施例还提供了一种基于车载毫米波雷达的道路边界检测系统1,包括:On the other hand, referring to FIG. 6 , an embodiment of the present invention further provides a road boundary detection system 1 based on a vehicle-mounted millimeter-wave radar, including:

目标点获取模块10,被配置为获取车载毫米波雷达输出的当前帧所有的静态目标点;The target point acquisition module 10 is configured to acquire all static target points of the current frame output by the vehicle-mounted millimeter-wave radar;

初步筛选模块20,被配置为根据初始筛选范围从所有的静态坐标点中筛选出多个稳定静态目标点,初始筛选范围根据雷达探测能力确定;The preliminary screening module 20 is configured to screen out a plurality of stable static target points from all the static coordinate points according to the initial screening range, and the initial screening range is determined according to the radar detection capability;

备选起始点获取模块30,被配置为将位于预设筛选范围内的稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行;The alternative starting point acquisition module 30 is configured to use a stable static target point located within a preset screening range as an alternative boundary starting point; wherein, taking the radar center axis as the axis of symmetry, within the range of the field of view of the vehicle-mounted millimeter-wave radar Two search ranges are set inside, and the set two search ranges are used as preset screening ranges. Each search range takes the sector radius boundary corresponding to the field of view of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other boundary is the same as the The sector radius boundaries are parallel;

初始起始点确定模块40,被配置为对备选边界起始点按照与雷达原点之间的距离由小到大进行排序,将距离最小的备选边界起始点作为边界起始点;The initial starting point determination module 40 is configured to sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point;

查找模块50,被配置为以边界起始点为起始点,在预设的查找范围内从稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止;The search module 50 is configured to take the boundary starting point as the starting point, search for the boundary point from the stable static target points within the preset search range, and use the found boundary point as a new starting point to continue to search for the next boundary point, until the next boundary point cannot be found;

曲线拟合模块60,被配置为利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数;The curve fitting module 60 is configured to use the least squares method to perform quadratic curve fitting on all boundary points obtained by searching to obtain quadratic curve coefficients of the current frame;

以及,滤波处理模块70,被配置为利用卡尔曼滤波算法对当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。And, the filtering processing module 70 is configured to perform filtering processing on the quadratic curve coefficient of the current frame and the quadratic curve coefficient of the previous frame by using the Kalman filtering algorithm, so as to obtain the target quadratic curve coefficient and the target road boundary curve.

上述中基于车载毫米波雷达的道路边界检测系统各模块的具体细节已经在对应的基于车载毫米波雷达的道路边界检测方法中进行了详细的描述,因此此处不再赘述。The specific details of each module of the above-mentioned vehicle-mounted millimeter-wave radar-based road boundary detection system have been described in detail in the corresponding vehicle-mounted millimeter-wave radar-based road boundary detection method, so they will not be repeated here.

再一方面,本发明实施例还提供了一种电子设备,包括:处理器和存储器,存储器上存储有计算机可读指令,计算机可读指令被处理器执行时实现如上述实施例所述的基于车载毫米波雷达的道路边界检测方法。In another aspect, an embodiment of the present invention also provides an electronic device, including: a processor and a memory, where computer-readable instructions are stored in the memory, and when the computer-readable instructions are executed by the processor, the based Road boundary detection method for vehicle-mounted millimeter-wave radar.

具体地,上述存储器和处理器能够为通用的存储器和处理器,这里不做具体限定,当处理器运行存储器存储的计算机可读指令时,能够执行上述实施例所述的基于车载毫米波雷达的道路边界检测方法。Specifically, the above-mentioned memory and processor can be a general-purpose memory and processor, which is not specifically limited here. When the processor executes the computer-readable instructions stored in the memory, it can execute the vehicle-mounted millimeter-wave radar-based method described in the above embodiment. Road boundary detection method.

又一方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现如上述实施例所述的基于车载毫米波雷达的道路边界检测方法。In yet another aspect, an embodiment of the present invention also provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the road boundary detection based on the vehicle-mounted millimeter-wave radar as described in the above-mentioned embodiment is implemented method.

本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:闪存盘、只读存储器(read-only memory,ROM)、随机存取器(randomaccessmemory,RAM)、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing relevant hardware through a program, and the program can be stored in a computer-readable storage medium, and the storage medium can include: Flash disk, read-only memory (read-only memory, ROM), random access memory (random access memory, RAM), magnetic disk or optical disk, etc.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the scope of the present invention. within the scope of protection.

Claims (10)

1.一种基于车载毫米波雷达的道路边界检测方法,其特征在于,包括:1. a road boundary detection method based on vehicle-mounted millimeter wave radar, is characterized in that, comprises: 获取车载毫米波雷达输出的当前帧所有的静态目标点;Obtain all static target points of the current frame output by the vehicle millimeter-wave radar; 根据初始筛选范围从所有的所述静态坐标点中筛选出多个稳定静态目标点,所述初始筛选范围根据雷达探测能力确定;Screen a plurality of stable static target points from all the static coordinate points according to the initial screening range, and the initial screening range is determined according to the radar detection capability; 将位于预设筛选范围内的所述稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为所述预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行;The stable static target point located in the preset screening range is used as the starting point of the alternative boundary; wherein, taking the central axis of the radar as the symmetry axis, two search ranges are set within the field of view of the vehicle-mounted millimeter-wave radar, and the The two set search ranges are used as the preset screening ranges, and each search range takes the sector radius boundary corresponding to the field of view range of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other boundary is parallel to the sector radius boundary; 对所述备选边界起始点按照与雷达原点之间的距离由小到大进行排序,将距离最小的所述备选边界起始点作为边界起始点;Sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point; 以所述边界起始点为起始点,在预设的查找范围内从所述稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止;Take the boundary starting point as the starting point, search for the boundary point from the stable static target point within the preset search range, and use the found boundary point as the new starting point to continue searching for the next boundary point until the search fails. until the next boundary point; 利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数;Use the least squares method to perform quadratic curve fitting on all the boundary points obtained by the search to obtain the quadratic curve coefficient of the current frame; 利用卡尔曼滤波算法对所述当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。The quadratic curve coefficients of the current frame and the quadratic curve coefficients of the previous frame are filtered by using the Kalman filtering algorithm, so as to obtain the target quadratic curve coefficients and the target road boundary curve. 2.根据权利要求1所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,还包括:2. the road boundary detection method based on vehicle-mounted millimeter wave radar according to claim 1, is characterized in that, also comprises: 如果在预设的查找范围内查找到多个边界点,则求取每个边界点与雷达原点之间的距离;If multiple boundary points are found within the preset search range, the distance between each boundary point and the origin of the radar is calculated; 将所有边界点按照求得的距离由小到大的顺序进行排序,并将距离最小值对应的边界点作为新的起始点。Sort all the boundary points in ascending order of the obtained distance, and use the boundary point corresponding to the minimum distance as the new starting point. 3.根据权利要求2所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,边界点与雷达原点之间的距离计算公式表示为:3. the road boundary detection method based on vehicle-mounted millimeter wave radar according to claim 2, is characterized in that, the distance calculation formula between boundary point and radar origin is expressed as: d2=αPx 2+βPy 2---(1)d 2 =αP x 2 +βP y 2 ---(1) 式中:d为计算得到的边界点与雷达原点之间的距离,α为横向的比例因子,β为纵向的比例因子,Px为边界点的横坐标,Py为边界点的纵坐标。In the formula: d is the distance between the calculated boundary point and the origin of the radar, α is the horizontal scale factor, β is the vertical scale factor, P x is the abscissa of the boundary point, and P y is the ordinate of the boundary point. 4.根据权利要求1所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,还包括:4. the road boundary detection method based on vehicle-mounted millimeter wave radar according to claim 1, is characterized in that, also comprises: 如果查找得到的所有边界点满足设定的数量阈值和长度阈值,则转入利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数的步骤;If all the boundary points obtained by the search meet the set quantity threshold and length threshold, then transfer to the step of using the least squares method to perform quadratic curve fitting on all the boundary points obtained by the search to obtain the quadratic curve coefficient of the current frame; 否则,去除当前的所述边界起始点,并将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点。Otherwise, the current boundary starting point is removed, and the candidate boundary starting point located next to the current boundary starting point in the sorting is used as a new boundary starting point. 5.根据权利要求4所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,所述去除当前的所述边界起始点,并将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点的步骤包括:5. The road boundary detection method based on the vehicle-mounted millimeter-wave radar according to claim 4, wherein, the current boundary starting point is removed, and the current boundary starting point is located next to the current boundary starting point in the sorting. The step of using the alternative boundary start point of the bit as a new boundary start point includes: 统计剩余的所述备选边界起始点的数量;count the number of the remaining candidate boundary starting points; 如果所述剩余的所述备选边界起始点的数量大于零,则将排序中位于当前的所述边界起始点的下一位的所述备选边界起始点作为新的边界起始点。If the number of the remaining candidate boundary starting points is greater than zero, the candidate boundary starting point located next to the current boundary starting point in the sorting is used as a new boundary starting point. 6.根据权利要求1所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,在雷达坐标系下,通过设置横向阈值与纵向阈值得到所述初始筛选范围;其中,所述横向阈值根据真实场景下道路宽度、雷达探测能力得到;所述纵向阈值根据雷达探测能力确定。6 . The road boundary detection method based on vehicle-mounted millimeter-wave radar according to claim 1 , wherein, in the radar coordinate system, the initial screening range is obtained by setting a horizontal threshold and a vertical threshold; wherein, the horizontal threshold Obtained according to the road width and radar detection capability in the real scene; the longitudinal threshold is determined according to the radar detection capability. 7.根据权利要求1所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数,包括:7. the road boundary detection method based on vehicle-mounted millimeter-wave radar according to claim 1, is characterized in that, utilize least squares method to carry out quadratic curve fitting to obtain the quadratic curve coefficient of current frame to all boundary points that search obtains, include: 获取上一帧查找得到的所有边界点;Get all the boundary points found in the previous frame; 将当前帧查找得到的所有边界点和上一帧查找得到的所有边界点进行融合;Fusion of all boundary points found in the current frame and all boundary points found in the previous frame; 利用最小二乘法对融合得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数。Use the least squares method to perform quadratic curve fitting on all the fused boundary points to obtain the quadratic curve coefficients of the current frame. 8.根据权利要求1-7任一项所述的基于车载毫米波雷达的道路边界检测方法,其特征在于,起始点的纵坐标越大,对应的预设的查找范围的横向宽度和纵向宽度越大。8. The road boundary detection method based on the vehicle-mounted millimeter-wave radar according to any one of claims 1-7, characterized in that, the larger the ordinate of the starting point is, the greater the horizontal width and the vertical width of the corresponding preset search range bigger. 9.一种基于车载毫米波雷达的道路边界检测系统,其特征在于,包括:9. A road boundary detection system based on a vehicle-mounted millimeter-wave radar, characterized in that, comprising: 目标点获取模块,被配置为获取车载毫米波雷达输出的当前帧所有的静态目标点;The target point acquisition module is configured to acquire all static target points of the current frame output by the vehicle-mounted millimeter-wave radar; 初步筛选模块,被配置为根据初始筛选范围从所有的所述静态坐标点中筛选出多个稳定静态目标点,所述初始筛选范围根据雷达探测能力确定;a preliminary screening module, configured to screen out a plurality of stable static target points from all the static coordinate points according to an initial screening range, the initial screening range being determined according to the radar detection capability; 备选起始点获取模块,被配置为将位于预设筛选范围内的所述稳定静态目标点作为备选边界起始点;其中,以雷达中心轴线为对称轴,在车载毫米波雷达的视场角范围内设定两个查找范围,将设定的两个查找范围作为所述预设筛选范围,每个查找范围以车载毫米波雷达的视场角范围对应的扇形半径边界为其中一个边界,另一个边界与扇形半径边界平行;The alternative starting point acquisition module is configured to use the stable static target point located within the preset screening range as the alternative boundary starting point; wherein, taking the radar center axis as the axis of symmetry, at the field of view of the vehicle-mounted millimeter-wave radar Two search ranges are set within the range, and the two set search ranges are used as the preset screening ranges, and each search range takes the sector radius boundary corresponding to the field of view angle range of the vehicle-mounted millimeter-wave radar as one of the boundaries, and the other a boundary parallel to the sector radius boundary; 初始起始点确定模块,被配置为对所述备选边界起始点按照与雷达原点之间的距离由小到大进行排序,将距离最小的所述备选边界起始点作为边界起始点;The initial starting point determination module is configured to sort the candidate boundary starting points according to the distance from the radar origin from small to large, and use the candidate boundary starting point with the smallest distance as the boundary starting point; 查找模块,被配置为以所述边界起始点为起始点,在预设的查找范围内从所述稳定静态目标点中查找边界点,并将找到的边界点作为新的起始点继续查找下一个边界点,直至查找不到下一个边界点为止;The search module is configured to take the boundary starting point as a starting point, search for a boundary point from the stable static target point within a preset search range, and use the found boundary point as a new starting point to continue searching for the next Boundary point until the next boundary point cannot be found; 曲线拟合模块,被配置为利用最小二乘法对查找得到的所有边界点进行二次曲线拟合得到当前帧的二次曲线系数;a curve fitting module, configured to perform quadratic curve fitting on all boundary points obtained by the least squares method to obtain quadratic curve coefficients of the current frame; 以及,滤波处理模块,被配置为利用卡尔曼滤波算法对所述当前帧的二次曲线系数和上一帧的二次曲线系数进行滤波处理,以得到目标二次曲线系数和目标道路边界曲线。And, the filtering processing module is configured to perform filtering processing on the quadratic curve coefficients of the current frame and the quadratic curve coefficients of the previous frame by using the Kalman filtering algorithm, so as to obtain the target quadratic curve coefficients and the target road boundary curve. 10.一种电子设备,其特征在于,包括:处理器和存储器,所述存储器上存储有计算机可读指令,所述计算机可读指令被所述处理器执行时实现如权利要求1至8中任一项所述的基于车载毫米波雷达的道路边界检测方法。10. An electronic device, characterized in that it comprises: a processor and a memory, wherein computer-readable instructions are stored on the memory, and when the computer-readable instructions are executed by the processor, implementation as in claims 1 to 8 Any one of the road boundary detection methods based on vehicle-mounted millimeter-wave radar.
CN202210072589.0A 2022-01-21 2022-01-21 Road boundary detection method and system based on vehicle-mounted millimeter-wave radar Pending CN114779235A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210072589.0A CN114779235A (en) 2022-01-21 2022-01-21 Road boundary detection method and system based on vehicle-mounted millimeter-wave radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210072589.0A CN114779235A (en) 2022-01-21 2022-01-21 Road boundary detection method and system based on vehicle-mounted millimeter-wave radar

Publications (1)

Publication Number Publication Date
CN114779235A true CN114779235A (en) 2022-07-22

Family

ID=82422728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210072589.0A Pending CN114779235A (en) 2022-01-21 2022-01-21 Road boundary detection method and system based on vehicle-mounted millimeter-wave radar

Country Status (1)

Country Link
CN (1) CN114779235A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115840227A (en) * 2023-02-27 2023-03-24 福思(杭州)智能科技有限公司 Road edge detection method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842231A (en) * 2016-11-08 2017-06-13 长安大学 A kind of road edge identification and tracking
WO2019140950A1 (en) * 2018-01-16 2019-07-25 华为技术有限公司 Vehicle positioning method and apparatus
CN112597839A (en) * 2020-12-14 2021-04-02 上海宏景智驾信息科技有限公司 Road boundary detection method based on vehicle-mounted millimeter wave radar

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842231A (en) * 2016-11-08 2017-06-13 长安大学 A kind of road edge identification and tracking
WO2019140950A1 (en) * 2018-01-16 2019-07-25 华为技术有限公司 Vehicle positioning method and apparatus
CN112597839A (en) * 2020-12-14 2021-04-02 上海宏景智驾信息科技有限公司 Road boundary detection method based on vehicle-mounted millimeter wave radar

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
于春和;: "基于激光雷达的道路边界检测", 电子器件, no. 03, 15 June 2008 (2008-06-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115840227A (en) * 2023-02-27 2023-03-24 福思(杭州)智能科技有限公司 Road edge detection method and device

Similar Documents

Publication Publication Date Title
CN111712731B (en) Target detection method, target detection system and movable platform
CN109087510B (en) Traffic monitoring method and device
CN112441022B (en) Lane center line determining method and device
CN113189975B (en) Method for determining position data and/or movement data of a vehicle
CN111693972A (en) Vehicle position and speed estimation method based on binocular sequence images
CN112781599B (en) Method for determining the position of a vehicle
CN114882701B (en) Parking space detection method and device, electronic equipment and machine readable storage medium
CN111007880B (en) Extended target tracking method based on automobile radar
CN115423019A (en) A density-based fuzzy clustering method and device
JP6038422B1 (en) Vehicle determination device, vehicle determination method, and vehicle determination program
CN113537411B (en) Improved fuzzy clustering method based on millimeter wave radar
CN113496199B (en) Histogram-based L-shaped detection of target objects
CN107991677A (en) A kind of pedestrian detection method
CN114973195A (en) Vehicle tracking method, device and system based on multi-information fusion
CN107480592B (en) Multi-lane detection method and tracking method
CN110413942B (en) Lane line equation screening method and screening module thereof
CN114779235A (en) Road boundary detection method and system based on vehicle-mounted millimeter-wave radar
CN108693517B (en) Vehicle positioning method and device and radar
CN116626630B (en) Object classification method and device, electronic equipment and storage medium
KR20230032628A (en) Method and apparatus for processing sensor information, and recording medium for recording program performing the method
CN117490710A (en) Parameter adjusting method and device for target fusion and target fusion system
WO2021035395A1 (en) Speed determining method and device, and storage medium
JP2018018215A (en) Object feature point detector
CN111830469A (en) A classification, identification and filtering method for noise data of vehicle-mounted millimeter-wave radar
CN118962682B (en) Target course angle fitting processing method of 4D millimeter wave Lei Dadian cloud based on data fusion

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