CN112082564B - 一种交互式地图匹配方法及系统 - Google Patents

一种交互式地图匹配方法及系统 Download PDF

Info

Publication number
CN112082564B
CN112082564B CN202010805405.8A CN202010805405A CN112082564B CN 112082564 B CN112082564 B CN 112082564B CN 202010805405 A CN202010805405 A CN 202010805405A CN 112082564 B CN112082564 B CN 112082564B
Authority
CN
China
Prior art keywords
matching
point
coordinate point
coordinate
map matching
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
Application number
CN202010805405.8A
Other languages
English (en)
Other versions
CN112082564A (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.)
Dongguan University of Technology
Original Assignee
Dongguan University of Technology
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 Dongguan University of Technology filed Critical Dongguan University of Technology
Priority to CN202010805405.8A priority Critical patent/CN112082564B/zh
Publication of CN112082564A publication Critical patent/CN112082564A/zh
Application granted granted Critical
Publication of CN112082564B publication Critical patent/CN112082564B/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/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

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种交互式地图匹配方法,包括:获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上;执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。本发明提供的交互式地图匹配方法,通过将人工计算和机器算法结合,提高了地图匹配结果的精度,同时结合特定的查询方法,有效减少人工标注,降低了人工作业的成本。

Description

一种交互式地图匹配方法及系统
技术领域
本发明涉及道路网络匹配及定位技术领域,尤其涉及一种交互式地图匹配方法、系统、终端设备和计算机可读存储介质。
背景技术
目前,许多基于GPS的车辆位置数据的手机基站定位系统在使用时都需要地图匹配,通过将GPS点与道路网络对齐来对其进行预处理。地图匹配的难点在于原始轨迹数据通常不包含运动对象的实际路径,并且切合实际的行走轨迹很难测试或评估任何的地图匹配算法。现有技术提供的地图匹配方法,往往存在以下问题:第一,通过人工手动收集地面实况花费成本高且效率低;第二,在采样点标注时,依赖人工手动将每个采样点标注到合适候选路段,不仅过程繁琐,操作复杂,且标注结果的精度无法得到保障;第三,现有的地图匹配系统在查询时无法利用注释者反馈及时更新,导致查询时经常出现无效查询,进一步增加了匹配作业的工作量,浪费时间。
发明内容
本发明的目的在于提供一种交互式地图匹配方法,通过将人工计算和机器算法结合,提高了地图匹配结果的精度,同时结合特定的查询方法,有效减少人工标注,降低了人工作业的成本。
为了克服上述现有技术中的缺陷,本发明实施例提供了一种交互式地图匹配方法,包括:
获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上;
执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。
进一步地,所述将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上,采用基于隐马尔可夫的地图匹配模型,包括:
计算坐标点pi匹配路段的发射概率,采用公式如下:
Figure GDA0003734158320000021
其中,δ是测量噪声的标准差,Pr(ei,j|pi)为坐标点与路段之间的最短垂直距离;
计算坐标点pi到pi+1的转变概率,采用公式如下:
Figure GDA0003734158320000022
其中,cdist(pi,pi+1)表示两点间的地球表面弧形距离,route(pi,pi+1)就表示沿着路段ei,j从点pi行驶到路段ei,j+1到达点pi+1的行驶距离;
根据所述发射概率和所述转变概率进行地图坐标匹配。
进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于维特比算法的交互式地图匹配算法,具体为:
计算坐标点pi的状态序列,采用公式如下:
Figure GDA0003734158320000023
Figure GDA0003734158320000024
其中,C(i,j)表示状态序列pi=(e1,j1,e2,j2,…,ei,ji)以ei,j为最终状态时的最高概率。
处理最后一个观察值时,递归终止。然后通过以下公式(也称为后向公式),从最后一个隐藏状态(每个步骤中的最高概率C(i,j))的最后一个隐藏状态反向检索产生C(i,j)的最佳状态序列:
Figure GDA0003734158320000031
进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于可疑坐标点的路程选择算法,具体为:
根据香农熵和每个坐标点匹配到不同路段的发射概率计算所述坐标点的不确定性,具体公式如下:
Figure GDA0003734158320000032
筛选出不确定性最小的坐标点作为下一个待匹配目标点。
进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于置信度的算法,具体为:
将坐标点的发射概率和转变概率的乘积作为置信度,采用公式如下:
Figure GDA0003734158320000033
进一步地,所述置信度的算法还包括动态置信度的选择方法,具体为:
在当前坐标点被匹配时,判断所述坐标点的前一个匹配点是否影响后续匹配的不确定性,若否,则继续进行下一个坐标点匹配;若是,对所述当前坐标点重新匹配。
进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于稳定性的算法,具体为:
计算影响当前坐标点匹配结果的其他所有坐标点的稳定性,采用公式如下:
S(pi)=|Dpi|;
其中,Dpi是一个集合,包含了所有对点pi没有影响的点,S(pi)为该集合的绝对值;
筛选S(pi)值最小的可疑点,将所述可疑点反馈给标注者进行再次考核。
本发明实施例还提供一种交互式地图匹配系统,包括:
初始坐标匹配单元,用于获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上;
可疑点选择单元,用于执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
匹配结果更新单元,用于响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。
本发明实施例还提供了一种计算机终端设备,包括:
一个或多个处理器;
存储器,与所述处理器耦接,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如上所述的交互式地图匹配方法。
本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行实现如上所述的交互式地图匹配方法。
相对于现有技术,本发明实施例至少存在如下有益效果:
(1)框架以迭代的方式将人类的工作与算法结合在一起,以实现较高的地图匹配精度;
(2)本发明提供的地图匹配算法基于HIMM,该算法可以考虑任意数量的人的标注,其输入信息可以大大提高准确性。
(3)提出了不同的查询选择方法,以有效减少需要人工标注的工作量。与传统方法相比,HIMM查询选择策略最多可以减少44%的查询;
(4)使用现实世界真实数据集、综合数据集和合成轨迹数据集对系统进行实验并分析实验结果,显着地降低了人工的成本。
附图说明
图1是本发明某一实施例提供的交互式地图匹配方法的流程示意图;
图2是本发明某一实施例提供的交互式地图匹配系统的总体框架图;
图3是本发明某一实施例提供的交互式地图匹配模型的结构示意图;
图4是本发明某一实施例提供的基于可疑点选择策略的几种选择结果图;
图5是本发明某一实施例提供的交互式地图匹配系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应当理解,文中所使用的步骤编号仅是为了方便描述,不对作为对步骤执行先后顺序的限定。
应当理解,在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样,除非上下文清楚地指明其它情况,否则单数形式的“一”、“一个”及“该”意在包括复数形式。
术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。
第一方面:
请参阅图1,本发明某一实施例提供了一种交互式地图匹配方法,包括:
S10、获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上;
本步骤中,需要说明的是,系统利用地图匹配算法生成初始路径,并将该路径以及原始轨迹显示在数字地图上;其中,轨迹是由一组包含经度、纬度、时间戳的位置点所构成,但是这些点由于噪音测试和样例采集的偏差也会变得不精确。地图匹配的基本难点在于,原始轨迹数据(尤其是在被动收集位置信息时)通常不包含运动对象的实际路径。并且贴切实际的行走轨迹,就很难测试或评估任何的地图匹配算法。本步骤中,将坐标点匹配到数字地图上时,采用交互式的地图匹配方式,利用地图匹配的算法将轨迹中所有的点匹配到地图中对应的路段上,其中,路段指的是处于两个路口之间径直的一段路。
S20、执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
在本步骤中,会针对可疑度设置一阈值,针对匹配的结果,将可疑度与该阈值进行比较,若小于或等于该阈值,则认为该坐标点匹配成功;否则,将其筛选出来反馈至用户终端,用户再重新匹配,直至匹配完所有的坐标点;其中,可疑度指的是坐标点错误匹配的概率,坐标点错误匹配的概率越大,则可疑度越高,二者呈正相关关系。
S30、响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。
本步骤中,在匹配完成之后,用户端会对坐标点的匹配结果进行查询,以更新匹配结果。因为系统会在每次人工标注之后继续更新显示路径,并利用注释者反馈的信息来改善匹配结果,在查询时可以直接以坐标点为关键字进行查询或者以相关注释为关键字进行查询。
本发明实施例提供的一种交互式地图匹配方法,通过将人工计算和机器算法结合,提高了地图匹配结果的精度,同时结合特定的查询方法,有效减少人工标注,降低了人工作业的成本。
第二方面:
请参阅图2,本发明某一实施例提供了一种交互式地图匹配系统的总体框架图,在一开始,需要取出需要地图匹配的轨迹,那么HIMM就用其地图匹配的算法将轨迹中所有的点匹配到地图中对应的路段上(路段指的是处于两个路口之间径直的一段路)然后HIMM会将匹配的结果绘制到数字地图上来让标注者直观的观察轨迹匹配的正确与否,如果标注者觉得轨迹匹配有问题,那么标注者会接收到一个来自疑点选择策略发送过来的最有可能错误匹配的坐标点,接着标注者就可以选择重新匹配该坐标点或者认为该坐标点没有问题。如此循环往复直到我们认为图中所有的点都匹配正确,此项工作就此完成。
第三方面:
本发明某一实施例提供的一种交互式地图匹配方法中,采用基于隐马尔可夫的地图匹配模型,首先任意点p有不同的可能性被匹配到地图上的路段,我们把坐标点p匹配到一条路段的概率成为发射概率,而每一条有可能与坐标点对应上的路段成为点p的一个隐状态,计算坐标点pi匹配路段的发射概率,采用公式如下:
Figure GDA0003734158320000071
其中,δ是测量噪声的标准差,pdist(ei,j|pi)为坐标点与路段之间的最短垂直距离;
在实际情况中,轨迹在行进间是有可能改变路段的,所以还涉及到一个转变概率的概念,当沿着路段ei,j行驶到路段ei,j+1的距离要小于pi与pi+1之间的地球表面弧形距离时,就意味着从pi到pi+1的转变概率是比较大的。计算坐标点pi到pi+1的转变概率,采用公式如下:
Figure GDA0003734158320000072
其中,cdist(pi,pi+1)表示两点间的地球表面弧形距离,route(pi,pi+1)就表示沿着路段ei,j从点pi行驶到路段ei,j+1到达点pi+1的行驶距离;
然后根据所述发射概率和所述转变概率进行地图坐标匹配,其中利用这两个概率进行交互式地图匹配的模型如图3所示。
本实施例中,进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于维特比算法的交互式地图匹配算法,具体为:
计算坐标点pi的状态序列,采用公式如下:
Figure GDA0003734158320000081
Figure GDA0003734158320000082
其中,C(i,j)表示状态序列pi=(e1,j1,e2,j2,…,ei,ji)以ei,j为最终状态时的最高概率。
处理最后一个观察值时,递归终止。然后通过以下公式(也称为后向公式),从最后一个隐藏状态(每个步骤中的最高概率C(i,j))的最后一个隐藏状态反向检索产生C(i,j)的最佳状态序列:
Figure GDA0003734158320000083
与维特比算法相似,我们的交互式地图匹配算法使用动态编程以递归的方式快速找到最佳状态序列。当算法计算局部最优概率C(i,j)时,对于每个隐藏状态ei,j作为第一次i观测的最终状态,它首先检查pi是否由注释者人工匹配的。如果是的话,算法会排除pi的所有候选隐藏状态,除了ei,j,由注释者选择(即<pi,ei,j>∈M’)通过修改发射概率Pr(ei,j|pi)为1,所有其他发射概率Pr(eik|pi)为0。否则,发射概率即为Pr(ei,j|pi)。这样,就保证了所有为0的C(i,k)(其中1≤k≤r,k≠j)与反向公式选择了包含对于特定C(i,j)的ei,j
本实施例中,还提供了基于路程的策略,即基于可疑坐标点的路程选择算法,具体为:
根据香农熵和每个坐标点匹配到不同路段的发射概率计算所述坐标点的不确定性,具体公式如下:
Figure GDA0003734158320000091
筛选出不确定性最小的坐标点作为下一个待匹配目标点。
需要说明的是,当要进行坐标点与路段之间的匹配时,首先想到的思路应该是基于距离的方法,但也会出现一些不确定的情况。例如图4(a)中的情况,p1非常地接近路段e1,并且p3非常地接近路段e7,所以p1就应该匹配到e1,p3就应该匹配到e7,这两点毋庸置疑。但考察到p2会发现该点与路段e2,e3,e6的距离几乎一样,这就为p2的匹配带来很多不确定性。在测量不确定性策略中,有一个概念叫香农熵。香农熵与概率分布有密切的关系。所以我们可以用每一个点匹配到不同路段的发射概率来计算一个点的不确定性。之后,再找出不确定性最小的坐标点来作为下一个要匹配的目标,这就是基于距离的可疑坐标点选择策略。
本实施例中,还提供了基于置信度的策略,即基于置信度的算法,具体为:
将坐标点的发射概率和转变概率的乘积作为置信度,采用公式如下:
Figure GDA0003734158320000092
需要说明的是,在基于距离的策略中,只考虑了点与路段之间的距离关系,而忽略了点与点之间的拓扑关系,一个错的点拓扑会导致匹配到的路段发生绕路。具体看图4(b),点p2到路段e2和e4的距离几乎一样,但对于点p3,p3是在路段e5和路段e3当中离e5是比较近的,那么这个时候在点p2处就会存在不确定性。如果p2错误地被匹配到e4,那么后面的轨迹就绕一段路(先经过e8再到e5)。所以为了将点的拓扑性也利用起来,定于了点匹配的置信度,这个置信度是所有点可能的发射概率和转变概率的乘积。可以看出来置信度不仅将点匹配到不同路段的概率考虑进来,并且将下一段可能的行驶路段概率也考虑进来了。这就使得点与点之间的拓扑更具有合理性。
进一步地,所述置信度的算法还包括动态置信度的选择方法,具体为:
在当前坐标点被匹配时,判断所述坐标点的前一个匹配点是否影响后续匹配的不确定性,若否,则继续进行下一个坐标点匹配;若是,对所述当前坐标点重新匹配。
需要说明的是,在应用置信度策略之前有一个前提就是在匹配一个点的时候,其他的点是没有被匹配到地图中的。有一种潜在的问题就是我们在匹配某个点的时候可能会受其他点的影响。考虑这样一种,如图4(c),点p1,p2,p3看似都可以被匹配到路段e1或e2,如果这时有一个点被错误地匹配了并且其不确定性还很低,意味着该点很大几率不再重新考核,那这个时候,图4(c)中不确定性大的三个点就会一一被考核,若p2被对应到e2,那么点p1和p3的不确定性就会降低,但实际上这两个点肯定是有一个不能被匹配到e1的,这会使得后续匹配出错。所以我们需要一种策略当一个点被匹配时还要综合考虑之前已经匹配过的点,并判断是否已匹配过的点会不会影响后续匹配的不确定性。在每次匹配过后,都要重新考量已匹配的点的不确定性。
进一步地,所述基于隐马尔可夫的地图匹配模型还包括,基于稳定性的算法,具体为:
计算影响当前坐标点匹配结果的其他所有坐标点的稳定性,采用公式如下:
S(pi)=|Dpi|;
其中,Dpi是一个集合,包含了所有对点pi没有影响的点,S(pi)为该集合的绝对值;
筛选S(pi)值最小的可疑点,将所述可疑点反馈给标注者进行再次考核。
其中,在实际中轨迹中的点是会相互影响的,而那些容易被其他点影响的坐标点在做匹配工作时,应该综合考量对其有影响的所有点的情况。具体如图4(d),点p1,p2,p3,p4被匹配到路段e1或e2的可能性几乎相同,当我们要对p1进行匹配时,可以看出其他三个点的有无对p1匹配到路段e1或e2的概率不会有什么影响。但对于点p5来说,若p6或p7其中一个去掉之后,就会对p5的匹配有很大的影响。这说明p5的稳定性更低。所以重新考察p5的可能性会远远高于考察p1的可能性。不确定度测量的另一个策略是稳定性。Dpi是一个集合,集合中包含了所有对点pi没有影响的点,可以看出S(pi)越高代表着对其有影响的点越少。稳定性策略即是要选出S(pi)值最小的可疑点反馈给标注着进行再考核。
第四方面:
请参阅图5,本发明某一实施例还提供了一种交互式地图匹配系统,包括:
初始坐标匹配单元01,用于获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上;
可疑点选择单元02,用于执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
匹配结果更新单元03,用于响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中所述查询选择请求携带修改后的坐标点和/或查询注释。
第五方面:
本发明某一实施例还提供了一种计算机终端设备,包括:一个或多个处理器;存储器,与所述处理器耦接,用于存储一个或多个程序;当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如上所述的交互式地图匹配方法。
处理器用于控制该计算机终端设备的整体操作,以完成上述的全自动用电量预测方法的全部或部分步骤。存储器用于存储各种类型的数据以支持在该计算机终端设备的操作,这些数据例如可以包括用于在该计算机终端设备上操作的任何应用程序或方法的指令,以及应用程序相关的数据。该存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,例如静态随机存取存储器(Static Random Access Memory,简称SRAM),电可擦除可编程只读存储器(Electrically Erasable Programmable Read-OnlyMemory,简称EEPROM),可擦除可编程只读存储器(Erasable Programmable Read-OnlyMemory,简称EPROM),可编程只读存储器(Programmable Read-Only Memory,简称PROM),只读存储器(Read-Only Memory,简称ROM),磁存储器,快闪存储器,磁盘或光盘。
计算机终端设备可以被一个或多个应用专用集成电路(Application Specific1ntegrated Circuit,简称AS1C)、数字信号处理器(Digital Signal Processor,简称DSP)、数字信号处理设备(Digital Signal Processing Device,简称DSPD)、可编程逻辑器件(Programmable Logic Device,简称PLD)、现场可编程门阵列(Field ProgrammableGate Array,简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现,用于执行如上述任一项实施例所述的交互式地图匹配方法,并达到如上述方法一致的技术效果。
第六方面:
本发明某一实施例还提供了一种包括程序指令的计算机可读存储介质,该程序指令被处理器执行时实现如上述任一项实施例所述的交互式地图匹配方法的步骤。例如,该计算机可读存储介质可以为上述包括程序指令的存储器,上述程序指令可由计算机终端设备的处理器执行以完成如上述任一项实施例所述的交互式地图匹配方法,并达到如上述方法一致的技术效果。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。

Claims (8)

1.一种交互式地图匹配方法,其特征在于,包括:
获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上,包括:基于维特比算法的交互式地图匹配算法,具体为:
计算坐标点pi的状态序列,采用公式如下:
Figure FDA0003773684590000011
Figure FDA0003773684590000012
其中,C(i,j)表示状态序列pi=(e1,j1,e2,j2,…,ei,ji)以ei,j为最终状态时的最高概率;
处理最后一个观察值时,递归终止,然后通过以下公式,从最后一个隐藏状态反向检索产生C(i,j)的最佳状态序列:
Figure FDA0003773684590000013
执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。
2.根据权利要求1所述的交互式地图匹配方法,其特征在于,所述将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上,采用基于隐马尔可夫的地图匹配模型,包括:
计算坐标点pi匹配路段的发射概率,采用公式如下:
Figure FDA0003773684590000021
其中,δ是测量噪声的标准差,pdist(ei,j,pi)为坐标点与路段之间的最短垂直距离;
计算坐标点pi到pi+1的转变概率,采用公式如下:
Figure FDA0003773684590000022
其中,cdist(pi,pi+1)表示两点间的地球表面弧形距离,route(pi,pi+1)就表示沿着路段ei,j从点pi行驶到路段ei,j+1到达点pi+1的行驶距离;
根据所述发射概率和所述转变概率进行地图坐标匹配。
3.根据权利要求2所述的交互式地图匹配方法,其特征在于,所述基于隐马尔可夫的地图匹配模型还包括,基于可疑坐标点的路程选择算法,具体为:
根据香农熵和每个坐标点匹配到不同路段的发射概率计算所述坐标点的不确定性,具体公式如下:
Figure FDA0003773684590000023
筛选出不确定性最小的坐标点作为下一个待匹配目标点。
4.根据权利要求2所述的交互式地图匹配方法,其特征在于,所述基于隐马尔可夫的地图匹配模型还包括,基于置信度的算法,具体为:
将坐标点的发射概率和转变概率的乘积作为置信度,采用公式如下:
Figure FDA0003773684590000024
5.根据权利要求4所述的交互式地图匹配方法,其特征在于,所述置信度的算法还包括动态置信度的选择方法,具体为:
在当前坐标点被匹配时,判断所述坐标点的前一个匹配点是否影响后续匹配的不确定性,若否,则继续进行下一个坐标点匹配;若是,对所述当前坐标点重新匹配。
6.一种交互式地图匹配系统,其特征在于,包括:
初始坐标匹配单元,用于获取用户终端的移动轨迹,并将所述移动轨迹中所有的坐标点匹配到数字地图的初始路径上,包括:基于维特比算法的交互式地图匹配算法,具体为:
计算坐标点pi的状态序列,采用公式如下:
Figure FDA0003773684590000031
Figure FDA0003773684590000032
其中,C(i,j)表示状态序列pi=(e1,j1,e2,j2,…,ei,ji)以ei,j为最终状态时的最高概率;
处理最后一个观察值时,递归终止,然后通过以下公式,从最后一个隐藏状态反向检索产生C(i,j)的最佳状态序列:
Figure FDA0003773684590000033
可疑点选择单元,用于执行可疑点选择策略,在匹配结果中筛选可疑度大于预设阈值的坐标点,反馈至所述用户终端;其中,所述可疑度为坐标点错误匹配的概率;
匹配结果更新单元,用于响应所述用户终端对某一坐标点的查询选择请求,更新匹配结果;其中,所述查询选择请求携带修改后的坐标点和/或查询注释。
7.一种计算机终端设备,其特征在于,包括:
一个或多个处理器;
存储器,与所述处理器耦接,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1至5任一项所述的交互式地图匹配方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行实现如权利要求1至5任一项所述的交互式地图匹配方法。
CN202010805405.8A 2020-08-11 2020-08-11 一种交互式地图匹配方法及系统 Active CN112082564B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010805405.8A CN112082564B (zh) 2020-08-11 2020-08-11 一种交互式地图匹配方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010805405.8A CN112082564B (zh) 2020-08-11 2020-08-11 一种交互式地图匹配方法及系统

Publications (2)

Publication Number Publication Date
CN112082564A CN112082564A (zh) 2020-12-15
CN112082564B true CN112082564B (zh) 2022-09-06

Family

ID=73729157

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010805405.8A Active CN112082564B (zh) 2020-08-11 2020-08-11 一种交互式地图匹配方法及系统

Country Status (1)

Country Link
CN (1) CN112082564B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113254562B (zh) * 2021-06-18 2021-09-28 长安大学 一种高效gps轨迹地图匹配方法
CN113506369A (zh) * 2021-07-13 2021-10-15 阿波罗智能技术(北京)有限公司 用于生成地图的方法、装置、电子设备和介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308047A (zh) * 2013-05-10 2013-09-18 杭州师范大学 一种具有人工在线修正功能的移动机器人定位方法
CN110427360A (zh) * 2019-06-28 2019-11-08 东软集团股份有限公司 轨迹数据的处理方法、处理装置、处理系统及计算机程序产品
CN110717350A (zh) * 2018-07-11 2020-01-21 沈阳美行科技有限公司 一种行车轨迹的校正方法及校正装置
CN111121791A (zh) * 2019-11-29 2020-05-08 上饶市中科院云计算中心大数据研究院 隐马尔可夫模型在地图匹配中的优化方法及gps定位方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10415984B2 (en) * 2017-12-29 2019-09-17 Uber Technologies, Inc. Measuring the accuracy of map matched trajectories

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308047A (zh) * 2013-05-10 2013-09-18 杭州师范大学 一种具有人工在线修正功能的移动机器人定位方法
CN110717350A (zh) * 2018-07-11 2020-01-21 沈阳美行科技有限公司 一种行车轨迹的校正方法及校正装置
CN110427360A (zh) * 2019-06-28 2019-11-08 东软集团股份有限公司 轨迹数据的处理方法、处理装置、处理系统及计算机程序产品
CN111121791A (zh) * 2019-11-29 2020-05-08 上饶市中科院云计算中心大数据研究院 隐马尔可夫模型在地图匹配中的优化方法及gps定位方法

Also Published As

Publication number Publication date
CN112082564A (zh) 2020-12-15

Similar Documents

Publication Publication Date Title
US10323948B2 (en) GPS data repair
CN112082564B (zh) 一种交互式地图匹配方法及系统
EP3961489A1 (en) Method and apparatus for identifying updated road, device and computer storage medium
Chawathe Segment-based map matching
CN109974718B (zh) 地图匹配方法、装置、设备和介质
JP2022501696A (ja) 品質スコアの勾配を使用した反復的なタンパク質構造予測
Koller et al. Fast hidden Markov model map-matching for sparse and noisy trajectories
US11205079B2 (en) Determining position data
CN110006439B (zh) 地图轨迹数据的匹配方法、装置、服务器及存储介质
CN109376761A (zh) 一种地址标识及其经纬度的挖掘方法及装置
TW201700955A (zh) 路徑規劃方法及裝置
CN111242148A (zh) 用于地图绑定的人工智能系统和方法
CN105205196A (zh) 用于生成路网的方法和系统
US11972578B2 (en) Method and system for object tracking using online training
CN107256230A (zh) 一种基于多元化地理信息点的融合方法
Millard-Ball et al. Map-matching poor-quality GPS data in urban environments: the pgMapMatch package
Luo et al. Incremental route inference from low-sampling GPS data: An opportunistic approach to online map matching
CN111881908A (zh) 目标检测模型的修正方法、检测方法、装置、设备及介质
CN114547223A (zh) 轨迹预测方法、轨迹预测模型的训练方法及装置
CN105516392A (zh) Ip地址信息库更新方法及装置
CN116026352A (zh) 定位信号的处理方法、装置、介质及电子设备
CN111078760B (zh) 货源搜索方法、装置、设备及存储介质
Dogramadzi et al. Accelerated map matching for GPS trajectories
CN115265555A (zh) 基于隐马尔科夫的多噪声感知的地图匹配校正方法及系统
CN111782748B (zh) 地图检索方法、信息点poi语义向量的计算方法和装置

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