CN109459045A - 一种针对低频gps轨迹的改进交互式投票匹配方法 - Google Patents

一种针对低频gps轨迹的改进交互式投票匹配方法 Download PDF

Info

Publication number
CN109459045A
CN109459045A CN201811147858.5A CN201811147858A CN109459045A CN 109459045 A CN109459045 A CN 109459045A CN 201811147858 A CN201811147858 A CN 201811147858A CN 109459045 A CN109459045 A CN 109459045A
Authority
CN
China
Prior art keywords
candidate
point
gps
weight
path
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
CN201811147858.5A
Other languages
English (en)
Other versions
CN109459045B (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201811147858.5A priority Critical patent/CN109459045B/zh
Publication of CN109459045A publication Critical patent/CN109459045A/zh
Application granted granted Critical
Publication of CN109459045B publication Critical patent/CN109459045B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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
    • G01C21/30Map- or contour-matching
    • G01C21/32Structuring or formatting of map data

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种针对低频GPS轨迹的改进交互式投票地图路径匹配方法,本发明不仅考虑距离特征、道路的拓扑结构以及路段的限速,还考虑了每个GPS的点实时移动方向和速度,以提高匹配准确率;另外,本发明中还加入了滤波器,通过约束条件去除候选噪声路段,以提高匹配效率。本发明具有匹配准确率高、效率高和鲁棒性强的优点。

Description

一种针对低频GPS轨迹的改进交互式投票匹配方法
技术领域
本算发明属于智能交通中的地图匹配技术领域,涉及一种地图匹配算法,具体涉及一种处理低频GPS轨迹的方法。
背景技术
随着内置GPS传感器的电子移动设备的普及,每天会产生大量的行驶轨迹。但是,由于卫星定位***的局限性,尤其是在卫星能见度有限、高楼林立的城市化地区,卫星信号被遮挡和折射,导致定位数据的不准确和丢失。所以我们通过电子设备采集的定位数据是不精确的。此外,由于设备电量、存储和网络传输带宽的局限,实际收集到的大多数轨迹数据都是低采样率的(例如:采样间隔为1分钟或以上)。携带GPS设备的移动对象可能在一分钟内行驶了相当长的距离,移动对象在两个GPS点之间通过的路段会有多种可能性。所以,在利用这些轨迹数据进行分析和应用之前,必须要对他们进行预处理过程,即地图匹配。
在过去的近二十年里,已有各种地图匹配算法相继提出,根据采用方法的不同,可分为:简单的几何和拓扑算法,涉及模糊逻辑、卡尔曼滤波、隐马尔科夫、条件随机场、神经网络、遗传算法、蚁群优化和机器学习等的高级算法。现有的这些地图匹配算法中大多数都专注于处理高采样率的GPS轨迹。相对而言,低采样率GPS轨迹的地图匹配比高采样率轨迹的匹配准确性降低,尤其是在城市复杂的路网中。
发明内容
本发明针对现有的地图匹配算法无法兼顾匹配准确率和效率的问题,提出了一种针对低频GPS轨迹的改进交互式投票匹配方法。本方法不仅考虑距离特征、道路的拓扑结构以及路段的限速,还考虑了每个GPS的点实时移动方向和速度,以提高算法的匹配准确率;另外,算法中还加入了一个滤波器,通过约束条件去除候选噪声路段,以提高算法的匹配效率。
一种针对低频GPS轨迹的改进交互式投票匹配方法,实现具体包括以下步骤:
步骤一:以一个GPS点为圆心,在路网中搜索半径范围内的所有路段作为候选路径;
步骤二:通过三个约束条件过滤错误的候选路径;
约束一:定义方向差阈值,如果车辆实时运动方向与对应路段方向之间的方向差值保留对应的候选点,并使用方向分析函数对候选路段分配权重,否则,当时,直接将对应的候选节点滤除,后续的处理过程将不再考虑此候选点;
约束二:考虑通过分析车辆平均速度与最短路径的权重限速的关系,将速度超过合理范围的对应候选路段过滤;定义参数α,1<α<2,当时,,将视为错误的匹配路径,其中为两个连续候选点之间的车辆平均速度,为他们之间的最短路径的权重限速;如果采样点i-1的所有候选点到pi的候选点全为错误匹配路径,则将GPS点i的候选点所对应的候选路段作为噪声路段过滤;
约束三:设置一个超速极限值将部分不符合条件的候选路段过滤;定义参数,设vi为GPS点pi的实时速度,vj是对应候选路段的限速,当vi>(1+γ)·vj时,我们将对应的候选路段rj作为噪声路段滤除;
步骤三:根据距离特征、道路的拓扑结构、路段的限速和GPS点的实时移动方向,利用分析函数分配候选点之间的权重,并建立一个具有权重的候选图;然后通过权重候选图建立一个反映两个相邻GPS点的候选点之间转移可能性的静态得分矩阵;
步骤四:加权分析在步骤三的基础上,考虑全局GPS轨迹点之间的相互影响关系,且相互影响的强弱程度取决于GPS点之间的距离,最后建立n个加权得分矩阵;n表示GPS点个数;
步骤五:在得到加权矩阵后,对于每个候选点,通过最大权重概率方法获得一条经过它的局部最优路径,统计在所有局部路径中每个GPS点对应候选点出现的次数,最后每个GPS点对应计数次数最多的候选点组成全局最优路径。
本发明的有益效果:
1.本方法在候选点选择阶段对路网建立基于Rtree的新型索引,提高了搜索候选路径的速度。
2.本方法通过加入基于GPS实时方向的分析函数。通过考虑汽车方向与道路方向的关系,提高了匹配准确率。
3.本方法中加入了滤波器,通过约束条件去除候选路段中的部分噪声,提高方法匹配准确率和效率。
本发明具有匹配准确率高、效率高和鲁棒性好的优点。
附图说明
图1为本发明框架图;
图2为本发明候选点选择示意图;
图3为本发明查找局部最优路径图。
具体实施方式
本发明体系结构如图1所示。它由以下几部分组成:候选点准备、位置环境分析、加权分析和交互式投票。
1.候选点准备
针对GPS轨迹P=(pi|i=1,2,3…N),选取在路网中以GPS点为圆心,半径r范围内的所有路段作为候选路段k表示GPS点pi的第k个候选路段,选取在路段上与GPS点最近的点作为候选点。如图2所示,GPS点p1拥有一条候选路段r3,对应有一个候选点GPS点p2拥有四条候选路段,路段r2和r3对应的候选点落在路段范围内,而路段r1和r4对应的候选点落在路段对应的起始节点。
2.位置环境分析
位置环境分析包括两部分:(1)约束分析;(2)时空分析。
(1)约束分析
约束一:定义方向差阈值,如果车辆实时运动方向与对应路段方向之间的方向差值保留对应的候选点,并使用方向分析函数对候选路段分配权重,否则,当时,直接将对应的候选节点滤除,后续的处理过程将不再考虑此候选点。阈值的设置可根据经验或统计数据进行设置。
约束二:现实中车辆在路网中的运动速度都是在一个合理的范围内的。根据这一事实,我们可考虑通过分析车辆平均速度与最短路径的权重限速的关系,将速度超过合理范围的对应候选路段过滤。定义参数α(1<α<2),当时(1<α<2),将视为错误的匹配路径,其中为两个连续候选点之间的车辆平均速度,为他们之间的最短路径的权重限速。如果采样点i-1的所有候选点到pi的候选点全为错误匹配路径,则将GPS点i的候选点所对应的候选路段作为噪声路段过滤。
约束三:现实中,车辆一般被允许超速15%左右,根据这一事实我们还可以设置一个超速极限值将部分不符合条件的候选路段过滤。定义参数,设vi为GPS点pi的实时速度,vj是对应候选路段的限速,当vi>(1+γ)·vj时,我们将对应的候选路段rj作为噪声路段滤除。
(2)时空分析
时空分析是考虑几何拓扑和道路限速信息,通过数学函数式来确定候选点的权重的分析过程。观测概率是由GPS与候选路段的距离和GPS实时方向与候选路段的夹角决定,状态转移概率是由道路的拓扑结构决定。此外,算法还加入了速度分析函数,通过道路限速信息寻找与GPS速度相似的道路。
一般而言,测量误差(距离误差和方向误差)的分布满足高斯分布N(μ,σ2),所以计算观测概率其中是GPS点pi的候选点,与pi的欧几里得距离,是GPS实时方向和对应候选路段之间的角度差。
状态转移概率是评估两个连续候选点之间的最短路径和直线路径的相似性,算法中计算状态转移概率为其中d(i-1,i)是两个相邻GPS采样点pi-1,pi之间的欧几里得距离,w(i-1,s),(i,t)是两个GPS点对应的两个候选之间的最短路径距离。
我们定义一个时间分析函数,考虑结合采样速度与路段的限速来计算权重,反映车辆的平均行驶速度与路段限速的相似性,倾向于帮助我们选择行驶速度与限速最为相似的路段。算法中将时间分析函数定义为 其中是候选点到候选点的平均速度,是最短路径的权重限速。
这一部分最后输出一个结合空间分析和时间分析后的权重候选图。节点是候选点集,边是两个相邻候选点之间最短路径的集合,节点和边所有的权重值都是基于位置和道路分析的结果。
3.加权分析
加权分析分析通过考虑全局GPS轨迹点之间的相互影响关系,在位置环境分析的基础上进行权重影响建模,即对静态矩阵的加权,它由两部分组成:(1)建立静态得分矩阵;(2)加权影响建模。
(1)建立静态得分矩阵
位置环境分析产生的静态矩阵表示为M=diag{M2,M3,M4,…,Mn},其中ai和ai-1分别代表第个和
i-1个GPS点的候选点个数权重影响建模。
(2)加权影响建模
考虑到两个GPS点之间距离影响成反比的关系,我们定义了一个反比权重函数对各个GPS点之间的影响进行分析。最后,每一个GPS点都有一个对应的距离加权矩阵:
4.交互式投票
交互式投票是通过加权分析得到的权重候选图构建全局最优路径的过程,它由两部分组成:(1)寻找局部路径;(2)全局投票。
(1)寻找局部最优路径
每一个候选点必定有一条经过它的路径,并称之为的局部最优路径。通过枚举,整个候选图中共有条局部最优路径,其中ai表示GPS点pi候选点的个数。例如:图3所示的一条轨迹的候选图中,我们假设的权重为-∞,这也就意味着整条路径必须经过每个候选点的累计权重计算为对于每个候选点,计算其累计权重后,记录得到该候选点累计权重的前一个候选点最后,计算完所有采样点对应的候选点的累计权重之后,由最后一个采样点的候选节点的累计权重得到一个最大累计权重值,即 图3中,最大累计权重值它来自于p4的三个候选点。最终,可通过全局最大权重值,以及每个最大权重值记录的其前一采样点对应的候选点的位置,得到权重最大的局部最优路径,即路径:
(2)全局投票
每找到一条局部最优路径,我们会对路径中出现的候选点进行计数。完成所有的局部路径查找后,我们拥有了一组拥有计数值的候选点。全局投票就是将每个GPS点中计数值最大的候选点取出来组成全局最优路径。

Claims (1)

1.一种针对低频GPS轨迹的改进交互式投票匹配方法,其特征在于,该方法具体包括以下步骤:
步骤一:以一个GPS点为圆心,在路网中搜索半径r范围内的所有路段作为候选路径;
步骤二:通过三个约束条件过滤错误的候选路径;
约束一:定义方向差阈值θ,如果车辆实时运动方向与对应路段方向之间的方向差值保留对应的候选点,并使用方向分析函数对候选路段分配权重,否则,当时,直接将对应的候选节点滤除,后续的处理过程将不再考虑此候选点;
约束二:考虑通过分析车辆平均速度与最短路径的权重限速的关系,将速度超过合理范围的对应候选路段过滤;定义参数α,1<α<2,当时,,将视为错误的匹配路径,其中为两个连续候选点之间的车辆平均速度,为他们之间的最短路径的权重限速;如果采样点pi-1的所有候选点到pi的候选点全为错误匹配路径,则将GPS点pi的候选点所对应的候选路段作为噪声路段过滤;
约束三:设置一个超速极限值将部分不符合条件的候选路段过滤;定义参数γ,设vi为GPS点pi的实时速度,vj是对应候选路段的限速,当vi>(1+γ)·vj时,我们将对应的候选路段rj作为噪声路段滤除;
步骤三:根据距离特征、道路的拓扑结构、路段的限速和GPS点的实时移动方向,利用分析函数分配候选点之间的权重,并建立一个具有权重的候选图;然后通过权重候选图建立一个反映两个相邻GPS点的候选点之间转移可能性的静态得分矩阵;
步骤四:加权分析在步骤三的基础上,考虑全局GPS轨迹点之间的相互影响关系,且相互影响的强弱程度取决于GPS点之间的距离,最后建立n个加权得分矩阵;n表示GPS点个数;
步骤五:在得到加权矩阵后,对于每个候选点,通过最大权重概率方法获得一条经过它的局部最优路径,统计在所有局部路径中每个GPS点对应候选点出现的次数,最后每个GPS点对应计数次数最多的候选点组成全局最优路径。
CN201811147858.5A 2018-09-29 2018-09-29 一种针对低频gps轨迹的改进交互式投票匹配方法 Active CN109459045B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811147858.5A CN109459045B (zh) 2018-09-29 2018-09-29 一种针对低频gps轨迹的改进交互式投票匹配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811147858.5A CN109459045B (zh) 2018-09-29 2018-09-29 一种针对低频gps轨迹的改进交互式投票匹配方法

Publications (2)

Publication Number Publication Date
CN109459045A true CN109459045A (zh) 2019-03-12
CN109459045B CN109459045B (zh) 2020-10-09

Family

ID=65607140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811147858.5A Active CN109459045B (zh) 2018-09-29 2018-09-29 一种针对低频gps轨迹的改进交互式投票匹配方法

Country Status (1)

Country Link
CN (1) CN109459045B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110727749A (zh) * 2019-09-30 2020-01-24 杭州电子科技大学 一种自适应候选路段搜索方法
CN110736471A (zh) * 2019-09-24 2020-01-31 武汉大学 基于隐马尔可夫模型的低频浮动车轨迹数据路网匹配方法
CN110986965A (zh) * 2019-11-28 2020-04-10 武汉大学 基于隐马尔可夫模型的低频浮动车轨迹数据路网匹配方法
CN111027783A (zh) * 2019-12-27 2020-04-17 江苏欣网视讯软件技术有限公司 一种用于路径选取的全局动态交互决策算法
CN111307164A (zh) * 2020-02-25 2020-06-19 西北大学 一种低采样率轨迹地图匹配方法
CN111694032A (zh) * 2020-05-09 2020-09-22 杭州电子科技大学 一种基于聚类的大规模轨迹数据的快速地图匹配方法
CN111982141A (zh) * 2020-07-31 2020-11-24 长安大学 一种面向低频车辆轨迹数据进行路径推断的方法、设备及存储介质
CN112988929A (zh) * 2021-02-18 2021-06-18 同济大学 基于全局合作式投票的地图匹配方法、***、介质及终端
CN113295173A (zh) * 2021-05-24 2021-08-24 安徽师范大学 环形路段的地图匹配方法
CN114879234A (zh) * 2021-10-14 2022-08-09 电子科技大学 一种复杂gps轨迹中的重要地点挖掘方法及装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20240003705A1 (en) * 2022-07-01 2024-01-04 Here Global B.V. Method and apparatus for providing a path-based map matcher

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120296560A1 (en) * 2011-05-19 2012-11-22 Microsoft Corporation Inferring a Behavioral State of a Vehicle
CN103235848A (zh) * 2013-04-15 2013-08-07 中国科学院软件研究所 一种基于简化路网模型的轻量级路网匹配方法
CN107103775A (zh) * 2017-05-18 2017-08-29 西安理工大学 一种基于群智计算的道路质量检测方法
CN107358629A (zh) * 2017-07-07 2017-11-17 北京大学深圳研究生院 一种基于目标识别的室内建图与定位方法
CN107911786A (zh) * 2017-10-24 2018-04-13 星际空间(天津)科技发展有限公司 一种基于路网校正的混合室内定位方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120296560A1 (en) * 2011-05-19 2012-11-22 Microsoft Corporation Inferring a Behavioral State of a Vehicle
CN103235848A (zh) * 2013-04-15 2013-08-07 中国科学院软件研究所 一种基于简化路网模型的轻量级路网匹配方法
CN107103775A (zh) * 2017-05-18 2017-08-29 西安理工大学 一种基于群智计算的道路质量检测方法
CN107358629A (zh) * 2017-07-07 2017-11-17 北京大学深圳研究生院 一种基于目标识别的室内建图与定位方法
CN107911786A (zh) * 2017-10-24 2018-04-13 星际空间(天津)科技发展有限公司 一种基于路网校正的混合室内定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YAYING ZHANG EL: ""An Advanced Interactive-Voting Based Map Matching Algorithm for Low-sampling-rate GPS data"", 《IEEE 15TH INTERNATIONAL CONFERENCE》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110736471B (zh) * 2019-09-24 2021-08-03 武汉大学 基于隐马尔可夫模型的低频浮动车轨迹数据路网匹配方法
CN110736471A (zh) * 2019-09-24 2020-01-31 武汉大学 基于隐马尔可夫模型的低频浮动车轨迹数据路网匹配方法
CN110727749A (zh) * 2019-09-30 2020-01-24 杭州电子科技大学 一种自适应候选路段搜索方法
CN110986965A (zh) * 2019-11-28 2020-04-10 武汉大学 基于隐马尔可夫模型的低频浮动车轨迹数据路网匹配方法
CN111027783A (zh) * 2019-12-27 2020-04-17 江苏欣网视讯软件技术有限公司 一种用于路径选取的全局动态交互决策算法
CN111307164A (zh) * 2020-02-25 2020-06-19 西北大学 一种低采样率轨迹地图匹配方法
CN111307164B (zh) * 2020-02-25 2022-12-02 西北大学 一种低采样率轨迹地图匹配方法
CN111694032A (zh) * 2020-05-09 2020-09-22 杭州电子科技大学 一种基于聚类的大规模轨迹数据的快速地图匹配方法
CN111982141A (zh) * 2020-07-31 2020-11-24 长安大学 一种面向低频车辆轨迹数据进行路径推断的方法、设备及存储介质
CN111982141B (zh) * 2020-07-31 2022-09-13 长安大学 一种面向低频车辆轨迹数据进行路径推断的方法、设备及存储介质
CN112988929A (zh) * 2021-02-18 2021-06-18 同济大学 基于全局合作式投票的地图匹配方法、***、介质及终端
CN113295173A (zh) * 2021-05-24 2021-08-24 安徽师范大学 环形路段的地图匹配方法
CN113295173B (zh) * 2021-05-24 2023-08-29 安徽师范大学 环形路段的地图匹配方法
CN114879234A (zh) * 2021-10-14 2022-08-09 电子科技大学 一种复杂gps轨迹中的重要地点挖掘方法及装置

Also Published As

Publication number Publication date
CN109459045B (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN109459045A (zh) 一种针对低频gps轨迹的改进交互式投票匹配方法
CN106969764B (zh) 一种道路匹配方法、装置及车载地图采集***
CN102102992B (zh) 基于多级网络划分的匹配道路初筛方法及地图匹配***
CN106574975B (zh) 使用***信号的轨迹匹配
CN106023587B (zh) 基于多信息融合的轨迹数据路网精确匹配方法
CN106781478B (zh) 一种基于lte信令数据的轨迹跟踪方法
CN107766808A (zh) 道路网络空间中车辆对象移动轨迹聚类的方法及***
CN106714110A (zh) 一种Wi‑Fi位置指纹地图自动构建方法及***
CN105241465B (zh) 一种道路更新的方法
CN105957342A (zh) 基于众包时空大数据的车道级道路测图方法及***
CN114005280A (zh) 一种基于不确定性估计的车辆轨迹预测方法
CN105260988A (zh) 一种高精地图数据的处理方法和装置
CN110095127A (zh) 一种基于分段的隐马尔可夫模型地图匹配方法
JP2007333385A (ja) 不動物位置記録装置
Davidson et al. Application of particle filters to a map-matching algorithm
CN112101527B (zh) 识别车道变化的方法和装置、电子设备和存储介质
CN108253976A (zh) 一种充分借助车辆航向的三阶段在线地图匹配算法
CN111429744A (zh) 一种路况信息融合分析的方法、装置及存储介质
CN113743469A (zh) 一种融合多源数据及综合多维指标的自动驾驶决策方法
CN110189362A (zh) 基于多分支自编码对抗网络的高效目标跟踪方法
Zhang et al. An advanced interactive-voting based map matching algorithm for low-sampling-rate GPS data
CN113985406A (zh) 一种海上雷达目标航迹拼接方法
CN111562605B (zh) 自适应的gps错误观测值识别方法
CN113295173A (zh) 环形路段的地图匹配方法
CN116975184A (zh) 一种车辆轨迹数据处理方法、装置、设备和可读存储介质

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190312

Assignee: Nanjing Qishengyun Information Technology Co.,Ltd.

Assignor: HANGZHOU DIANZI University

Contract record no.: X2021330000844

Denomination of invention: An improved interactive voting matching method for low frequency GPS trajectory

Granted publication date: 20201009

License type: Common License

Record date: 20211225

EE01 Entry into force of recordation of patent licensing contract