基于点密度与卡尔曼滤波的路面标识线提取方法.pdf

第 40 卷第 6 期2023 年 11 月中 国 科 学 院 大 学 学 报Journal of University of Chinese Academy of SciencesVol.40NovemberNo.62023国家自然科学基金(61771456)资助 通信作者,E-mail:文章编号:2095-6134(2023)06-0800-10基于点密度与卡尔曼滤波的路面标识线提取方法李晓宇1,2,周梅1,王金虎1,姚强强1(1 中国科学院空天信息创新研究院 中国科学院定量遥感信息技术重点实验室,北京 100094;2 中国科学院大学光电学院,北京 100049)(2022 年 2 月 17 日收稿;2022 年 3 月 28 日收修改稿)Li X Y,Zhou M,Wang J H,et al.Road marking detection and extraction method based on neighborhood density and Kalman filterJ.Journal of University of Chinese Academy of Sciences,2023,40(6):800-809.DOI:10.7523/j.ucas.2022.024.摘 要 提出一种基于邻域点密度和卡尔曼滤波的车载激光点云路面标识线提取与补全方法。该方法由 3 个部分组成:1)利用时间信息将路面点云数据分割为单条扫描线;2)根据邻域点密度自动生成卷积核以提取标识线轮廓点;3)结合最小二乘算法拟合标识线轮廓,并利用卡尔曼滤波算法填补缺损标识线。经实验验证,最终标识线提取结果的平均中心点位置偏差为 0.04 m,平均方向偏差为 0.04,平均完整度为 99.69%。该方法降低了点云密度分布不均对路面标识线检测精度的影响,可有效提高标识线检测的完整度。关键词 车载激光扫描;点云;卷积;路面标识线;卡尔曼滤波 中图分类号:TP75 文献标志码:A DOI:10.7523/j.ucas.2022.024Road marking detection and extraction method based onneighborhood density and Kalman filter LI Xiaoyu1,2,ZHOU Mei1,WANG Jinhu1,YAO Qiangqiang1(1 CAS Key Laboratory of Quantitative Remote Sensing Information Technology,Aerospace Information Research Institute,Chinese Academyof Sciences,Beijing 100094,China;2 School of Optoelectronics,University of Chinese Academy of Sciences,Beijing 100049,China)Abstract This paper presents a methodology for detection and extraction of road marking based on neighborhood point density and Kalman filter from mobile laser scanning in three steps:1)Segmenting road point clouds into scan lines;2)Generating convolution kernel based on point density in neighborhood for extracting road marking contour points;3)Fitting contour line using least squares algorithm and completing omitted road markings using Kalman filter.A quantitative validation shows that the average deviation of center points and orientation are 0.04 m and 0.04,respectively,and the average completeness of results are 99.69%.The proposed method reduces the influence of unevenly distributed points on road marking extraction and improves the overall extraction completeness.Keywords mobile laser scanning;point clouds;convolution;road markings;Kalman filter 第 6 期李晓宇,等:基于点密度与卡尔曼滤波的路面标识线提取方法 路面信息的精确检测是自动驾驶技术的重要组成部分1,其中线型交通标识线(包括车道边缘线和车道分界线)的检测可对车辆行为起到引导、规范和禁止作用2。车载激光扫描系统是一种集成全球卫星导航系统、惯性导航系统、激光扫描仪和数码相机等多种传感器的移动三维激光扫描系统(mobile laser scanning,MLS)3,该系统不受光照限制,可高效获取路面高精度稠密三维点云数据,为线型交通标识线的检测提供了可靠数据基础。当前基于点云数据的路面交通标识线提取方法可分为两类:将点云投影为特征图像再提取的方法;基于离散点云的方法。第 1 类方法是根据点云回波强度或空间分布特征将点云投影为特征图像,再结合阈值分割4、语义分割5和神经网络模型6等方法提取路面交通标识线。该类方法提取的路面交通标识线往往缺乏高程信息7且将点云投影为图像的过程容易造成精度损失。第 2 类方法直接基于离散点云提取路面交通标识线,Yu 等8和 Kim 等9提出基于点云回波强度的多阈值分割方法,通过沿行车轨迹垂直方向分割点云并结合 Otsu10方法自适应设定强度阈值提取路面交通标识线,该方法通过对点云数据分段设定阈值,减少了扫描距离对路面交通标识线提取精度的影响,但其中点云分段方法的鲁棒性不佳,容易漏提取强度较低的交通标识线。尹辉增和方莉娜11对点云回波强度进行直方图统计,根据统计结果划定强度阈值提取交通标识线点,并结合交通标识线的几何、语义特征细化提取结果。邹晓亮等12根据交通标识线点的回波强度分布范围划定强度阈值,并利用扫描角和扫描距离等信息判定点的交通标识线类型。上述两种方 法自动化程度较低,难以识别受遮挡或磨损的路面交通标识线。Yan 等13提出基于扫描线的路面交通标识线提取方法,通过检测每条扫描线内相邻两点间的强度变化确定交通标识线点,该方法精度和自动化程度较高,但提取精度易受噪声影响,此外提取结果的完整性因点云缺失或交通标识线的磨损而降低。针对上述问题,提出一种基于邻域点密度和卡尔曼滤波14的路面线型交通标识线提取与补全方法。首先,为避免入射角和扫描距离变化造成的点云密度不均现象对提取结果精度的影响,以一阶微分算子为卷积模板,根据单条扫描线内每个点邻域内的点密度分布自适应生成对应尺寸卷积核,通过对邻域点的强度值进行卷积以提取交通标识线轮廓点。随后,在轮廓点提取结果基础上利用最小二乘算法拟合线型交通标识线轮廓,并结合卡尔曼滤波算法补全点云数据缺失或标识线磨损造成的轮廓提取结果不完整的区域,以提高线型交通标识线提取结果的完整度。1 方法原理 本文提出的路面线型交通标识线提取算法总体流程如图 1 所示,在预处理模块,根据扫描点的时间信息将路面点云数据分割为多条扫描线。随后在轮廓点提取模块,计算扫描线内每个点对应的卷积核,对点的强度值进行自适应卷积;根据每条扫描线卷积结果生成的阈值和强度变化规律确定交通标识线轮廓点位置,利用空间密度聚类法(density-based spatial clustering of applications with noise,DBSCAN)15剔除轮廓点提取结果中的噪声。最后在轮廓拟合与填补模块中,基于轮廓点提取结果和最小二乘法拟合标识线轮廓,并结合 图 1 路面标识线提取与填补算法总体流程Fig.1 Overall workflow of the proposed road marking extraction and completion algorithm108中国科学院大学学报第 40 卷卡尔曼滤波算法预测缺损区域的轮廓位置,最终实现线型交通标识线的完整提取。1.1 路面交通标识线轮廓点提取 图 2(a)为某城区路面点云数据,经预处理后的扫描线分割结果如图 2(b)所示,其中每条扫描线被随机赋色;图 2(c)中红色线为路面中某条扫描线数据,矩形框内为该扫描线中经过路面标识线的部分,其对应的强度分布曲线如图 2(d)所示,强度曲线中路面交通标识线点回波强度(图 2(d)矩形框内点)明显高于沥青路面点的回波强度,在二者衔接处点的回波强度会出现明显的上升或下降。本文以此为依据,提取交通标识线轮廓点(图 2(d)中红色点)。1.1.1 自适应卷积 使用以一阶微分算子为模板的卷积核,对每条扫描线点及其邻域内点的强度值进行卷积,卷积结果中数值绝对值较大的极值点对应交通标识线轮廓点。由于一阶微分算子对噪声敏感,在卷积之前使用高斯滤波器抑制噪声。如图 3(a)中路面点云数据所示,一条扫描线分别经过宽度相同的车行道边缘线(区域)和车道分界线(区域),由于扫描距离与入射角度的影响,区域中点分布较稀疏,而区域中点分布更为密集。若使用单一尺寸的卷积核提取轮廓点,容易降低提 取结果的位置精度。为验证不同尺寸卷积核对轮廓点提取精度的影响,本文采用尺寸为 15、115 的卷积核分别对区域和区域中扫描线点进行卷积,区域和区域的卷积结果如图 3(b)、3(c)所示,在区域中 15 卷积核卷积结果的极值点(图 3(b)中三角形)更准确对应原始强度曲线的标识线轮廓点位置,然而在区域中 115 卷积核卷积结果的极值点(图 3(c)中五角星)更准确对应原始强度曲线的标识线轮廓点位置。由此可得,在提取轮廓点位置的过程中,扫描线内不同的点密度适用于不同尺寸的卷积核,在密度分布不均的扫描线内使用单一尺寸卷积核进行卷积,将使得轮廓点定位精度下降。因此,本文提出一种根据邻域点密度自适应确定卷积核尺寸的方法,简称为自适应卷积。以单条扫描线的卷积过程为例,具体流程如下所示:1)单条扫描线内有 n 个点 p1,p2,pn,以其中第 i 个点 pi为圆心、距离 r 为半径,搜索此邻域内的所有点 pi1,pi1,pin 共 m 个。为保证轮廓点提取结果精度,本文设定搜索半径 r 长度等于扫描线中标识线宽度。2)确定邻域内每个点的相对位置xj=j-m2(j 0,m).(1)图 2 路面点云数据与单条扫描线强度分布特征Fig.2 Point cloud data of road surface and intensity distribution of a scanline208第 6 期李晓宇,等:基于点密度与卡尔曼滤波的路面标识线提取方法图 3 路面中某条包含交通标识线的扫描线及其对应的强度卷积结果Fig.3 A scan line containing multiple road markings and its corresponding intensity convolution result其中 xj为邻域内第 j 个点在所有邻域点中的相对位置。3)将邻域内每个点的相对位置代入下式,计算其对应的初始卷积核数值kj=-xj212 e-xj2212xj-m2,m2().(2)式(2)为高斯一阶导函数,其中 kj 表示邻域内 m个点中第 j 个点对应的初始卷积核值。如图 4 所示,本文设定扫描线中点 pi邻域内的 m 个点的范围约等于高斯一阶导函数的(-31,31)范围,因此高斯一阶导函数中的方差 1=m6。4)为避免卷积核尺寸变化对结果的影响,将邻域内点对应的初始卷积核数值进行归一化,得到点 pi邻域内 m 个点对应的最终卷积核 k0,k1,kmkj=kjm/20|kj|,j (0,m).(3)式中 kj为 kj 归一化后的结果。5)对点 pi邻域内 m 个点的强度值进行卷积,得到点 pi对应的卷积结果 cici=mj=0Ij kj.(4)图 4 扫描线中某点 pi处根据高斯一阶导函数模板生成卷积核的示意图Fig.4 Schematic diagram of convolution kernel at a point piin a scan line according to the Gaussianfirst-order derivative function template式中 Ij为邻域内第 j 个点的强度值。6)重复上述步骤直至得到扫描线内 n 个点的卷积结果 c1,c2,cn。1.1.2 轮廓点提取 路面交通标识线轮廓点在卷积结果中表现为数值绝对值较大的极值点,因此本文通过设定阈值在卷积结果中提取可能的交通标识线轮廓点,具体方法如下:1)计算扫描线卷积结果 c1,c2,cn 的均值 、方差 2。308中国科学院大学学报第 40 卷2)由于扫描线卷积结果近似满足正态分布,根据 3 准则,本文设定阈值 T=+a 2,其中常数 a 与扫描线内标识线点和沥青路面点的数量相关。经多次试验发现,当扫描线内标识线点数远小于路面点数时,常数 a 为 3 的提取效果最佳。此外可根据比例关系适当调整。3)若点 pi对应的卷积结果 ci大于阈值 T,则点 pi为可能的交通标识线轮廓点。经上述阈值提取的点为强度发生明显变化的点,然而该类点中除交通标识线轮廓点,还包含如路面裂缝或路面拼接处强度发生明显变化的点。因此,本文通过交通标识线轮廓点的强度变化规律和空间距离信息进一步优化提取结果。如图 2(c)、2(d)所示,每个交通标识线在扫描线内均存在一对强度的上升点和下降点(图 2(d)中红色点),在卷积结果中则对应为一个极小值和一个极大值,且一对强度变化点之间的欧式距离应近似等于交通标识线的宽度。经上述规则优化后,提取结果中的点即可基本确定为路面交通标识线轮廓点。最后,为保证提取结果的精确性,本文选择 DBSCAN 算法进一步滤除提取结果中的噪声点。1.2 路面交通标识线轮廓拟合与填补 由于路面交通标识线受磨损区域与沥青路面的强度对比度低、物体遮挡造成的路面点云数据缺失等原因,路面交通标识线提取结果往往存在缺失。为提高提取结果完整性,本文在路面交通标识线轮廓点提取结果基础上利用最小二乘算法 拟合路面中线型交通标识线轮廓,并结合卡尔曼滤波算法预测缺失区域的标识线轮廓位置,进而提取出缺失区域的路面交通标识线点。以图 5(a)的路面点云数据为例,路面中的交通标识线包含 2 条车行道边缘线和 2 条车道分界线,按图中方向自上而下分为 4 条车道线 Line1Line4。图 5(b)为交通标识线的轮廓点提取结果,其中 Line1、Line2 和 Line4 中存在遮挡或漏提取造成的轮廓点缺失。经轮廓拟合后(图 5(c)中黑色矩形框所示),Line1 和 Line4 中的缺失部分可以直接通过拟合步骤进行填补,而 Line2 中缺失的虚线标识(图 5(c)中 Mloss所在位置)仍存在缺失。因此本文使用卡尔曼滤波算法进行进一步填补,具体流程如下:1)确定车道线中缺失区域的位置 Mloss,沿道路方向分别向前、后搜寻完整标识线点集,其中前向搜索的 p 个标识线为 Mfront=Mf1,Mf2,Mfp,后向搜索的 q 个标识线为 Mback=Mb1,Mb2,Mbq,在图 5(c)中,p 和 q 分别为 3 和 2。2)以前向为例,计算所有前向搜寻的完整标识线点集中相邻 2 个标识线点集之间的中心点间距 d=d1,2,d2,3,dp-1,p 和方向向量 l=l1,2,l2,3,lp-1,p。其中,dp-1,p为前向搜寻的第 p 个标识线点集 Mfp和第 p-1 个标识线点集Mfp-1之间的中心点间距,lp-1,p为第 p 个标识线点集 Mfp的中心点和第 p-1 个标识线点集 Mfp-1的中心点之间的方向向量。图 5 路面交通标识线拟合与填补方法示意图Fig.5 Schematic diagram of road marking fitting and completion408第 6 期李晓宇,等:基于点密度与卡尔曼滤波的路面标识线提取方法 3)将上述计算获得的中心点间距和方向向量作为观测值 Z=dl(),从距离缺失标识线区域最远的观测值位置处开始(如图 5(c)中 Mf3),向缺失区域(如图 5(c)中 Mloss所在位置)进行迭代,预测 Mf1与 Mloss之间的中心点间距 d和方向向量 l(,具体计算公式如下:x-t=x-t-1.(5)P-t=FPt-1FT+Q.(6)Kt=P-tHTHP-tHT+R()-1.(7)x t=x-t+KtZt-Hx-t().(8)Pt=(I-KtH)P-t.(9)其中,式(5)为状态预测公式,x-t表示迭代至第 t次的预测值;式(6)为误差矩阵的预测公式,P-t表示迭代至第 t 次时由 t-1 次的误差矩阵 Pt-1预测得到的误差矩阵,F 为状态转移矩阵,Q 为预测噪声协方差矩阵;式(7)为卡尔曼增益计算公式,Kt表示第 t 次的卡尔曼增益,H 为观测矩阵,R 为测量的噪声协方差矩阵;式(8)为状态矫正公式,结合观测值对预测结果进行矫正,其中 Zt表示第t 次的观测值;式(9)为误差矩阵更新公式,I 表示单位矩阵。4)根据完整标识线点集 Mf1的矩形轮廓角点位置和预测得到的 Mf1与 Mloss之间的中心点间距 d和方向向量 l,计算 Mloss的轮廓角点位置。5)重复步骤 2)到 4)对后向搜寻的标识线点集进行相同操作,获取 Mloss后向预测得到的轮廓角点位置。6)按前向和后向迭代预测次数对前向、后向预测的角点位置进行加权,获得最终的缺失区域标识线轮廓位置。2 结果与分析 2.1 实验数据 为验证本文提出的线型标识线提取与补全方法的可行性,采用 RIEGL VMX-2HA 双扫描仪测图系统采集的北京城区道路点云数据进行实验(实验数据由中测瑞格公司提供),并通过人工去除部分路面信息后作为原始实验数据(图 6(c)所示)。实验数据的路面全长约为 500 m,宽 15 m。数据中点密度由路中央向道路两侧逐渐递减,在路中央处点密度最高可达 10 000 pts/m2(pts 为点的个数),而在道路两侧点密度约为 3 000 pts/m2;在单条扫描线内,点密度在路中央处可达 120pts/m,在道路两侧仅为 20 pts/m。路面包括车道分界线(虚线)和车道边缘线(长实线),同时还存在多个裂缝(如图 6(c)中区域)和道路拼接造成的标识线中断(如图 6(c)中区域);在该路 图 6 实验数据Fig.6 Test data508中国科学院大学学报第 40 卷面区域,车道分界线(图 6(c)中区域)与车道边缘线宽度(图 6(c)中区域)相同,但车道边缘线处点云分布更为稀疏。2.2 实验结果分析与精度评价2.2.1 标识线轮廓点提取结果 图 7(a)为实验的路面点云数据,使用大小为15、115 的固定尺寸卷积核和本文提出的基于邻域点密度的自适应尺寸卷积核提取标识线轮廓点,提取结果分别如图 7(b)、7(c)、7(d)所示,并通过轮廓点提取内部标识线点(内部标识线点提取效果如图 7(e)所示),以人工提取的标识线点为参考数据进行精度评价,以准确率 P、召回率 R和 F1-Measure 作为评价指标,间接评价单一尺寸卷积核与自适应卷积核提取轮廓点的精度差异(如表 1 所示)。P=TPTP+FP.(10)R=TPTP+FN.(11)F1-Measure=2 P RP+R.(12)式中,TP为正确提取的标识线点数,FP为错误提取的标识线点数,FN为未被提取的标识线点数。如图 7 与表 1 所示,在使用 15 尺寸卷积核的轮廓点提取结果中,由于小尺寸的卷积核对噪声敏感,因此路中央点分布密集区域存在大量噪声点,其准确率仅为 63.14%,路面上侧的车道边缘线点由于分布稀疏而被漏提,召回率(即被正确分类的点)仅为 55.28%;当卷积核尺寸增大为1 15 时,噪 声 点 的 减 少 使 得 准 确 率 提 升 为 表 1 基于不同尺寸卷积核的路面交通标识线轮廓点提取精度Table 1 The extraction accuracy of road marking pointsbased on convolutional kernels of different size%卷积核尺寸RPF1-Measure1555.2863.1458.9511541.0397.2457.71自适应93.1099.1596.03图 7 实验数据的标识线轮廓点提取结果和标识线点提取效果示意图Fig.7 Contour points extraction result(blue points in the figure are contour points)and the extracting effect of road marking point608第 6 期李晓宇,等:基于点密度与卡尔曼滤波的路面标识线提取方法97.24%,然而尺寸过大的卷积核使得道路两侧的车道 边 缘 线 漏 提 现 象 更 明 显,召 回 率 仅 为41.03%;使用自适应卷积的轮廓点提取结果中,在保证准确率的同时,轮廓点提取结果的完整性明显提高,召回率为 93.10%,综合评价指标 F1-Measure 可达到 96.03%。2.2.2 标识线轮廓拟合与填补结果 在轮廓点提取结果之上,本文对标识线轮廓进行了拟合与填补,如图 8 所示,为验证其效果,自上而下分将标识线分为 Line1Line4,以标识线中心点偏差 L(绝对误差)、标识线方向偏差 d(绝对误差)和标识线拟合完整度 (即标识线拟合长度与真实长度的相对误差)作为精度评价指标,对轮廓拟合(如图 8(a)和填补结果(如图 8(b)进行精度评价,结果如表 2 所示。如表 2 所示,标识线 Line1 Line4 的轮廓拟合结果中,平均的中心点偏差约为 0.04 m,拟合方向与真实方向之间平均的方向偏差为 0.03,由于 Line1 与 Line4 的缺失区域可通过拟合直接填补,其标识线拟合完整度达到 99%以上,而Line2 与 Line3 中的拟合完整度仅为 91.14%和90.73%。经卡尔曼滤波方法填补后,标识线Line1Line4 中平均的中心点偏差为 0.04 m,平均的方向偏差为 0.04,无明显变化,而 Line2 与Line3 的拟合完整度提升为 99.28%和 99.68%,总体平均拟合完整度从填补前的 95.42%提升至99.69%。为进一步验证填补方法的准确性,将缺失区域的填补结果(图 8(b)中区域 所示)和该区域的真实路面标识线进行比较,如表 3 所示,在区域的填补结果中,平均的中心点偏差和方向偏差分别为 0.03 m 和 0.21,平均拟合完整度为 99.37%。上述结果表明,经过轮廓拟 合和卡尔曼滤波方法填补后的标识线轮廓精度较高,且能够有效提取缺损区域的标识线轮廓,减少数据缺失对提取结果完整性的影响。上述实验通过对数据中的全部扫描线进行处理以提取轮廓点,进而对标识线轮廓进行拟合与填补。但当数据量较大时,为保证该方法运算效率,可以对扫描线进行抽稀,由于扫描线排列紧密,因此数据抽稀后的轮廓点提取结果对后续轮廓拟合与填补精度影响较小。如图 9 所示,本文方法运算时间随着数据量的减少而逐渐下降,与对所有扫描线进行处理的时间相比,每间隔 1 条和 2 条扫描线进行数据抽稀后,其运算时间明显表 2 标识线轮廓拟合与填补后的精度Table 2 The accuracy of road marking fitting and completion results方法交通标识线中心点偏差L/m方向偏差d/()拟合完整度/%轮廓拟合Line10.040.0099.90Line20.030.0891.14Line30.020.0590.73Line40.050.0099.91轮廓填补Line10.040.0099.90Line20.030.0899.28Line30.020.0899.68Line40.050.0099.91表 3 标识线缺失区域的填补精度Table 3 The accuracy of road marking completionresults in missing areas缺失区域中心点偏差L/m方向偏差d/()拟合完整度/%区域0.040.3499.58区域0.020.1698.77区域0.020.0399.54区域0.030.3099.60蓝色点为标识线轮廓点,区域为需要填补的区域。图 8 轮廓拟合与填补结果Fig.8 Contour fitting and completion results708中国科学院大学学报第 40 卷图 9 不同数据量的运算时间对比Fig.9 The comparison of computingtime for different data sizes降低;当间隔 2 条以上扫描线抽稀数据后,其运算时间仍随数据量减少而降低,但时间变化不明显。因此,在实际应用中为保证运算效率,可根据数据量和提取精度的需求对数据进行抽稀。2.2.3 标识线提取方法对比分析 本文提出的方法与已有方法的对比如表 4 所示。在提取路面标识线前,Cheng 等16和 Ma等6提出的方法需要在将路面点云数据栅格化 为二维图像,并结合图像增强和去噪算法提高图像对比度以提取多种路面标识,提取结果的 F1-Measure 为 94.00%和 92.43%;Yang 等17提出的方法需要将路面点云数据进行体素化,再结合机器学习分类器进行线型标识线的提取,提取结果的 F1-Measure 为95.10%;Yang 等18提出的方法需预先将路面数据分块,且不同的分块策略极大影响标识线提取结果的精度,鲁棒性不强,文中多种路面标识提取结果的 F1-Measure 为 95.70%;相比于上述方法,本文提出的方法仅需要对点云进行扫描线分割,实现简单,且对后续步骤的提取精度无影响。Yan 等13提出的线型标识线提取方法的 F1-Measure 为 94.00%,该方法虽与本文具有相似之处,但本文方法还考虑了点密度分布不均、标识线磨损等因素,适用于多种道路环境。此外,上述方法未考虑缺损标识线的填补问题,且提取的标识线结果多为像素或点集的形式,而本文提出了缺损标识线填补策略,经填补后的最终标识线提取结果为多个独立的对象实体。综上所述,本文方法具有较强的鲁棒性,可在不同质量的道路环境下进行标识线的完整提取,具有一定的工程应用价值。表 4 本文方法与已有方法对比Table 4 Comparison of the proposed method with existing approaches方法预处理数据类型原理缺损填补输出Cheng 等16点云栅格化、图像增强图像阈值分割无像素Ma 等6点云栅格化、高通滤波去噪图像神经网络无像素Yan 等13扫描线分割点云扫描线内强度变化无点Yang 等17体素化点云局部特征描述子与随机森林分类器无点Yang 等18数据分块点云阈值分割无点本文方法扫描线分割点云扫描线内强度变化有对象3 结论 为降低点云密度不均匀对路面交通标识线提取精度的影响,本文提出一种新的标识线提取及补全算法。该方法根据邻域点密度自动生成对应尺寸卷积核,通过卷积实现了点云强度变化位置(即标识线轮廓点)的检测;针对点云数据缺失和标识线磨损造成的提取结果不完整的问题,本文在标识线轮廓点基础上结合最小二乘方法拟合标识线轮廓,并结合卡尔曼滤波算法填补缺损标识线信息。经验证,本文提出的标识线轮廓点提取方法的准确率可达 99.15%,经卡尔曼滤波填补后,文章提取的标识线轮廓的平均拟合完整度可达 99.69%,与真实标识线的中心点位置和方向相比,拟合结果中平均的中心点偏差和方向偏差仅为 0.04 m 和 0.04。综上所述,本文提出的方法有效改善了路面标识线磨损、车辆遮挡、点云密度分布不均等因素对路面标识线提取结果精度的影响,实现了线型路面标识线的快速、精确提取。目前本文仅讨论了线型标识线的拟合与填补方法,如何基于标识线轮廓点实现多种类型的路面标识线拟合与填补有待进一步研究。808第 6 期李晓宇,等:基于点密度与卡尔曼滤波的路面标识线提取方法参考文献 1 Tan H C,Zhou Y,Zhu Y,et al.A novel curve lane detection based on improved river flow and RANSAC 17th International IEEE Conference on Intelligent Transportation Systems(itsc).October 8-11,2014,Qingdao,China.IEEE,2014:133-138.DOI:10.1109/ITSC.2014.6957679.2 方莉娜,黄志文,罗海峰,等.车载激光扫描数据中实线型交通标线提取J.测绘学报,2019,48(8):960-974.DOI:10.11947/j.AGCS.2019.20180579.3 牟春霖.车载激光雷达在城市轨道交通勘测中的应用J.铁道勘察,2020,46(6):44-49.DOI:10.19630/ki.tdkc.202007140004.4 Liu L R,Ma H,Chen S Y,et al.Image-translation-based road marking extraction from mobile laser point cloudsJ.IEEE Access,2020,8:64297-64309.DOI:10.1109/ACCESS.2020.2985413.5 Yang B S,Fang L N,Li Q Q,et al.Automated extraction of road markings from mobile Lidar point clouds J.Photogrammetric Engineering&Remote Sensing,2012,78(4):331-338.DOI:10.14358/PERS.78.4.331.6 Ma L F,Li Y,Li J,et al.Capsule-based networks for road marking extraction and classification from mobile LiDAR point cloudsJ.IEEE Transactions on Intelligent Transportation Systems,2021,22(4):1981-1995.DOI:10.1109/TITS.2020.2990120.7 闫利,李赞.车载激光点云道路标线提取方法J.遥感信息,2018,33(1):1-6.DOI:10.3969/j.issn.1000-3177.2018.01.001.8 Yu Y T,Li J,Guan H Y,et al.Learning hierarchical features for automated extraction of road markings from 3-D mobile LiDAR point cloudsJ.IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2015,8(2):709-726.DOI:10.1109/JSTARS.2014.2347276.9 Kim H,Liu B B,Myung H.Road-feature extraction using point cloud and 3D LiDAR sensor for vehicle localizationC 2017 14th International Conference on Ubiquitous Robots and Ambient Intelligence(URAI).June 28-July 1,2017,Jeju,Korea(South).IEEE,2017:891-892.DOI:10.1109/URAI.2017.7992858.10 Otsu N.A threshold selection method from gray-level histogramsJ.IEEE Transactions on Systems,Man,and Cybernetics,1979,9(1):62-66.DOI:10.1109/TSMC.1979.4310076.11 尹辉增,方莉娜.智能交通系统中道路标识线自动提取方法 J.科技 资讯,2012,10(26):167-168.DOI:10.16661/ki.1672-3791.2012.26.182.12 邹晓亮,缪剑,郭锐增,等.移动车载激光点云的道路标线自动识别与提取J.测绘与空间地理信息,2012,35(9):5-8.DOI:10.3969/j.issn.1672-5867.2012.09.002.13 Yan L,Liu H,Tan J X,et al.Scan line based road marking extraction from mobile LiDAR point clouds J.Sensors(Basel,Switzerland),2016,16(6):903.DOI:10.3390/s16060903.14 Kalman R E.A new approach to linear filtering and prediction problemsJ.Journal of Basic Engineering,1960,82(1):35-45.DOI:10.1115/1.3662552.15 Ester M,Kriegel H P,Sander J,et al.A density-based algorithm for discovering clusters in large spatial databases with noise C Proceedings of the Second International Conference on Knowledge Discovery and Data Mining.August 2-4,1996,Portland Oregon,America.AAAI Press,1996:226-231.16 Cheng M,Zhang H C,Wang C,et al.Extraction and classification of road markings using mobile laser scanning point cloudsJ.IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2017,10(3):1182-1196.DOI:10.1109/JSTARS.2016.2606507.17 Yang B S,Liu Y,Dong Z,et al.3D local feature BKD to extract road