一种道路面及其标志线的自动提取方法

1.一种道路面的自动提取方法，其特征在于，包括以下步骤：
步骤1，利用统计滤波法从获得的采集数据中剔除离群点，并基于获得的轨迹数据和点云的高程信息剔除非地面点，得到候选点集合；
步骤2，沿着车辆行驶方向，将轨迹划分成具有一定宽度的轨迹片段，然后根据候选点集合在轨迹上的垂足点，分别将候选点集合分割成不同的伪扫描线，并将伪扫描线内候选点的坐标系统转换为基于伪扫描线的局部坐标系统；
步骤3，在分割的伪扫描线中，采用移动最小二乘直线拟合获取道路点的拟合直线，并通过计算点到直线的距离，来检测路缘石与道路面之间的高度突变特征，实现道路面提取；
在步骤2中，首先沿着车辆行驶方向构建轨迹直线，然后利用预设的宽度阈值 将该直线划分为等宽的连续片段集合；再分别将候选点集合投影在轨迹直线上形成垂足点集合，最后将不同片段内的垂足点集合所对应的候选点集合分割成不同的伪扫描线。
2.根据权利要求1所述的道路面的自动提取方法，其特征在于，在步骤一中，构建k‑d树索引每个点云的邻域，计算各点到它所有临近点的平均距离 ；根据高斯分布原理的均值和标准差 ，若 ，则该点被视为离群点。
3.根据权利要求1所述的道路面的自动提取方法，其特征在于，在步骤一中，根据轨迹高程 、轨迹高度 和点云高程 ，判断每个点的 的绝对值是否小
于预设的高度阈值，若绝对值小于预设高度阈值则该点被视为候选点保留，反之则视为非地面点剔除。
4.根据权利要求1所述的道路面的自动提取方法，其特征在于，在步骤2中，将每条伪扫描线内候选点的坐标系统转换为，以该伪扫描线对应的轨迹片段中心为原点的局部右手坐标系，并根据极角 将点云按逆序排序。
5.根据权利要求1所述的道路面的自动提取方法，其特征在于，在步骤3中，在每条伪扫描线内，定义轨迹正下方的一个道路点为种子点，沿着种子点到伪扫描线两端，分别移动窗口逐点检测路缘石。
6.根据权利要求5所述的道路面的自动提取方法，其特征在于，首先基于最小二乘直线拟合依次拟合移动窗口内的道路点，然后计算候选点到拟合直线的距离 ，如果距离 大于预设距离阈值 ，则该点被视为非道路点；当连续的两个候选点被判定为非道路点时，移动窗口停止检测。
7.一种道路面标志线的自动提取方法，其特征在于，包括以下步骤：
步骤1，利用统计滤波法从获得的采集数据中剔除离群点，并基于获得的轨迹数据和点云的高程信息剔除非地面点，得到候选点集合；
步骤2，沿着车辆行驶方向，将轨迹划分成具有一定宽度的轨迹片段，然后根据候选点集合在轨迹上的垂足点，分别将候选点集合分割成不同的伪扫描线，并将伪扫描线内候选点的坐标系统转换为基于伪扫描线的局部坐标系统；
步骤3，在分割的伪扫描线中，采用移动最小二乘直线拟合获取道路点的拟合直线，并通过计算点到直线的距离，来检测路缘石与道路面之间的高度突变特征，实现道路面提取；
步骤4，在提取的道路面中，采用基于局部密度的自适应窗口中值滤波器平滑道路面的强度信息；
步骤5，在平滑后的道路面中，采用多约束方法检测道路和标志线之间的边缘点，实现标志线提取；
步骤6，在提取的标志线中，基于标志线的结构特征和维度特征，优化道路标志线；
在步骤2中，首先沿着车辆行驶方向构建轨迹直线，然后利用预设的宽度阈值 将该直线划分为等宽的连续片段集合；再分别将候选点集合投影在轨迹直线上形成垂足点集合，最后将不同片段内的垂足点集合所对应的候选点集合分割成不同的伪扫描线。
8.根据权利要求7所述的道路面标志线的自动提取方法，其特征在于，在步骤4中，在每条伪扫描线内，构建k‑d树索引每个道路点的邻域并统计最近邻点的数量，基于不同的局部密度自动选择不同尺寸的滤波窗口；在滤波窗口中，将窗口的强度中值替代窗口中心的强度值。
9.根据权利要求7所述的道路面标志线的自动提取方法，其特征在于，在步骤5中，首先在每条平滑后的伪扫描线内，计算每个点云的强度梯度；然后根据平滑后的道路面的强度统计直方图，计算标志线强度特征；通过结合强度梯度和标志线强度特征检测道路与标志线之间的边缘点，并直接提取标志线。

一种道路面及其标志线的自动提取方法\n技术领域\n[0001] 本发明涉及数据处理技术领域，尤其涉及一种通用的道路面及其标志线的自动提取方法。\n背景技术\n[0002] 道路面与标志线是城市空间数据库的重要元素，在高精度定位、环境感知、路径规划以及仿真实验中发挥着重要作用，是智慧城市和智慧交通发展所必需的空间数据。\n[0003] 随着城市数字化进程加快，传统测绘手段已无法满足现代空间数据的快速获取与更新。移动激光雷达技术作为一种主动遥感技术，通过发射激光脉冲和接收回波信号，能在行驶过程中快速地获取高精度的道路场景数据。目前，商业领域可用于获取高精度数据的移动激光雷达系统分为两大类：单线移动激光系统(简称单线系统)和多线移动雷达系统(简称多线系统)。\n[0004] 单线系统采集到的点云具有精度高、内部一致性好和结构简单等特点，如图2中所示。多线系统能采集比单线系统更稠密的点云，但多线系统点云精度较低、内部一致性极其复杂且存在更多的噪声点，如图3中所示。由于两种类型的点云具有不同的内部结构，国内外的相关解译算法主要针对单一类型的数据，尤其是单线系统点云。研究一种通用的移动激光雷达数据道路面与标志线自动提取方法，这对推动移动激光雷达技术的大规模普及，和加快智慧城市、智慧交通的发展具有重要意义。\n[0005] 现有移动激光雷达数据的道路面提取方法主要基于道路面的几何特征和材料特征，包括平面特征、局部高程和路沿检测等。Yuan(2008)等基于最大熵理论的模糊聚类方法将点云聚类成扫描线，根据坡度信息提取近似的水平直线作为路面点云。刘如飞(2015)运用近似平面约束法、有序最小二乘坡度估计法和多尺度窗口迭代分析法提取道路面种子点，然后基于局部坡度滤波提取路面点。方莉娜(2013)利用时间戳将点云分割成不同的扫描线，并根据局部高程、密度和坡度等特征检测路沿，最后通过聚类、去噪和优化提取道路面。\n[0006] 现有移动激光雷达数据的道路面标志线提取方法主要分为两类：(1)基于强度特征图像的标志线提取；(2)直接从3D点云中提取标志线。Yang(2012)利用反距加权插值算法将点云投影成地理参考强度特征图像，然后利用高差信息分割地理参考特征图像，最后通过语义信息提取标志线的轮廓。Soilán(2017)将点云投影成二值图像，然后根据像素分布和几何特征计算道路标线的特征，最后通过前馈网络实现标志线分类。Jung(2019)采用拉普拉斯方程生成栅格图像，然后用期望最大化算法对标志线进行提取。虽然利用成熟的图像处理技术能快速地提取标志线，但3D点云转换为2D图像的过程可能会导致精度损失。为了避免精度损失，部分研究员致力于从3D点云中直接提取标志线。Chen(2009)等沿着轨迹将点云分割成扫描线，然后利用一个基于轮廓线的强度分析方法从扫描线中直接提取道路标线点。Yan(2016)等根据扫描角度将点云分离成扫描线，然后利用强度梯度检测标志线的边缘点。\n[0007] 虽然现在有上述诸多道路面提取方法和道路面标志线提取方法，但是目前仍然存在以下难以克服的问题：\n[0008] (1)道路面与标志线提取方法大多数还是基于单线系统数据。由于两种系统的点云具有不同的结构特征，因此现有的提取算法并不能通用于不同系统的点云。且由于现有移动激光雷达采集技术的限制，多线系统不同的激光束在不同的时间点采集数据，导致无法利用时间戳、聚类和扫描角度等方式将离散的多线系统点云组织形成扫描线结构；\n[0009] (2)标志线的提取主要还是基于强度特征图像，但在3D点云转换为2D图像的过程中，会导致精度损失。从激光雷达数据中提取标志线主要仍然是基于标志线的反射强度，但点云的强度信息容易受到扫描距离、扫描角度和材料性质等因素的影响，且多线系统点云的强度信息更为复杂。\n发明内容\n[0010] 本发明提供了一种道路面自动提取方法，以解决现有道路面自动提取方法过于依赖单线系统数据而无法通用的问题。\n[0011] 本方案中的道路面自动提取方法，包括以下步骤：\n[0012] 步骤1，利用统计滤波法从获得的采集数据中剔除离群点，并基于获得的轨迹数据和点云的高程信息剔除非地面点，得到候选点集合；\n[0013] 步骤2，沿着车辆行驶方向，将轨迹划分成具有一定宽度的轨迹片段，然后根据候选点集合在轨迹上的垂足点，分别将候选点集合分割成不同的伪扫描线，并将伪扫描线内候选点的坐标系统转换为基于伪扫描线的局部坐标系统；\n[0014] 步骤3，在分割的伪扫描线中，采用移动最小二乘直线拟合获取道路点的拟合直线，并通过计算点到直线的距离，来检测路缘石与道路面之间的高度突变特征，实现道路面提取。\n[0015] 本方案的优点在于：\n[0016] 先通过统计滤波法和高程信息去除掉离群点和非地面点，得到候选点集合，减少需要处理的点云的体积，同时为提取的精确性做了前期准备。然后，根据垂足点将候选点集合分割成不同的伪扫描线，并将通过现有技术获得的三维坐标点都转换为以伪扫描线为基础建立的坐标体系当中，能够有效屏蔽掉现有采集技术中的不足，实现了将不同结构的点云数据组织成扫描线结构的目的。使无论是多线系统还是单线系统采集到的数据，都能够被本方法进行精确处理，这也是本方法能够具有通用性的关键点之一。\n[0017] 此外，通过种子点的定义和路缘石的检测，利用了道路面的几何结构，更加准确直观地完成了道路面的自动提取。\n[0018] 本方案有效解决了现有方法不具通用性的缺点，能够很好的适应各种现有技术采集提取到三维数据完成道路面的自动提取。\n[0019] 进一步，在步骤一中，构建k‑d树索引每个点云的邻域，计算各点到它所有临近点的平均距离 根据高斯分布原理的均值μ和标准差σ，若 则该点被视为离群点。\n[0020] 通过这个判断方法，快速进行离群点筛选和剔除。\n[0021] 进一步，在步骤一中，根据轨迹高程ZPOS、轨迹高度HPOS和点云高程Zpt，判断每个点的ZPOS‑Zpt‑HPOS的绝对值是否小于预设的高度阈值，若绝对值小于预设高度阈值则该点被视为候选点保留，反之则视为非地面点剔除。\n[0022] 通过这个判断方法，快速确定候选点集合。\n[0023] 进一步，在步骤2中，首先沿着车辆行驶方向构建轨迹直线，然后利用预设的宽度阈值wth将该直线划分为等宽的连续片段集合；再分别将候选点集合投影在轨迹直线上形成垂足点集合，最后将不同片段内的垂足点集合所对应的候选点集合分割成不同的伪扫描线。\n[0024] 通过此方法，在每个候选点集合片段中形成伪扫描线，便于后面的三维点的坐标变换。\n[0025] 进一步，在步骤2中，将每条伪扫描线内候选点的坐标系统转换为，以该伪扫描线对应的轨迹片段中心为原点的局部右手坐标系，并根据极角θi将点云按逆序排序。\n[0026] 后续的道路面与标志线提取，均基于此局部坐标系统。\n[0027] 进一步，在步骤3中，在每条伪扫描线内，定义轨迹正下方的一个道路点为种子点，沿着种子点到伪扫描线两端，分别移动窗口逐点检测路缘石。\n[0028] 先定义相当于位于道路中间的道路上的三维点为种子点，然后再分别向两侧扫描检测路缘石，能够缩短找到路缘石的时间。\n[0029] 进一步，首先基于最小二乘直线拟合依次拟合移动窗口内的道路点，然后计算候选点到拟合直线的距离d，如果距离d大于预设距离阈值dth，则该点被视为非道路点；当连续的两个候选点被判定为非道路点时，移动窗口停止检测。\n[0030] 在没有路缘石的情况下通过找到非道路点来区分道路边缘界限，方便对道路面的自动提取。\n[0031] 基于上述方法，本发明还提供了一种道路面标志线的自动提取方法，以解决现有标志线提取方法过于依赖强度特征图像精准度不高的问题。\n[0032] 本方案中的道路面标志线自动提取方法，包括以下步骤：\n[0033] 步骤1，利用统计滤波法从获得的采集数据中剔除离群点，并基于获得的轨迹数据和点云的高程信息剔除非地面点，得到候选点集合；\n[0034] 步骤2，沿着车辆行驶方向，将轨迹划分成具有一定宽度的轨迹片段，然后根据候选点集合在轨迹上的垂足点，分别将候选点集合分割成不同的伪扫描线，并将伪扫描线内候选点的坐标系统转换为基于伪扫描线的局部坐标系统；\n[0035] 步骤3，在分割的伪扫描线中，采用移动最小二乘直线拟合获取道路点的拟合直线，并通过计算点到直线的距离，来检测路缘石与道路面之间的高度突变特征，实现道路面提取；\n[0036] 步骤4，在提取的道路面中，采用基于局部密度的自适应窗口中值滤波器平滑道路面的强度信息；\n[0037] 步骤5，在平滑后的道路面中，采用多约束方法检测道路和标志线之间的边缘点，实现标志线提取；\n[0038] 步骤6，在提取的标志线中，基于标志线的结构特征和维度特征，优化道路标志线。\n[0039] 本方案的优点在于：现有技术中，Yan(2016)等使用强度梯度检测单线系统点云中的标志线边缘，由于多线系统点云的复杂性，标志线内部的点云强度值，波动非常剧烈且不具有规律性，如图5所示。使用一对非零正负梯度值的方法可能会将噪点错误地检测为标志线边缘点。Wang(2018)等在提取的道路面上，根据强度的分布特征使用直通滤波器提取标志线，但其不能提取低强度的标志线点云，尤其是标志线边缘点云。\n[0040] 而本方案为了处理不同系统的点云，且能有效提取低强度的标志线边缘点。在提取的道路面上，先是通过基于局部密度的自适应窗口中值滤波器平滑道路面的强度信息，有效抑制了扫描距离、扫描角度和材料性质对点云强度值的影响，增强了算法的稳健性；然后通过多约束方法来检测道路面和标志线之间的边缘点，进而达到提取直接标志线的目的。在提取标志线的同时，基于标志线的结构特征和维度特征，优化道路标志线，剔除提取成果中的非标志噪点，使标志线提取得更加精准。本方案克服了现有标志线提取过于依赖强度特征图像的问题，有效避免了将3D点云转换为2D图像的过程中造成精度损失。且能有效从不同系统的点云中提取标志线，实现了对不同系统的通用化处理。\n[0041] 进一步，在步骤4中，在每条伪扫描线内，构建k‑d树索引每个道路点的邻域并统计最近邻点的数量，基于不同的局部密度自动选择不同尺寸的滤波窗口；在滤波窗口中，将窗口的强度中值替代窗口中心的强度值。从而达到平滑强度噪声的目的。\n[0042] 沿着种子点到道路面边缘，由于道路点的局部密度逐渐变得稀疏，使用一个固定尺寸的滤波窗口，可能会导致稀疏区域的标志边缘点的强度值被过滤，还可能导致稠密区域的标志线平滑效果不足。为了达到有效地抑制道路面的强度噪声，同时保护标志边缘点的强度值，而本方案采用了一种基于局部密度的自适应窗口中值滤波器，根据不同的局部密度，自动选择不同尺寸的窗口。这样，达到最佳的道路面平滑效果，使得后续处理中计算强度梯度和建立强度统计直方图变得更加稳健可靠。\n[0043] 进一步，在步骤5中，首先在每条平滑后的伪扫描线内，计算每个点云的强度梯度；\n然后根据平滑后的道路面的强度统计直方图，计算标志线强度特征；通过结合强度梯度和标志线强度特征检测道路与标志线之间的边缘点，并直接提取标志线。\n[0044] 由于多线系统点云的复杂性，即使道路面经过平滑处理，仍有小部分噪点没有被彻底地抑制，如果只使用强度梯度检测标志线边缘点，这些噪点很可能会被错误地检测为标志线边缘点。为了准确地检测伪扫描线中标志线的边缘点，本方案基于强度梯度和强度值统计直方图的多约束边缘检测方法，能够克服掉过度依赖强度特征图像的缺点，实现了对不同类型点云的道路面和标志线边缘精确检测，并直接提取标志线。\n[0045] 在步骤6中，首先基于欧式聚类法将点云分割成单独的类，根据标志线的结构特征剔除被错误提取的非标志点；然后构建k‑d树计算每个点云的维度特征，进一步剔除非标志点。\n[0046] 通过这些方法进一步精确标志线，提高标志线自动提取的精准度。\n[0047] 综上，本发明具有如下优点：\n[0048] (1)基于垂足点使点云分割成了伪扫描线，实现了对不同类型的数据进行有效地组织；(2)伪扫描线是一种简单的线性数据结构，处理的点云数据量更少，且能运用并行运算，处理效率得到了加大提高；(3)利用自适应窗口中值滤波器平滑道路面强度信息，克服了强度噪声影响，有效地提升了该发明的抗噪能力；(4)基于强度梯度和强度值统计直方图的检测策略，能有效地从不同类型的数据中直接提取标志线。\n[0049] 本发明中，获得的采集数据可以是移动激光雷达数据、也可以是其他采集技术采集出来的三维数据，路面数据的获取仪器不局限于汽车、推车或摩托车等移动载体，也可以是无人机、直升机等现有的其他采集设备。本发明中的道路面自动提取方法具有通用性，适用于现有的各种采集设备和采集场景。同样的，提取标志线的方法也不局限于移动采集设备，也可以是其他现有的采集设备，该方法同样具有通用性。\n附图说明\n[0050] 图1是本发明实施例一的流程图；\n[0051] 图2是本发明实施例一的单线系统数据的示意图；\n[0052] 图3是本发明实施例一的多线系统数据的示意图；\n[0053] 图4是本发明实施例一中的基于伪扫描线的局部坐标系统示意图；\n[0054] 图5是本发明实施例一平滑前局部的标志线点云强度值示意图；\n[0055] 图6是本发明实施例一平滑后局部的标志线点云强度值示意图；\n[0056] 图7是本发明实施例一中的非标志点示意图。\n具体实施方式\n[0057] 下面通过具体实施方式进一步详细说明：\n[0058] 实施例一\n[0059] 本实施例提供的技术方案能够使用已有技术手段利用计算机语言实现道路面和标志线的自动提取，道路面及其标志线的自动提取方法流程如图1所示，包括以下步骤：\n[0060] 步骤1：从现有的采集设备中获得道路采集数据，可以是从移动激光雷达采集系统中获取移动激光雷达数据(可以是单系统数据，也可以是多系统数据)，也可以是其他现有采集设备/系统中获得的采集数据，利用统计滤波法剔除离群点，并根据GNSS(所有国家的卫星)轨迹数据和高程信息从激光雷达数据中选取道路面候选点；\n[0061] 本步骤具体实施如下：\n[0062] 步骤1.1，构建k‑d树索引每个点云的邻域，并计算各点到它所有临近点的平均距离 根据高斯分布原理的均值μ和标准差σ，建立判别方程式：μ+σ，若 则该点被视为离群点，给予剔除。\n[0063] 步骤1.2，根据轨迹高程ZPOS、轨迹高度HPOS和点云高程Zpt，建立判别方程式：ZPOS‑Zpt‑HPOS，若判别式的绝对值小于高度阈值，则该点被视为候选点，给予保留。高度阈值根据经验取值，一般设置为0.6m。\n[0064] 步骤2：利用GNSS轨迹数据和点云在轨迹上的垂足点，将候选点分割成伪扫描线；\n[0065] 本步骤具体实施如下：\n[0066] 步骤2.1，首先沿着车辆行驶方向构建轨迹直线，然后利用宽度阈值wth将该直线划分为宽度约0.05m的片段集合并进行编号。再分别将候选点集合投影在直线上，根据垂足点坐落的片段编号，将垂足点对应的候选点分割为不同的伪扫描线。\n[0067] 步骤2.1，为了简化后续处理，如图4所示，将每条伪扫描线内候选点的GNSS坐标系统转换为，以片段中心为原点的局部右手坐标系统，并根据极角θi将点云按逆序排序。后续的道路面与标志线提取，均基于此局部坐标系统。\n[0068] 步骤3：利用道路面与路缘石(或排水沟)之间的高度突变，实现道路面提取；\n[0069] 为了保护道路面，路缘石(或排水沟)被修建在道路边界，道路面低于路缘石且高于排水沟。因此，在伪扫描线的横截面中，道路面与连接物之间存在一定的高度差异。基于道路面的几何结构，通过检测路缘石即可提取完整的道路面。\n[0070] 本步骤具体实施如下：\n[0071] 步骤3.1，在每条伪扫描线内，定义轨迹正下方的一个道路点为种子点。\n[0072] 步骤3.2，沿着种子点到伪扫描线两端，移动窗口逐点检测路缘石。首先基于最小二乘直线拟合依次拟合移动窗口内的道路点，见公式(1)，然后计算候选点到拟合直线的距离d，见公式(2)。如果距离d大于距离阈值dth，则该点被视为非道路点。当连续的两个候选点被判定为非道路点时，移动窗口停止检测。距离阈值根据经验取值，一般设置为0.04m。\n[0073]\n[0074] 公式(1)中，k为拟合直线的斜率；b为拟合直线的截距xi和yi是移动窗口内第i个点的2D坐标，i＝1,2,...,n；n为移动窗口的尺寸，根据经验取值，一般设置为20；\n[0075]\n[0076] 公式(2)中，A、B和C为拟合直线的回归系数，xp和yp为候选点的2D坐标。\n[0077] 步骤4：采用基于密度的自适应窗口中值滤波器平滑道路面的强度信息；\n[0078] 沿着种子点到道路面边缘，由于道路点的局部密度逐渐变得稀疏，使用一个固定的滤波窗口，可能会导致稀疏区域的标志边缘点的强度值被过滤，还可能导致稠密区域的标志线平滑效果不足。为达到最佳的平滑效果，我们采用了一种基于局部密度的自适应窗口中值滤波器。\n[0079] 本步骤具体实施如下：\n[0080] 步骤4.1，在伪扫描线内，构建k‑d树索引每个道路点的邻域并统计最近邻点的数量，基于不同的局部密度选择不同尺寸的窗口，见公式(3)：\n[0081]\n[0082] 公式(3)中，Sw为中值滤波窗口的尺寸；Smax为最大窗口尺寸，根据经验取值，一般设置为7；ζ为偶系数，根据经验取值，一般设置为2；pn为局部点云的数量，根据经验取值，N1和N2一般设置为5和10。\n[0083] 步骤4.2，在滤波窗口中，将窗口的强度中值替代窗口中心的强度值，从而达到平滑强度噪声的目的。在密度稀疏区域，选择小尺寸滤波窗口；在密度稠密区域，选择大尺寸滤波窗口。这样，根据不同的局部密度，自动选择不同的窗口尺寸，即可有效地抑制道路面的强度噪声，又保护标志边缘点的强度值，使得后续处理中计算强度梯度和构建强度统计直方图变得更加稳健。如图5和图6所示，伪扫描线内点云的强度值经自适应窗口中值滤波后，强度噪声被有效抑制，同时标志线边缘点的强度值被保护。\n[0084] 步骤5：采用多约束方法检测道路和标志线之间的边缘；\n[0085] 经过基于局部密度的自适应窗口中值滤波器处理后，标志点的强度值变得较为平滑。如图6所示，由于多线系统点云的复杂性，即使道路面经过平滑处理，仍有小部分噪点没有被彻底地抑制，实用强度梯度的方法依旧可能会将这些噪点错误地检测为标志线边缘点。为了准确地检测伪扫描线中标志线的边缘点，我们提出了一种基于强度梯度和强度统计直方图的多约束检测方法。\n[0086] 步骤5.1，在平滑后的伪扫描线中，计算各道路点的强度梯度，见公式(4)，并根据平滑后的道路面的强度值构建强度统计直方图。\n[0087] Gi＝Ii‑Ii‑k   (4)\n[0088] 公式(4)中，Gi为第i个点的强度梯度；Ii和Ii‑k为第i和第i‑k个点的平滑后的强度值，i＝1,2,...,n；k为梯度的尺寸，根据经验取值，一般设置为3；\n[0089] 步骤5.2，沿着伪扫描线从左到右，进入标志线的边缘被定义为左边缘，离开标志线的边缘被定义为右边缘。基于强度梯度和强度统计直方图的多约束方法，标志边缘点被逐点检测，见公式(5)：\n[0090]\n[0091] 公式(5)中，pth为非零正梯度阈值；nth为非零负梯度阈值；ξth为从强度统计直方图提取的强度阈值，其表征了大多数标志线的强度值；由于不同激光雷达系统采集的数据，其强度值量级具有非常大的差异性，pth、nth和ξth应根据采集数据的强度特征进行设置。\n[0092] 在每条伪扫描线中，左边缘的点为标志线点，其强度梯度大于非零正值，且强度值接近大多数标志线的强度值；右边缘的点是非标志线点，其强度梯度小于非零负值，且强度值小于大多数标志线的强度值。通过检测成对的左右边缘点，并提取左右边缘之间的点云即为标志线点，因此通过逐伪扫描线检测，所有的标志线能够被直接提取。\n[0093] 在该方案中我们使用了多个约束条件，即强度提取和标志线强度特征。强度梯度能有效检测强度值突变的点云，但容易受到噪点的干扰。为了准确探测真实的标志线边缘点，还使用了标志线强度特征。由于数据的复杂性，标志线的强度值是不统一的。一个恒定的强度阈值，无法适应不同的系统数据，因此我根据统计学原理，计算不同系统数据标志线的强度特征，提取道路点与标志线之间的边缘点强度特征。通过结合强度梯度和标志线强度特征精确地检测标志线边缘点，避免了标志线内部噪点的干扰，普适于具有不同强度值量级的标志线。\n[0094] 步骤6：基于标志线的结构特征和维度特征剔除非标志点；\n[0095] 由于道路面和标志线强度的相似性，提取的标志线中存在小部分非标志点。这些噪点表现为两种形式：(1)标志线外的离散段，如图7中#1所示；(2)标志线内的线性段，如图\n7中#2所示。\n[0096] 本步骤具体实施如下：\n[0097] 步骤6.1，基于欧式距离聚类，将标志点分割成单独的标志类。标志线能被分割成多条伪扫描线，而离散段通常只能被分割成少量的伪扫描线。根据该结构特征，如果标志类被分割成伪扫描线的数量小于数量阈值 则该标志类被判断为噪点，给予剔除。计算该数量阈值，见公式(6)：\n[0098]\n[0099] 公式(6)中， 为伪扫描线的数量阈值；wth为上文的伪扫描线宽度阈值；lmin为标志线的最短长度，沿着行驶方向，中国最短的标志线长度为0.2m。\n[0100] 步骤6.2，在邻域范围内，真实的标志点构成二维平面，而噪点构成一维线段。构建k‑d树索引每个点云的邻域，并计算该点的线性度，见公式(7)：如果该点的线性度大于特征阈值lth，则该点被判定为噪点，给予剔除。根据经验取值，lth一般设置为0.95。\n[0101]\n[0102] 公式(7)中，linearity为线性度；λ1、λ2和λ3是协方差矩阵的三个特征值，且λ1≥λ2≥λ3≥0。\n[0103] 本实施例针对不同类型的移动激光雷达数据(即：单线系统点云和多线系统点云)，公开了一种通用的道路面与标志线自动提取方法。利用统计滤波法、GNSS轨迹数据和点云的高程选取道路面候选点，减少点云体积；根据轨迹直线和点云的垂足点分割候选点成不同的伪扫描线，实现点云的有效组织；根据道路面和路缘石(或排水沟)的几何特征，实现道路面提取；采用基于密度的自适应窗口中值滤波器平滑道路面的强度信息，增强算法的稳健性；通过强度梯度和强度统计直方图的多约束方法，避免了扫描距离和扫描角度对点云强度值的影响，实现标志线直接提取；基于标志线的结构特征和维度特征，剔除标志线中的非标志点，最终实现了道路面与标志线的自动提取。\n[0104] 实施例二\n[0105] 本实施例与实施例一的区别在于，采用北斗卫星导航系统获取轨迹数据和点云数据。\n[0106] 实施例三\n[0107] 本实施例与实施例一的区别在于，在实现道路面的提取步骤中，对道路面与道路面边缘的高度差进行定义。道路面边缘包括但不限于，道旁树、地砖、盲道、路灯等。\n[0108] 每个高度差对应至少一种道路面预设宽度范围。在道路面提取过程中，通过对道路面与道路面边缘的高度差计算，提取该高度差对应的道路面预设宽度范围，将该道路面预设宽度范围与已经测得的道路面边缘之间的距离进行比较，在比较时，通过同一段道路的至少三个点进行比较。若道路面边缘之间的距离小于道路面预设宽度范围，则判断该道路面是否被道路边缘遮挡，如果判断是被遮挡，则将前后一定长度内的道路边缘点拟合成两条光滑曲线，再通过聚类和追踪的方法连接这两条曲线，使局部的道路面边缘形成一条完整曲线，并自动通过预设密度范围对路面的缺口进行补偿，将补偿后的道路面作为被提取的道路面。\n[0109] 本实施例把路旁常见的道旁树、地砖、盲道、路灯等相关要素作为边界标示，结合具体区域的参考数据(如路灯距离路肩的距离)，来进行特定环境的路面提取。还可以解决路旁停车遮挡路肩等问题。\n[0110] 实施例四\n[0111] 本实施例与实施例一的区别在于，在实现标志线的提取步骤中，对标志线的类型进行定义。标志线的类型包括但不限于，矩形标志、箭头标志等。\n[0112] 每个标志线至少对应一种预定义标志线类型。在提取标志线的过程中，通过欧式距离聚类的方法，将提取的标志线点聚类成单独的标志类。首先，将标志线不同边缘的点拟合成直线；然后，计算每条拟合边缘的趋向，并计算依次两相邻拟合边缘之间的夹角。如该标志线任意两条相邻边缘的趋向和夹角，呈近似垂直状态，则判定其为矩阵标志。此外，将该矩形标志的重心坐标，与道路面预设网格坐标范围比较，根据重心坐落的区域，将该标志分类为道路中心线、道路边界线或车道分界线等。如该标志线的任意两条相邻边缘的趋向和夹角，存在非垂直状态，则判定其为箭头标志。此外，结合语义知识该标志线被进一步分类，若存在圆弧边缘，则分类为掉头线；若存在钝角转折边缘，则分类为转弯线等，根据其他几何特征箭头标志被分类但不局限于掉头线、转弯线、出口线或方向线等。\n[0113] 本实施例把道路面中常见的矩形、箭头标志等相关要素作为分类表示，结合具体几何特征，来进行特定类型的标志线提取。还可以解决标志线功能分类等问题。\n[0114] 以上说明书中的具体实施方式等记载可以用于解释权利要求的内容。