CN113405560A - 车辆定位和路径规划统一建模方法 - Google Patents

车辆定位和路径规划统一建模方法 Download PDF

Info

Publication number
CN113405560A
CN113405560A CN202110588879.6A CN202110588879A CN113405560A CN 113405560 A CN113405560 A CN 113405560A CN 202110588879 A CN202110588879 A CN 202110588879A CN 113405560 A CN113405560 A CN 113405560A
Authority
CN
China
Prior art keywords
laser
coordinate system
gps
point
positioning
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.)
Granted
Application number
CN202110588879.6A
Other languages
English (en)
Other versions
CN113405560B (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.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
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 Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN202110588879.6A priority Critical patent/CN113405560B/zh
Publication of CN113405560A publication Critical patent/CN113405560A/zh
Application granted granted Critical
Publication of CN113405560B publication Critical patent/CN113405560B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/34Route searching; Route guidance
    • G01C21/3446Details of route searching algorithms, e.g. Dijkstra, A*, arc-flags, using precalculated routes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/28Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
    • 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
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/86Combinations of lidar systems with systems other than lidar, radar or sonar, e.g. with direction finders
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Abstract

本发明提供了一种车辆定位和路径规划统一建模方法,将定位与路径规划创新性的进行统一建模,以路径规划的轨迹作为一种约束,计算GPS定位结果与激光雷达定位结果到约束轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果,利用激光定位数据修正GPS定位结果,得到车辆的精确定位信息。与现有技术相比,本发明将定位与路径规划关联起来,统一建模,融合GPS定位结果与激光雷达定位结果,解决了GPS在信号差的区域定位结果不准和激光雷达定位存在累计误差的问题,可以为车辆循迹导航提供精确的定位精度。

Description

车辆定位和路径规划统一建模方法
技术领域
本发明属于定位技术领域,具体涉及一种车辆定位和路径规划统一建模方法。
背景技术
近年来,随着汽车智能化的程度越来越高,车辆的定位问题与路径规划问题成为各方研 究的热点。目前主流的车辆定位方法主要有基于GPS的、基于视觉的和基于激光的定位。当 前基于GPS的车辆定位技术日趋成熟,然而其定位误差在10米左右,定位精度低,并且在 某些GPS信号弱或者接受不到的地方,其定位精度更是难以得到保证,单纯的GPS已经不能 满足车辆的定位需求。随着激光雷达传感器的广泛应用,基于激光的车辆定位技术得到快速 发展,激光雷达的优点是测量精确,能够比较精准的提供角度和距离信息,可以达到<1°的 角度精度以及厘米级别的测距精度,扫描范围广,但存在由于累计误差而导致的定位偏差问 题。
发明内容
针对现有技术的不足,本发明提供一种车辆定位和路径规划统一建模方法,提高定位精 度。
本发明所采用的技术方案如下:
一种车辆定位和路径规划统一建模方法,包括以下步骤:
S1、处理激光点云数据,得到二维栅格地图;
S2、通过标定点云与GPS数据之间的映射关系,将激光坐标系中的车辆轨迹信息和二 维栅格地图映射到GPS坐标系下;
S3、基于A星算法在二维栅格地图中确定路径规划的轨迹结果;
S4、以路径规划的轨迹作为约束轨迹,计算GPS定位结果与激光雷达定位结果到约束 轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。
优选地,步骤S1包括:
S11、对激光点云数据进行预处理,剔除噪声点、离群点,同时过滤掉场景中路面点云 信息,只保留场景中除路面信息之外的固有静态障碍信息;
S12、将预处理后的三维激光点云数据进行XOY平面的投影处理,然后将整个XOY平面均匀地划分为一定尺寸的栅格,规定存在点云的栅格为占据状态,无点云存在的栅格为自 由状态,构建二维栅格地图。
优选地,步骤S11包括:
S111、采用直通滤波器保留距离车辆一定距离内的点云数据,即:若激光点云数据中某 一点的坐标(xi,yi,zi)满足:xi∈(X1,X2)∩yi∈(Y1,Y2)∩zi∈(Z1,Z2),则该点被保留;其 中,Xi、Yi和Zi分别为激光点云坐标系沿X、Y和Z轴方向设定的阈值;
S112、采用体素网格滤波器对激光点云数据进行滤波降采样处理,将点云数据划分为多 个体素网格,计算每一个体素网格的重心,以重心点代表体素网格中的所有点;
S113、将z坐标小于一定高度的点作为候选地面点集P',并对这些点按照高度z从小 到大排序,求出平均高度zave,再计算高度小于zave+zthreshold的点的法线,Zthreshold表示预设阈 值,将与z轴的夹角大于预设角度的点保留,剩下的所有点作为最终地面点集Pground,从而 剔除最终地面点集Pground
优选地,步骤S2包括:
S21、在激光点云扫描范围内选择一标志物G;
S22、采集标志物G的GPS坐标位置信息(Gx,Gy,Gz)T
S23、提取标志物G在激光坐标系下的激光坐标位置(x,y,z)T
S24、根据最小二乘问题求出旋转矩阵R和平移矩阵T。
优选地,步骤S23包括:
S231、剔除激光点云中的离群点与噪声点;
S232、利用RANSAC算法拟合得到标志物立杆G的点集,通过计算其密度公式得到立杆 的质心坐标位置(x,y,z)T,以此坐标位置作为立杆G在激光坐标系下的激光坐标位置(x,y,z)T
优选地,步骤S24包括:
S241、标志物的坐标之间存在以下关系:
Figure BDA0003088702420000021
式中,R是激光坐标系到GPS坐标系下的旋转矩阵,T是激光坐标系到GPS坐标系下的 平移矩阵;
S242、将Z方向的坐标做近似0的处理,上式转换如下:
Figure BDA0003088702420000031
式中,[x,y]T是标志物G在GPS坐标系下的坐标位置,[Gx,Gy]T是标志物G在激光坐标系下的坐标位置;
S243、近似处理后的式子可变为:
Figure BDA0003088702420000032
转换为方程组形式:
Figure BDA0003088702420000033
进而变为矩阵形式:
Figure BDA0003088702420000034
即可转换为最小二乘问题:Ax=b,求解得出方程的解x=(ATA)-1ATb即求解出旋转矩 阵R和平移矩阵T。
优选地,步骤S3包括:
S31、在二维栅格地图中设定路径的起点与终点,确定起点与终点在二维栅格地图中的 位置信息;
S32、根据A*算法的启发函数选择最优的方向展开搜索;
S33、计算欧几里得距离:
Figure BDA0003088702420000035
式中,(xn,yn)代表当前节点的坐标,(xd,yd)代表目标节点的坐标;
S34、结合A*算法使用OPEN列表和CLOSED列表来进行搜索,在栅格地图中生成由起点 到终点的路径规划轨迹结果。
优选地,步骤S4包括:
S41、将当前车辆位置采集点记做点P,实时获取车辆位置采集点P的GPS位置信息[x,y]T,并将其转换到UTM坐标系下:UP=(xP,yP)=trans(x,y);(xP,yP)为位置采集点P在UTM坐标系下的GPS位置信息,trans()为将GPS坐标系转换为UTM坐标系下的转换关系;
同时,根据旋转矩阵R和平移矩阵T以及GPS坐标系到UTM坐标系下的转换关系trans(),得出车辆位置采集点P从激光坐标系下的位置[Gx,Gy]T转换到UTM坐标系下的位置:UG=(xG,yG)=trans(R*[Gx,Gy]T+T);(xG,yG)为位置采集点P在UTM坐标系下的激光雷达位 置信息;
根据变换关系将路径规划轨的迹转换到UTM坐标系下;
S42、分别计算GPS定位结果与激光雷达定位结果到规划轨迹之间的横向距离;
S43、以路径规划的轨迹作为约束轨迹,GPS定位结果与激光雷达定位结果到约束轨迹 之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。
优选地,步骤S42包括:
S421、路径规划轨迹由直线段构成,由两点式直线方程计算出规划轨迹的直线方程: l:aix+biy+ci=0;
S422、根据车辆位置采集点当前时刻在UTM坐标系下的GPS位置信息UP=(xP,yP)及激光 雷达位置信息UG=(xG,yG),分别求解出横向距离:
Figure BDA0003088702420000041
优选地,步骤S43包括:
S431、以当前位置xk和N个历史位置共同组成的轨迹作为当前状态: Xk=[xk xk-1… xk-N],每个位置均包括横纵两个方向的信息:xk=[xk X xk Y]T;基于运 动学匀变速公式:Xk=2Xk-1-Xk-2,得到状态转移矩阵
Figure BDA0003088702420000042
预测状态 协方差矩阵
Figure BDA0003088702420000043
其中Q为GPS定位精度;
S432、融合当前车辆的GPS定位结果UP和激光雷达定位结果UG以及横向距离D1、D2,建 立观测值z:λ1UP2UG;当D1<D2时,取
Figure BDA0003088702420000044
当D2<D1时, 取
Figure BDA0003088702420000045
S433、基于运动学匀变速运动公式预测下一时刻的位置,融合车辆的GPS定位结果UP和 激光雷达定位结果,以λ1UP2UG作为观测值,以GPS定位结果与激光雷达定位结果到路径规划轨迹之间的横向距离作为卡尔曼滤波器中的约束,由kalman滤波公式得到融合后的定位结果[xp X xp Y],该融合的定位结果为UTM坐标系下的结果。
本发明的有益效果为:本发明将定位与路径规划创新性的进行统一建模,以路径规划的 轨迹作为一种约束,计算GPS定位结果与激光雷达定位结果到约束轨迹之间的横向距离作为 卡尔曼滤波器中的约束,融合GPS定位与激光定位结果,利用激光定位数据修正GPS定位结 果,得到车辆的精确定位信息。与现有技术相比,本发明将定位与路径规划相互关联,统一 建模,融合GPS定位结果与激光雷达定位结果,解决了GPS在信号差的区域定位结果不准和 激光雷达定位存在累计误差的问题,可以为车辆循迹导航提供精确的定位精度。
附图说明
图1为本发明一实施例的车载设备示意图。
图2为本发明一实施例的坐标标定示意图。
图3为本发明一实施例的车辆定位和路径规划统一建模方法流程图。
图中:1-车辆,2-电源,3-工控机,4-激光雷达点云采集设备,5-GPS接收器,6-显示器,7-差分GPS基站。
具体实施方式
下面将结合附图对本发明作进一步的说明:
本发明提供一种车辆定位和路径规划统一建模方法,首先通过激光雷达采集到的三维点 云数据进行相关滤波处理和俯视投影得到二维栅格地图;然后通过标定点云局部坐标与GPS 数据的映射关系将激光坐标系中的轨迹信息和栅格地图映射到GPS坐标系下;最后基于A星 算法在二维栅格地图中确定路径规划的轨迹结果,以路径规划的轨迹作为约束轨迹,计算 GPS定位结果与激光定位结果到约束轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合 GPS定位与激光定位结果,利用激光定位数据修正GPS定位结果,得到车辆在全局坐标系下 的精确定位信息。该方法利用卡尔曼滤波器融合GPS、激光雷达与约束轨迹三者之间的结果, 解决了GPS在信号差的区域定位结果不准和激光雷达定位存在累计误差的问题,可以为车辆 的循迹导航提供精确的定位精度。
本发明实施例的车辆定位和路径规划统一建模方法,如图3所示,包括地图构建、坐标 标定和车辆定位三个阶段,其中车辆定位包括GPS定位、激光雷达定位和融合定位三个部分。 为了实现本发明的方法,所需的硬件如图1所示,它包括:
数据采集模块,具体包括车辆1及车辆1上的激光雷达点云采集设备4,GPS信号接收 器5;其中激光雷达点云采集设备用于采集车辆当前时刻下的点云信息,点云中包含有标定 过程所用的标志物-立杆。GPS信号接收器用于接收来自差分GPS基站7所发出的信号,实 现当前车辆在GPS坐标系下的定位。
数据传输模块,包括将激光雷达采集的点云信息和GPS接收器接收到的报文信息传输到 工控机3。
数据处理模块,包括对点云信息与坐标信息在工控机3中进行信息同步处理,并提前对 采集的点云信息进行建图处理,构建栅格地图,规划路径轨迹,完成对定位结果的融合。
此外,电源2为上述设备供电,显示器6实时显示规划路径与当前定位结果。
本发明实施例中,车辆定位和路径规划统一建模方法,具体包括以下步骤:
S1、激光点云数据处理得到二维栅格地图。
S11、激光点云预处理,在预处理中,需要将噪声点、离群点剔除,同时还需要过滤掉 场景中路面点云信息,只保留场景中除路面信息之外的固有静态障碍信息,如建筑物、树木 等。
采用直通滤波器剔除掉距离车辆较远处的点云数据,假设原始激光点云中某一点的坐标, 若该坐标满足:xi∈(X1,X2)∩yi∈(Y1,Y2)∩zi∈(Z1,Z2),则该点被保留,否则将该点剔除掉。 其中,Xi、Yi和Zi分别为激光点云坐标系沿X、Y和Z轴方向设定的阈值。
采用体素网格滤波器对激光原始点云数据进行滤波降采样处理,通过对点云的获取得到 激光原始点云数据,对原始点云数据网格化处理,将点云数据划分为5cm×5cm的体素网格, 然后通过计算每一个体素网格的重心,以重心点来代表该体素网格中的所有点。
利用条件筛选法检测出地面点集Pground。将z小于一定高度的点作为候选地面点集P', 然后对这些点按照高度z从小到大排序,求出平均高度zave,再计算高度小于zave+zthreshold的 点的法线,Zthreshold表示预设阈值,将与z轴的夹角大于10°的点保留,剩下的所有点作为 最终地面点集Pground,剔除最终地面点集Pground
S12、将预处理后的三维激光点云数据进行XOY平面的投影处理。然后将整个XOY平面 均匀地划分为一定尺寸的栅格,栅格尺寸根据实际情况人为设定,规定存在点云的栅格为占 据状态“1”,无点云存在的栅格为自由状态“0”,构建二维栅格地图。
S2、通过标定点云与RTK/INS数据的映射关系将激光坐标系中的轨迹信息和栅格地图映 射到GPS坐标系下。
S21、为了标定点云与RTK/INS数据的映射关系,在激光点云扫描范围内选择一个标志 物,如图2所示,本实施例选用立杆。
S22、运用GPS接收机、组合惯导系统以及RTK基站采集该立杆的GPS坐标位置信息(Gx,Gy,Gz)T
S23、提取立杆G在激光坐标系下的激光坐标位置(x,y,z)T
首先,在获取得到的激光扫描数据中,通过对点云的获取得到激光原始点云数据,剔除 激光点云中的离群点与噪声点。最后,利用RANSAC(Random Sample Consensus)算法对标 志物立杆G稳健估计其位置。立杆G可以近似的看作一个细小的圆柱体,圆柱方程可以表示 为:
Figure BDA0003088702420000071
通过RANSAC算法拟 合得到的标志物立杆G的点集,通过计算其密度公式f(x,y),可以得到圆柱体的质心坐标位 置(x,y,z)T,以此坐标位置作为立杆G在激光坐标系下的激光坐标位置(x,y,z)T
S24、根据最小二乘问题求出旋转矩阵R和平移矩阵T。
通过对标志物立杆G获取其GPS坐标系下的坐标(Gx,Gy,Gz)T,那么在不同坐标系下, 立杆的坐标之间存在以下关系:
Figure BDA0003088702420000072
其中R是激光坐标系到GPS坐标系下的旋 转矩阵,T是激光坐标系到GPS坐标系下的平移矩阵。
考虑到车辆在正常行驶过程中车辆紧贴地面行驶,且立杆紧靠地面,可将其Z方向的坐 标做近似0的处理,因此上式可转换如下:
Figure BDA0003088702420000073
[x,y]T是立杆G在GPS坐标系下的坐标位置,[Gx,Gy]T是立杆G在激光坐标系下的坐标位置。在二维情况下,旋转矩阵由两个未知数构成,即
Figure BDA0003088702420000074
平移向量可表示为
Figure BDA0003088702420000075
则上式可变为:
Figure BDA0003088702420000076
转换为方程组形式:
Figure BDA0003088702420000077
可变为矩阵形式:
Figure BDA0003088702420000078
即可转换为最小二乘问题:Ax=b,求解得出方程的解 x=(ATA)-1ATb。即可求解出旋转矩阵R和平移矩阵T。
S3、基于A星算法在二维栅格地图中确定路径规划的轨迹结果。
S31、在栅格地图设定路径的起点与终点,确定起点与终点在栅格地图中的位置信息。
S32、根据启发函数来选择最优的方向展开搜索。A*算法的启发函数表达式为: f(n)=g(n)+h(n),其中,f(n)表示在当前节点n时,从初始节点到达目标节点的总代价函数;g(n)表示从初始节点到达当前节点n实际代价值的大小;h(n)表示从当前节点n到达目标节点的最小估计代价值的大小。
S33、计算欧几里得距离
Figure BDA0003088702420000081
其中(xn,yn)代表当前节点的坐标; (xd,yd)代表目标节点的坐标。
S34、结合A*算法使用OPEN列表和CLOSED列表来进行搜索,在栅格地图中生成由起点 到终点的路径规划轨迹结果。
S4、以路径规划的轨迹作为约束轨迹,计算GPS定位结果与激光定位结果到约束轨迹之 间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。
S41、坐标转换。为了便于计算横向距离,将GPS坐标系与激光雷达坐标系转换到UTM (Universal Transverse Mercator)坐标系,即通用横轴墨卡托投影系统。把当前车辆测试采集点位置记做点P,利用GPS采集设备实时获取到车辆测试点P的GPS位置信息[x,y]T,将其转换到UTM坐标系底下UP=(xP,yP)=trans(x,y),(xP,yP)为测试采集点P的UTM坐标位置 信息,trans()为将GPS坐标系转换为UTM坐标系下的转换关系。同时,根据旋转矩阵R和平移矩阵T,以及GPS坐标系到UTM坐标系下的转换关系trans(),可以得出车辆测试采集点P从激光坐标系下的位置[Gx,Gy]T转换到UTM坐标系下的位置 UG=(xG,yG)=trans(R*[Gx,Gy]T+T)。同时根据变换关系将路径规划轨迹转换到UTM坐标系下。
S42、通过GPS定位与激光雷达定位计算定位结果到规划轨迹之间的横向距离。路径规 划轨迹由直线段构成,由两点式直线方程可以计算出规划轨迹的直线方程l:aix+biy+ci=0, 根据上一步骤的车辆采集点当前时刻在UTM坐标系下的GPS位置信息UP=(xP,yP)和激光雷达 位置信息UG=(xG,yG)可以求解出横向距离
Figure BDA0003088702420000082
S43、以路径规划的轨迹作为约束轨迹,上一步骤所述的GPS定位结果与激光定位结果 到约束轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。卡 尔曼滤波器融合定位具体为:
以当前位置和N个历史位置,由N+1个位置组成的轨迹作为当前状态(本实例中,取N=2):Xk=[xk xk-1 … xk-N],每个位置包括横纵两个方向的信息:xk=[xk X xk Y]T, 基于运动学匀变速公式:Xk=2Xk-1-Xk-2,从而可以得到状态转移矩阵
Figure BDA0003088702420000091
预测状态协方差矩阵
Figure BDA0003088702420000092
其中Q为GPS定位精度。这里取Q=10。
将当前车辆在GPS定位下的结果UP和激光定位下的结果UG,以及横向距离D1、D2,把λ1UP2UG作为观测值z,当D1<D2时,取
Figure BDA0003088702420000093
当D2<D1时, 取
Figure BDA0003088702420000094
为了求解出观测矩阵H,根据关系有观测值 z=H[xk xk-1xk-2]T;其中H为观测矩阵,求解得
Figure BDA0003088702420000095
观测噪声协方差矩阵
Figure BDA0003088702420000096
其中δ1 2,δ2 2,...,δn 2为多次测量的各测量误差数据的方差。
根据Kalman滤波公式:
Figure BDA0003088702420000097
其中,Xk'为融合定位后当前车辆位置的结果,误差精度可达到厘米级,满足当前车辆 定位需求。
基于运动学匀变速运动公式得到的预测数据和上述步骤中融合车辆在GPS定位下的结果 和激光定位下的结果λ1UP2UG作为观测值,以GPS定位结果与激光定位结果到路径规划 轨迹之间的横向距离作为卡尔曼滤波器中的约束,由Kalman滤波公式得到融合后的定位结 果[xp X xp Y],该融合的结果为UTM坐标系下的结果,再根据坐标转换关系: (x,y)=trans-1(xp X,xp Y),得到在全局坐标系下精确地融合定位结果。由于解决了GPS在信 号差的区域定位结果不准和激光雷达定位存在累计误差的问题,本发明设备简单,定位结果 鲁棒性高,可以为车辆循迹导航提供精确的定位精度,适用于复杂环境下的循迹车辆的导航 定位,尤其适用于无人驾驶中智能车的定位。
本领域的技术人员容易理解,以上仅为本发明的较佳实施例而已,并不用以限制本发明, 凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保 护范围之内。

Claims (10)

1.一种车辆定位和路径规划统一建模方法,其特征在于,包括以下步骤:
S1、处理激光点云数据,得到二维栅格地图;
S2、通过标定点云与GPS数据之间的映射关系,将激光坐标系中的车辆轨迹信息和二维栅格地图映射到GPS坐标系下;
S3、基于A星算法在二维栅格地图中确定路径规划的轨迹结果;
S4、以路径规划的轨迹作为约束轨迹,计算GPS定位结果与激光雷达定位结果到约束轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。
2.根据权利要求1所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S1包括:
S11、对激光点云数据进行预处理,剔除噪声点、离群点,同时过滤掉场景中路面点云信息,只保留场景中除路面信息之外的固有静态障碍信息;
S12、将预处理后的三维激光点云数据进行XOY平面的投影处理,然后将整个XOY平面均匀地划分为一定尺寸的栅格,规定存在点云的栅格为占据状态,无点云存在的栅格为自由状态,构建二维栅格地图。
3.根据权利要求2所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S11包括:
S111、采用直通滤波器保留距离车辆一定距离内的点云数据,即:若激光点云数据中某一点的坐标(xi,yi,zi)满足:xi∈(X1,X2)∩yi∈(Y1,Y2)∩zi∈(Z1,Z2),则该点被保留;其中,Xi、Yi和Zi分别为激光点云坐标系沿X、Y和Z轴方向设定的阈值;
S112、采用体素网格滤波器对激光点云数据进行滤波降采样处理,将点云数据划分为多个体素网格,计算每一个体素网格的重心,以重心点代表体素网格中的所有点;
S113、将z坐标小于一定高度的点作为候选地面点集P',并对这些点按照高度z从小到大排序,求出平均高度zave,再计算高度小于zave+zthreshold的点的法线,Zthreshold表示预设阈值,将与z轴的夹角大于预设角度的点保留,剩下的所有点作为最终地面点集Pground,从而剔除最终地面点集Pground
4.根据权利要求1所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S2包括:
S21、在激光点云扫描范围内选择一标志物G;
S22、采集标志物G的GPS坐标位置信息(Gx,Gy,Gz)T
S23、提取标志物G在激光坐标系下的激光坐标位置(x,y,z)T
S24、根据最小二乘问题求出旋转矩阵R和平移矩阵T。
5.根据权利要求4所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S23包括:
S231、剔除激光点云中的离群点与噪声点;
S232、利用RANSAC算法拟合得到标志物立杆G的点集,通过计算其密度公式得到立杆的质心坐标位置(x,y,z)T,以此坐标位置作为立杆G在激光坐标系下的激光坐标位置(x,y,z)T
6.根据权利要求4所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S24包括:
S241、标志物的坐标之间存在以下关系:
Figure FDA0003088702410000021
式中,R是激光坐标系到GPS坐标系下的旋转矩阵,T是激光坐标系到GPS坐标系下的平移矩阵;
S242、将Z方向的坐标做近似0的处理,上式转换如下:
Figure FDA0003088702410000022
式中,[x,y]T是标志物G在GPS坐标系下的坐标位置,[Gx,Gy]T是标志物G在激光坐标系下的坐标位置;
S243、近似处理后的式子可变为:
Figure FDA0003088702410000023
转换为方程组形式:
Figure FDA0003088702410000024
进而变为矩阵形式:
Figure FDA0003088702410000025
即可转换为最小二乘问题:Ax=b,求解得出方程的解x=(ATA)-1ATb即求解出旋转矩阵R和平移矩阵T。
7.根据权利要求1所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S3包括:
S31、在二维栅格地图中设定路径的起点与终点,确定起点与终点在二维栅格地图中的位置信息;
S32、根据A*算法的启发函数选择最优的方向展开搜索;
S33、计算欧几里得距离:
Figure FDA0003088702410000031
式中,(xn,yn)代表当前节点的坐标,(xd,yd)代表目标节点的坐标;
S34、结合A*算法使用OPEN列表和CLOSED列表来进行搜索,在栅格地图中生成由起点到终点的路径规划轨迹结果。
8.根据权利要求1所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S4包括:
S41、将当前车辆位置采集点记做点P,实时获取车辆位置采集点P的GPS位置信息[x,y]T,并将其转换到UTM坐标系下:UP=(xP,yP)=trans(x,y);(xP,yP)为位置采集点P在UTM坐标系下的GPS位置信息,trans()为将GPS坐标系转换为UTM坐标系下的转换关系;
同时,根据旋转矩阵R和平移矩阵T以及GPS坐标系到UTM坐标系下的转换关系trans(),得出车辆位置采集点P从激光坐标系下的位置[Gx,Gy]T转换到UTM坐标系下的位置:UG=(xG,yG)=trans(R*[Gx,Gy]T+T);(xG,yG)为位置采集点P在UTM坐标系下的激光雷达位置信息;
根据变换关系将路径规划轨的迹转换到UTM坐标系下;
S42、分别计算GPS定位结果与激光雷达定位结果到规划轨迹之间的横向距离;
S43、以路径规划的轨迹作为约束轨迹,GPS定位结果与激光雷达定位结果到约束轨迹之间的横向距离作为卡尔曼滤波器中的约束,融合GPS定位与激光定位结果。
9.根据权利要求8所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S42包括:
S421、路径规划轨迹由直线段构成,由两点式直线方程计算出规划轨迹的直线方程:l:aix+biy+ci=0;
S422、根据车辆位置采集点当前时刻在UTM坐标系下的GPS位置信息UP=(xP,yP)及激光雷达位置信息UG=(xG,yG),分别求解出横向距离:
Figure FDA0003088702410000032
10.根据权利要求9所述的车辆定位和路径规划统一建模方法,其特征在于,步骤S43包括:
S431、以当前位置xk和N个历史位置共同组成的轨迹作为当前状态:Xk=[xk xk-1 …xk-N],每个位置均包括横纵两个方向的信息:xk=[xk X xk Y]T;基于运动学匀变速公式:Xk=2Xk-1-Xk-2,得到状态转移矩阵
Figure FDA0003088702410000041
预测状态协方差矩阵
Figure FDA0003088702410000042
其中Q为GPS定位精度;
S432、融合当前车辆的GPS定位结果UP和激光雷达定位结果UG以及横向距离D1、D2,建立观测值z:λ1UP2UG;当D1<D2时,取
Figure FDA0003088702410000043
当D2<D1时,取
Figure FDA0003088702410000044
S433、基于运动学匀变速运动公式预测下一时刻的位置,融合车辆的GPS定位结果UP和激光雷达定位结果,以λ1UP2UG作为观测值,以GPS定位结果与激光雷达定位结果到路径规划轨迹之间的横向距离作为卡尔曼滤波器中的约束,由kalman滤波公式得到融合后的定位结果[xp X xp Y],该融合的定位结果为UTM坐标系下的结果。
CN202110588879.6A 2021-05-28 2021-05-28 车辆定位和路径规划统一建模方法 Active CN113405560B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110588879.6A CN113405560B (zh) 2021-05-28 2021-05-28 车辆定位和路径规划统一建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110588879.6A CN113405560B (zh) 2021-05-28 2021-05-28 车辆定位和路径规划统一建模方法

Publications (2)

Publication Number Publication Date
CN113405560A true CN113405560A (zh) 2021-09-17
CN113405560B CN113405560B (zh) 2023-03-24

Family

ID=77674883

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110588879.6A Active CN113405560B (zh) 2021-05-28 2021-05-28 车辆定位和路径规划统一建模方法

Country Status (1)

Country Link
CN (1) CN113405560B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115877429A (zh) * 2023-02-07 2023-03-31 安徽蔚来智驾科技有限公司 自动驾驶车辆的定位方法、装置、存储介质及车辆
CN116907512A (zh) * 2023-09-14 2023-10-20 山东博昂信息科技有限公司 一种室外无人清扫机器的组合打点清扫方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090043504A1 (en) * 2007-05-31 2009-02-12 Amrit Bandyopadhyay System and method for locating, tracking, and/or monitoring the status of personnel and/or assets both indoors and outdoors
US20110208429A1 (en) * 2010-02-24 2011-08-25 Microsoft Corporation Route Computation Based on Route-Oriented Vehicle Trajectories
CN106772516A (zh) * 2016-12-01 2017-05-31 湖南大学 一种基于模糊理论的复合定位新方法
US20180058861A1 (en) * 2016-08-26 2018-03-01 Here Global B.V. Automatic localization geometry detection
CN108519615A (zh) * 2018-04-19 2018-09-11 河南科技学院 基于组合导航和特征点匹配的移动机器人自主导航方法
CN111337030A (zh) * 2020-05-19 2020-06-26 北京数字绿土科技有限公司 一种基于背负式的激光雷达扫描系统、导航定位方法
CN111457929A (zh) * 2019-12-31 2020-07-28 南京工大数控科技有限公司 一种基于地理信息系统的物流车辆自主路径规划与导航方法
US20200256686A1 (en) * 2019-02-13 2020-08-13 The Boeing Company Methods and apparatus for determining a vehicle path
CN111928860A (zh) * 2020-07-24 2020-11-13 上海交通大学烟台信息技术研究院 一种基于三维曲面定位能力的自主车辆主动定位方法
CN112082545A (zh) * 2020-07-29 2020-12-15 武汉威图传视科技有限公司 一种基于imu和激光雷达的地图生成方法、装置及系统
WO2021022615A1 (zh) * 2019-08-02 2021-02-11 深圳大学 机器人探索路径生成方法、计算机设备和存储介质
CN112762957A (zh) * 2020-12-29 2021-05-07 西北工业大学 一种基于多传感器融合的环境建模及路径规划方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090043504A1 (en) * 2007-05-31 2009-02-12 Amrit Bandyopadhyay System and method for locating, tracking, and/or monitoring the status of personnel and/or assets both indoors and outdoors
US20110208429A1 (en) * 2010-02-24 2011-08-25 Microsoft Corporation Route Computation Based on Route-Oriented Vehicle Trajectories
US20180058861A1 (en) * 2016-08-26 2018-03-01 Here Global B.V. Automatic localization geometry detection
CN106772516A (zh) * 2016-12-01 2017-05-31 湖南大学 一种基于模糊理论的复合定位新方法
CN108519615A (zh) * 2018-04-19 2018-09-11 河南科技学院 基于组合导航和特征点匹配的移动机器人自主导航方法
US20200256686A1 (en) * 2019-02-13 2020-08-13 The Boeing Company Methods and apparatus for determining a vehicle path
WO2021022615A1 (zh) * 2019-08-02 2021-02-11 深圳大学 机器人探索路径生成方法、计算机设备和存储介质
CN111457929A (zh) * 2019-12-31 2020-07-28 南京工大数控科技有限公司 一种基于地理信息系统的物流车辆自主路径规划与导航方法
CN111337030A (zh) * 2020-05-19 2020-06-26 北京数字绿土科技有限公司 一种基于背负式的激光雷达扫描系统、导航定位方法
CN111928860A (zh) * 2020-07-24 2020-11-13 上海交通大学烟台信息技术研究院 一种基于三维曲面定位能力的自主车辆主动定位方法
CN112082545A (zh) * 2020-07-29 2020-12-15 武汉威图传视科技有限公司 一种基于imu和激光雷达的地图生成方法、装置及系统
CN112762957A (zh) * 2020-12-29 2021-05-07 西北工业大学 一种基于多传感器融合的环境建模及路径规划方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENG XIAO-DONG 等: "New method for UAV online path planning", 《2013 IEEE INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING, COMMUNICATION AND COMPUTING (ICSPCC 2013)》 *
叶锦华 等: "多信息融合定位的自动导向车鲁棒导引控制", 《电机与控制学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115877429A (zh) * 2023-02-07 2023-03-31 安徽蔚来智驾科技有限公司 自动驾驶车辆的定位方法、装置、存储介质及车辆
CN116907512A (zh) * 2023-09-14 2023-10-20 山东博昂信息科技有限公司 一种室外无人清扫机器的组合打点清扫方法
CN116907512B (zh) * 2023-09-14 2023-12-22 山东博昂信息科技有限公司 一种室外无人清扫机器的组合打点清扫方法

Also Published As

Publication number Publication date
CN113405560B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
Khan et al. A comparative survey of lidar-slam and lidar based sensor technologies
CN108225302B (zh) 一种石化工厂巡检机器人定位系统和方法
CN110243358A (zh) 多源融合的无人车室内外定位方法及系统
CN109282808B (zh) 用于桥梁三维巡航检测的无人机与多传感器融合定位方法
US6728608B2 (en) System and method for the creation of a terrain density model
CN108519615A (zh) 基于组合导航和特征点匹配的移动机器人自主导航方法
US20100191461A1 (en) System and method of lane path estimation using sensor fusion
CN110702091B (zh) 一种沿地铁轨道移动机器人的高精度定位方法
EP3617749B1 (en) Method and arrangement for sourcing of location information, generating and updating maps representing the location
CN113405560B (zh) 车辆定位和路径规划统一建模方法
CN102831646A (zh) 一种基于扫描激光的大尺度三维地形建模方法
Pfaff et al. Towards mapping of cities
CN112444246B (zh) 高精度的数字孪生场景中的激光融合定位方法
CN108759834A (zh) 一种基于全局视觉的定位方法
CN114565674B (zh) 自动驾驶车辆城市结构化场景纯视觉定位方法及装置
CN110412596A (zh) 一种基于图像信息和激光点云的机器人定位方法
Bai et al. A sensor fusion framework using multiple particle filters for video-based navigation
CN113791074A (zh) 一种基于多传感器融合的无人机桥梁裂缝巡检系统及方法
CN113822944B (zh) 一种外参标定方法、装置、电子设备及存储介质
CN114077249B (zh) 一种作业方法、作业设备、装置、存储介质
CN110441760A (zh) 一种基于先验地形图的大范围海底地形图拓展构图方法
Klein et al. LiDAR and INS fusion in periods of GPS outages for mobile laser scanning mapping systems
CN103777196B (zh) 基于地理信息的地面目标距离单站测量方法及其测量系统
CN115166721B (zh) 路侧感知设备中雷达与gnss信息标定融合方法及装置
KR102506411B1 (ko) 차량의 위치 및 자세 추정 방법, 장치 및 이를 위한 기록매체

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant