基于众包时空大数据的车道级道路测图方法及系统的制作方法

基于众包时空大数据的车道级道路测图方法及系统的制作方法
【专利摘要】本发明提供一种基于众包时空大数据的车道级道路测图方法，包括建立轨迹向量的相似度评估模型，基于融合经验知识的生长聚类方法进行轨迹优选，构建高斯约束混合模型，并使用EM算法求解模型参数；探测车道信息，得到道路路段车道数量初次探测结果；基于道路建设规则，对初次探测结果进行修正；根据修正后的车道数量，根据相邻情况对车道中心线进行修正。本发明降低了获取城市精细道路信息的成本，且探测方法简单、容易实现。
【专利说明】
基于众包时空大数据的车道级道路测图方法及系统
技术领域
[0001] 本发明设及众包时空大数据的高精度车道级道路测图，属于地理信息系统与智能 交通研究领域。
【背景技术】
[0002] 高精度道路地图是未来自动驾驶的血液和灵魂，而车道级精细道路信息是构建高 精度道路地图的关键组成部分。目前，研究者提出基于视觉方法从高分辨率遥感影像中提 取车道线，亦或是利用高精度激光点云数据获取道路细节信息，W及从测量车采集的大量 高精度GPS轨迹数据中提取道路路面信息（道路边界线、车道数量、车道中屯、线KRogerset al. (1999)是最早尝试利用时空DGPS(差分全球定位系统)轨迹数据提取道路中屯、线W及车 道边界线的研究者之一。随后在Rogers等人研究的基础上，利用时空GI^轨迹数据获取道路 信息逐渐发展成为一种端对端的模式。运种端对端模式的道路信息获取可W总结为如下几 个过程:首先对DGPS轨迹数据进行优化，然后将DGPS轨迹数据与现有地图数据匹配，样条曲 线拟合道路中屯、线，最后通过聚类方法提取车道信息W及细化交叉口的几何结构。John Krumm提出一种脱离原始地图的道路信息获取模式，该模式首先采用轨迹分类和融合方法 从大量DGPS轨迹数据中提取道路级别信息，然后利用高斯混合模型从归属于每一条路段的 大量轨迹数据中提取车道信息。然而，运些获取车道级精细道路信息的途径都存在数据采 集成本高、采集时间长、更新速度慢、数据处理复杂等缺点。
[0003] 随着传感器技术、无线通信和网络技术的飞速发展，"人人都是传感器"，人们的出 行会产生大量时空轨迹大数据，蕴含着丰富的精细道路信息和人类行为活动信息。轨迹数 据的采集逐渐由专业部口测量车或者专业人员采集演变为由非专业人±自由自愿的记录 其出行轨迹的形式，数据的采集开始转变为众包模式。众包模式下的车载轨迹数据(众包大 数据)无疑是目前可W提供车道级道路信息提取的最佳数据源。与现有的出租车数据相比， 由众包模式采集的车载轨迹数据属于大数据(大数据是指无法在可承受的时间范围内用常 规软件工具进行捕捉、管理和处理的数据集合，是需要新处理模式才能具有更强的决策力、 桐察发现力和流程优化能力来适应海量、高增长率和多样化的信息资产）。目前国内学者唐 炉亮等人(2015,2016)提出利用低精度GI^轨迹数据提取城市车道级道路信息，包括车道数 量、车道转向、车道中屯、线，然而如何利用众包大数据，开展车道级精细道路测图是全世界 科学家们面临的难题。

【发明内容】

[0004] 本发明在W上研究的基础上，提出了一种基于众包时空大数据的高精度车道级道 路测图（优质轨迹数据滤选和高精度道路信息提取)的新技术方案。
[0005] 本发明技术方案提供一种基于众包时空大数据的车道级道路测图方法，包括W下 步骤，
[0006] 步骤1，建立轨迹向量的相似度评估模型，设Va和Vb是两个不同的轨迹向量，所述相 似度评价模型如下，
[0007]
[000引其中，'>'叫康示向量之间的相似度值，e为自然底数，…和CO 2分别表示距离因 子diffHd和角度因子diffeab的权重值，且W 1+?2 = 1;距离因子diffHd和角度因子diff0￡ib分 别表示向量Va和Vb的距离差异和角度差异；
[0009] 步骤2,基于融合经验知识的生长聚类方法进行轨迹优选，包括根据已有的高精度 GPS轨迹数据与同步的低精度GI^轨迹数据，确定相似度评价模型的权重值CO 1和《2,提取轨 迹优选的先验知识，基于众包轨迹数据之间的相似度采用生长聚类方式进行数据优选；
[0010] 步骤3,构建高斯约束混合模型，并使用算法求解模型参数;所述高斯约束混合 模型定义如下，
[0011]
[0012] 其中，P(X)表示为高斯约束混合模型的综合概率值，X表示待计算样本值，在进行 车道计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值;k是高斯成分的 数量，每一个高斯成分对应一个车道；CO J是第j个高斯成分的权重，对应车道的交通流量； 参数yr-化是每一个高斯成分中轨迹的平均值，等于每个车道的中屯、线，iij表示参数 内的任意一个值，j = 1，2，…，k; 0是每一个高斯成分中轨迹的标准差；
[0013] 所述高斯约束混合模型中高斯成分的数量k获取方式为，计算结构风险模型的值， W结构风险模型值最小为原则确定k;
[0014] 步骤4,根据步骤3所得结果，探测车道信息，得到道路路段车道数量初次探测结 果;实现方式如下，
[0015] 将处于同一条路段上的所有轨迹作为一个提取单元，设给定一组从交叉口 Intersectioni到交叉口 Intersections的轨迹集合AT，从轨迹集合aT的一端开始，构建移动 矩形窗口，其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中屯、线，移动矩形窗口的 宽边则垂直于当前覆盖所有轨迹的中屯、线，矩形窗口长边的中屯、线垂直于其覆盖的所有轨 迹数据的中屯、线，矩形窗口宽边的中线与其覆盖的轨迹数据的中屯、线重合；
[0016] 将移动矩形窗口从轨迹集合AT的一端开始，按照矩形窗口的长边长度开始平移， 依次利用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中屯、 线，包括根据移动矩形窗口，将移动矩形窗口内的所有轨迹点投影到矩形窗口的长边中屯、 线上，得到投影后的轨迹数据集X=(xi，x2,…，XN)，t = l，2,3，…，N，其中，xt表示投影后第t 个轨迹点的纵坐标值，N为参加投影轨迹点的个数;将轨迹数据集X代入高斯约束混合模型， 提取矩形窗口内路段的车道数量和车道中屯、线；假设从交叉口 Intersectiom到交叉口 Intersections的轨迹集合AT，矩形窗口一共进行了 1次平移，每一次平移确定的车道数量记 为Nlanef，f= 1，2，…，1，作为道路路段车道数量初次探测结果；
[0017] 步骤5,根据步骤4获取的道路路段车道数量初次探测结果，基于道路建设规则，对 初次探测结果进行修正；
[0018] 步骤6,根据步骤5得到的修正后的车道数量，根据相邻情况对车道中屯、线进行修 正。
[0019] 而且，步骤5中对初次探测结果进行修正，实现方式如下，
[0020] 第一步，对于步骤4中某次平移确定的车道数量Nlanef，比较Nlanef4和^曰]16:、 Nlanef+2，如果Nlanef和Nlanef+l是不同的，贝lJ用Nlanef替换Nlanef+l，f=l,2,…，l-2;
[0021] 第二步，根据Nlanef的值和分布对第一步的结果分类，设存在S个类，记为Cg = < Nlg，ncg〉，Nlg是类簇Cg的车道数量，ncg是Nlaneg中属于Cg，g = 1，2，...，S的总数量；
[00剖第立步，比较Cg+i和Cg，如果Nlg+i不同于Nlg，且ncg+i<cv，令Cg的Nlg替换Cg+i的Nlg+i， g = l，2,…，S，完成车道数量结果的最终优化，其中，CV是预设的阔值。
[0023] 而且，步骤6对车道中屯、线进行修正，实现方式如下，
[0024] 设某一段路段La的车道数量被修正，如果La相邻路段化和Lc同时满足与La具有相 同的车道数量，且修正前与修正后的车道数量并未发生变化，那么就将化与Lc的车道中屯、 线连接，得到La最终的车道中屯、线；如果La相邻路段化或者Lc的车道数量修正前与修正后 也发生变化，那么就根据La修正前提取的车道中屯、线位置，计算基于La修正前车道中屯、线 位置推算La路段的道路中屯、线，按照La修正后车道宽度和车道数量重新确定La修正后的车 道中屯、线位置。
[0025] 本发明提供一种基于众包时空大数据的车道级道路测图系统，包括W下模块，
[0026] 第一模块，用于建立轨迹向量的相似度评估模型，设Va和Vb是两个不同的轨迹向 量，所述相似度评价模型如下，
[0027]
[002引其中，明"。如表示向量之间的相似度值，e为自然底数，…和CO 2分别表示距离因 子diffHd和角度因子diffeab的权重值，且W 1+?2 = 1;距离因子diffHd和角度因子diff0￡ib分 别表示向量Va和Vb的距离差异和角度差异；
[0029] 第二模块，用于基于融合经验知识的生长聚类方法进行轨迹优选，包括根据已有 的高精度GPS轨迹数据与同步的低精度GPS轨迹数据，确定相似度评价模型的权重值CO 1和 ?2，提取轨迹优选的先验知识，基于众包轨迹数据之间的相似度采用生长聚类方式进行数 据优选；
[0030] 第=模块，用于构建高斯约束混合模型，并使用算法求解模型参数;所述高斯约 束混合模型定义如下，
[0031]
[0032] 其中，P(X)表示为高斯约束混合模型的综合概率值，X表示待计算样本值，在进行 车道计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值;k是高斯成分的 数量，每一个高斯成分对应一个车道；CO J是第j个高斯成分的权重，对应车道的交通流量； 参数yr-化是每一个高斯成分中轨迹的平均值，等于每个车道的中屯、线，iij表示参数 内的任意一个值，j = 1，2，…，k; 0是每一个高斯成分中轨迹的标准差；
[0033] 所述高斯约束混合模型中高斯成分的数量k获取方式为，计算结构风险模型的值， W结构风险模型值最小为原则确定k;
[0034] 第四模块，用于根据第=模块所得结果，探测车道信息，得到道路路段车道数量初 次探测结果;实现方式如下，
[0035] 将处于同一条路段上的所有轨迹作为一个提取单元，设给定一组从交叉口 Intersectioni到交叉口 Intersections的轨迹集合AT，从轨迹集合aT的一端开始，构建移动 矩形窗口，其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中屯、线，移动矩形窗口的 宽边则垂直于当前覆盖所有轨迹的中屯、线，矩形窗口长边的中屯、线垂直于其覆盖的所有轨 迹数据的中屯、线，矩形窗口宽边的中线与其覆盖的轨迹数据的中屯、线重合；
[0036] 将移动矩形窗口从轨迹集合AT的一端开始，按照矩形窗口的长边长度开始平移， 依次利用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中屯、 线，包括根据移动矩形窗口，将移动矩形窗口内的所有轨迹点投影到矩形窗口的长边中屯、 线上，得到投影后的轨迹数据集X=(xi，x2,…，XN)，t = l，2,3，…，N，其中，xt表示投影后第t 个轨迹点的纵坐标值，N为参加投影轨迹点的个数;将轨迹数据集X代入高斯约束混合模型， 提取矩形窗口内路段的车道数量和车道中屯、线；假设从交叉口 Intersectiom到交叉口 Intersections的轨迹集合AT，矩形窗口一共进行了 1次平移，每一次平移确定的车道数量记 为Nlanef，f= 1，2，…，1，作为道路路段车道数量初次探测结果；
[0037] 第五模块，用于根据第四模块获取的道路路段车道数量初次探测结果，基于道路 建设规则，对初次探测结果进行修正；
[0038] 第六模块，用于根据第五模块得到的修正后的车道数量，根据相邻情况对车道中 屯、线进行修正。
[0039] 而且，第五模块中对初次探测结果进行修正，实现方式如下，
[0040] 第一步，对于第四模块中某次平移确定的车道数量Nlanef,比较Nlanew和Nlanef、 Nlanef+2，如果Nlanef和Nlanef+l是不同的，贝lJ用Nlanef替换Nlanef+l，f=l,2,…，l-2;
[0041] 第二步，根据Nlanef的值和分布对第一步的结果分类，设存在S个类，记为Cg = < Nlg，ncg〉，Nlg是类簇Cg的车道数量，ncg是Nlaneg中属于Cg，g = 1，2，...，S的总数量；
[00创第立步，比较Cg+i和Cg，如果Nlg+i不同于Nlg，且ncg+i<cv，令Cg的Nlg替换Cg+i的Nlg+i， g = l，2,…，S，完成车道数量结果的最终优化，其中，CV是预设的阔值。
[0043] 而且，第六模块对车道中屯、线进行修正，实现方式如下，
[0044] 设某一段路段La的车道数量被修正，如果La相邻路段化和Lc同时满足与La具有相 同的车道数量，且修正前与修正后的车道数量并未发生变化，那么就将化与Lc的车道中屯、 线连接，得到La最终的车道中屯、线；如果La相邻路段化或者Lc的车道数量修正前与修正后 也发生变化，那么就根据La修正前提取的车道中屯、线位置，计算基于La修正前车道中屯、线 位置推算La路段的道路中屯、线，按照La修正后车道宽度和车道数量重新确定La修正后的车 道中屯、线位置。
[0045] 本发明构筑了一种众包时空大数据的高精度车道级道路测图技术方案，降低了获 取城市精细道路信息的成本，且探测方法简单、容易实现。本发明所提供技术方案包括:首 先，通过对比高精度车载GPS轨迹与其同步低精度GPS轨迹数据之间的空间相似度，采用基 于经验知识的生长聚类方法，从众包轨迹数据中挑选出定位精度相对较高的数据;其次，构 建垂直于轨迹数据的移动窗口；然后，采用优化后的高斯混合模型方法对处于道路路段的 所有轨迹进行纵向探测，获取探测窗口内的车道数量;进一步，利用道路建设规则，也即，同 一条路段只有在接近交叉口位置会出现增设车道的情况，而路段中间部分车道数量通常保 持不变，提出车道数量优化策略，对车道数量信息进行修正;最后，利用修正后的车道数量 信息，对已提取的车道中屯、线进行修正，完成相应路段车道级道路信息的提取。本发明得到 的车道数量判断正确率为85%，车道中屯、线的定位精度在0.35m左右，降低了获取城市精细 道路信息的成本，且探测方法简单、容易实现。
【附图说明】：
[0046] 图1是本发明实施例的方法流程图；
[0047] 图2是本发明实施例的轨迹向量相似度示意图；
[0048] 图3是本发明实施例的基于经验知识的生长聚类示意图；
[0049] 图4是本发明实施例的基于经验知识的生长聚类方法轨迹优选结果示意图，其中 图4a为众包轨迹数据实验区域示意图，图4b轨迹优选结果示意图；
[0050] 图5是本发明实施例的高斯混合模型及车道中屯、线位置探测示意图，其中图5a为 高斯混合模型优化结果，图化为车道中屯、线位置探测结果；
[0051] 图6是本发明实施例的构建矩形窗口探测车道信息示意图；
[0052] 图7是本发明实施例的车道数量优化示意图；
[0053] 图8是本发明实施例的车道数量探测结果示意图；
[0054] 图9是本发明实施例的车道数量优化结果示意图。
【具体实施方式】
[0055] W下结合附图和实施例详细说明本发明技术方案。
[0056] 本发明提供一种众包时空大数据的高精度车道级道路测图方法，参见图1，实施例 包括如下步骤：
[0057] 步骤1，建立轨迹向量的相似度评估模型。根据车辆行驶特点，一般驾驶员按照驾 驶规则会沿着车道中屯、线行驶。因此，在不考虑短暂的变车道行为条件下，可W真实刻画驾 驶员驾驶轨迹的高精度的轨迹数据通常集中在车道中屯、线附近，且相邻轨迹之间距离小于 车道宽度，轨迹之间的航向角夹角趋近于0°左右。为了评估运种属于同一条车道中屯、线周 围的轨迹之间的相似度，本发明建立了一种新的轨迹向量相似度评估模型，该相似度模型 从轨迹向量之间的垂直距离和夹角两方面进行相似度度量，其中轨迹向量是指一条轨迹中 由每一个轨迹点及其航向角构成的单位轨迹向量，也即运些轨迹向量的模都相同，且可W 随意定义。两个轨迹向量的差异从方向和距离两方面度量。如图2所示，N表示正北方向，Va< (Xa，ya)，（Xa+l，ya+l)>和Vb<(Xb，yb)，（Xb + l，yb + l)>是两个不同的轨迹向量，矢量Va和Vb的方位 角分别是9a和0b，相似度评价模型如下：
[0化引
公式7
[0059]其中，tw'%,.,,康示向量之间的相似度值，且^啼。'当相似度值为1时表示 两个轨迹向量完全相同，相似度值为0时表示两个轨迹向量完全不相似;e为自然底数；CO 1 和O 2分别表示距离因子(diffHd)和角度因子(diffeab)的权重值，且《 1+? 2 = 1 和 diffeab分别表示向量Va和Vb的距离差异和角度差异。公式8定义两个轨迹向量的距离差异 diffHd，公式9定义了两个轨迹向量的角度差异differ:
[0060] 公式 8
[0061 ] 公式 9
[0062] 公式8中的DiSconfine由车道的宽度决定，是一个常量，用来约束相同车道上靠近每 个车道中屯、线的GPS轨迹的相似度。Hdab向量Va起点到Vb起点的垂直距离，计算公式如公式 1 1所示;Hdba是向量Vb起点到向量Va起点的垂直距离，计算公式如公式1 2所示；A 0是向量Va 与向量Vb的航向角度差值，其中矢量Va和Vb的方位角分别是03和0b，计算公式分别如下：
[0063]
[0064] 公式 11
[00化] 公式12
[0066] 步骤2,基于融合经验知识的生长聚类方法进行轨迹优选。高精度的轨迹点其轨迹 向量之间的相似度值也相对较高。在进行轨迹优选时，首先需要已有的高精度GI^轨迹数据 与其同步低精度GPS轨迹数据，提取轨迹优选的先验知识。实施例中，高精度DGPS轨迹和同 步的低精度轨迹的精度分别为0.5m和10-15m，采样频率均为1S。按照步骤一建立的相似度 估算模型计算低精度轨迹和DGPS轨迹的相似度，进行相似度评价模型权重值确定及经验知 识获取。
[0067] 1)相似度评价模型权重值确定
[0068] 相似度不仅用于DGPS和低精度GPS数据中提取先验知识，而且还用于众源数据聚 类和优选，因此本发明进一步提出使用垂直距离、角度差异与测量误差的相关性来估计 和W 2的值。
[0069] 对于GPS轨迹集合T = <Tracei ,化日。日2,... ,!"races〉，其包含S条GPS轨迹：Tracei, Traces，…，Traces，轨迹集合T的同步高精度DGPS轨迹表示为DT = <Dti，Dt2，…，Dts〉，Dti， 化2，…，Dts是集合DT内的轨迹数据。轨迹集合DT和轨迹集合T的精度分别是0.5m和10-15m。 假设第i条GPS轨迹化acei = <pi，P2，. ? ?，Pn>，Pi，P2，. ? ?，Pn分别表示轨迹化acei的轨迹点，共有 n个轨迹点；Dti = <巧1，巧2，...，巧。〉，巧1，巧2，...，巧n则为轨迹hacei的同步高精度轨迹Dti 的轨迹点；Tracei G T，Dti GDT，i = 1，2，…，S。由轨迹Tracei和其同步高精度轨迹化i内的轨 迹点，构成的轨迹向量表示为：Tvi = <Vl，V2，…，Vn-1>，VI，V2，…，Vn-I分别表示为由轨迹 Tracei的轨迹点Pi与P2,…，Pn-I与Pn构成的轨迹向量;Dvi = <;rvi,;rv2,…，rvn-i〉，:rvi,;rv2,…， rvn-i分别表示为由轨迹化acei的同步高精度轨迹化i轨迹点巧1与巧2，'''，巧n-l与巧n构成的 轨迹向量。Tvi和Dvi的距离和角度分别表示为:Di = <山,d2，…，山-1〉，di, d2，…，dn-1表示轨迹 向量集合Tvi和Dvi内相对应的轨迹向量Vi ,V2,…，Vn-I到rvi,;rv2,…，rvn-i的距离;Ai = <ai, B2，…，an-l>，ai，曰2，…，an-l表示轨迹向量集合Tvi和Dvi内相对应的轨迹向量Vl , V2，…，Vn-I与 rvi，;rv2，? ? ?，rvn-i的角度差异，i = l，2，'''，s。轨迹集合T内所有的轨迹数据的定位误差可W 表示为Ei = <ei, ￡2,…，En〉，其中Ei表示轨迹集合T内的轨迹化acei所有轨迹点与其对应的同 步高精度轨迹化1的所有轨迹点的空间距离，Ej是集合Ei内的任意一个误差值，其中Ej= Ipj- rpj I，Pj为化acei内任意一个轨迹点，rpj为Pj相对应的高精度轨迹点，i = l，2,…，s，j = l, 2，…，n。相似度评价模型中，权值…和CO 2的计算公式分别如下(掀和Me分别表示Di和Ei的 相关系数，可W基于协方差矩阵估计Al和Ei的值）：
[0070]
公式 13
[007"1] (〇2二1-"1 公式 14
[0072] 2)先验知识提取
[0073] 第一步:按照步骤1提出的向量相似度评价模型及权重值的确定，计算高精度DGPS 与其同步低精度GPS之间的相似度，同时通过对比低精度GI^轨迹数据的位置和其同步高精 度DGPS数据的位置，计算低精度GI^数据的测量误差；
[0074] 第二步，根据第一步获得的相似度W及其对应的GPS测量误差，构建归属于同一个 GPS轨迹数据的属性描述对，也即（相似度值，GPS测量误差)是某一个轨迹的属性描述对。 [007引第立步，从GPS数据的属性描述对内，按照相似度阔值，也即从相似度值为0.5开 始，选择出相似度值大于0.5,大于0.6,大于0.7…，大于0.9的所有GPS数据，并统计属于每 一个相似度阔值内所有GPS数据的GPS现慢误差，并计算运些GPS现慢误差的平均值W及满 足运些阔值的数据的个数占总体数据的比例。
[0076] 第四步，将第=步设定的相似度阔值（如：相似度阔值大于0.5,相似度阔值大于 0.6…定义为:了34，11 = 1，2,3,4,5;将满足运些阔值的6?5现慢误差的平均值定义为134，11 =1，2，3，4，5;将满足阔值的GPS测量误差占总体数据的比例定义为:Perh，h=l，2,3,4,5; 完成先验知识获取。
[0077] 其中，S化表示满足TSh的数据集，S化GT，T是用于经验提取的实验数据的总体数据 集，S化的百分比计算公式如下：
[0078]
公式 15
[0079] 化n表示STh的百分比，N(S化)和N(T)是分别是STh和T轨迹点的数量，S化的计算公 式为：
[0080] 公式化
[0081 ] TEh是Tsh的测量误差，Ee是STh中所有轨迹点误差的和，h=l，2,…，5，Tsh，Teh， 化n记为RS化=<Tsh，TEh，Perh〉，RS化为先验知识集合，作为生长聚类方法的先验知识。
[0082] 根据提取的先验知识，对相似度评价模型中的距离权值和角度权值进行计算，然 后基于众包轨迹数据之间的相似度采用生长聚类方法进行数据优选，其中先验知识内的Te 是聚类的阔值，Per被用于从整个轨迹簇中选取高精度的数据的比例，参见图3(图3中Vs表 不种子向量，Vsn表不与种子向量进行相似度计算的任意一个向量，Clusteri ,Clusters， Clusters表示通过聚类后得到的几个轨迹向量类簇；图3中（d)部分内的'Selected data' 则表示最终被选取的轨迹数据），生长聚类方法的主要步骤如下：
[0083] 第一步:初始化所有的轨迹向量，标记为未聚类;初始化当前类簇的编号(CCl),记 CCl- 1 ；
[0084] 第二步：如果存在未聚类的轨迹向量，则从剩余的未聚类轨迹向量中随机选择一 个轨迹向量作为种子轨迹向量Vs,种子轨迹向量的聚类标记即为化(Vs),且化(Vs)=CCl,如 图3中(a)部分，进入第=步;如果所有轨迹都已经被聚类，如图3中（C)部分，则进入第五步；
[0085] 第；步:捜索Vs临近的轨迹向量，记为Vsn:如果满足Sim( Vs，Vsn)〉Te，将Vs和Vsn融合 为一个轨迹簇，Vsn的聚类标记记为化(Vsn),且化(Vsn) =CCl,进入第四步；如果找不到与当 前种子向量Vs满足相似度阔值的轨迹向量时，令CCl = CCl+1，返回第二步；（Sim(vs，vsn)是Vs 和Vsn的相似度值，Te是相似度阔值）
[0086] 第四步:令轨迹向量Vsn作为种子轨迹Vs,返回第S步，如图3中(b)部分；
[0087] 第五步:根据第一步到第四步的运行结果，也即最终的聚类类簇，计算所有类簇内 轨迹点数量占参与聚类轨迹点总数的比例，然后对所有的轨迹簇按照其轨迹点数占总体轨 迹点数比值大小倒序排列，其中化r表示数据优选的选择度，从第一个类簇的轨迹点比例数 开始累积求和直到满足Per,运些参与累积求和的类簇被作为高精度数据优选选取，如图3 中（d)部分，选取的实验区域众包轨迹数据如图4a，轨迹优选的结果如图4b。
[0088] 步骤3,构建高斯约束混合模型，并使用BI算法求解模型参数。在已有的高斯约束 混合算法求取车道信息的算法基础上，本发明对其进行了优化。高斯约束混合模型定义如 下：
[00例 公式17
[0090] 其中，P(X)表示为高斯约束混合模型的综合概率值，X表示待计算样本值(在进行 车道计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值）；k是高斯成分的 数量，也即约束高斯混合模型中高斯峰的个数，代表车道的数量，每一个高斯成分对应一个 车道；《 r k是每个成分的权重，对应每个车道的交通流量，其中权重值为正且被标准化， 即W j是第j个高斯成分的权重，W j〉〇, j = l，2,…，k, ?1+?2+…Wk=I;参数化…化是每一 个高斯成分中轨迹的平均值，等于每个车道的中屯、线，表示yr-iik参数内的任意一个值，j =1，2,…，k;〇是每一个高斯成分中轨迹的标准差，并且由于每个车道与邻近车道的宽度通 常相同，因此把O设为一个常数，通常设定为1.75(根据国内道路建设标准，车道宽度一般为 3.75m左右，在具体实施过程中，本领域技术人员可W按照所选区域道路建设标准进行重新 设定）。使用EM法求解模型参数:0产> (W产>，y产> yml) J二1，2，…，k，其中m是迭代次数。目 前如何利用EM算法求解高斯混合模型的各个参数已经有很多成熟的方法，具体实施过程 中，技术人员可W参考现有方法，本发明不予寶述。
[0091] 高斯约束混合模型的关键是获得高斯成分的数量，也即计算每一个k值对应下，结 构风险模型的值，然后从其中选择出结构风险模型值最小时对应的k作为其车道数量。结构 风险模型的构建方法如下所示：
[0092]
公式 18
[0093] k=min(Rsrm(p(xi I 0k)))公式 19
[0094] L(xi,p(xi| 0k))=-log(p(xi| 0k))公式20
[00对公式帥的Rsrm(P(Xi|00)是结构风险模型，1(町口佔|目0)是用于评估适合度的经 验风险模型，J(P(Xi|00)是正则项，用于标示模型复杂度，也即表示为:jTSW(P(Xi|0k))，A〉O 是正则参数，P(X1 I目0表示样本值Xi在模型参数目k条件下的高斯概率值，其中模型参数目k可 表示为:目k( Wk,化，0)，n表示样本的个数，i = l，n;计算公式如下：
[0096] 公式 21 公式22
[0097]
[0098] 其中Dw是优化后的轨迹在道路路面的平铺宽度，如图5,图5a内长虚线标示的 ' IStGaussian component '表示高斯混合模型的第一个成分，短虚线标示的' 2ndGaussian component'表示高斯混合模型的第二个成分;其中位于图化中的参数则分别对应为图 5a内第一个高斯成分的均值及第二个高斯成分的均值，也即为第一个车道的车道中屯、线位 置和第二个车道的车道中屯、线位置。（如何获取轨迹在道路路面的平铺宽度目前已经有很 多方法提出，具体实施过程中，本领域技术人员可自行选取，本发明不予寶述），其中k表示 可能存在的车道数量，按照目前国内道路建设标准，城市车道数量一般包括两车道，=车 道，四车道，五车道，也即k = 2,3,4,5。Ayk是k对应的高斯成分中两个临近高斯峰的平均值 的变化值，其变化值也反映了探测车道宽度，j = l，2,…，k。Ayk的计算方法如公式6,公式 6内K、ri、丫 1功基于最大后验估计的EM算法的超参数;X表示样本值化进行车道计算时，X表 示判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值）；n表示样本数据的个数，CO W表 示第j+1个高斯成分的权重值，j = l，2,…，k-l，k为高斯成分数，k = 2,3,4,5;现有研究中已 有非常成熟的参数推荐，具体实施时可W参考现有方法案例中给出的数值进行计算，具体 不再寶述。
[0099] 步骤4,根据步骤3所述的优化后的高斯混合模型方法，完成高精度车道级道路测 图，实现探测车道信息，得到道路路段车道数量初次探测结果。在具体实施高精度车道级道 路测图过程中，利用现有方法将处于同一条路段上的所有轨迹作为一个提取单元(具体如 何将同一条路段上的所有轨迹归为一个提取单元，目前已经有很多成熟方法，具体实施过 程中本领域技术人员可W参见现有方法，本发明不予寶述）。
[0100] 假设给定一组从交叉口 Intersectioni到交叉口 Intersections的轨迹集合AT，从 轨迹集合AT的一端开始，构建移动矩形窗口，如图6所示。移动矩形窗口的长度和宽度分别 为rh和rw(矩形窗口的长和宽推荐定义为IOm和30m，具体实施时本领域技术人员可自行预 设取值），其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中屯、线，移动矩形窗口的宽 边则垂直于当前覆盖所有轨迹的中屯、线，矩形窗口长边的中屯、线垂直于其覆盖的所有轨迹 数据的中屯、线，矩形窗口宽边的中线与其覆盖的轨迹数据的中屯、线重合，如何得到一段轨 迹数据的中屯、线，目前已经有很多方法，具体实施过程中本领域技术人员可W参见现有方 法，本发明不予寶述。
[0101 ]实施例进一步提供的优选实现方式如下：
[0102]将移动矩形窗口从轨迹集合AT的一端开始，按照矩形窗口的长边长度开始平移， 依次利用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中屯、 线。具体方法包括:根据构建好的矩形窗口，将矩形窗口内的所有轨迹点投影到矩形窗口的 长边中屯、线上，得到投影后的轨迹数据集X=(xi,X2，，，，，XN)，t = l ，N，其中，Xt表示 投影后第t个轨迹点的纵坐标值，N为参加投影轨迹点的个数。按照步骤3所述，将轨迹数据 集X代入相应计算公式（即公式17所示高斯约束混合模型），提取矩形窗口内路段的车道数 量和车道中屯、线。假设从交叉口 Intersectioni至Ij交叉口 Intersections的轨迹集合AT，矩形 窗口一共进行了 1次平移，每一次平移确定的车道数量记为Nlanef，f = 1，2，…，1，作为道路 路段车道数量初次探测结果。
[0103]步骤5,根据步骤4获取的道路路段车道数量初次探测结果，基于道路建设规则，提 出车道数量优化策略，对初次探测结果进行修正，如图7所示，有路段化Kl、化的2、化的3，其 中化和化存在增设车道区域。大多数情况下，在两个交叉口之间的路段总是在靠近交 叉口附近处会出现增设车道，而路段中间部分车道数量通常保持不变，因此本发明提出一 种优化车道数量提取结果的方法，具体方法如下：
[0104] 第一步：对于步骤4中某次平移确定的车道数量Nlanef，比较Nlanef4和^曰116!、 Nlanef+2，如果Nlanef和Nlanef+l是不同的，贝lJ用Nlanef替换Nlanef+l，f=l,2,…，l-2。
[0105] 第二步：根据Nlanef的值和分布对第一步的结果分类，例如如果Nlanee，Nlanee+i, Nlanee+2，…，Nlanee+c的值相同则被分为一类，其中e<l，e+(3<l。假设存在S个类，记为Cg = < Nlg，ncg〉，Nlg是类簇Cg的车道数量，ncg是Nlaneg中属于Cg，g = 1，2，...，S的总数量。
[0106] 第立步：比较Cg+i和Cg,如果Nlg+i不同于Nlg,且ncg+i<cv(cv是一个依赖于道路修建 规则的阔值，其修改规则主要体现在城市道路在接近交叉口时，出现增设车道的缓冲区范 围，例如：按照当前道路设计规定，车道增设的长度为50m左右，也即从位于交叉口处的路 段，其新增车道的长度为50m，因此当车道数量判断的分割段设为10米，本发明推荐将CV设 为5，具体实施时本领域技术人员可自行预设取值），令Cg的Nlg替换Cg+i的Nlg+i，g= 1，2，…， S，完成车道数量结果的最终优化。车道数量探测结果如图8,车道数量探测优化如图9(其中 图8的横坐标表示移动窗口进行滑动探测过程中的移动数量，纵坐标则表示每一次滑动探 测过程的车道数量探测结果）。
[0107] 步骤6,根据步骤5得到的修正后的车道数量，对其相对应的车道中屯、线进行修正。 当某一段路段的车道数量被修正后，其对应的车道中屯、线则采用相邻原则也被修正。具体 方法包括：
[0108] 假设某一段路段La的车道数量被修正，那么寻找La前后相邻路段，如果La相邻路 段化和Lc同时满足与La具有相同的车道数量，且运些路段修正前与修正后的车道数量并未 发生变化，那么就将化与Lc的车道中屯、线连接，得到La最终的车道中屯、线;如果La相邻路段 化或者Lc的车道数量修正前与修正后也发生变化，那么就根据La修正前提取的车道中屯、线 位置，计算基于La修正前车道中屯、线位置推算La路段的道路中屯、线，（如何根据La修正前车 道中屯、线位置推算La路段的道路中屯、线，目前已有成熟方法，具体不在寶述)按照La修正后 车道宽度和车道数量重新定义La修正后的车道中屯、线位置，也即从根据La路段的道路中屯、 线开始依次按照车道数量和车道宽度等距离平行于道路中屯、线得到La修正后每一个车道 的车道中屯、线位置。
[0109] 基于本发明，可W方便地从GPS轨迹数据中获取待城市道路的车道信息，为未来智 能导航及无人驾驶提供基础路网数据。
[0110] 具体实施时，本发明所提供方法可基于软件技术实现自动运行流程，也可采用模 块化方式实现相应系统。
[0111] 本发明提供一种基于众包时空大数据的车道级道路测图系统，包括W下模块，
[0112] 第一模块，用于建立轨迹向量的相似度评估模型，设Va和Vb是两个不同的轨迹向 量，所述相似度评价模型如下，
[0113]
[0114] 其中，W'"如。,,康示向量之间的相似度值，e为自然底数，…和CO 2分别表示距离因 子diffHd和角度因子diffeab的权重值，且W 1+?2 = 1;距离因子diffHd和角度因子diff0￡ib分 别表示向量Va和Vb的距离差异和角度差异；
[0115] 第二模块，用于基于融合经验知识的生长聚类方法进行轨迹优选，包括根据已有 的高精度GPS轨迹数据与同步的低精度GPS轨迹数据，确定相似度评价模型的权重值CO 1和 ?2，提取轨迹优选的先验知识，基于众包轨迹数据之间的相似度采用生长聚类方式进行数 据优选；
[0116] 第=模块，用于构建高斯约束混合模型，并使用算法求解模型参数;所述高斯约 束混合模型定义如下，
[0117]
[0118] 其中，P(X)表示为高斯约束混合模型的综合概率值，X表示待计算样本值，在进行 车道计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值;k是高斯成分的 数量，每一个高斯成分对应一个车道；Oj是第j个高斯成分的权重，对应车道的交通流量； 参数是每一个高斯成分中轨迹的平均值，等于每个车道的中屯、线，iij表示参数 内的任意一个值，j = 1，2，…，k; 0是每一个高斯成分中轨迹的标准差；
[0119] 所述高斯约束混合模型中高斯成分的数量k获取方式为，计算结构风险模型的值， W结构风险模型值最小为原则确定k;
[0120] 第四模块，用于根据第=模块所得结果，探测车道信息，得到道路路段车道数量初 次探测结果;实现方式如下，
[0121] 将处于同一条路段上的所有轨迹作为一个提取单元，设给定一组从交叉口 Intersectioni到交叉口 Intersections的轨迹集合AT，从轨迹集合aT的一端开始，构建移动 矩形窗口，其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中屯、线，移动矩形窗口的 宽边则垂直于当前覆盖所有轨迹的中屯、线，矩形窗口长边的中屯、线垂直于其覆盖的所有轨 迹数据的中屯、线，矩形窗口宽边的中线与其覆盖的轨迹数据的中屯、线重合；
[0122] 将移动矩形窗口从轨迹集合AT的一端开始，按照矩形窗口的长边长度开始平移， 依次利用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中屯、 线，包括根据移动矩形窗口，将移动矩形窗口内的所有轨迹点投影到矩形窗口的长边中屯、 线上，得到投影后的轨迹数据集X=(xi，x2,…，XN)，t = l，2,3，…，N，其中，xt表示投影后第t 个轨迹点的纵坐标值，N为参加投影轨迹点的个数;将轨迹数据集X代入高斯约束混合模型， 提取矩形窗口内路段的车道数量和车道中屯、线；假设从交叉口 Intersectiom到交叉口 Intersections的轨迹集合AT，矩形窗口一共进行了 1次平移，每一次平移确定的车道数量记 为Nlanef，f= 1，2，…，1，作为道路路段车道数量初次探测结果；
[0123] 第五模块，用于根据第四模块获取的道路路段车道数量初次探测结果，基于道路 建设规则，对初次探测结果进行修正；
[0124] 第六模块，用于根据第五模块得到的修正后的车道数量，根据相邻情况对车道中 屯、线进行修正。
[0125] 各模块具体实现可参见相应步骤，本发明不予寶述。
[0126] 本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领 域的技术人员可W对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替 代，但并不会偏离本发明的精神或超越所附权利要求书所定义的范围。
【主权项】
1. 一种基于众包时空大数据的车道级道路测图方法，其特征在于:包括以下步骤， 步骤1，建立轨迹向量的相似度评估模型，设vdPvb是两个不同的轨迹向量，所述相似度 评价模型如下，其中，表示向量之间的相似度值，e为自然底数，ω :和ω 2分别表示距离因子 diffHd和角度因子diffeab的权重值，且ω # ω 2= 1;距离因子diffHd和角度因子diffeab分别 表不向量Va和Vb的距尚差异和角度差异； 步骤2,基于融合经验知识的生长聚类方法进行轨迹优选，包括根据已有的高精度GPS 轨迹数据与同步的低精度GPS轨迹数据，确定相似度评价模型的权重值ω :和ω 2，提取轨迹 优选的先验知识，基于众包轨迹数据之间的相似度采用生长聚类方式进行数据优选； 步骤3，构建高斯约束混合模型，并使用ΕΜ算法求解模型参数;所述高斯约束混合模型 定义如下，其中，P(x)表示为高斯约束混合模型的综合概率值，X表示待计算样本值，在进行车道 计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值;k是高斯成分的数量， 每一个高斯成分对应一个车道；ω」是第j个高斯成分的权重，对应车道的交通流量；参数 ydk是每一个高斯成分中轨迹的平均值，等于每个车道的中心线，表示ydk参数内的 任意一个值，j = l，2,…，k;〇是每一个高斯成分中轨迹的标准差； 所述高斯约束混合模型中高斯成分的数量k获取方式为，计算结构风险模型的值，以结 构风险模型值最小为原则确定k; 步骤4,根据步骤3所得结果，探测车道信息，得到道路路段车道数量初次探测结果;实 现方式如下， 将处于同一条路段上的所有轨迹作为一个提取单元，设给定一组从交叉口 Intersectiom到交叉口 Intersectiom的轨迹集合Ατ，从轨迹集合Ατ的一端开始，构建移动 矩形窗口，其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中心线，移动矩形窗口的 宽边则垂直于当前覆盖所有轨迹的中心线，矩形窗口长边的中心线垂直于其覆盖的所有轨 迹数据的中心线，矩形窗口宽边的中线与其覆盖的轨迹数据的中心线重合； 将移动矩形窗口从轨迹集合Ατ的一端开始，按照矩形窗口的长边长度开始平移，依次利 用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中心线，包括 根据移动矩形窗口，将移动矩形窗口内的所有轨迹点投影到矩形窗口的长边中心线上，得 到投影后的轨迹数据集乂=(^，1 2，一，別）4 = 1，2,3，一，1其中，^表示投影后第七个轨迹 点的纵坐标值，N为参加投影轨迹点的个数;将轨迹数据集X代入高斯约束混合模型，提取矩 形窗口内路段的车道数量和车道中心线；假设从交叉口 Intersect ioru到交叉口 Intersection的轨迹集合Ατ，矩形窗口一共进行了 1次平移，每一次平移确定的车道数量记 为Nlanef，f = 1，2，…，1，作为道路路段车道数量初次探测结果； 步骤5,根据步骤4获取的道路路段车道数量初次探测结果，基于道路建设规则，对初次 探测结果进行修正； 步骤6,根据步骤5得到的修正后的车道数量，根据相邻情况对车道中心线进行修正。2. 根据权利要求1所述基于众包时空大数据的车道级道路测图方法，其特征在于:步骤 5中对初次探测结果进行修正，实现方式如下， 第一步，对于步骤4中某次平移确定的车道数量Nlanef，比较Nlanef + jPNlanef、 Nlanef+2,如果Nlanef和Nlanef+i是不同的，贝lj用Nlanef替换Nlanef+i，f = 1，2，…，1-2; 第二步，根据Nlanef的值和分布对第一步的结果分类，设存在s个类，记为Cg =〈Nig，ncg >，Nig是类簇Cg的车道数量，ncg是Nlaneg中属于Cg，g = 1，2，…，s的总数量； 第三步，比较 Cg+jPCg，如果 Nlg+1 不同于 Nlg，且 ncg+1〈cV，令 Q^Nlg 替换 Cg+1 的 Nlg+1，g=l， 2，···，s，完成车道数量结果的最终优化，其中，cv是预设的阈值。3. 根据权利要求1或2所述基于众包时空大数据的车道级道路测图方法，其特征在于： 步骤6对车道中心线进行修正，实现方式如下， 设某一段路段La的车道数量被修正，如果La相邻路段Lb和Lc同时满足与La具有相同的 车道数量，且修正前与修正后的车道数量并未发生变化，那么就将Lb与Lc的车道中心线连 接，得到La最终的车道中心线；如果La相邻路段Lb或者Lc的车道数量修正前与修正后也发 生变化，那么就根据La修正前提取的车道中心线位置，计算基于La修正前车道中心线位置 推算La路段的道路中心线，按照La修正后车道宽度和车道数量重新确定La修正后的车道中 心线位置。4. 一种基于众包时空大数据的车道级道路测图系统，其特征在于:包括以下模块， 第一模块，用于建立轨迹向量的相似度评估模型，设vjPvb是两个不同的轨迹向量，所 述相似度评价模型如下，其中表示向量之间的相似度值，e为自然底数，ω^ω2分别表示距离因子 diffHd和角度因子diffeab的权重值，且ω # ω 2= 1;距离因子diffHd和角度因子diffeab分别 表不向量Va和Vb的距尚差异和角度差异； 第二模块，用于基于融合经验知识的生长聚类方法进行轨迹优选，包括根据已有的高 精度GPS轨迹数据与同步的低精度GPS轨迹数据，确定相似度评价模型的权重值ω :和ω 2，提 取轨迹优选的先验知识，基于众包轨迹数据之间的相似度采用生长聚类方式进行数据优 选； 第三模块，用于构建高斯约束混合模型，并使用ΕΜ算法求解模型参数;所述高斯约束混 合模型定义如下，其中，P(x)表示为高斯约束混合模型的综合概率值，X表示待计算样本值，在进行车道 计算时，X代表判断窗口内轨迹点在其纵剖面上垂直投影的纵坐标值;k是高斯成分的数量， 每一个高斯成分对应一个车道；ω」是第j个高斯成分的权重，对应车道的交通流量；参数 ydk是每一个高斯成分中轨迹的平均值，等于每个车道的中心线，表示ydk参数内的 任意一个值，j = l，2,…，k;〇是每一个高斯成分中轨迹的标准差； 所述高斯约束混合模型中高斯成分的数量k获取方式为，计算结构风险模型的值，以结 构风险模型值最小为原则确定k; 第四模块，用于根据第三模块所得结果，探测车道信息，得到道路路段车道数量初次探 测结果;实现方式如下， 将处于同一条路段上的所有轨迹作为一个提取单元，设给定一组从交叉口 Intersectiom到交叉口 Intersectiom的轨迹集合Ατ，从轨迹集合Ατ的一端开始，构建移动 矩形窗口，其中移动矩形窗口的长边平行于当前覆盖所有轨迹的中心线，移动矩形窗口的 宽边则垂直于当前覆盖所有轨迹的中心线，矩形窗口长边的中心线垂直于其覆盖的所有轨 迹数据的中心线，矩形窗口宽边的中线与其覆盖的轨迹数据的中心线重合； 将移动矩形窗口从轨迹集合Ατ的一端开始，按照矩形窗口的长边长度开始平移，依次利 用高斯约束混合模型来探测每一个矩形窗口内覆盖的路段的车道数量及车道中心线，包括 根据移动矩形窗口，将移动矩形窗口内的所有轨迹点投影到矩形窗口的长边中心线上，得 到投影后的轨迹数据集乂=(^，1 2，一，別）4 = 1，2,3，一，1其中，^表示投影后第七个轨迹 点的纵坐标值，N为参加投影轨迹点的个数;将轨迹数据集X代入高斯约束混合模型，提取矩 形窗口内路段的车道数量和车道中心线；假设从交叉口 Intersect ioru到交叉口 Intersection的轨迹集合Ατ，矩形窗口一共进行了 1次平移，每一次平移确定的车道数量记 为Nlanef，f = 1，2，…，1，作为道路路段车道数量初次探测结果； 第五模块，用于根据第四模块获取的道路路段车道数量初次探测结果，基于道路建设 规则，对初次探测结果进行修正； 第六模块，用于根据第五模块得到的修正后的车道数量，根据相邻情况对车道中心线 进行修正。5. 根据权利要求4所述基于众包时空大数据的车道级道路测图系统，其特征在于:第五 模块中对初次探测结果进行修正，实现方式如下， 第一步，对于第四模块中某次平移确定的车道数量Nlanef，比较Nlanef+jPNlanef、 犯已116￡+2，如果犯3116￡和犯3116￡+1是不同的，贝1|用犯3116￡替换犯3116￡+1，€=1，2，"，1-2; 第二步，根据Nlanef的值和分布对第一步的结果分类，设存在s个类，记为Cg =〈Nig，ncg >，Nig是类簇Cg的车道数量，ncg是Nlaneg中属于Cg，g = 1，2，…，s的总数量； 第三步，比较 Cg+jPCg，如果 Nlg+1 不同于 Nlg，且 ncg+1〈cV，令 Q^Nlg 替换 Cg+1 的 Nlg+1，g=l， 2，···，s，完成车道数量结果的最终优化，其中，cv是预设的阈值。6. 根据权利要求4或5所述基于众包时空大数据的车道级道路测图系统，其特征在于： 第六模块对车道中心线进行修正，实现方式如下， 设某一段路段La的车道数量被修正，如果La相邻路段Lb和Lc同时满足与La具有相同的 车道数量，且修正前与修正后的车道数量并未发生变化，那么就将Lb与Lc的车道中心线连 接，得到La最终的车道中心线；如果La相邻路段Lb或者Lc的车道数量修正前与修正后也发 生变化，那么就根据La修正前提取的车道中心线位置，计算基于La修正前车道中心线位置 推算La路段的道路中心线，按照La修正后车道宽度和车道数量重新确定La修正后的车道中 心线位置。
【文档编号】G08G1/01GK105957342SQ201610370700
【公开日】2016年9月21日
【申请日】2016年5月30日
【发明人】唐炉亮, 杨雪, 李清泉
【申请人】武汉大学