一种基于决策树的激光雷达波形信号分类方法与流程

本公开涉及图像处理技术领域，具体而言，涉及一种基于决策树的激光雷达波形信号分类方法。

背景技术：

与传统的微波雷达相比，激光雷达具有高精度、高分辨率、高探测灵敏度、高保密性、体积小，重量轻便于机载、舰载的特点。此外，由于工作机制的不同，激光雷达所发出激光脉冲相比于传统的微波雷达信号具有更强的抗干扰能力、获取的数据信息更加丰富，这导致其具有更高的探测和识别能力。尤其是新型全波形采样激光雷达的发展，使得复杂目标的探测及识别成为可能。然而，激光雷达海量回波信号的处理已经成为限制新型激光雷达技术发展的主要技术瓶颈。

在某些实际应用中，我们只对部分目标感兴趣并不需要对测量范围内的全部点云数据进行处理。以林业部门为例，其只关心测量区域的植被信息，而建筑物、公路等信息属于干扰信息，在处理数据时就可以只筛选有用信息处理，减少冗余计算，加快数据处理速度。在对数据进行筛选时，就需要对数据进行分类。

决策树是一种由多个节点构成的类似于流程图的树状结构。大多数决策树的根节点位于树的最顶层，均是采用自上向下的方法。节点包括内部节点和树叶节点，其中内部节点表示对一个属性上的测试，测试输出由分枝表示；树叶节点位于决策树的最底部，用于存放类标号。决策树这种自上而下的构造结构，使得训练集被各内部节点属性阈值的限定递归地分成更小的子集。决策树简单易用、不需要对数据有任何的先验假设、计算速度较快、结果解释性强、稳健性强，这些优势使得决策树被广泛应用于分类技术中。

当前，关于激光雷达的数据分类技术绝大部分是以激光雷达点云数据为基础，在此之上利用算法进行地形地物分类。点云数据是指激光雷达获取的激光脚点数据，其中包含了目标点的空间位置信息(x,y,z)、物体反射的强度信息、回波次数信息等。

由于激光原始信号到点云数据生成过程较为复杂，点云数据生成需要大量时间及运算资源，且分类效果严重依赖于点云数据质量。因此，现有的基于激光点云数据的分类方案存在数据处理过程繁琐，分类时间长，计算量大，难以满足实际需要等问题。

公开内容

本公开的目的在于提供一种基于决策树的激光雷达波形信号分类方法，能够解决上述提到的至少一个技术问题。具体方案如下：

根据本公开的具体实施方式，本公开提供一种基于决策树的激光雷达波形信号分类方法，包括：

通过波形采样雷达回波信号获取条纹图数据，基于所述条纹图数据提取特征值；

将地形地物类别分别编码为1、2、3、4四个类别，分别表示平原、丘陵、建筑物和树木；

基于所述特征值计算所述条纹图数据的拘束条件，其中所述拘束条件包括连通域数量、霍夫直线峰值、峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、高程差值、强度均值和面积；

根据所述地形地物的特征值以及所述拘束条件，构建决策树；

通过所述决策树对所述条纹图数据进行分析，识别出所述条纹图数据中的地形地物类别。

可选的，所述根据所述地形地物的特征值以及所述拘束条件，构建决策树，包括：

获取降噪后的所述地形地物的特征；

计算所述连通域数量，当所述连通域数量为1时，识别出所述地形地物为平原或丘陵；

当所述连通域数量大于1时，识别出所述地形地物为建筑物或树木。

可选的，所述当所述连通域数量为1时，识别出所述地形地物为平原或丘陵，包括：

当所述连通域数量为1时，获取所述霍夫直线峰值；

当所述霍夫直线峰值大于第一阈值时，识别出所述地形地物为平原；

当所述霍夫直线峰值小于所述第一阈值时，识别出所述地形地物为丘陵。

可选的，所述当所述连通域数量大于1时，识别出所述地形地物为建筑物或树木，包括：

当所述连通域数量大于1时，获取所述高程差值；

当所述高程差值大于第二阈值时，进一步确定所述地形地物为建筑物或树木；

当所述高程差值小于第二阈值时，进一步确定所述地形地物为平原或丘陵。

可选的，所述当所述高程差值大于第二阈值时，进一步确定所述地形地物为建筑物或树木，包括：

当所述高程差值大于第二阈值时，获取所述条纹图数据的峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积；

根据所述峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积，确定所述地形地物为建筑物或树木。

可选的，所述根据所述峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积，确定所述地形地物为建筑物或树木，包括：

所述强度均值大于70时，确定所述地形地物为树木，否则，确定所述地形地物为建筑物。

可选的，所述强度均值表达式如下：

其中中，i表示灰度级，i＝0,1,2……l-1；l表示灰度级的种类数目，回波信号的灰度级数目为256；ni表示具有灰度级i的总的像素个数；n表示图像中总的像素数目。

可选的，所述圆形度表达式如下：

其中，s表示目标区域的面积，l表示目标区域的周长，圆形度c越大，目标的圆形度越好。

可选的，所述矩形度的计算公式为：

其中，s表示目标区域的面积；smer表示最小外接矩形的面积，矩形度r越大，目标越接近矩形，最大值为1。

可选的，所述长宽比的计算公式为：

其中，w为短轴长度，l为长轴的长度，长宽比k值越大表明物体越细长。

本公开实施例的上述方案与现有技术相比，至少具有以下有益效果：本公开基于决策树的激光雷达波形信号分类方法，无需对条纹图像进行点云数据的转化，可以直接进行图像识别，通过分类后得到的建筑类、坡地类以及平地类目标绝大部分被正确分类标记出来，分类效果良好。而植被等其它类目标主要是针对那些回波信号缺失比较多的原始条纹回波信号进行识别，也能够准确的识别出树木或建筑，整个图像识别方法简单、高效。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1示出了根据本公开实施例的基于决策树的激光雷达波形信号分类方法总体流程图；

图2示出了根据本公开实施例的基于决策树的激光雷达波形信号分类方法流程图；

图3示出了根据本公开实施例的传统机载激光雷达探测系统工作原理图；

图4示出了根据本公开实施例的新体制激光雷达中条纹管成像机制原理图；

图5示出了根据本公开实施例的周长示意图；

图6示出了根据本公开实施例的高程差示意图；

图7示出了根据本公开实施例的选择属性进行节点分裂流程图；

图8示出了根据本公开实施例的决策树训练流程图；

图9示出了根据本公开实施例的决策树结构；

图10示出了根据本公开的实施例的电子设备连接结构示意图。

具体实施方式

为了使本公开的目的、技术方案和优点更加清楚，下面将结合附图对本公开作进一步地详细描述，显然，所描述的实施例仅仅是本公开一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本公开保护的范围。

在本公开实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本公开。在本公开实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义，“多种”一般包含至少两种。

应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。

应当理解，尽管在本公开实施例中可能采用术语第一、第二、第三等来描述……，但这些……不应限于这些术语。这些术语仅用来将……区分开。例如，在不脱离本公开实施例范围的情况下，第一……也可以被称为第二……，类似地，第二……也可以被称为第一……。

取决于语境，如在此所使用的词语“如果”、“若”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地，取决于语境，短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的商品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种商品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的商品或者装置中还存在另外的相同要素。

基于条纹管的机载激光雷达采用波形采样技术，可以根据其原始回波信号条纹图形态学特征直接分类，减少了将从所有原始回波信号转化为点云数据图这一繁琐过程。此外，基于条纹管的机载激光雷达每个回波信号的数据量可以达到单点探测的1000倍，获取的信息更加丰富。通过分析各类地形地物的属性特征，构建决策树，可以提高分类的速度和准确度。所以，本公开旨在提出一种基于激光条纹原始回波信号的能够快速准确分类的点云数据分类技术。

本公开用基于条纹管的新体制激光雷达的原始数据分类，总体流程如图1所示。

下面结合附图详细说明本公开的可选实施例。

实施例1

如图2所示，根据本公开的具体实施方式，本公开提供一种基于决策树的激光雷达波形信号分类方法，具体包括如下方法步骤：

步骤s102：通过波形采样雷达回波信号获取条纹图数据，基于所述条纹图数据提取特征值。

如图3所示，图示为传统机载激光雷达探测系统工作原理，传统的机载单点扫描激光雷达测量系统主要包括以下几个组成部分：(1)激光测距单元(runningunit)，光学机械扫描单元(opto-mechanicalscanner)，控制记录单元(control-monitoringandrecordingunits)，差分全球定位系统dgps(differentialglobepositioningsystem)和惯性测量单元imu(inertialmeasurementunit)等。当前大部分激光雷达使用脉冲激光器向目标发射光脉冲，到达目标后会有部分激光脉冲被反射回激光器，通过计算从发射到返回的时间差t即可测出激光器到目标的距离d，如公式(式中c表示光速)：

为了得到关于目标地形地物的更多特征细节，实现更加准确的分类。本公开采用基于条纹管的新体制激光雷达，条纹管的工作机理如图4所示。

激光脉冲发射到物体表面之后，返回部分光信号。这部分光信号经过狭缝整形为线形光束，再由后面的光学聚焦系统聚焦于条纹管的光电阴极上。光电阴极发生光电效应将光子转变为光电子，其瞬时发射密度与该时刻的脉冲强度成正比，所以光电阴极发射的电子与入射光信号在时空结构上等效。之后，光电子脉冲进入偏转系统，偏转系统会对不同时刻的电子进行线性地偏转，使不同时刻的电子在荧光屏上按照时间顺序展开。

步骤s104：将地形地物类别分别编码为1、2、3、4四个类别，分别表示平原、丘陵、建筑物和树木。

步骤s106：基于所述特征值计算所述条纹图数据的拘束条件，其中所述拘束条件包括连通域数量、霍夫直线峰值、峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、高程差值、强度均值和面积。

决策树是一个分类预测模型，为了构建决策树，需要对待分类的信号进行特征分析，确定各约束条件的阈值。

在本发明中，我们选取条纹图像的连通域数量、霍夫直线峰值、峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、高程差值、强度均值、面积作为地形地物的分类特征。

1、霍夫直线峰值

霍夫变换(hough)直线检测是一种图像识别和特征提取技术，它通过一种投票算法检测具有直线形状的物体。该变换通过累加器统计复合直线形状的集合，一般地，有几条直线就会出现几个峰值，峰值的大小也代表了该部分直线的形态大小。

2、周长

周长是目标区域的边界长度，如图5所示。

3、面积

计算面积常用的方法是统计目标区域内及边界像素的个数，由于二值的图像中目标值均为1，背景值均为0，所以通过对图像强度值进行累加即可求得，计算公式为：

4、强度均值表达式如下：

其中中，i表示灰度级，i＝0,1,2……l-1；l表示灰度级的种类数目，回波信号的灰度级数目为256；ni表示具有灰度级i的总的像素个数；n表示图像中总的像素数目。

5、圆形度表达式如下：

其中，s表示目标区域的面积，l表示目标区域的周长，圆形度c越大，目标的圆形度越好。

6、矩形度的计算公式为：

其中，s表示目标区域的面积；smer表示最小外接矩形的面积，矩形度r越大，目标越接近矩形，最大值为1。

7、长宽比的计算公式为：

其中，w为短轴长度，l为长轴的长度，长宽比k值越大表明物体越细长。

8、高程差值：

如图6所示，连通域1和连通域2之间的差值即为两个连通域之间的高程差值信息。

分别选取平原、坡地、建筑物、树木四种地形地物进行上述特征的分析，结果如下表。从表中可以看出，这四个类型的地形地物可以通过这些特征区分开来。

步骤s108：根据所述地形地物的特征值以及所述拘束条件，构建决策树。

在训练样本集的属性集中，计算出每个属性的信息熵，并求得该属性的信息增益率，依照信息增益率最大的标准逐层挑选出最佳属性对结点内数据进行分裂，直至所在同一节点内的样本数据属于同一类或者属性集没有新的属性可用于判断，就停止生长。其中最核心的属性选择集节点分裂流程图如图7所示。

在实际建树过程中，样本集中的属性会出现离散、连续、缺失等各种情况，针对每种情况需要具体问题具体分析。

用于创建决策树模型的程序，其中最核心的是信息熵理论算法，通过循环挑选最佳分裂属性以及属性对应的最佳分裂值最为树分支的判别条件，直至满足不能再对训练样本再划分等止决策树的生长。

实验中对已经提取得到的特征数据集中的每一特征属性进行信息增益率计算，逐次挑选出最优属性及对应最佳分裂值进行决策树的生长。整个处理流程如图8所示，其中的挑选典型目标的原始条纹回波信号与提取其特征数据已在第三章中完成。将已提取到的特征数据整形得到适用于决策树程序的输入参数，进而训练目标分类模型。遵循挑选目标条纹回波信号→提取特征数据→训练决策树模型。

其中，traindata和testdata为从挑选出来的原始条纹训练样本集和测试集提取的特征数据矩阵；a为特征属性的种类；dim为特征属性的编号；split_loc为最佳属性的分裂值。整个处理流程为从挑选好的建筑、坡地以及平地这三类目标原始条纹回波信号中提取得到目标面积、平均强度、线条数目等特征，整形得到traindata与testdata。作为输入矩阵开始训练决策树分类模型，依次计算出最佳分裂属性以及对应属性的最佳分裂值，逐次对节点进行分裂，直至满足叶结点条件停止分裂。

具体的，如图9所示，所述根据所述地形地物的特征值以及所述拘束条件，构建决策树，包括：

获取降噪后的所述地形地物的特征；

计算所述连通域数量，当所述连通域数量为1时，识别出所述地形地物为平原或丘陵；

当所述连通域数量大于1时，识别出所述地形地物为建筑物或树木。

可选的，所述当所述连通域数量为1时，识别出所述地形地物为平原或丘陵，包括：

当所述连通域数量为1时，获取所述霍夫直线峰值；

当所述霍夫直线峰值大于第一阈值时，识别出所述地形地物为平原；当所述霍夫直线峰值小于所述第一阈值时，识别出所述地形地物为丘陵。其中，所示第一阈值范围为800-1000，优选900。

可选的，所述当所述连通域数量大于1时，识别出所述地形地物为建筑物或树木，包括：

当所述连通域数量大于1时，获取所述高程差值；

当所述高程差值大于第二阈值时，进一步确定所述地形地物为建筑物或树木；当所述高程差值小于第二阈值时，进一步确定所述地形地物为平原或丘陵。其中第二阈值范围为0.5-10米。

可选的，所述当所述高程差值大于第二阈值时，进一步确定所述地形地物为建筑物或树木，包括：

当所述高程差值大于第二阈值时，获取所述条纹图数据的峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积；

根据所述峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积，确定所述地形地物为建筑物或树木。

可选的，所述根据所述峰值所占比例、峰值的θ、峰值的ρ、矩形度、圆形度、长宽比、强度均值和面积，确定所述地形地物为建筑物或树木，包括：

所述强度均值大于70时，确定所述地形地物为树木，否则，确定所述地形地物为建筑物。

步骤s110：通过所述决策树对所述条纹图数据进行分析，识别出所述条纹图数据中的地形地物类别。

本公开基于决策树的激光雷达波形信号分类方法，无需对条纹图像进行点云数据的转化，可以直接进行图像识别，通过分类后得到的建筑类、坡地类以及平地类目标绝大部分被正确分类标记出来，分类效果良好。而植被等其它类目标主要是针对那些回波信号缺失比较多的原始条纹回波信号进行识别，也能够准确的识别出树木或建筑，整个图像识别方法简单、高效。

实施例2

如图10所示，本实施例提供一种电子设备，所述电子设备，包括：至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；其中，

所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行如上实施例所述的方法步骤。

实施例3

本公开实施例提供了一种非易失性计算机存储介质，所述计算机存储介质存储有计算机可执行指令，该计算机可执行指令可执行如上实施例所述的方法步骤。

实施例4

下面参考图10，其示出了适于用来实现本公开实施例的电子设备的结构示意图。本公开实施例中的终端设备可以包括但不限于诸如移动电话、笔记本电脑、数字广播接收器、pda(个人数字助理)、pad(平板电脑)、pmp(便携式多媒体播放器)、车载终端(例如车载导航终端)等等的移动终端以及诸如数字tv、台式计算机等等的固定终端。图10示出的电子设备仅仅是一个示例，不应对本公开实施例的功能和使用范围带来任何限制。

如图10所示，电子设备可以包括处理装置(例如中央处理器、图形处理器等)1001，其可以根据存储在只读存储器(rom)1002中的程序或者从存储装置1008加载到随机访问存储器(ram)1003中的程序而执行各种适当的动作和处理。在ram1003中，还存储有电子设备操作所需的各种程序和数据。处理装置1001、rom1002以及ram1003通过总线1005彼此相连。输入/输出(i/o)接口1005也连接至总线1005。

通常，以下装置可以连接至i/o接口1005：包括例如触摸屏、触摸板、键盘、鼠标、摄像头、麦克风、加速度计、陀螺仪等的输入装置1006；包括例如液晶显示器(lcd)、扬声器、振动器等的输出装置1005；包括例如磁带、硬盘等的存储装置1008；以及通信装置1005。通信装置1005可以允许电子设备与其他设备进行无线或有线通信以交换数据。虽然图10示出了具有各种装置的电子设备，但是应理解的是，并不要求实施或具备所有示出的装置。可以替代地实施或具备更多或更少的装置。

特别地，根据本公开的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本公开的实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信装置1005从网络上被下载和安装，或者从存储装置1008被安装，或者从rom1002被安装。在该计算机程序被处理装置1001执行时，执行本公开实施例的方法中限定的上述功能。

需要说明的是，本公开上述的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑磁盘只读存储器(cd-rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本公开中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本公开中，计算机可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读信号介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：电线、光缆、rf(射频)等等，或者上述的任意合适的组合。

上述计算机可读介质可以是上述电子设备中所包含的；也可以是单独存在，而未装配入该电子设备中。

可以以一种或多种程序设计语言或其组合来编写用于执行本公开的操作的计算机程序代码，上述程序设计语言包括面向对象的程序设计语言—诸如java、smalltalk、c++，还包括常规的过程式程序设计语言—诸如“c”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(lan)或广域网(wan)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。

附图中的流程图和框图，图示了按照本公开各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

描述于本公开实施例中所涉及到的单元可以通过软件的方式实现，也可以通过硬件的方式来实现。其中，单元的名称在某种情况下并不构成对该单元本身的限定。