CN111459953A - 多源数据融合的地图匹配更新方法和系统 - Google Patents

Abstract

本发明公开了多源数据融合的地图匹配更新方法，包括以下步骤：对源地图的道路信息矢量化处理生成矢量化数据源；对基础地图处理生成多层次空间索引分区；将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型，使得所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地图数据；将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和修复生成新基础地图完成更新。本发明还公开了多源数据融合的地图匹配更新系统。本发明多源数据融合的地图匹配更新方法和系统，实现了对多源数据的处理，并完成了将多源数据更新到离线地图，实现了高效的、准确率高的、对数据属性依赖较小的地图匹配更新。

Description

多源数据融合的地图匹配更新方法和系统

技术领域

本发明涉及地理信息技术领域，具体涉及多源数据融合的地图匹配更新方法和系统。

背景技术

复杂离线地图精确匹配和更新任务最初源自实际应用需求，要求实现将GPS导航地图中 的道路与PGIS中的道路进行匹配，并对GPS导航地图中的独有的道路(PGIS由于道路更新 周期长等原因，无法及时绘制新建道路)更新到PGIS平台的基础地图中，以支撑公安和消 防应急救援的导航规划等。然而由于数据来源多样，导致数据具有测绘标准不统一、相对误 差大、个别地图数据属性维度低等特点，使得离线地图的道路匹配难度较大；同时在地图更 新时，将新道路更新到既有的地图中时保证其拓扑关系正确也是一个难题。因此需要设计一 种高效的、准确率高的、对数据属性依赖较小的地图匹配更新算法以满足大规模的离线复杂 地图的匹配和更新。

发明内容

本发明所要解决的技术问题是现有技术中由于数据来源多样，导致数据具有测绘标准不 统一、相对误差大、个别地图数据属性维度低等特点，使得离线地图的道路匹配难度较大， 目的在于提供多源数据融合的地图匹配更新方法和系统，解决上述问题。

本发明通过下述技术方案实现：

多源数据融合的地图匹配更新方法，包括以下步骤：S1：对源地图的道路信息矢量化处 理生成矢量化数据源；对基础地图处理生成多层次空间索引分区；S2：将所述矢量化数据源 和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型，使得所述矢量化数据源匹 配到所述多层次空间索引分区中生成预更新基础地图数据；其中所述矢量化数据源为所述路 网匹配模型的观察状态，所述多层次空间索引分区为所述路网匹配模型的隐含状态；S3：将 预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和修复生成新基础地图完成更 新。

本发明应用时，源地图一般为外部来源的数据地图，例如互联网GPS导航地图，而基础 地图为需要根据源地图进行更新的地图，例如公安地理信息系统(PGIS)的道路底图。将源地 图的道路信息矢量化处理可以便于后续数据输入，无论源地图的来源为何种，都可以将这种 来源化为同样的数据格式，而基础地图处理生成的多层次空间索引分区中，一般每一条道路 为最底层的空间，多条道路形成上一层空间，多个区域形成更上一层空间，以此类推直至形 成多层次空间索引分区。在基础地图和源地图的匹配中，所述矢量化数据源中每条轨迹匹配 所述多层次空间索引分区中道路的计算方法都是独立且相同的，所以可以采用并行计算的方 式进行运算。

在本申请中通过模型进行道路匹配时是已知GPS坐标点序列，即源地图的数据，计算其 在道路网络中的真实经过路径的过程，即基础地图的数据。我们假设GPS坐标点序列为观察 状态，真实经过路径为隐含状态，基于隐马尔科夫(HMM)的路网匹配模型则是在给定一系 列观察序列的前提下，寻找最有可能产生这个观察序列的隐含状态序列，也就是实际经过的 道路。然后需要检查缺失道路信息和拓扑关系信息并完成更新。本发明通过上述步骤，实现 了对多源数据的处理，并完成了将多源数据更新到离线地图，实现了高效的、准确率高的、 对数据属性依赖较小的地图匹配更新。

进一步的，步骤S1包括以下子步骤：将所述源地图的每条道路提取成独立道路轨迹，并 生成包含每条道路原始属性项的数据矩阵作为所述矢量化数据源；对所述基础地图建立基于 R-Tree的空间索引分区作为所述多层次空间索引分区；所述多层次空间索引分区中，叶节点 为基础地图中的道路，任意一个父节点为该父节点所有子节点的集合直至根节点。

进一步的，步骤S2包括以下子步骤：将所述矢量化数据源和所述多层次空间索引分区输 入基于隐马尔科夫的路网匹配模型后获取所述路网匹配模型的观察概率和转移概率；根据所 述观察概率和所述转移概率寻找并匹配与所述矢量化数据源最大似然的所述多层次空间索引 分区中的道路。

进一步的，所述路网匹配模型的观察概率采用下式获取：

式中，p(o_t,i|c_t,i)为所述观察概率，o_t,i为所述矢量化数据源中的离散坐标点序列，c_t,i为 离散坐标点在候选路段上的投影点，d_t,i为观察距离，u为距离均值，σ_z为距离标准差；

所述路网匹配模型的转移概率采用下式获取：

式中，p(d_t,θ_t)为所述转移概率，d_t为所述矢量化数据源中轨迹点之间距离与其在候选道 路上的投影点间距离之差，根据d_t近似符合指数分布的特征，拟合了β作为转移概率的系数， θ_t为所述矢量化数据源中相邻两个轨迹点构成的向量与其在候选路段上的投影点构成的向量 之间的夹角。

进一步的，步骤S2还包括以下子步骤：

当所述矢量化数据源无法准确匹配于所述多层次空间索引分区中的道路时，在保证候选 道路连通的情况下弹性扩大缓冲半径或跳过矢量化数据源中该轨迹点；

扩大缓冲半径或跳过的轨迹点的累计长度不超过轨迹总长度的5％。

进一步的，步骤S3包括以下子步骤：

将所述矢量化数据源中转移概率和/或观察概率为零的道路定义为缺失道路；

用所述多层次空间索引分区中的道路匹配所述矢量化数据源，并标记所述矢量化数据源 中无法被所述多层次空间索引分区中的道路匹配到的轨迹作为缺失道路；

用所述矢量化数据源检测所述多层次空间索引分区中的道路，并标记所述矢量化数据源 中无法匹配所述多层次空间索引分区中的道路的轨迹作为缺失道路；

进一步的，步骤S3还包括以下子步骤：

获取所述缺失道路在所述矢量化数据源中的关联道路，并根据下式获取所述缺失道路和 所述关联道路的夹角：

θ＝min(arccos(v₁,v₂),π-arccos(v₁,v₂))

式中，v₁为缺失道路在关联点最近两个顶点坐标组成的向量，v₂为关联道路中离关联点 最相邻的两个坐标点组成的向量，θ为所述缺失道路和所述关联道路的夹角；

根据缺失道路与关联道路的位置关系和θ的值将缺失道路分为已相交、延长相交和泛平行 三种类型；

其中已相交为所述缺失道路与所述关联道路已有明确交点；

延长相交为所述缺失道路与所述关联道路无明确交点，但是所述关联道路与所述缺失道 路的延长线上存在交点，且

泛平行为

进一步的，步骤S3还包括以下子步骤：

当所述缺失道路的类型为已相交时，在交点处将多余部分截断删除完成更新；

当所述缺失道路的类型为延长相交时，将所述缺失道路沿端点方向延长，并在交点处截 断，删去多余部分完成更新；

当所述缺失道路的类型为泛平行时，获取所述缺失道路中首个拐角

的点作为接入点； 将所述缺失道路的端点在所述关联道路上的投影点连接所述接入点生成替换道路，将所述缺 失道路的端点到所述接入点的道路作为待替换道路；将所述替换道路替换所述待替换道路完 成更新。

多源数据融合的地图匹配更新系统，包括：

ETL单元：用于对源地图的道路信息矢量化处理生成矢量化数据源；对基础地图处理生 成多层次空间索引分区；

模型单元：用于将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的 路网匹配模型，使得所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地 图数据；其中所述矢量化数据源为所述路网匹配模型的观察状态，所述多层次空间索引分区 为所述路网匹配模型的隐含状态；

检测修复单元：用于将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和 修复生成新基础地图完成更新。

进一步的，所述ETL单元将所述源地图的每条道路提取成独立道路轨迹，并生成包含每 条道路原始属性项的数据矩阵作为所述矢量化数据源；

所述ETL单元对所述基础地图建立基于R-Tree的空间索引分区作为所述多层次空间索引 分区；所述多层次空间索引分区中，叶节点为基础地图中的道路，任意一个父节点为该父节 点所有子节点的集合直至根节点。

所述模型单元将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路 网匹配模型后获取所述路网匹配模型的观察概率和转移概率；所述模型单元根据所述观察概 率和所述转移概率寻找并匹配与所述矢量化数据源最大似然的所述多层次空间索引分区中的 道路。

本发明与现有技术相比，具有如下的优点和有益效果：

1、本发明多源数据融合的地图匹配更新方法和系统，实现了对多源数据的处理，并完成 了将多源数据更新到离线地图，实现了高效的、准确率高的、对数据属性依赖较小的地图匹 配更新；

2、本发明多源数据融合的地图匹配更新方法和系统，对传统的HMM模型的转移概率进 行了改进，使其适用于低维地图数据，从而不依赖速度、方向、限速等高维数据；

3、本发明多源数据融合的地图匹配更新方法和系统，在匹配时，对GPS轨迹异常点进 行了容差处理，使其在较大位置偏差的情况下依然具有较好的鲁棒性。地图更新则是在地图 匹配的基础上，定义源地图中无法成功匹配且满足一定条件的道路为缺失道路，本发明创造 的基于拓扑和几何修正的道路更新方法，能够简单高效地为基础地图提供高质量道路更新功 能；

4、地图匹配和更新中广泛采用了R-Tree空间索引技术、向量化计算以及多任务并行计 算等方法，在保证了匹配更新准确率的前提下，显著提升了运行效率。解决了多源地图融合 更新的准确率低、运行速度慢的难题。

附图说明

此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不 构成对本发明实施例的限定。在附图中：

图1为本发明系统结构示意图；

图2为本发明实施例示意图；

图3为本发明实施例示意图；

图4为本发明实施例示意图；

图5为本发明实施例示意图；

图6为本发明实施例示意图；

图7为本发明实施例示意图；

图8为本发明实施例示意图；

图9为本发明实施例示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明 作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本 发明的限定。

实施例

如图1所示，本发明多源数据融合的地图匹配更新方法，包括以下步骤：S1：对源地图 的道路信息矢量化处理生成矢量化数据源；对基础地图处理生成多层次空间索引分区；S2： 将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型，使得 所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地图数据；其中所述矢 量化数据源为所述路网匹配模型的观察状态，所述多层次空间索引分区为所述路网匹配模型 的隐含状态；S3：将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和修复生 成新基础地图完成更新。

本实施例实施时，源地图一般为外部来源的数据地图，例如互联网GPS导航地图，而基 础地图为需要根据源地图进行更新的地图，例如公安地理信息系统(PGIS)的道路底图。将源 地图的道路信息矢量化处理可以便于后续数据输入，无论源地图的来源为何种，都可以将这 种来源化为同样的数据格式，而基础地图处理生成的多层次空间索引分区中，一般每一条道 路为最底层的空间，多条道路形成上一层空间，多个区域形成更上一层空间，以此类推直至 形成多层次空间索引分区。在基础地图和源地图的处理中，由于二者完全独立，所以可以采 用并行计算的方式进行运算。

在本申请中通过模型进行道路匹配时是已知GPS坐标点序列，即源地图的数据，计算其 在道路网络中的真实经过路径的过程，即基础地图的数据。我们假设GPS坐标点序列为观察 状态，真实经过路径为隐含状态，基于隐马尔科夫(HMM)的路网匹配模型则是在给定一系 列观察序列的前提下，寻找最有可能产生这个观察序列的隐含状态序列，也就是实际经过的 道路。然后需要检查缺失道路信息和拓扑关系信息并完成更新。本发明通过上述步骤，实现 了对多源数据的处理，并完成了将多源数据更新到离线地图，实现了高效的、准确率高的、 对数据属性依赖较小的地图匹配更新。

为了进一步的说明本实施例的工作过程，步骤S1包括以下子步骤：将所述源地图的每条 道路提取成独立道路轨迹，并生成包含每条道路原始属性项的数据矩阵作为所述矢量化数据 源；对所述基础地图建立基于R-Tree的空间索引分区作为所述多层次空间索引分区；所述多 层次空间索引分区中，叶节点为基础地图中的道路，任意一个父节点为该父节点所有子节点 的集合直至根节点。

在本实施例中，将源地图的每条道路提取成独立道路轨迹，生成包含原始属性项的 DataFrame数据格式。对数据采用向量化计算方法，能够以极小的空间和时间成本上产生衍 生变量(如道路坐标点轨迹序列、起始点、几何长度、道路等级等)，每一行表示一条道路的 属性数据。通过预处理后，每条轨迹可以独立匹配，能够避免单条轨迹匹配错误影响到整个 网络匹配准确率，同时独立轨迹的匹配能够利用多任务并行计算以提高匹配效率。

同时，本实施例中为基础地图建立了基于R-Tree的空间索引分区，以提升道路匹配时的 空间搜索效率。如图2所示，黑色实线的几何形状是道路，实线矩形R8～R20是给每条道路 添加最小限定箱，虚线矩形R3～R7则是包含最近的几个区域的更高一层的限定区域，同理， 再次迭代生成更高层次的限定区域如虚线矩形R1和R2，直到根节点。采用R-Tree进行空间 搜索时无需穷举整个地图，只需要按照遍历少数几个叶节点即可完成搜索，其时间复杂度为

其中N为道路总数，M为最大子节点个数。

为了进一步的说明本实施例的工作过程，步骤S2包括以下子步骤：将所述矢量化数据源 和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型后获取所述路网匹配模型的 观察概率和转移概率；根据所述观察概率和所述转移概率寻找并匹配与所述矢量化数据源最 大似然的所述多层次空间索引分区中的道路。

基于ETL产生的DataFrame格式源地图数据和进行R-Tree空间索引的基础地图，我们建 立了地图匹配并行计算框架。在源地图数据中每一条数据都是独立的一条道路，因此我们把 地图匹配看成若干条独立的源地图道路在基础地图中进行道路匹配任务。基于此，利用计算 机多核CPU进行并行计算，最终将所有匹配结果汇总即完成了匹配任务。在一台具有32核 心CPU的计算机上进行实验发现多任务并行计算能够比单任务计算提升约30倍的效率。

为了进一步的说明本实施例的工作过程，所述路网匹配模型的观察概率采用下式获取：

式中，p(o_t,i|c_t,i)为所述观察概率，o_t,i为所述矢量化数据源中的离散坐标点序列，c_t,i为 离散坐标点在候选路段上的投影点，d_t,i为观察距离，u为距离均值，σ_z为距离标准差；

所述路网匹配模型的转移概率采用下式获取：

式中，p(d_t,θ_t)为所述转移概率，d_t为所述矢量化数据源中轨迹点之间距离与其在候选道 路上的投影点间距离之差，根据d_t近似符合指数分布的特征，拟合了β作为转移概率的系数， θ_t为所述矢量化数据源中相邻两个轨迹点构成的向量与其在候选路段上的投影点构成的向量 之间的夹角。

本实施例实施时，将源地图中的道路的坐标点序列视作GPS坐标序列，基础地图中的道 路视作真实道路，则可以通过HMM解码完成道路的匹配，图3是基于HMM的道路匹配模型，图3中左图灰色实线为基础地图中的道路网络，空心圆点o₁～o₄为原地图中的某条道路的 离散坐标点序列，灰色圆形区域后每个坐标点后候选区域，实心圆点为坐标点序列在基础地 图中候选路段上的投影点，基于隐马尔科夫的匹配模型就是根据观察概率和转移概率寻找最 大似然的真实道路，即左图中红色实线道路。图3中右图是采用维特比解码最优道路的示意 图，其时间复杂度为O(N·D²)，其中N表示源地图一条道路的GPS轨迹点个数，D表示每 个坐标点的候选道路个数，相对于穷举的时间复杂度O(D^N)，维特比算法更够显著提升道路 匹配的运行效率。在本实施例中，u取0，σ_z根据经验以及数据统计分析得出为7。

为了进一步的说明本实施例的工作过程，步骤S2还包括以下子步骤：

当所述矢量化数据源无法准确匹配于所述多层次空间索引分区中的道路时，在保证候选 道路连通的情况下弹性扩大缓冲半径或跳过矢量化数据源中该轨迹点；

扩大缓冲半径或跳过的轨迹点的累计长度不超过轨迹总长度的5％。

本实施例实施时，GPS定位设备在车辆行驶中容易受到干扰到定位变差较大，特别是车 辆在山区和隧道中定位误差有时甚至达到数百米之大。我们在道路匹配中设置了容差距离来 避免个别异常值对匹配准确率的影响。我们规定在保证候选道路连通的情况下，若GPS轨迹 点无法搜索到获选路段，则可以弹性扩大缓冲半径或跳过该GPS轨迹点，如图4所示，虚线 内的GPS点将被作为容差点。同时，我们规定扩大缓冲半径或跳过的轨迹点的累计长度不超 过该轨迹总长度的5％。

为了进一步的说明本实施例的工作过程，步骤S3包括以下子步骤：

将所述矢量化数据源中转移概率和/或观察概率为零的道路定义为缺失道路；

用所述多层次空间索引分区中的道路匹配所述矢量化数据源，并标记所述矢量化数据源 中无法被所述多层次空间索引分区中的道路匹配到的轨迹作为缺失道路；

用所述矢量化数据源检测所述多层次空间索引分区中的道路，并标记所述矢量化数据源 中无法匹配所述多层次空间索引分区中的道路的轨迹作为缺失道路；

本实施例实施时，定义缺失道路为源地图中无法成功匹配的道路，有两种情况会导致匹 配失败:观察概率为零，即道路坐标点无法在缓冲半径内好到候选路段，在实际匹配中，我们 允许道路中5％的坐标点无法找到候选路段，以排除异常值对匹配的影响；转移概率为零,转移 概率为零存在三种情况：(a)道路中相邻坐标点在候选路段上的投影点之间无法建立联通道 路；(b)相邻坐标之间的欧氏距离与其在候选路段上的路径距离之差大于阈值；(c)道路中 相邻坐标点之间的夹角与其候选路段的夹角大于阈值。

满足上述条件之一则为缺失道路，我们在实际检测中会采用交叉匹配来最终确定缺失道 路，即先用基础地图中的道路去匹配源地图，标记源地图中无法被匹配到的道路，以筛选源 地图中多条位置相近的平行道路在基础地图中只对应一条道路的情况；然后再用源地图中的 道路去检测基础地图中的缺失道路，以此保证缺失道路的查全率。

为了进一步的说明本实施例的工作过程，步骤S3还包括以下子步骤：

获取所述缺失道路在所述矢量化数据源中的关联道路，并根据下式获取所述缺失道路和 所述关联道路的夹角：

θ＝min(arccos(v₁,v₂),π-arccos(v₁,v₂))

式中，v₁为缺失道路在关联点最近两个顶点坐标组成的向量，v₂为关联道路中离关联点 最相邻的两个坐标点组成的向量，θ为所述缺失道路和所述关联道路的夹角；

根据缺失道路与关联道路的位置关系和θ的值将缺失道路分为已相交、延长相交和泛平行 三种类型；

其中已相交为所述缺失道路与所述关联道路已有明确交点；

延长相交为所述缺失道路与所述关联道路无明确交点，但是所述关联道路与所述缺失道 路的延长线上存在交点，且

泛平行为

本实施例实施时，对基础地图更新时，不仅要依据缺失道路本身的坐标位置，同时还要 考虑到其与其他道路的拓扑关系。对于某条缺失道路，我们需要获取其在源地图中相关联道 路的拓扑关系，并在基础地图中连接到到对应的关联道路上。

我们通过缺失道路与关联道路的坐标点来计算其夹角，如图5所示，其中v₁为缺失道路 在关联点最近两个顶点坐标组成的向量，v₂为关联道路中离关联点最相邻的两个坐标点组成 的向量。我们定义缺失道路与关联道路的夹角为

θ＝min(arccos(v₁,v₂),π-arccos(v₁,v₂))

根据θ的值以及缺失道路与关联道路是否有交点，可以将缺失道路更新分为三种拓扑关系： 已相交、延长相交、泛平行。如图6所示红色曲线表示源地图中的道路，蓝色曲线表示基础 地图中的道路，图6中a(1)表示已相交，即缺失道路与关联道路已有明确交点p_i；图6中b(1) 表示延长相交，即缺失道路与关联道路无明确交点，但是在缺失道路的延长线上存在交点， 同时我们规定

以避免延长线交点距离投影点太远；图6中c(1)表示泛平行，即缺失道 路末端和关联道路几乎平行，我们规定

为了进一步的说明本实施例的工作过程，步骤S3还包括以下子步骤：

当所述缺失道路的类型为已相交时，在交点处将多余部分截断删除完成更新；

当所述缺失道路的类型为延长相交时，将所述缺失道路沿端点方向延长，并在交点处截 断，删去多余部分完成更新；

当所述缺失道路的类型为泛平行时，获取所述缺失道路中首个拐角

的点作为接入点； 将所述缺失道路的端点在所述关联道路上的投影点连接所述接入点生成替换道路，将所述缺 失道路的端点到所述接入点的道路作为待替换道路；将所述替换道路替换所述待替换道路完 成更新。

本实施例实施时，道路更新本质上是将缺失的道路添加到基础地图中，但是在添加缺失 道路的时候需要考虑其与关联道路的接边问题。根据缺失道路拓扑分类，我们分别制定了3 种更新方法。

(1)已相交的更新方法。

如图6所示，缺失道路与关联道路已有交点，则只需要在交点处将多余部分截断删除即 可。如图6中a(1)所示，将缺失道路在交点处截断成线段l₁,l₂，删除l₁即可。截断操作使用 shapely工具的split方法即可，更新却是道路后如a(2)所示。

(2)延长相交的更新方法。

缺失道路与关联道路虽无交点，但是夹角

则需要缺失道路沿端点方向延长，并在 交点处截断，删去多余部分即可。如图7所示，我们将缺失道路沿着

方向延长a倍(在我 们模型中a＝5)，得到新的新的坐标点(x_*,y_*)加入到缺失道路l₁中，该问题则转换成已相交问 题并更新。新的坐标点的计算公式为：

(3)泛平行的更新方法

在地图中，道路一般表现形式为L＝linestring(c1,c2,c3,c4,…,ci,…cn)，其中pi＝(x_i,y_i)。缺失道路末端与关联道路泛平行时，无法在关联点附近找到交点，但是我们可以计算满足拐角

的最小第一个点i(若所有点的都不满足，则取最后一个点)，该交计算 公式为：

更新方法如图8(图中仅为更新示意图，实际道路中p₃p₅的长度远大于p₃p₄的长度，因此 实更新效果不会有如图中的较大变形)所示：(a)计算缺失道路中首个拐角

的点p₅；(b) 缺失道路端点p₃在关联道路上的投影点p₄；(c)用线段p₄p₅替换缺失道路的p₃p₅段，并完成 更新。

更新后的效果如图9所示，该图为某县的基础地图整体和局部更新效果(图中蓝色为基 础地图中的道路，红色为更新的缺失道路)，我们实现发现整体匹配更新准确率约为93％，运 行时间约为13分钟。

如图1所示，多源数据融合的地图匹配更新系统，包括：

ETL单元：用于对源地图的道路信息矢量化处理生成矢量化数据源；对基础地图处理生 成多层次空间索引分区；

模型单元：用于将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的 路网匹配模型，使得所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地 图数据；其中所述矢量化数据源为所述路网匹配模型的观察状态，所述多层次空间索引分区 为所述路网匹配模型的隐含状态；

检测修复单元：用于将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和 修复生成新基础地图完成更新。

为了进一步的说明本实施例的工作过程，所述ETL单元将所述源地图的每条道路提取成 独立道路轨迹，并生成包含每条道路原始属性项的数据矩阵作为所述矢量化数据源；

所述ETL单元对所述基础地图建立基于R-Tree的空间索引分区作为所述多层次空间索引 分区；所述多层次空间索引分区中，叶节点为基础地图中的道路，任意一个父节点为该父节 点所有子节点的集合直至根节点。

所述模型单元将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路 网匹配模型后获取所述路网匹配模型的观察概率和转移概率；所述模型单元根据所述观察概 率和所述转移概率寻找并匹配与所述矢量化数据源最大似然的所述多层次空间索引分区中的 道路。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说 明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护 范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。

Claims (10)

1.多源数据融合的地图匹配更新方法，其特征在于，包括以下步骤：

S1：对源地图的道路信息矢量化处理生成矢量化数据源；对基础地图处理生成多层次空间索引分区；

S2：将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型，使得所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地图数据；其中所述矢量化数据源为所述路网匹配模型的观察状态，所述多层次空间索引分区为所述路网匹配模型的隐含状态；

S3：将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和修复生成新基础地图完成更新。

2.根据权利要求1所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S1包括以下子步骤：

将所述源地图的每条道路提取成独立道路轨迹，并生成包含每条道路原始属性项的数据矩阵作为所述矢量化数据源；

对所述基础地图建立基于R-Tree的空间索引分区作为所述多层次空间索引分区；所述多层次空间索引分区中，叶节点为基础地图中的道路，任意一个父节点为该父节点所有子节点的集合直至根节点。

3.根据权利要求1所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S2包括以下子步骤：

将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型后获取所述路网匹配模型的观察概率和转移概率；

根据所述观察概率和所述转移概率寻找并匹配与所述矢量化数据源最大似然的所述多层次空间索引分区中的道路。

4.根据权利要求3所述的多源数据融合的地图匹配更新方法，其特征在于，所述路网匹配模型的观察概率采用下式获取：

式中，p(o_t,i|c_t,i)为所述观察概率，o_t,i为所述矢量化数据源中的离散坐标点序列，c_t,i为离散坐标点在候选路段上的投影点，d_t,i为观察距离，u为距离均值，σ_z为距离标准差；

所述路网匹配模型的转移概率采用下式获取：

式中，p(d_t,θ_t)为所述转移概率，d_t为所述矢量化数据源中轨迹点之间距离与其在候选道路上的投影点间距离之差，根据d_t近似符合指数分布的特征，拟合了β作为转移概率的系数，θ_t为所述矢量化数据源中相邻两个轨迹点构成的向量与其在候选路段上的投影点构成的向量之间的夹角。

5.根据权利要求3所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S2还包括以下子步骤：

当所述矢量化数据源无法准确匹配于所述多层次空间索引分区中的道路时，在保证候选道路连通的情况下弹性扩大缓冲半径或跳过矢量化数据源中该轨迹点；

扩大缓冲半径或跳过的轨迹点的累计长度不超过轨迹总长度的5％。

6.根据权利要求3所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S3包括以下子步骤：

将所述矢量化数据源中转移概率和/或观察概率为零的道路定义为缺失道路；

用所述多层次空间索引分区中的道路匹配所述矢量化数据源，并标记所述矢量化数据源中无法被所述多层次空间索引分区中的道路匹配到的轨迹作为缺失道路；

用所述矢量化数据源检测所述多层次空间索引分区中的道路，并标记所述矢量化数据源中无法匹配所述多层次空间索引分区中的道路的轨迹作为缺失道路。

7.根据权利要求6所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S3还包括以下子步骤：

获取所述缺失道路在所述矢量化数据源中的关联道路，并根据下式获取所述缺失道路和所述关联道路的夹角：

θ＝min(arccos(v₁,v₂),π-arc cos(v₁,v₂))

式中，v₁为缺失道路在关联点最近两个顶点坐标组成的向量，v₂为关联道路中离关联点最相邻的两个坐标点组成的向量，θ为所述缺失道路和所述关联道路的夹角；

根据缺失道路与关联道路的位置关系和θ的值将缺失道路分为已相交、延长相交和泛平行三种类型；

其中已相交为所述缺失道路与所述关联道路已有明确交点；

延长相交为所述缺失道路与所述关联道路无明确交点，但是所述关联道路与所述缺失道路的延长线上存在交点，且

泛平行为

8.根据权利要求7所述的多源数据融合的地图匹配更新方法，其特征在于，步骤S3还包括以下子步骤：

当所述缺失道路的类型为已相交时，在交点处将多余部分截断删除完成更新；

当所述缺失道路的类型为延长相交时，将所述缺失道路沿端点方向延长，并在交点处截断，删去多余部分完成更新；

当所述缺失道路的类型为泛平行时，获取所述缺失道路中首个拐角

的点作为接入点；将所述缺失道路的端点在所述关联道路上的投影点连接所述接入点生成替换道路，将所述缺失道路的端点到所述接入点的道路作为待替换道路；将所述替换道路替换所述待替换道路完成更新。

9.使用权利要求1～8中任意一项多源数据融合的地图匹配更新方法的系统，其特征在于，包括：

ETL单元：用于对源地图的道路信息矢量化处理生成矢量化数据源；对基础地图处理生成多层次空间索引分区；

模型单元：用于将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型，使得所述矢量化数据源匹配到所述多层次空间索引分区中生成预更新基础地图数据；其中所述矢量化数据源为所述路网匹配模型的观察状态，所述多层次空间索引分区为所述路网匹配模型的隐含状态；

检测修复单元：用于将预更新基础地图数据中缺失道路信息和拓扑关系信息进行检测和修复生成新基础地图完成更新。

10.根据权利要求9所述的多源数据融合的地图匹配更新系统，其特征在于，所述ETL单元将所述源地图的每条道路提取成独立道路轨迹，并生成包含每条道路原始属性项的数据矩阵作为所述矢量化数据源；

所述ETL单元对所述基础地图建立基于R-Tree的空间索引分区作为所述多层次空间索引分区；所述多层次空间索引分区中，叶节点为基础地图中的道路，任意一个父节点为该父节点所有子节点的集合直至根节点。

所述模型单元将所述矢量化数据源和所述多层次空间索引分区输入基于隐马尔科夫的路网匹配模型后获取所述路网匹配模型的观察概率和转移概率；所述模型单元根据所述观察概率和所述转移概率寻找并匹配与所述矢量化数据源最大似然的所述多层次空间索引分区中的道路。

Images