基于移动激光点云的交通标志牌提取 - 图文

基于移动激光点云的交通标志牌特征提取

瓮升霞，陈一平1*

(厦门大学信息科学与技术学院,福建 厦门,361005)

摘要：针对有效采集到的道路两旁高密度点云数据，本文提出了一种基于移动激光测量技术进行户外交通标志牌的自动检测方法。首先阐述基于交通标志牌反射率的检测方法，利用交通标志牌在点云中的反射强度与其他地物目标区别较为明显的特性进行目标提取。再结合基于交通标志牌的形状特征检测方法，通过提取杆状物并分析杆状物中非杆状部分的维度特征，提取出具有面状维度特征的目标物作为所要检测的交通标志牌。最后综合检测结果，有效检测到点云数据中的交通标志牌，且召回率优于常用的基于反射率的方法，实验表明该方法能够取得好的提取结果。

关键词：移动激光扫描；交通标志牌检测；主成分分析

交通标志牌作为一种重要的交通基础设施，可以为驾驶员和行人提供道路指引信息，规范驾驶员行为的同时也起到了便利交通、保障安全的作用。鉴于此，交通部门定期对交通标志牌的测绘数据进行更新和维护。传统的测绘技术效率低，且受天气地形等因素影响。采用现代移动测量技术，移动激光车载系统沿道路行驶，可有效获得高精度点云数据和动态影像，极大的削减了人工成本，减少了人为误差。对于获得的点云数据，人工检测显得尤为繁琐，自动化提取越来越成为未来的发展方向，但是需要高效的检测率才能真正应用于实际场景中。

基于视频图像的交通标志牌检测已经进行了十多年的研究[1,2,3]，图像的采集需要较好的天气和光照条件对目标物的位置进行正确的识别。然而，点云数据获取即使在恶劣的天气条件下也能完成。与在交通工程中普遍应用的人工视觉系统相比，激光扫描对点云的采集更加准确，密度也更高[4,5]，例如，RIEGL VMX-450移动激光扫描系统每平方米可以采集到高达7000-8000多个点的三维点云数据。移动激光雷达系统可以完全捕获道路环境，包括真实的地理空间信息和回波反射能量，其中回波反射可以反映物体表面的材料特性。因此，利用移动激光雷达技术对城市道路目标进行采集是非常理想的。

近年来，基于移动激光雷达技术的道路目标提取也取得了一定的进展。Yokoyama等[6]

提出了一种基于拉普拉斯算子平滑和主成分分析的杆状物提取算法。在该方法中，拉普拉斯算子用来对数据进行平滑处理，以消除数据中的噪声点和数据点的分布偏差。而主成分分析过程则用于对平滑后的数据分段进行特征分析，并从中检测柱状结构目标。Shi Pu[7]等利用形状先验知识和拓扑结构来识别道路结构。董震[8]等通过从原始点云中生成多尺度超体素，在进行体素分割后根据分割区域的显著性进行邻域聚类，最后结合聚类区域的几何特性来判断目标所属类别。由于交通标识牌一般由高反射材料制成，具有清晰可见性，因此点云数据中的反射强度信息可以有效的将标志牌与道路其他目标区别开来，从而广泛的应用于交通标

[9,10]

志牌检测中。另外，基于模板驱动的方法也用于检测交通标志牌[11,12]，利用交通标志牌的高反射强度、对称性等特征来实现检测。

交通标志牌的检测依然存在着一些待解决的难点，主要在于：交通标志牌往往会被其他目标遮挡，如树枝，路灯等，导致检测难度增大；标志牌长期在户外，难免会有损毁老化， 收稿日期：2015-12-18 录用日期：2016-02-23

基金项目：重点实验室科技基金项目（9140C800401150C80339） *通信作者：ypchenhk@gmail.com

致使它的反射率达不到标准水平；杆部可能倾斜变形；在采集数据的过程中，一些离扫描设备比较远的标志牌在点云中的密度较低，影响检测效果。因此，一个鲁棒的检测算法应该克服这些难点，得到精确的检测结果。目前基于三维点云进行交通标志牌提取的研究主要分为基于分割和特征识别的方法和基于反射强度的方法，前者主要受限于分割算法，由于自然场景的复杂性，重叠在一起的目标很难被分割开。因此，对标志牌检测的算法性能产生一定的影响。另一方面，基于反射强度的方法也有其缺陷，由于该方法主要通过标志牌在点云中的高强度值与其他道路目标区分开来，如果交通牌因老化而褪色，或是数据采集时由标志牌的反面扫描，会导致点云中标志牌的反射强度缺乏明显的区分度，基于反射强度的检测方法就会发生漏检。

基于以上分析，本文综合运用两种方法进行交通标志牌检测，既可以弥补基于反射率的检测算法在交通牌的强度值与其他地物目标未明显区分时导致漏检的缺陷，又可以对两种方法的检测结果相互验证。

1方法

本文提出了基于反射强度和形状特征相结合进行交通标志牌检测的方法。基于反射强度的检测方法利用交通标志牌在点云中的高强度特性进行检测，并在后续聚类过程中定义规则以获得含有交通标志牌的聚类。基于形状特征的检测方法首先从点云中识别杆状物，通过分析杆状物中非杆状部分的维度特征，提取出具有面状维度特征的目标物作为所要检测的交通标志牌。最后将两次检测结果进行叠加和合并形成最终的检测成果。工作流程如下：

点云点云反射强度过滤地面滤除聚类聚类杆状物检测点云数量滤波海拔滤波高度滤波提取具有面状维度特征的目标物标志牌检测结果标志牌检测结果交通标志牌最终检测结果 图 1 交通标志牌检测流程图

Fig.1 The flowchart of traffic sign detection

1.1基于反射率的检测

1.1.1 反射强度滤波

扫描得到的点云数据一般包含空间坐标、反射强度以及颜色信息，其反射强度受扫描距离、扫描角度、目标材质及传播介质等因素影响。为了行车安全，交通标志牌使用高反射率的材料进行喷刷，反射率的值必须达到国家安全生产标准的规定值。理论上，激光回波强度与地物目标的反射率成正比关系，即反射率越大，点云反射强度值越高。因此在点云中，交通标志牌的反射强度比其他目标都要高。尽管点云的反射强度也会受大气衰减、激光入射角度等多种因素的影响，还是能够区分介质属性区别较明显的目标。通过将反射强度的阈值设定为超出某一值，道路场景中的大部分非目标物都将被滤除[10]。阈值的大小需要经过实验进行调整以获得最佳滤波效果。滤波后的点云仍包括除交通标志牌之外的其他高强度目标，如汽车牌照、路锥等，需要进行进一步处理。 1.1.2 欧式距离聚类

聚类的目的在于将离散、无组织的数据划分成一系列具有一定空间拓扑结构的三维语义目标，使每个目标中的数据“相似”或“相近”，而不同目标中的数据差异尽可能大。欧式距离聚类是一种根据相邻点之间的欧式距离关系进行划分的方法，其原理是：对于含有n个数据点的集合X?{x1,x2,??,xn}，随机选取1个点pi，并标记为已聚类点，将pi为球心，dc为半径的球形区域内的所有未聚类点标记为已聚类状态，并且给以相同的类别标识Li。接着，对于每一个被标记的点，重复上述过程，直到Li已无法再吸收新的点，在数据集合的未标记点中再次进行随机选取。当数据集中所有的点都已经被标记，该递归聚类过程结束。

初始聚类之后，部分聚类中包含感兴趣的目标，为了判断每个聚类 Ci是否包含交通标志牌，本文中定义了以下三种规则： 1) 点云数量滤波

结合交通标志牌的尺寸限制以及车载移动激光扫描点云的密度， Ci应该至少包含70个点，少于这个数量的聚类不足以构成交通标志牌。 2) 海拔滤波

依据交通标志牌的高度标准， Ci的z轴最小值与地面点的距离至少应为2m。这一规则可以过滤掉海拔较低的聚类目标，如汽车牌照，路锥，临时标志牌等。 3) 高度滤波

依据交通标志牌尺寸标准， Ci的z轴最大值与最小值之差至少应为0.4m。

根据上述的三个规则，不符合条件的聚类将被滤除。点云中剩下的聚类目标即为交通标志牌目标物。

1. 2 基于形状的检测

基于反射率的方法会遗漏因老化或被损坏等因素致使反射率变低的交通标志牌，以及背面相对于扫描仪的交通标志牌。因此，结合基于形状的方法，以确保检测的完整性。首先，采用基于体素的向上生长方法将地面点从原始三维场景中滤除。其次，采用欧氏距离聚类方法对非地面点进行聚类，分离出空间上相互独立的目标聚类，从这些聚类中识别出杆状物目标。最后，采用主成分分析（PCA）方法对杆状目标的非杆状部分点的维度特征进行分析，

从而检测到交通标志牌目标。 1.2.1 地面点去除

车载移动激光扫描系统采集到的被测场景的三维点云数据中包含大量、高密度的地面点。由于我们的检测目标为非地面点，这些地面点无疑会增加检测算法的时间和空间复杂度，因此对点云数据进行地面点预处理将有效提高后续的检测效率。采用基于体素向上增长的方法进行地面点滤除，步骤如下：

1)将整个三维点云场景在XY平面内按一定的宽度垂直切分成一系列局部点云块其中，Nb为生成的局部点云块的总数。这些局部点云块将被单独Blocki，i ? 1, 2,?, Nb，地进行地面点滤除处理。

2)对于每个局部点云块Blocki，i ? 1, 2,?, Nb，利用八叉树（Octree）索引结构划分，按一定的宽度将其划分成一系列空间上连续的点云体素Vj, j? 1, 2,?, Nv，其中，Nv为生成的点云体素的总数。

3)点云体素Vj, j? 1, 2,?, Nv的上9个邻域定义为与Vj相邻接并且在空间上位于Vj上方的9个点云体素。接下来，向上生长过程以点云体素Vj的上9个邻域（L1, L2, ?, L9）

为起点，按同样的方式，分别沿着L1, L2, ?, L9的上9个邻域继续向上生长。以此类推，当所有已生长到的点云体素再无上9个邻域可生长时，该递归式向上生长过程结束。从已生长的区域中计算出具有最大局部高度值的点云体素Vh，即：该生长区域的至高点。

4)定义高度Hg为地面高度临界值。如果Vh小于Hg，那么，点云体素Vh将被标记为地面体素。同时，滤除包含点云体素Vh的局部点云块中的所有点。如果点云体素Vh的高度值大于Hg，点云体素Vh将被标记为非地面体素。同时，保留包含点云体素Vh的局部点云块中的所有点。

1.2.2基于比例的杆状物检测

首先采用欧式距离聚类方法，将非地面点划分成一系列三维目标聚类。欧氏聚类的具体过程已在方法1.1.2中介绍。

对于初始分割得到的每个聚类Ci，为判断Ci是否为杆状物，对每个聚类Ci作如下处理： 1)根据z坐标的值将Ci中的点分为4等分，从下至上编号为0~3。接着，取编号为1的部分，等距离分割成N个水平截面，并将每一个点都映射到相应的平面上（如图2）。

2)这样，对于每一个平面，易得包含所有点的最小外接圆。计算其最小外接圆的半径长度，并记录它与相邻平面外接圆半径的差值rdif。规定若Ci为杆状物，相邻半径的差值应不超过dneigh。当一次迭代完成后，计算rdif超过了dneigh的比例x?Nrdif?dneigh/N，如果x超

过了预设阈值，则认为聚类Ci不是杆状物。

选择编号为1的部分作为判断杆状物的依据，既避免了草丛等物体附着在杆状物上所造成的干扰，也避免涉及到杆状物上部不同的非杆状部分，如路灯的灯头，树木的树冠等。

为了进一步对杆状物进行区分，对于杆状物聚类Ci作如下处理:

1)取Ci中编号为2、3的部分，等距离分割为水平截面，并将每一个点都映射到相应的平面上。

2)从Ci编号为2的部分开始，从下至上遍历每一个平面，拟合其平面上所映射的点的最小外接圆并计算其半径。当圆的半径长度大于某一个设定的最小值时停止。这时，在这个平面以上的部分即为杆状物的非杆状部分。