CN101986102A - 一种城市地理信息系统中电子地图匹配方法 - Google Patents
一种城市地理信息系统中电子地图匹配方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000004891 communication Methods 0.000 claims description 4
- 101100247669 Quaranfil virus (isolate QrfV/Tick/Afghanistan/EG_T_377/1968) PB1 gene Proteins 0.000 description 6
- 101150025928 Segment-1 gene Proteins 0.000 description 6
- 101100242902 Thogoto virus (isolate SiAr 126) Segment 1 gene Proteins 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 5
- 101100242890 Quaranfil virus (isolate QrfV/Tick/Afghanistan/EG_T_377/1968) PA gene Proteins 0.000 description 4
- 101150027881 Segment-3 gene Proteins 0.000 description 4
- 101100242891 Thogoto virus (isolate SiAr 126) Segment 3 gene Proteins 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 101100242901 Quaranfil virus (isolate QrfV/Tick/Afghanistan/EG_T_377/1968) PB2 gene Proteins 0.000 description 3
- 101150082826 Segment-2 gene Proteins 0.000 description 3
- 101100194052 Thogoto virus (isolate SiAr 126) Segment 2 gene Proteins 0.000 description 3
- 101150091797 Segment-4 gene Proteins 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
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)中的所述第二计算公式具体为:
步骤(10)中的所述第三计算公式具体为:
其中,为行驶方向与第二参数segment方向夹角的归一化值,θh为行驶方向与第二参数segment方向夹角|(θseg-θG)|,θseg为第二参数segment的斜率,θG为行驶方向的斜率,G为行驶方向,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为例进行说明,其中,第二计算公式为:
其中,δsegh为第h个第二参数segment的匹配权重值,为GPS点距第二参数segment垂直距离的归一化值,dh为GPS点到第二参数segmnet的垂直距离,ΔGPS为GPS点的平均误差,这里设定为15米。
109:判断行驶方向和当前第二参数segment的方向是否大于等于第五阈值,如果是,删除当前第二参数segment;如果否,执行步骤110;
其中,第五阈值的设定根据实际应用中的具体应用情况进行设定,具体实现时本发明实施例对此不做限制,本发明实施例以第五阈值为60度为例进行说明。
110:判断GPS点到当前第二参数segment的距离是否小于第三阈值,如果是,根据GPS点到当前第二参数segment的距离、行驶方向和第三计算公式获取当前路段的匹配权重值;如果否,获取GPS点的第二参数候选集合,对匹配权重值进行顺序排列,根据GPS点在第二参数segment上投影得到匹配位置;
其中,第三计算公式为:
其中,为行驶方向与第二参数segment方向夹角的归一化值,θh为行驶方向与第二参数segment方向夹角|(θseg-θG)|,θseg为第二参数segment的斜率,θG为行驶方向的斜率,G为行驶方向,kθ为行驶方向权值系数,kd为距离权值系数,且满足kθ+kd=1。根据实际的应用情况可以对kθ、kd值进行调整,本发明实施例将两个参数设置为:kd=0.6,kθ=04。
其中,参见表2,该顺序排列可以为从大到小的顺序排列或从小到大的顺序排列,或其它类型的排序方式,具体实现时,本发明实施例对此不做限制,本发明实施例以从大到小的顺序排列为例进行说明。
表2匹配权重排列顺序表
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)中的所述第二计算公式具体为:
其中,δsegh为第h个第二参数segment的匹配权重值,为GPS点距第二参数segment垂直距离的归一化值,dh为GPS点到第二参数segmnet的垂直距离,ΔGPS为GPS点的平均误差。
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)
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)
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 | 杭州北海信息系统有限公司 | 电子地图自动定位系统方法 |
-
2010
- 2010-10-14 CN CN2010105065254A patent/CN101986102B/zh not_active Expired - Fee Related
Patent Citations (10)
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)
Title |
---|
《中国博士学位论文全文数据库基础科学辑》 20080615 柳林 移动终端导航系统中地图匹配技术的研究与应用 A008-9 1-7 , 第6期 2 * |
Cited By (18)
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 |