CN101986102A - 一种城市地理信息系统中电子地图匹配方法 - Google Patents

一种城市地理信息系统中电子地图匹配方法 Download PDF

Info

Publication number
CN101986102A
CN101986102A CN 201010506525 CN201010506525A CN101986102A CN 101986102 A CN101986102 A CN 101986102A CN 201010506525 CN201010506525 CN 201010506525 CN 201010506525 A CN201010506525 A CN 201010506525A CN 101986102 A CN101986102 A CN 101986102A
Authority
CN
China
Prior art keywords
parameter
matched
segment
road section
executing
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
CN 201010506525
Other languages
English (en)
Other versions
CN101986102B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN2010105065254A priority Critical patent/CN101986102B/zh
Publication of CN101986102A publication Critical patent/CN101986102A/zh
Application granted granted Critical
Publication of CN101986102B publication Critical patent/CN101986102B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

本发明公开了一种城市地理信息系统中电子地图匹配方法,涉及地理信息系统领域,通过对电子地图的预处理,可使在电子地图匹配中需要参与计算的路段减少到预处理之前需要参与计算的总路段的50%以内;通过基于正方形网格的地图存储结构,使参与计算的路段进一步减少到总路段的1/N(N为网格的数量);进一步地,使得电子地图匹配的速度可以提高到1/2N,准确度可提高到90%以上。

Description

一种城市地理信息系统中电子地图匹配方法
技术领域
本发明涉及地理信息系统领域,特别涉及一种城市地理信息系统中电子地图匹配方法。
背景技术
随着城市地理信息系统的发展,对于电子地图匹配的要求越来越高。由于现有电子地图本身地图格式的限制及在绘制过程中的缺陷,使对电子地图的预处理成为在电子地图匹配之前的一项必要的工作。现有的电子地图处理方法主要集中在电子地图压缩、特征提取等方面,没有关注电子地图本身点线关系的处理。另一方面,现有电子地图匹配算法过于复杂、效率不高以及精度不高,不适宜在实际工程项目中使用。
发明内容
为了解决上述问题,减少计算复杂度、提高效率以及提高精度,本发明提供了一种城市地理信息系统中电子地图匹配方法,所述方法包括以下步骤:
(1)定义第一参数line和第二参数segment,并对第一参数line和第二参数segment分别编写ID,从Tab文件中读取初始数据,将初始数据保存在数据库中,对于lineID相同或邻接的第一参数line,如果有连续的第二参数segment斜率差小于等于第一阈值,则按照第一预设规则对连续的第二参数segment进行合并;
(2)将一个第一参数line中的任意一个顶点作为第一顶点,判断在其他第一参数line中是否查找到第一顶点,如果是,两个第一参数line邻接,流程结束;如果否,执行步骤(3);
(3)在其他第一参数line中查找距离第一顶点最近的第二顶点,根据第二预设规则判断在其他第一参数line中是否查找到第二顶点,如果是,第一顶点和第二顶点相同,两个第一参数line邻接,流程结束;如果否,建立第一参数line的拓扑结构;
(4)根据预先设定的车辆的行驶路线来获取道路网络范围,并将所述道路网络范围划分成m横n列的正方形网格,为每个所述正方形网格建立缓冲区,按照预设格式建立索引文件,为每辆车辆分别初始化已匹配路段列表和待匹配路段列表;
(5)接收车辆的GPS点的当前位置、时间、速度和行驶方向,判断速度是否为0,如果是,重新执行本步骤;如果否,执行步骤(6);
(6)根据GPS点的当前位置判断GPS点是否在所述道路网络范围内,如果是,根据第一计算公式获取GPS点所在的正方形网格的序号;如果否,根据获取到的正方形网格的序号,获取正方形网格内的所有第二参数segment;
(7)分别计算GPS点到每个第二参数segment的垂直距离,判断所述垂直距离是否大于第三阈值,如果是,抛弃相应的第二参数segment;如果否,执行步骤(8);
(8)判断速度是否小于第四阈值,如果是,根据第二计算公式获取当前路段的匹配权重值;如果否,执行步骤(9);
(9)判断行驶方向和当前第二参数segment的方向是否大于等于第五阈值,如果是,删除当前第二参数segment;如果否,执行步骤(10);
(10)判断GPS点到当前第二参数segment的距离是否小于所述第三阈值,如果是,根据GPS点到当前第二参数segment的距离、行驶方向和第三计算公式获取当前路段的匹配权重值;如果否,获取GPS点的第二参数候选集合,对匹配权重值进行顺序排列,根据GPS点在第二参数segment上的投影得到匹配位置;
(11)判断第二参数候选集合中的记录数是否为1,如果是,得到GPS点的匹配位置,执行步骤(12);如果否,执行步骤(15);
(12)判断所述待匹配路段列表中是否没有记录,如果是,将唯一匹配的第二参数segment加入到所述已匹配路段列表中,重新执行步骤(5);如果否,执行步骤(13);
(13)判断所述待匹配路段列表中最后一条记录与候选第二参数segment是否存在联通性,如果是,执行步骤(14);如果否,在所述待匹配路段列表中选匹配权值最小的第二参数segment作为匹配的第二参数,将匹配的第二参数segment放入所述已匹配路段列表中;
(14)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将所述待匹配路段列表中符合联通条件的第二参数segment作为匹配的第二参数加入到所述已匹配路段列表中,并在所述待匹配路段列表中,删除符合联通条件的第二参数segment;如果否,根据未匹配第二参数segment的匹配权值,选择所述匹配权值最小的第二参数segment加入到所述已匹配路段列表中,在所述待匹配路段列表中删除所述匹配权值最小的第二参数segment,将不满足联通条件的第二参数segment删除,重新执行步骤(5);
(15)判断候选第二参数segment和所述已匹配路段列表中的最后一条第二参数segment是否存在联通性,如果是,执行步骤(16);如果否,执行步骤(17);
(16)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将候选第二参数segment记录加入所述已匹配路段列表中,重新执行步骤(5);如果否,执行步骤(17);
(17)判断所述待匹配路段列表中的记录数是否为0,如果是,将候选第二参数segment记录加入所述待匹配路段列表中,重新执行步骤(5);如果否,执行步骤(18);
(18)判断所述待匹配路段列表中的记录与候选第二参数segment是否存在联通性,如果是,执行步骤(19);如果否,将所述待匹配路段列表中权值最小的一条加入到所述已匹配路段列表中,再删除所述待匹配路段列表中的相应记录,将剩余的候选第二参数segment加入到所述待匹配路段列表中;
(19)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将符合条件的候选第二参数segment加入到所述已匹配路段列表中,并将所述待匹配路段列表清空;如果否,将符合条件的候选第二参数segment保存在所述待匹配路段列表中,将不满足联通关系的第二参数segment删除,重新执行步骤(5)。
步骤(1)中的所述第一预设规则具体为:
选第一个segment的起点为合并后的segment的起点,最后一个segment的终点为合并后的segment的终点。
步骤(3)中的所述第二预设规则具体为:
第一顶点和第二顶点之间的距离小于第二阈值,则所述第一顶点和所述第二顶点相同。
步骤(4)中的所述预设格式为:
------------
N1,Count1,line001,line005,linexxx…
N2,Count2,line002,line 004,lineyyy…
......
Ni,Countj,line 002,line010,linezzz…
--------
其中,Ni为正方形网格的序号,N为正方形网格的数量,Countj为第j个正方形网格中第一参数line的数量,linexxx、lineyyy和linezzz为第一参数line的标识ID,每个正方形网格索引缓冲区内的第一参数数据。
步骤(6)中的所述第一计算公式具体为:
Ni=int[(y-y0)/CellSize]×n+int[(x0-x)/CellSize]+1
其中,(x0,y0)为GPS点的初始点坐标,(x,y)为GPS点的坐标,CellSize为正方形网格的边长。
步骤(8)中的所述第二计算公式具体为:
δ seg h = 100 d ‾ h
其中,δsegh为第h个第二参数segment的匹配权重值,
Figure BDA0000028171470000042
为GPS点距第二参数segment垂直距离的归一化值,
Figure BDA0000028171470000043
dh为GPS点到第二参数segmnet的垂直距离,ΔGPS为GPS点的平均误差。
步骤(10)中的所述第三计算公式具体为:
δ s eg h = 100 ( k d d ‾ h + k θ θ h ‾ )
其中,
Figure BDA0000028171470000045
为行驶方向与第二参数segment方向夹角的归一化值,θh为行驶方向与第二参数segment方向夹角|(θsegG)|,θseg为第二参数segment的斜率,θG为行驶方向的斜率,G为行驶方向,
Figure BDA0000028171470000046
kθ为行驶方向权值系数,kd为距离权值系数,且满足kθ+kd=1。
本发明提供的技术方案的有益效果是:
本发明提供了一种城市地理信息系统中电子地图匹配方法,通过对电子地图的预处理,可使在电子地图匹配中需要参与计算的路段减少到预处理之前需要参与计算的总路段的50%以内;通过基于正方形网格的地图存储结构,使参与计算的路段进一步减少到总路段的1/N(N为网格的数量);进一步地,使得电子地图匹配的速度可以提高到1/2N,准确度可提高到90%以上。
附图说明
图1为本发明提供的城市地理信息系统中电子地图匹配方法的流程图;
图2为本发明提供的索引结构示意图;
图3为本发明提供的索引结构的另一示意图;
图4为本发明提供的正方形网格CELL的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
为了解决上述问题,减少计算复杂度、提高效率以及提高精度,本发明提供了一种城市地理信息系统中电子地图匹配方法,在电子地图匹配过程中,需要电子地图上的点、线以及拓扑结构的相关信息,如路段的邻接信息。但是由于电子地图格式本身及电子地图在绘制过程中的问题,无法直接通过读取电子地图文件来得到这些信息,或者因为得到的电子地图拓扑会使电子地图匹配的计算量过大,需要对读取之后的信息进行必要的处理,使得到的拓扑结构能使电子地图匹配高效的运行。电子地图匹配是将车辆定位轨迹与电子地图中的道路网络联系起来,本发明实施例提出基于正方形网格的电子地图匹配方法,能够快速、准确的将车辆轨迹匹配到电子地图上。该方法主要包括以下内容,参见图1,详见下文描述:
101:定义第一参数line和第二参数segment,并对第一参数line和第二参数segment分别编写ID,从Tab文件中读取初始数据,将初始数据保存在数据库中,对于第一参数lineID相同或邻接的第一参数line,如果有连续的第二参数segment斜率差小于等于第一阈值,则按照第一预设规则对连续的第二参数segment进行合并;
其中,第一参数line表示一段道路,一个电子地图中有多个第一参数line,建立一个第一参数line表,每一段道路line都有自身的ID,例如:line1的ID为1、line2的ID为2、line3的ID为3等,具体实现时,还可以采用其他的编码方式,本发明实施例对此不做限制。第二参数segment是第一参数line的一部分,即第一参数line由至少一个第二参数segment组成,每一个第二参数segment都有自身的ID,例如:line1有3个segment,分别为segment1、segment2和segment3,则segment1的ID为10、segment2的ID为11和segment3的ID为13。具体实现时,可以根据实际应用的情况来划分line,得到符合实际数量要求的第二参数segment,具体实现时,本发明实施例对此不做限制。电子地图的经纬度坐标信息存储在第二参数segment这一层。其中,Tab文件是预先得到的,初始数据是一条一条的记录,每一条记录对应于一个第二参数Segment。
其中,第一阈值的设定根据实际应用中具体应用情况进行设定,具体实现时本发明实施例对此不做限制,本发明实施例以第一阈值为0.01为例进行说明。第一预设规则具体为:选第一个segment的起点为合并后的segment的起点,最后一个segment的终点为合并后的segment的终点。合并后的segment的斜率由合并后的起点和终点之间的斜率计算出来,合并后的segment的长度是所有被合并的segment的长度之和。现实中两个邻接的l ine是有公共顶点的,而在电子地图中有时会出现两个邻接的line没有公共顶点的情况,这时就需要对邻接的line在电子地图中没有公共顶点的问题进行处理,处理过程如下:
102:将一个第一参数l ine中的任意一个顶点作为第一顶点,判断在其他第一参数line中是否查找到第一顶点,如果是,两个第一参数line邻接,流程结束;如果否,执行步骤103;
103:在其他第一参数line中查找距离第一顶点最近的第二顶点,根据第二预设规则判断在其他第一参数line中是否查找到第二顶点,如果是,第一顶点和第二顶点相同,两个第一参数line邻接,流程结束;如果否,建立第一参数line的拓扑结构;
其中,第二预设规则具体为:第一顶点和第二顶点之间的距离小于第二阈值,则第一顶点和第二顶点相同,即两个第一参数line邻接。第二阈值的具体取值根据实际应用中的情况来确定,本发明实施例以3.6m为例进行说明,具体实现时,本发明实施例对此不做限制。
在line表中增加字段来表示每个line的连通关系,参见表1,给出了拓扑结构的索引表。
表1索引表
  字段   内容描述
  PrelID_1   进入该line的第一条line
  PrelID_2   进入该line的第二条line
  PrelID_3   进入该line的第三条line
  PrelID_5   进入该line的第五条line
  NextlID_1   驶出该line的第一条line
  NextlID_2   驶出该line的第二条line
  NextlID_3   驶出该line的第三条line
  NextlID_4   驶出该line的第四条line
  NextlID_5   驶出该line的第五条line
图2和图3中给出了索引表对应的结构示意图,通过图2和图3可以清楚的获知当前line的连接情况。
104:根据预先设定的车辆的行驶路线来获取道路网络范围,并将道路网络范围划分成m横n列的正方形网格,为每个正方形网格建立缓冲区,按照预设格式建立索引文件,为每辆车分别初始化已匹配路段列表和待匹配路段列表;
其中,本发明实施例建立基于正方形网格CELL的索引文件,索引文件的预设格式如下所示:
------------
N1,Count1,line001,line005,linexxx…
N2,Count2,line002,line 004,lineyyy…
......
Ni,Countj,line002,line010,linezzz…
--------
其中,参见图4,Ni为正方形网格的序号,N为正方形网格的数量,Countj为第j个正方形网格中第一参数line的数量,line xxx、lineyyy和linezzz为第一参数的标识ID,每个正方形网格索引缓冲区内的第一参数数据。图4中横纵坐标之间为道路网络范围,将道路网络范围划分为3*3的正方形网格,图中斜线部分表示缓冲区,图中共有9个网格,其中,每个正方形网格的边长为CellSize,图层中的第二参数segment的平均长度大约为250m,为减少每个正方形网格中的第二参数segment数量,建立平均长度4倍或更多倍数的正方形网格,本发明实施例优选边长为1000米×1000米正方形网格,缓冲区设置为正方形网格周围40米的区域。其中,已匹配路段列表用AML(Already Matched List)表示,待匹配路段列表用PML(Pending Matched List)表示,待匹配路段列表PML中的每个元素是车辆当前待匹配路段集合车辆匹配历史中的每个元素是单个车辆及其所有待匹配路段的集合。
105:接收车辆的GPS点的当前位置、时间、速度和行驶方向,判断速度是否为0,如果是,重新执行本步骤;如果否,执行步骤106;
其中,(x,y)为GPS点的当前位置的坐标。
106:根据GPS点的当前位置判断GPS点是否在道路网络范围内,如果是,根据第一计算公式获取GPS点所在的正方形网格的序号;如果否,根据获取到的正方形网格的序号,获取正方形网格内的所有第二参数segment;
其中,第一计算公式:
Ni=int[(y-y0)/CellSize]×n+int[(x0-x)/CellSize]+1        (1)
其中,(x0,y0)为GPS点的初始点坐标,(x,y)为GPS点的坐标,CellSize为正方形网格的边长。
107:分别计算GPS点到每个第二参数segment的垂直距离,判断垂直距离是否大于第三阈值,如果是,抛弃相应的第二参数segment;如果否,执行步骤108;
其中,第三阈值的设定根据实际应用中具体应用情况进行设定,具体实现时本发明实施例对此不做限制,本发明实施例以第三阈值为50米为例进行说明。
108:判断速度是否小于第四阈值,如果是,根据第二计算公式获取当前路段的匹配权重值;如果否,执行步骤109;
其中,第四阈值的设定根据实际应用中的具体应用情况进行设定,具体实现时本发明实施例对此不做限制,本发明实施例以第四阈值为5km/h为例进行说明,其中,第二计算公式为:
δ seg h = 100 d ‾ h - - - ( 2 )
其中,δsegh为第h个第二参数segment的匹配权重值,为GPS点距第二参数segment垂直距离的归一化值,
Figure BDA0000028171470000083
dh为GPS点到第二参数segmnet的垂直距离,ΔGPS为GPS点的平均误差,这里设定为15米。
109:判断行驶方向和当前第二参数segment的方向是否大于等于第五阈值,如果是,删除当前第二参数segment;如果否,执行步骤110;
其中,第五阈值的设定根据实际应用中的具体应用情况进行设定,具体实现时本发明实施例对此不做限制,本发明实施例以第五阈值为60度为例进行说明。
110:判断GPS点到当前第二参数segment的距离是否小于第三阈值,如果是,根据GPS点到当前第二参数segment的距离、行驶方向和第三计算公式获取当前路段的匹配权重值;如果否,获取GPS点的第二参数候选集合,对匹配权重值进行顺序排列,根据GPS点在第二参数segment上投影得到匹配位置;
其中,第三计算公式为:
δ seg h = 100 ( k d d ‾ h + k θ θ h ‾ ) - - - ( 3 )
其中,
Figure BDA0000028171470000092
为行驶方向与第二参数segment方向夹角的归一化值,θh为行驶方向与第二参数segment方向夹角|(θsegG)|,θseg为第二参数segment的斜率,θG为行驶方向的斜率,G为行驶方向,
Figure BDA0000028171470000093
kθ为行驶方向权值系数,kd为距离权值系数,且满足kθ+kd=1。根据实际的应用情况可以对kθ、kd值进行调整,本发明实施例将两个参数设置为:kd=0.6,kθ=04。
其中,参见表2,该顺序排列可以为从大到小的顺序排列或从小到大的顺序排列,或其它类型的排序方式,具体实现时,本发明实施例对此不做限制,本发明实施例以从大到小的顺序排列为例进行说明。
表2匹配权重排列顺序表
Figure BDA0000028171470000094
111:判断第二参数候选集合中的记录数是否为1,如果是,得到GPS点的匹配位置,执行步骤112;如果否,执行步骤115;
112:判断待匹配路段列表PML中是否没有记录,如果是,将唯一匹配的第二参数segment加入到已匹配路段列表AML中,重新执行步骤105;如果否,执行步骤113;
113:判断待匹配路段列表PML中最后一条记录与候选第二参数segment是否存在联通性,如果是,执行步骤114;如果否,在待匹配路段列表PML中选匹配权值最小的第二参数segment作为匹配的第二参数segment,将匹配的第二参数segment放入已匹配路段列表AML中;
114:判断符合联通条件的第二参数segment的数量是否等于1,如果是,将待匹配路段列表PML中符合联通条件的第二参数segment作为匹配的第二参数segment加入到已匹配路段列表AML中,并在待匹配路段列表PML中,删除符合联通条件的第二参数segment;如果否,根据未匹配第二参数segment的匹配权值,选择匹配权值最小的第二参数segment加入到已匹配路段列表AML中,在待匹配路段列表PML中删除匹配权值最小的第二参数segment,将不满足联通条件的第二参数segment删除,重新执行步骤105;
其中,若匹配权值最小的第二参数segment有多个时,则将该多个第二参数segment都加入到已匹配路段列表AML中,在待匹配路段列表PML中删除该多个匹配权值最小的第二参数segment。
115:判断候选第二参数segment和已匹配路段列表AML中的最后一条第二参数segment是否存在联通性,如果是,执行步骤116;如果否,执行步骤117;
116:判断符合联通条件的第二参数segment的数量是否等于1,如果是,将候选第二参数segment记录加入已匹配路段列表AML中,重新执行步骤105;如果否,执行步骤117;
117:判断待匹配路段列表PML中的记录数是否为0,如果是,将候选第二参数segment记录加入待匹配路段列表PML中,重新执行步骤105;如果否,执行步骤118;
118:判断待匹配路段列表PML中的记录与候选第二参数segment是否存在联通性,如果是,执行步骤119;如果否,将待匹配路段列表PML中权值最小的一条加入到已匹配路段列表AML中,再删除待匹配路段列表PML中的记录后将剩余的候选第二参数segment加入到待匹配路段列表PML中;
119:判断符合联通条件的第二参数segment的数量是否等于1,如果是,将符合联通条件的候选第二参数segment加入到已匹配路段列表AML中,并将待匹配路段列表PML清空;如果否,将符合联通条件的候选第二参数segment保存在待匹配路段列表PML中,将不满足联通关系的第二参数segment删除,重新执行步骤105。
例如:待匹配路段列表PML中存放的当前待匹配路段的候选segment有4个,分别为segment1、segment2、segment3和segment4,其中,权值最小的一条为segment1,则将segment1放入已匹配路段列表AML中,再删除待匹配路段列表PML中的segment1后,将剩余的候选segment2、segment3和segment4加入到待匹配路段列表PML中。
综上所述,本发明实施例提供了一种城市地理信息系统中电子地图匹配方法,通过对电子地图的预处理,可使在电子地图匹配中需要参与计算的路段减少到预处理之前需要参与计算的总路段的50%以内;通过基于正方形网格的地图存储结构,使参与计算的路段进一步减少到总路段的1/N(N为正方形网格的数量);进一步地,使得电子地图匹配的速度可以提高到1/2N,准确度可提高到90%以上。
本发明实施例以第一阈值小于等于0.01、第二阈值小于3.6m、第三阈值为50米、第四阈值为5km/h和第五阈值为60度为例进行了浮动车的实验验证,经过实验验证得到电子地图匹配速度可以提高到1/2N,准确度可提高到90%以上,验证了本发明实施例的可行性,满足了实际应用中的需要。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种城市地理信息系统中电子地图匹配方法,其特征在于,所述方法包括以下步骤:
(1)定义第一参数line和第二参数segment,并对第一参数line和第二参数segment分别编写ID,从Tab文件中读取初始数据,将初始数据保存在数据库中,对于lineID相同或邻接的第一参数line,如果有连续的第二参数segment斜率差小于等于第一阈值,则按照第一预设规则对连续的第二参数segment进行合并;
(2)将一个第一参数line中的任意一个顶点作为第一顶点,判断在其他第一参数line中是否查找到第一顶点,如果是,两个第一参数line邻接,流程结束;如果否,执行步骤(3);
(3)在其他第一参数line中查找距离第一顶点最近的第二顶点,根据第二预设规则判断在其他第一参数line中是否查找到第二顶点,如果是,第一顶点和第二顶点相同,两个第一参数line邻接,流程结束;如果否,建立第一参数line的拓扑结构;
(4)根据预先设定的车辆的行驶路线来获取道路网络范围,并将所述道路网络范围划分成m横n列的正方形网格,为每个所述正方形网格建立缓冲区,按照预设格式建立索引文件,为每辆车辆分别初始化已匹配路段列表和待匹配路段列表;
(5)接收车辆的GPS点的当前位置、时间、速度和行驶方向,判断速度是否为0,如果是,重新执行本步骤;如果否,执行步骤(6);
(6)根据GPS点的当前位置判断GPS点是否在所述道路网络范围内,如果是,根据第一计算公式获取GPS点所在的正方形网格的序号;如果否,根据获取到的正方形网格的序号,获取正方形网格内的所有第二参数segment;
(7)分别计算GPS点到每个第二参数segment的垂直距离,判断所述垂直距离是否大于第三阈值,如果是,抛弃相应的第二参数segment;如果否,执行步骤(8);
(8)判断速度是否小于第四阈值,如果是,根据第二计算公式获取当前路段的匹配权重值;如果否,执行步骤(9);
(9)判断行驶方向和当前第二参数segment的方向是否大于等于第五阈值,如果是,删除当前第二参数segment;如果否,执行步骤(10);
(10)判断GPS点到当前第二参数segment的距离是否小于所述第三阈值,如果是,根据GPS点到当前第二参数segment的距离、行驶方向和第三计算公式获取当前路段的匹配权重值;如果否,获取GPS点的第二参数候选集合,对匹配权重值进行顺序排列,根据GPS点在第二参数segment上的投影得到匹配位置;
(11)判断第二参数候选集合中的记录数是否为1,如果是,得到GPS点的匹配位置,执行步骤(12);如果否,执行步骤(15);
(12)判断所述待匹配路段列表中是否没有记录,如果是,将唯一匹配的第二参数segment加入到所述已匹配路段列表中,重新执行步骤(5);如果否,执行步骤(13);
(13)判断所述待匹配路段列表中最后一条记录与候选第二参数segment是否存在联通性,如果是,执行步骤(14);如果否,在所述待匹配路段列表中选匹配权值最小的第二参数segment作为匹配的第二参数,将匹配的第二参数segment放入所述已匹配路段列表中;
(14)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将所述待匹配路段列表中符合联通条件的第二参数segment作为匹配的第二参数加入到所述已匹配路段列表中,并在所述待匹配路段列表中,删除符合联通条件的第二参数segment;如果否,根据未匹配第二参数segment的匹配权值,选择所述匹配权值最小的第二参数segment加入到所述已匹配路段列表中,在所述待匹配路段列表中删除所述匹配权值最小的第二参数segment,将不满足联通条件的第二参数segment删除,重新执行步骤(5);
(15)判断候选第二参数segment和所述已匹配路段列表中的最后一条第二参数segment是否存在联通性,如果是,执行步骤(16);如果否,执行步骤(17);
(16)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将候选第二参数segment记录加入所述已匹配路段列表中,重新执行步骤(5);如果否,执行步骤(17);
(17)判断所述待匹配路段列表中的记录数是否为0,如果是,将候选第二参数segment记录加入所述待匹配路段列表中,重新执行步骤(5);如果否,执行步骤(18);
(18)判断所述待匹配路段列表中的记录与候选第二参数segment是否存在联通性,如果是,执行步骤(19);如果否,将所述待匹配路段列表中权值最小的一条加入到所述已匹配路段列表中,再删除所述待匹配路段列表中的相应记录,将剩余的候选第二参数segment加入到所述待匹配路段列表中;
(19)判断符合联通条件的第二参数segment的数量是否等于1,如果是,将符合条件的候选第二参数segment加入到所述已匹配路段列表中,并将所述待匹配路段列表清空;如果否,将符合条件的候选第二参数segment保存在所述待匹配路段列表中,将不满足联通关系的第二参数segment删除,重新执行步骤(5)。
2.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(1)中的所述第一预设规则具体为:
选第一个segment的起点为合并后的segment的起点,最后一个segment的终点为合并后的segment的终点。
3.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(3)中的所述第二预设规则具体为:
第一顶点和第二顶点之间的距离小于第二阈值,则所述第一顶点和所述第二顶点相同。
4.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(4)中的所述预设格式为:
------------
N1,Count1,line001,line005,linexxx…
N2,Count2,line002,line 004,lineyyy…
......
Ni,Countj,line 002,line010,linezzz…
 --------
其中,Ni为正方形网格的序号,N为正方形网格的数量,Countj为第j个正方形网格中第一参数line的数量,linexxx、lineyyy和linezzz为第一参数line的标识ID,每个正方形网格索引缓冲区内的第一参数数据。
5.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(6)中的所述第一计算公式具体为:
Ni=int[(y-y0)/CellSize]×n+int[(x0-x)/CellSize]+1
其中,(x0,y0)为GPS点的初始点坐标,(x,y)为GPS点的坐标,CellSize为正方形网格的边长。
6.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(8)中的所述第二计算公式具体为:
δ seg h = 100 d ‾ h
其中,δsegh为第h个第二参数segment的匹配权重值,为GPS点距第二参数segment垂直距离的归一化值,dh为GPS点到第二参数segmnet的垂直距离,ΔGPS为GPS点的平均误差。
7.根据权利要求1所述的城市地理信息系统中电子地图匹配方法,其特征在于,步骤(10)中的所述第三计算公式具体为:
δ seg h = 100 ( k d d ‾ h + k θ θ h ‾ )
其中,
Figure FDA0000028171460000045
为行驶方向与第二参数segment方向夹角的归一化值,θh为行驶方向与第二参数segment方向夹角|(θsegG)|,θseg为第二参数segment的斜率,θG为行驶方向的斜率,G为行驶方向,
Figure FDA0000028171460000046
kθ为行驶方向权值系数,kd为距离权值系数,且满足kθ+kd=1。
CN2010105065254A 2010-10-14 2010-10-14 一种城市地理信息系统中电子地图匹配方法 Expired - Fee Related CN101986102B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010105065254A CN101986102B (zh) 2010-10-14 2010-10-14 一种城市地理信息系统中电子地图匹配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010105065254A CN101986102B (zh) 2010-10-14 2010-10-14 一种城市地理信息系统中电子地图匹配方法

Publications (2)

Publication Number Publication Date
CN101986102A true CN101986102A (zh) 2011-03-16
CN101986102B CN101986102B (zh) 2012-01-11

Family

ID=43710467

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010105065254A Expired - Fee Related CN101986102B (zh) 2010-10-14 2010-10-14 一种城市地理信息系统中电子地图匹配方法

Country Status (1)

Country Link
CN (1) CN101986102B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634352A (zh) * 2015-03-02 2015-05-20 吉林大学 一种基于浮动车移动轨迹与电子地图融合的道路匹配方法
CN104833362A (zh) * 2015-05-18 2015-08-12 山东省计算中心(国家超级计算济南中心) 基于百度地图的线路固定车辆定位及行驶路径生成方法
CN105243863A (zh) * 2015-08-13 2016-01-13 南京理工大学 一种车辆实时测速监控系统
CN106225791A (zh) * 2016-08-03 2016-12-14 福建工程学院 一种基于网格划分的gps定位与道路匹配方法
CN106469505A (zh) * 2015-08-20 2017-03-01 方正国际软件(北京)有限公司 一种浮动车轨迹纠偏方法和装置
CN106980029A (zh) * 2016-01-15 2017-07-25 厦门雅迅网络股份有限公司 车辆超速判断方法及其系统
CN107644533A (zh) * 2017-10-27 2018-01-30 上海云砥信息科技有限公司 基于移动网络数据的高速公路虚拟断面车流量监测方法
CN109754631A (zh) * 2017-11-01 2019-05-14 腾讯科技(深圳)有限公司 车辆位置追踪的方法、装置、计算机设备及存储介质
CN111609860A (zh) * 2020-03-30 2020-09-01 北京拓明科技有限公司 同轨用户识别方法及装置
CN112050820A (zh) * 2020-09-02 2020-12-08 平安科技(深圳)有限公司 道路匹配方法、装置、电子设备及可读存储介质
CN113032507A (zh) * 2021-03-22 2021-06-25 湖南汽车工程职业学院 一种基于扑拓结构的地理信息云存储系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2322445A (en) * 1997-02-21 1998-08-26 Pellon Limited A terrestrial navigation system and method
US20020065603A1 (en) * 2000-11-30 2002-05-30 Nissan Motor Co., Ltd. Vehicle position calculation apparatus and method
CN1908588A (zh) * 2006-08-17 2007-02-07 北京航空航天大学 一种基于小网格路网组织结构的快速地图匹配方法
US20080091347A1 (en) * 2006-10-16 2008-04-17 Eric Tashiro Map matching method and apparatus for navigation system
WO2008117787A1 (ja) * 2007-03-27 2008-10-02 Nec Corporation マップマッチングシステム、マップマッチング方法およびプログラム
KR20080094992A (ko) * 2007-04-23 2008-10-28 엘지전자 주식회사 맵 매칭 방법
CN101324440A (zh) * 2008-07-29 2008-12-17 光庭导航数据(武汉)有限公司 基于预测思想的地图匹配方法
CN101424540A (zh) * 2008-12-12 2009-05-06 清华大学 一种车辆导航中基于动态虚拟路口的地图匹配方法
CN101482419A (zh) * 2008-01-11 2009-07-15 上海邮电设计院有限公司 一种基于a-gps和3g网络的车辆动态导航业务系统
CN101673282A (zh) * 2009-10-15 2010-03-17 杭州北海信息系统有限公司 电子地图自动定位系统方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2322445A (en) * 1997-02-21 1998-08-26 Pellon Limited A terrestrial navigation system and method
US20020065603A1 (en) * 2000-11-30 2002-05-30 Nissan Motor Co., Ltd. Vehicle position calculation apparatus and method
CN1908588A (zh) * 2006-08-17 2007-02-07 北京航空航天大学 一种基于小网格路网组织结构的快速地图匹配方法
US20080091347A1 (en) * 2006-10-16 2008-04-17 Eric Tashiro Map matching method and apparatus for navigation system
WO2008117787A1 (ja) * 2007-03-27 2008-10-02 Nec Corporation マップマッチングシステム、マップマッチング方法およびプログラム
KR20080094992A (ko) * 2007-04-23 2008-10-28 엘지전자 주식회사 맵 매칭 방법
CN101482419A (zh) * 2008-01-11 2009-07-15 上海邮电设计院有限公司 一种基于a-gps和3g网络的车辆动态导航业务系统
CN101324440A (zh) * 2008-07-29 2008-12-17 光庭导航数据(武汉)有限公司 基于预测思想的地图匹配方法
CN101424540A (zh) * 2008-12-12 2009-05-06 清华大学 一种车辆导航中基于动态虚拟路口的地图匹配方法
CN101673282A (zh) * 2009-10-15 2010-03-17 杭州北海信息系统有限公司 电子地图自动定位系统方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《中国博士学位论文全文数据库基础科学辑》 20080615 柳林 移动终端导航系统中地图匹配技术的研究与应用 A008-9 1-7 , 第6期 2 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634352B (zh) * 2015-03-02 2015-11-11 吉林大学 一种基于浮动车移动轨迹与电子地图融合的道路匹配方法
CN104634352A (zh) * 2015-03-02 2015-05-20 吉林大学 一种基于浮动车移动轨迹与电子地图融合的道路匹配方法
CN104833362A (zh) * 2015-05-18 2015-08-12 山东省计算中心(国家超级计算济南中心) 基于百度地图的线路固定车辆定位及行驶路径生成方法
CN105243863A (zh) * 2015-08-13 2016-01-13 南京理工大学 一种车辆实时测速监控系统
CN106469505A (zh) * 2015-08-20 2017-03-01 方正国际软件(北京)有限公司 一种浮动车轨迹纠偏方法和装置
CN106980029A (zh) * 2016-01-15 2017-07-25 厦门雅迅网络股份有限公司 车辆超速判断方法及其系统
CN106980029B (zh) * 2016-01-15 2021-11-02 厦门雅迅网络股份有限公司 车辆超速判断方法及其系统
CN106225791B (zh) * 2016-08-03 2019-09-20 福建工程学院 一种基于网格划分的gps定位与道路匹配方法
CN106225791A (zh) * 2016-08-03 2016-12-14 福建工程学院 一种基于网格划分的gps定位与道路匹配方法
CN107644533B (zh) * 2017-10-27 2021-02-12 上海云砥信息科技有限公司 基于移动网络数据的高速公路虚拟断面车流量监测方法
CN107644533A (zh) * 2017-10-27 2018-01-30 上海云砥信息科技有限公司 基于移动网络数据的高速公路虚拟断面车流量监测方法
CN109754631A (zh) * 2017-11-01 2019-05-14 腾讯科技(深圳)有限公司 车辆位置追踪的方法、装置、计算机设备及存储介质
CN109754631B (zh) * 2017-11-01 2021-07-20 腾讯科技(深圳)有限公司 车辆位置追踪的方法、装置、计算机设备及存储介质
CN111609860A (zh) * 2020-03-30 2020-09-01 北京拓明科技有限公司 同轨用户识别方法及装置
CN111609860B (zh) * 2020-03-30 2022-02-22 北京拓明科技有限公司 同轨用户识别方法及装置
CN112050820A (zh) * 2020-09-02 2020-12-08 平安科技(深圳)有限公司 道路匹配方法、装置、电子设备及可读存储介质
CN112050820B (zh) * 2020-09-02 2024-05-07 平安科技(深圳)有限公司 道路匹配方法、装置、电子设备及可读存储介质
CN113032507A (zh) * 2021-03-22 2021-06-25 湖南汽车工程职业学院 一种基于扑拓结构的地理信息云存储系统

Also Published As

Publication number Publication date
CN101986102B (zh) 2012-01-11

Similar Documents

Publication Publication Date Title
CN101986102B (zh) 一种城市地理信息系统中电子地图匹配方法
CN111209291B (zh) 一种利用众包感知地图更新高精度地图的方法及系统
CN104062671B (zh) 曲率约束的gnss浮动车地图匹配方法及装置
CN107195180B (zh) 一种基于电警数据的交通出行轨迹提取方法和装置
CN104575075B (zh) 一种基于北斗的城市路网车辆坐标校正方法及装置
CN110427360A (zh) 轨迹数据的处理方法、处理装置、处理系统及计算机程序产品
CN102968900B (zh) 一种对rfid交通数据进行处理的方法
CN104567906A (zh) 一种基于北斗的城市路网车辆路径规划方法及装置
CN108286979A (zh) 一种获取高精度导航路径数据的方法和装置以及导航系统、驾驶系统
CN111143440B (zh) 基于车辆轨迹的车辆停靠行为确定方法、确定系统及存储介质
CN108021686A (zh) 一种快速匹配电子地图中公交线路与道路网络的方法
CN112530158B (zh) 一种基于历史轨迹的路网补充方法
CN109035783A (zh) 一种基于公交gps轨迹的虚拟路网缺失路段自动识别方法
CN113989451A (zh) 高精地图构建方法、装置及电子设备
CN112947485A (zh) 路径规划方法及系统
CN104034337B (zh) 一种浮动车地理位置点的地图匹配方法及装置
CN109816982A (zh) 基于共享单车轨迹的虚拟路网非机动车道属性修正方法
CN111982135B (zh) 一种基于不同协议的地图格式之间的转换方法
CN110830915B (zh) 一种确定起点位置的方法及装置
CN117789066A (zh) 基于深度学习的无人机多视角影像钢轨线自动提取方法
CN111882913B (zh) 一种自动生成公交车辆线路走向的方法
CN101308205A (zh) 一种从轨道卫星定位数据中自动获取关键数据的方法
CN105717517B (zh) 一种车载北斗多模gnss高精度道路基础数据采集方法
CN117423000A (zh) 一种基于mars的道路纵断面线形智能重构的方法
CN111507503B (zh) 供水管网参数确定方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120111

Termination date: 20121014