CN113721254B - 一种基于道路指纹空间关联矩阵的车辆定位方法 - Google Patents
一种基于道路指纹空间关联矩阵的车辆定位方法 Download PDFInfo
- Publication number
- CN113721254B CN113721254B CN202110916802.7A CN202110916802A CN113721254B CN 113721254 B CN113721254 B CN 113721254B CN 202110916802 A CN202110916802 A CN 202110916802A CN 113721254 B CN113721254 B CN 113721254B
- Authority
- CN
- China
- Prior art keywords
- road
- point cloud
- fingerprint
- matrix
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 120
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000010276 construction Methods 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000013519 translation Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000010426 asphalt Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/86—Combinations of lidar systems with systems other than lidar, radar or sonar, e.g. with direction finders
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/28—Navigation; 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Automation & Control Theory (AREA)
- Traffic Control Systems (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开一种基于道路指纹空间关联矩阵的车辆定位方法,该方法包括道路指纹空间关联矩阵构建和定位两个阶段。在道路指纹空间关联矩阵构建阶段,首先对高精度地图数据进行道路指纹提取,然后计算每个指纹之间的空间距离构建空间关联矩阵,得到每个场景具有唯一性的空间特征。在定位阶段,检测激光雷达点云中的道路指纹,并计算当前道路指纹的空间关联矩阵;再计算当前空间关联矩阵与地图数据库中空间关联矩阵之间的匹配度,选取匹配度最大的位置节点;最后利用当前点云道路指纹的坐标信息与该节点地图中道路指纹所包含的位置信息进行匹配,从而计算当前车辆的位置。本发明利用道路指纹关联矩阵的唯一性进行表征,保证了车辆定位的精确度。
Description
技术领域
本发明属于车辆定位领域,具体涉及一种基于道路指纹空间关联矩阵的车辆定位方法。
背景技术
随着经济发展与社会进步,智能汽车已逐渐进入人们的视野并取得了实质性的发展。而其中,车辆的自定位是智能车最核心的问题之一。目前,GPS技术被广泛应用于车辆定位,但低成本GPS定位精度较低,而高精度的GPS因成本较高而无法普遍使用。
近年来,开始出现了基于相机和激光雷达的SLAM定位方法,该方法通过提取和匹配特征点来计算运动状态而实现定位,但该方法依然存在着误差累积问题,由于没有先验地图信息进行矫正,随着运动距离的增加,该方法的误差将越来越大。
发明内容
本发明的目的在于,提供一种基于道路指纹空间关联矩阵的车辆定位方法,消除SLAM定位的误差累积问题,实现精度更高的定位。
本发明所采用的技术方案如下:
一种基于道路指纹空间关联矩阵的车辆定位方法,包括以下步骤:
道路指纹关联矩阵构建:
(1)对高精度点云地图进行道路指纹信息提取,把道路指纹聚类成一个个类别,并计算出每个类别的聚类中心,用聚类中心代表指纹的空间位置;
(2)计算每个指纹中心之间的欧式距离,利用欧式距离构建道路指纹关联矩阵,同时保留每个指纹的位置信息;
车辆定位:
(3)利用GPS获取车辆的初始位置信息,基于初始位置信息获取高精度点云地图中距离当前位置最接近的多个节点并提取节点中所包含的道路指纹关联矩阵;
(4)对当前位置所获取的点云数据进行道路指纹信息提取,获取当前位置的道路指纹关联矩阵;将当前位置的道路指纹关联矩阵与节点的道路指纹关联矩阵进行匹配,获取与当前位置最接近的节点;
(5)将当前位置的道路指纹位置信息和与当前位置最接近的节点的地图点云进行位置匹配,实现车辆精确定位。
进一步的,高精度点云地图的构建包括以下步骤:
在特定路段上,利用激光雷达采集包含道路指纹信息的三维激光点云,同时利用高精度惯导获取车辆的位置数据;
将位置数据转化成与激光坐标系尺度一样的UTM坐标系,接着利用最小二乘法来优化两坐标系轨迹之间的残差,得到两者之间的映射关系,即旋转平移矩阵;最后利用得到的激光雷达和高精度惯导之间的映射关系,将激光点云数据投影到激光坐标系下,获得高精度点云地图。
进一步的,采用欧式聚类来提取指纹信息,以将离散、无组织且信息量较少的单一激光点云聚合成一个个具有一定空间拓扑结构且信息量丰富的三维语义目标;聚类完成后,会获得多个类别,即多个指纹;
将每个类别的点云数据投影至平面,然后计算该类别中所有点云数据的中心坐标,以计算得到的中心坐标作为该类别所对应的指纹的空间位置。
进一步的,利用欧式距离构建道路指纹关联矩阵包括以下步骤:
得到指纹的空间位置之后,计算各个指纹之间的欧式距离,计算公式如下:
式中,d表示欧式距离,(x1,y1)和(x2,y2)分别表示不同指纹的空间位置;
得到指纹之间的距离后,选择其中最大和最小的距离值dmax,dmin,然后两者相减得到间隔大小,进而将这大小分成N个子空间,每个子空间的大小为ΔI=(dmax-dmin)/N,最后可以形成[dmin+k*ΔI,dmin+(k+1)*ΔI]的区间,k取值0至N-1之间;
然后将各指纹之间的距离进行投影,计算落在每个区间内的个数,最后就形成了一个1xN维的道路指纹关联矩阵。
进一步的,选取路杆-路杆、路杆-交通标志牌、路杆-车道线、井盖或减速带这三种关联来构建关联矩阵,矩阵表示如下:
式中,H1表示路杆-路杆形成的关联矩阵,H2表示路杆-交通标志牌形成的关联矩阵,H3表示路杆-车道线、井盖或减速带形成的关联矩阵。
进一步的,对当前位置所获取的点云数据进行道路指纹信息提取时,只选择道路指纹信息,剔除其他冗余信息。
进一步的,将当前位置的道路指纹关联矩阵与节点的道路指纹关联矩阵进行匹配,获取与当前位置最接近的节点具体为:匹配采用的是余弦相似度,计算公式如下:
其中,表示当前位置的道路指纹关联矩阵的第i行,/>表示地图中节点的道路指纹关联矩阵的第i行,将每行计算的得分进行相加;
两个关联矩阵越相似,向量夹角越小,cosθ绝对值越大,选取得分最大的那个节点地图信息,实现节点匹配定位。
进一步的,道路指纹分为交通标志牌、车道线、减速带、井盖和路杆。
进一步的,交通标志牌的提取方法为:设置一个高度阈值和一个强度阈值,将符合高度阈值和强度阈值的点云提取出来,即可得到正对的交通标志牌;对于背对的交通标志牌,首先利用PCA算法遍历三维激光点云数据以获取整体点云集,其次利用该点云集构造协方差矩阵Pcov,最后对协方差矩阵Pcov进行奇异值分解得到特征值λ1,λ2,λ3,并根据分解结果分析点云维度特征:
当λ1>>λ2=λ3时,D=1,说明点集主要在1个维度上分散开来,则该点集具有线性维度特征;当λ1=λ2>>λ3时,D=2,说明点集主要在2个维度上分散开来,则该点集具有平面维度特征;当λ1=λ2=λ3时,D=3,说明点集主要在3个维度上分散开来,则该点集为离散结构;对于交通标志牌,其点云集具有平面维度特征,因此提取D=2的点集;
车道线、井盖与减速带的提取方法为:车道线是附着在地面点中的,故提取车道线前需要提取地面点云,即设定高度阈值,将包含地面点的点云保留下来;之后再利用RANSAC算法对地面点进行提取,提取出地面点后,设定一个阈值threshold,对包含地面点的点云进行如下操作;
式中,intensity表示包含地面点的点云发射强度值;
选取flag为1的点云,得到车道线;井盖采取同样的方法进行提取;
减速带的提取方法为:利用曲率进行,定义某帧点云中某一点的曲率为:
式中,|S|为点集S中的点云数量,S为与该点在同一激光线束上的点集,且S表示的点集均匀分布在该点的两侧;ri为该点云与激光坐标系原点间的欧式距离,rj为另一侧与激光坐标系原点间的欧式距离;设定一个阈值,选择大于阈值的点云,提取出减速带;
路杆的提取方法为:
路杆具有垂直特性,故将点云投影至平面进行处理;将点云投影到平面后,在二维空间对点进行聚类,聚类的目的是把相似点划分为一类,聚类后路杆包含的点被划分在一起,鉴于路杆竖直方向的半径变化均匀且高度较高,对聚类后的点进行横向切片处理,计算各切片的对角线长度,如果两两相邻切片的对角线长增长在预设范围内,则认为该聚类为路杆。
进一步的,将当前位置的道路指纹位置信息和与当前位置最接近的节点的地图点云进行位置匹配,实现车辆精确定位具体为:
获取与当前位置最接近的节点的地图点云集Q,以及当前位置经过道路指纹提取的目标点云集P,按照一定的约束条件,找到最邻近点对(pi,qi),采用ICP算法计算出最优匹配参数R和T,使得误差函数E最小,计算公式如下:
式中,Pi、Qi分别代表地图点云集与目标点云集中的对应点云点对,n表示对应的点云的总对数;
最终得到的R、T则为最优旋转平移变换矩阵,也表示当前车辆的位姿。
本发明的有益效果为:
本发明首先在离线状态下利用高精度惯导和激光雷达来构建地图数据库,地图之中包含了位置信息和道路指纹关联矩阵。在定位时,首先利用普通GPS传感器采集得到的位置信息进行初定位,得到相近的多个节点地图信息,然后利用当前点云数据得到的道路指纹关联矩阵与节点地图中的矩阵进行匹配,选择相似度得分最大的一个节点,完成节点级定位,最终利用最近节点中的位置信息与当前点云数据的位置信息进行匹配,最终获得车辆高精度的定位。本发明提出的定位方法精度高、不易受到光线干扰、鲁棒性强。
附图说明
图1是本发明一实施例的车载设备示意图。
图2是本发明一实施例的坐标标定示意图。
图3是本发明一实施例的方法流程示意图。
图中:1-激光雷达,2-惯导。
具体实施方式
下面将结合附图对本发明作进一步的说明:
本发明公开了一种基于道路指纹空间关联矩阵的车辆定位方法,该方法包括道路指纹空间关联矩阵构建和定位两个阶段。在道路指纹空间关联矩阵构建阶段,首先对高精度地图数据进行道路指纹提取,道路指纹包含车道线、交通标志牌、减速带、井盖、路杆等信息,提取出道路指纹后,再计算每个指纹之间的空间距离构建空间关联矩阵,得到每个场景具有唯一性的空间特征。在定位阶段,通过激光雷达采集激光雷达点云,检测激光雷达点云中的道路指纹,并计算当前道路指纹的空间关联矩阵。之后再计算当前空间关联矩阵与地图数据库中空间关联矩阵之间的匹配度,选取匹配度最大的位置节点。最后利用当前点云道路指纹的坐标信息与该节点地图中道路指纹所包含的位置信息进行匹配,从而计算当前车辆的位置。本发明利用道路指纹关联矩阵的唯一性进行表征,保证了车辆定位的精确度。
图3是本发明实施例提供的一种方法流程示意图,包括地图数据库构建和车辆定位两个阶段;其中,地图数据库构建阶段主要包括激光雷达和高精度惯导数据,利用惯导数据和提取的道路指纹来得到位置和空间关联矩阵,定位阶段主要包括普通GPS的初定位、矩阵匹配的节点级定位和位置匹配的车辆定位。首先利用GPS做一个初定位,选择当前车辆所在的一个范围,获得多个节点的关联矩阵,例如当前位置前后40m范围的地图节点,这样可以减少全局搜索的范围,加快计算速率;之后对单帧点云进行处理,通过计算曲率和聚类等操作把道路指纹信息提取出来并获得当前的关联矩阵,利用这个关联矩阵与初定位选取的多个道路指纹关联矩阵库进行匹配,选择相似度最高的节点。最后利用当前位置的指纹位置信息与地图位置信息进行匹配实现车辆定位。
具体包括以下步骤:
S1:数据采集;
首先在本发明实施例路段上,如图1所示,利用数据采集车上装载的激光雷达1实时扫描得到当前环境包含道路指纹的三维激光点云,同时,数据采集车上还装载有高精度惯导定位模块,利用该惯导2得到当前的车辆位置信息;其中定位模块精度为厘米级。
S2:数据处理;
利用高精度惯导和激光雷达采集一段轨迹的GPS和激光雷达点云信息,然后将GPS数据转化成与激光坐标系尺度一样的UTM坐标系,如图2所示,坐标系统一后,需要对两轨迹进行拟合,首先构建损失函数,再利用最小二乘法来优化两轨迹之间的残差,得到两者之间的映射关系,即旋转平移矩阵。之后利用得到的两传感器的映射关系,将每帧点云利用对应的旋转平移矩阵将局部坐标系下的点云投影到全局坐标,得到包含道路指纹的高精度点云地图。
S3:道路指纹库构建;
S31:得到高精度地图后,对高精度点云地图进行道路指纹信息的提取,把道路指纹聚类成一个个类别,并计算出每个类别的聚类中心,用聚类中心代表指纹的空间位置。可以采用欧式聚类来提取指纹信息,聚类的目的在于将离散、无组织且信息量较少的单一激光点云聚合成一个个具有一定空间拓扑结构的且信息量丰富的三维语义目标,即道路指纹。本实施例采用的DBSCAN算法,它不需要指定簇的个数,具有擅长找到离群点,聚类结果没有偏倚等优点。
聚类完成后,使用聚类中心来表示每个类别的位置。为了减小点云扫描不均匀的影响,可以将点云投影至平面中,然后计算投影后点云的中心坐标,其计算公式可以表示成如下形式:
其中,xn为该类别的某点数据,为计算的该类别中心坐标。
S32:计算每个指纹中心之间的欧式距离,利用这个欧式距离构建关联矩阵,同时保留每个指纹的位置信息,用于位置匹配。每个类别的欧式距离可以表示成如下形式:
得到距离后,就可以构建关联矩阵。空间关联矩阵是一个二维矩阵。本实施例采用路杆-路杆、路杆-交通标志牌、路杆-车道线、井盖和减速带(路杆-车道线、路杆-井盖、路杆-减速带视为一类)三种关联来构建矩阵。把这个矩阵保存在文本中,在定位时从文本读取用于匹配。关联矩阵可以表示成如下:
式中,H1表示路杆-路杆形成的关联矩阵,H2表示路杆-交通标志牌形成的关联矩阵,H3表示路杆-车道线、井盖和减速带形成的关联矩阵。除了保存关联矩阵,我们还保留了每个类别的位置坐标。关联矩阵是为了寻找到相似节点,位置坐标是为了匹配实现车辆定位。
S4:初定位
利用待定位车辆上的普通GPS接收机来提供一个初始位置信息,利用这个信息来选择地图中距离当前位置最接近的多个节点并提取节点中所包含的点云和空间关联矩阵。在后续中,不再使用GPS数据,GPS只用于初始化阶段。
S5:道路指纹提取
道路指纹提取分为交通标志牌、车道线、减速带、井盖和路杆等标志。每个标志具有不同的提取方法。对当前激光雷达采集的点云进行特征提取时,可以只选择道路指纹信息,剔除不需要的冗余的其他信息,这样不仅减少了数据量,也可以加快匹配时间。
S51:交通标志牌提取
交通标志牌是道路中常见的道路标志,其布置一般具有高度要求,因此可以先设置一定的高度,将一些地面点滤除,也把包含车牌等高反射强度的物体滤除。同时考虑到交通标志的高反射率,激光雷达扫描到交通标志牌后其反射回来的激光束强度很大,因此利用这种特性将正对的交通标志牌提取出来。而对背对的情况,利用PCA算计进行提取。
PCA算法首先遍历三维激光点云数据以获取整体点云集,其次利用该点云集构造协方差矩阵Pcov,最后对协方差矩阵Pcov进行奇异值分解得到特征值λ1,λ2,λ3,并根据分解结果分析点云维度特征。
当λ1>>λ2=λ3时,D=1,说明点集主要在1个维度上分散开来,则该点集具有线性维度特征;当λ1=λ2>>λ3时,D=2,说明点集主要在2个维度上分散开来,则该点集具有平面维度特征;当λ1=λ2=λ3时,D=3,说明点集主要在3个维度上分散开来,则该点集为离散结构。对于交通标志牌,其点云集具有平面维度特征,因此提取D=2的点集。
S52:车道线、井盖与减速带提取
车道线是附着在地面点中的,因此提取车道线前需要提取地面点云。在知道激光雷达放置高度的前提下,先设定高度阈值,将包含地面点的点云保留下来,之后再利用RANSAC算法对地面点进行提取,提取出地面点后,考虑到车道线、井盖等材质组成,其反射能力较沥青、水泥等材质有明显区别,其反射能力较强,可以利用强度信息提取车道线、井盖等点云。本实施例设定一个阈值threshold,对其进行如下操作;
式中,intensity表示包含地面点的点云的反射强度值。选取flag为1的点云,经过上述操作,车道线就被检测出来了。井盖采取同样的方法进行提取。
减速带的提取,利用曲率进行,定义某帧点云中某一点的曲率为:
式中,|S|为点集S中的点云数量,S为与该点在同一激光线束上的点集,且S表示的点集均匀分布在该点的两侧;ri为该点云与激光坐标系原点间的欧式距离,rj为另一侧与激光坐标系原点间的欧式距离。设定一个阈值,选择大于阈值的点云,这样减速带就被提取出来。
S53:路杆提取
对于路杆,考虑到其具有垂直特性,要将其投影至平面进行处理。将点云投影到平面后,在二维空间对点进行聚类,聚类的目的是把相似点划分为一类,聚类后路杆包含的点被划分在一起,鉴于路杆竖直方向的半径变化均匀且高度较高,对其进行横向切片处理,计算各切片的对角线长度,如果两两相邻切片的对角线长增长不大,则认为该聚类为路杆。
S6:关联矩阵匹配
在提取出指纹后,需要对其进行关联矩阵计算,获得待定位位置的关联矩阵。在得到当前位置的关联矩阵后,将该矩阵与地图中矩阵进行匹配,选择矩阵相似度最大的节点,并保存下该节点地图中的点云信息,实现节点定位。
S61:对于提取出来的道路指纹,首先计算每个指纹的质心,用质心来代替每个指纹的位置信息。然后计算每个指纹之间的欧式距离用于关联矩阵计算。
S62:得到指纹之间的距离后,先进行排序选择其中最大和最小的距离值dmax,dmin,然后两者相减得到间隔大小,然后将这大小分成N子空间,每个子空间为ΔI=(dmax-dmin)/N,我们最后可以形成[dmin+k*ΔI,dmin+(k+1)*ΔI](k取值0-N-1之间)的区间,然后将各指纹之间的距离进行投影,计算落在每个区间内各个数,最后就形成了一个1xN维的道路指纹关联矩阵。对路杆-路杆、路杆-交通标志牌、路杆-车道线、井盖和减速带三种类别矩阵进行上述投影操作,最终可以形成3xN的矩阵。
S6.2:得到当前关联矩阵后,用这个矩阵和筛选得到的地图节点矩阵进行匹配,匹配采用的是余弦相似度。
其中表示当前矩阵的第i行,其中/>表示地图中关联矩阵的第i行,最终得分是将三行计算的得分进行相加:
两个关联矩阵越相似,向量夹角越小,cosθ绝对值越大,选取得分score最大的那个节点地图信息,实现最后的节点级定位
S7:车辆定位
每帧激光点云获取的道路指纹与地图匹配采用的是ICP算法,ICP算法不必对点云进行分割和特征提取,其精度和收敛性都不错。假设源点云Pi为激光实时采集的经过道路指纹提取的点云,目标点云Qi为已经构建好的点云,则ICP算法公式为
通过迭代来搜索旋转矩阵R和平移向量T,当E小于给定阈值时,则认为ICP已配准完成,停止迭代,否则将通过计算得到的R、T重新获得新的变换点云集继续与目标点云集配准,直至满足迭代结束条件。当迭代结束,得到的R、T则为最优旋转平移变换矩阵,也表示当前车辆的位姿。
本领域的技术人员容易理解,以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,包括以下步骤:
道路指纹关联矩阵构建:
(1)对高精度点云地图进行道路指纹信息提取,以将离散、无组织且信息量较少的单一激光点云聚合成一个个具有一定空间拓扑结构且信息量丰富的三维语义目标,即道路指纹;把道路指纹聚类成一个个类别后,计算出每个类别的聚类中心,用聚类中心代表指纹的空间位置;
(2)计算每个指纹中心之间的欧式距离,利用欧式距离构建道路指纹关联矩阵,同时保留每个指纹的位置信息,结合道路指纹关联矩阵以及指纹的位置信息构建地图节点数据库;
车辆定位:
(3)利用GPS获取车辆的初始位置信息,基于初始位置信息获取高精度点云地图中距离当前位置最接近的多个节点并提取节点中所包含的道路指纹关联矩阵;
(4)对当前位置所获取的点云数据进行道路指纹信息提取,获取当前位置的道路指纹关联矩阵;将当前位置的道路指纹关联矩阵与节点的道路指纹关联矩阵进行匹配,获取与当前位置最接近的节点;
(5)将当前位置的道路指纹位置信息和与当前位置最接近的节点的地图点云进行位置匹配,实现车辆精确定位。
2.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,高精度点云地图的构建包括以下步骤:
在特定路段上,利用激光雷达采集包含道路指纹信息的三维激光点云,同时利用高精度惯导获取车辆的位置数据;
将位置数据转化成与激光坐标系尺度一样的UTM坐标系,接着利用最小二乘法来优化两坐标系轨迹之间的残差,得到两者之间的映射关系,即旋转平移矩阵;最后利用得到的激光雷达和高精度惯导之间的映射关系,将每帧点云利用对应的旋转平移矩阵将局部坐标系下的点云投影到全局坐标系下,获得包含道路指纹的高精度点云地图。
3.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,采用欧式聚类来提取指纹信息,以将离散、无组织且信息量较少的单一激光点云聚合成一个个具有一定空间拓扑结构且信息量丰富的三维语义目标;聚类完成后,会获得多个类别,即多个指纹;
将每个类别的点云数据投影至平面,然后计算该类别中所有点云数据的中心坐标,以计算得到的中心坐标作为该类别所对应的指纹的空间位置。
4.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,利用欧式距离构建道路指纹关联矩阵包括以下步骤:
得到指纹的空间位置之后,计算各个指纹之间的欧式距离,计算公式如下:
式中,d表示欧式距离,(x1,y1)和(x2,y2)分别表示不同指纹的空间位置;
得到指纹之间的距离后,选择其中最大和最小的距离值dmax,dmin,然后两者相减得到间隔大小,进而将这大小分成N个子空间,每个子空间的大小为ΔI=(dmax-dmin)/N,最后可以形成[dmin+k*ΔI,dmin+(k+1)*ΔI]的区间,k取值0至N-1之间;
然后将各指纹之间的距离进行投影,计算落在每个区间内的个数,最后就形成了一个1xN维的道路指纹关联矩阵。
5.根据权利要求4所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,选取路杆-路杆、路杆-交通标志牌、路杆-车道线、井盖或减速带这三种关联来构建关联矩阵,矩阵表示如下:
式中,H1表示路杆-路杆形成的关联矩阵,H2表示路杆-交通标志牌形成的关联矩阵,H3表示路杆-车道线、井盖或减速带形成的关联矩阵。
6.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,对当前位置所获取的点云数据进行道路指纹信息提取时,只选择道路指纹信息,剔除其他冗余信息。
7.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,将当前位置的道路指纹关联矩阵与节点的道路指纹关联矩阵进行匹配,获取与当前位置最接近的节点具体为:匹配采用的是余弦相似度,计算公式如下:
其中,表示当前位置的道路指纹关联矩阵的第i行,/>表示地图中节点的道路指纹关联矩阵的第i行,将每行计算的得分进行相加;
两个关联矩阵越相似,向量夹角越小,cosθ绝对值越大,选取得分最大的那个节点地图信息,实现节点匹配定位。
8.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,道路指纹分为交通标志牌、车道线、减速带、井盖和路杆。
9.根据权利要求8所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,交通标志牌的提取方法为:设置一个高度阈值和一个强度阈值,将符合高度阈值和强度阈值的点云提取出来,即可得到正对的交通标志牌;对于背对的交通标志牌,首先利用PCA算法遍历三维激光点云数据以获取整体点云集,其次利用该点云集构造协方差矩阵Pcov,最后对协方差矩阵Pcov进行奇异值分解得到特征值λ1,λ2,λ3,并根据分解结果分析点云维度特征:
当λ1>>λ2=λ3时,D=1,说明点集主要在1个维度上分散开来,则该点集具有线性维度特征;当λ1=λ2>>λ3时,D=2,说明点集主要在2个维度上分散开来,则该点集具有平面维度特征;当λ1=λ2=λ3时,D=3,说明点集主要在3个维度上分散开来,则该点集为离散结构;对于交通标志牌,其点云集具有平面维度特征,因此提取D=2的点集;
车道线与井盖的提取方法为:车道线是附着在地面点中的,故提取车道线前需要提取地面点云,即设定高度阈值,将包含地面点的点云保留下来;之后再利用RANSAC算法对地面点进行提取,提取出地面点后,设定一个阈值threshold,对包含地面点的点云进行如下操作;
式中,intensity表示包含地面点的点云发射强度值;
选取flag为1的点云,得到车道线;井盖采取同样的方法进行提取;
减速带的提取方法为:利用曲率进行,定义某帧点云中某一点的曲率为:
式中,|S|为点集S中的点云数量,S为与该点在同一激光线束上的点集,且S表示的点集均匀分布在该点的两侧;ri为该点云与激光坐标系原点间的欧式距离,rj为另一侧与激光坐标系原点间的欧式距离;设定一个阈值,选择大于阈值的点云,提取出减速带;
路杆的提取方法为:
路杆具有垂直特性,故将点云投影至平面进行处理;将点云投影到平面后,在二维空间对点进行聚类,聚类的目的是把相似点划分为一类,聚类后路杆包含的点被划分在一起,鉴于路杆竖直方向的半径变化均匀且高度较高,对聚类后的点进行横向切片处理,计算各切片的对角线长度,如果两两相邻切片的对角线长增长在预设范围内,则认为该聚类为路杆。
10.根据权利要求1所述的基于道路指纹空间关联矩阵的车辆定位方法,其特征在于,将当前位置的道路指纹位置信息和与当前位置最接近的节点的地图点云进行位置匹配,实现车辆精确定位具体为:
获取与当前位置最接近的节点的地图点云集Q,以及当前位置经过道路指纹提取的目标点云集P,按照一定的约束条件,找到最邻近点对(pi,qi),采用ICP算法计算出最优匹配参数R和T,使得误差函数E最小,计算公式如下:
式中,Pi、Qi分别代表地图点云集与目标点云集中的对应点云点对,n表示对应的点云的总对数;
最终得到的R、T则为最优旋转平移变换矩阵,也表示当前车辆的位姿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110916802.7A CN113721254B (zh) | 2021-08-11 | 2021-08-11 | 一种基于道路指纹空间关联矩阵的车辆定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110916802.7A CN113721254B (zh) | 2021-08-11 | 2021-08-11 | 一种基于道路指纹空间关联矩阵的车辆定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113721254A CN113721254A (zh) | 2021-11-30 |
CN113721254B true CN113721254B (zh) | 2023-10-17 |
Family
ID=78675419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110916802.7A Active CN113721254B (zh) | 2021-08-11 | 2021-08-11 | 一种基于道路指纹空间关联矩阵的车辆定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113721254B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11906647B2 (en) | 2022-04-29 | 2024-02-20 | Koko Home, Inc. | Person location determination using multipath |
CN115620264B (zh) * | 2022-12-02 | 2023-03-07 | 广汽埃安新能源汽车股份有限公司 | 车辆定位方法、装置、电子设备和计算机可读介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20170063002A (ko) * | 2015-11-30 | 2017-06-08 | 현대엠엔소프트 주식회사 | Las 데이터를 이용한 도로 지문 데이터 구축 시스템 및 그 방법 |
CN109099929A (zh) * | 2018-07-13 | 2018-12-28 | 武汉理工大学 | 基于场景指纹的智能车定位装置及方法 |
WO2019242174A1 (zh) * | 2018-06-21 | 2019-12-26 | 华南理工大学 | 基于激光雷达的建筑结构自动测量及3d模型生成方法 |
WO2020020701A1 (en) * | 2018-07-26 | 2020-01-30 | Signify Holding B.V. | Method for configuring a tracking system, tracking system, lighting system incorporating a tracking system and computer program |
CN111239715A (zh) * | 2020-01-13 | 2020-06-05 | 哈尔滨工业大学 | 一种联合灰色关联和神经网络的指纹定位方法 |
CN111882612A (zh) * | 2020-07-21 | 2020-11-03 | 武汉理工大学 | 一种基于三维激光检测车道线的车辆多尺度定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9369986B2 (en) * | 2012-10-24 | 2016-06-14 | Telefonaktiebolaget Lm Ericsson (Publ) | Wireless communication network for estimating the accuracy of fingerprinting positioning algorithms |
US11195065B2 (en) * | 2019-11-22 | 2021-12-07 | Samsung Electronics Co., Ltd. | System and method for joint image and lidar annotation and calibration |
-
2021
- 2021-08-11 CN CN202110916802.7A patent/CN113721254B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20170063002A (ko) * | 2015-11-30 | 2017-06-08 | 현대엠엔소프트 주식회사 | Las 데이터를 이용한 도로 지문 데이터 구축 시스템 및 그 방법 |
WO2019242174A1 (zh) * | 2018-06-21 | 2019-12-26 | 华南理工大学 | 基于激光雷达的建筑结构自动测量及3d模型生成方法 |
CN109099929A (zh) * | 2018-07-13 | 2018-12-28 | 武汉理工大学 | 基于场景指纹的智能车定位装置及方法 |
WO2020020701A1 (en) * | 2018-07-26 | 2020-01-30 | Signify Holding B.V. | Method for configuring a tracking system, tracking system, lighting system incorporating a tracking system and computer program |
CN111239715A (zh) * | 2020-01-13 | 2020-06-05 | 哈尔滨工业大学 | 一种联合灰色关联和神经网络的指纹定位方法 |
CN111882612A (zh) * | 2020-07-21 | 2020-11-03 | 武汉理工大学 | 一种基于三维激光检测车道线的车辆多尺度定位方法 |
Non-Patent Citations (10)
Title |
---|
"道路指纹"关键技术及其在智能车路系统中的应用;胡钊政;交通信息与安全;全文 * |
A novel method of WiFi fingerprint positioning using spatial multi-points matching;Wen Li;2016 International Conference on Indoor Positioning and Indoor Navigation (IPIN);全文 * |
Improved intelligent vehicle self-localization with integration of sparse visual map and high-speed pavement visual odometry;Huang, G (Huang, Gang);PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS PART D-JOURNAL OF AUTOMOBILE ENGINEERING;全文 * |
Nonminutiae Based Fingerprint Matching;Jaspreet Kour;2009 International Association of Computer Science and Information Technology - Spring Conference;全文 * |
Performance Analysis of Received Signal Strength Fingerprinting Based Distributed Location Estimation System for Indoor WLAN;Jain, VK (Jain, V. K.);WIRELESS PERSONAL COMMUNICATIONS;第70卷(第1期);全文 * |
Precise Location by Fingerprinting Road Segments with Variation of Broadcast Wave Reception;Kamiya, Y (Kamiya, Yuji);IEICE-INST ELECTRONICS INFORMATION COMMUNICATION ENGINEERS;全文 * |
一种基于K-L变换聚类的室内定位算法;金纯;邱灿;马金辉;陈光勇;;科学技术与工程(第20期);全文 * |
基于射频信号及信道状态信息的被动式目标定位方法研究;王举;中国博士学位论文全文数据库 信息科技辑;全文 * |
多类智能终端基于灰色关联度的WKNN定位效果分析;王威;地理空间信息;第19卷(第3期);全文 * |
顾及尺度变换的矢量地图道路网数据数字指纹嵌入点选取算法研究;陈维林;中国优秀硕士学位论文全文数据库 基础科学辑;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113721254A (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110119438B (zh) | 基于主动学习的机载LiDAR点云滤波方法 | |
CN108267141B (zh) | 道路点云数据处理系统 | |
CN113721254B (zh) | 一种基于道路指纹空间关联矩阵的车辆定位方法 | |
CN110542908A (zh) | 应用于智能驾驶车辆上的激光雷达动态物体感知方法 | |
CN113470090A (zh) | 基于sift-shot特征的多固态激光雷达外参标定方法 | |
CN116863357A (zh) | 一种无人机遥感堤坝影像定标与智能分割变化检测方法 | |
CN117078870A (zh) | 融合高精地图与激光稀疏点云的道路环境三维重建方法 | |
CN115867939A (zh) | 用于空中到地面配准的系统和方法 | |
CN110927765B (zh) | 激光雷达与卫星导航融合的目标在线定位方法 | |
Yao et al. | Automatic extraction of road markings from mobile laser-point cloud using intensity data | |
CN106709432B (zh) | 基于双目立体视觉的人头检测计数方法 | |
CN113838129B (zh) | 一种获得位姿信息的方法、装置以及系统 | |
CN113096181A (zh) | 设备位姿的确定方法、装置、存储介质及电子装置 | |
CN114387293A (zh) | 道路边缘检测方法、装置、电子设备及车辆 | |
CN117932333A (zh) | 一种顾及不同地形场景下的城市建筑高度提取方法 | |
CN117706577A (zh) | 基于激光雷达三维点云算法的船舶尺寸测量方法 | |
CN117053779A (zh) | 一种基于冗余关键帧去除的紧耦合激光slam方法及装置 | |
CN117629215A (zh) | 一种基于单线激光雷达点云配准的底盘充电回桩方法 | |
CN117029870A (zh) | 一种基于路面点云的激光里程计 | |
WO2018217193A1 (en) | Bayesian methodology for geospatial object/characteristic detection | |
CN116429126A (zh) | 一种动态环境下地面车辆点云地图快速重定位方法 | |
Tan et al. | Semantic segmentation of UAV lidar point clouds of a stack interchange with deep neural networks | |
Gu et al. | A spatial alignment method for UAV LiDAR strip adjustment in nonurban scenes | |
CN114063107A (zh) | 一种基于激光束的地面点云提取方法 | |
Wei et al. | Plane-based scan registration with moving vehicles exclusion |
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 |