一种基于激光测量点云的室内物体提取方法与流程

本发明属于地图制图技术领域，特别涉及一种基于激光测量点云的室内物体提取方法。

背景技术：

随着经济活动的快速发展，建筑物的数量大幅增加，人类的室内活动显著增长。与此同时，建筑物内部结构也呈现大型化、复杂化趋势。因此，人们对室内环境的了解和认知需求也随之不断提升，室内场景理解与识别引发了广泛的关注。室内物体提取与识别，作为室内场景理解的基础部分与主要组成部分，对于室内环境的认知与理解上具有重要意义。

通过点云实现室内物体的分割与分类一直以来都是研究的热点，目前的方法大致分为两类：一类是基于深度学习(或机器学习)的分类分割方法，该类方法虽然可以取得良好的效果，但其针对不同的室内环境需要不同的分类样本集，为分割与分类带来众多麻烦；另一类是假设所有的物体与地面保持垂直，不存在倒放、斜放的物体，该类方法大大限制了室内环境，从而导致在真实场景中分类精度不高。

运用传统方法对室内物体直接进行分类与分割，费时费力。随激光传感器技术发展，基于激光点云进行室内物体识别具有速度快、精度高等优势，成为室内场景理解与识别的研究热点。然而，通过室内三维点云进行室内物体分割与分类还存在若干问题，如由于室内设施相互遮挡导致的墙面数据大面积缺失，由于室内的窗户、玻璃等强反射面导致的噪声，通过深度学习方法的室内物体识别需要大量样本，通过提取平面为主的物体识别无法保证物体识别精度，以及目前方法无法摆脱曼哈顿空间假设等问题。这些问题都给基于点云的室内物体分割与分类带来了困难。

技术实现要素：

为了解决上述技术问题，本发明提出了一种在多楼层室内场景中，以含噪声、遮挡的点云数据为输入数据，实现点云的分割与分类，为之后的室内场景理解提供基础。

本发明所采用的技术方案为：

一种基于激光测量点云的室内物体提取方法，包括以下步骤：

步骤1，利用区域生长算法对室内点云进行类平面分割，得到分割后的类平面点簇；

步骤2，对每一个类平面点簇做邻域分析，并连接相邻类平面点簇，生成邻接图集；

步骤3，针对每个邻接图提取图中的节点属性与边属性，节点属性包括主截面标识、类平面中心点与主截面节点描述子，边属性包括主边标识与边描述子；主截面标识为邻接图中相邻节点最多的节点，将该节点定义为主截面节点，主边标识为连接主截面中心点与其他类平面中心点的边；

步骤4，根据提取的每个邻接图的节点属性与边属性，计算任意两个邻接图的相似性，根据相似性以超图分割的方法进行图分割和图匹配；

步骤5，对步骤4处理后的邻接图，沿主截面节点的法向量方向上下延伸，做主截面包围盒，将包围盒内未分类的点云加入主截面所在类别，完成室内物体的提取。

其中，步骤4具体包括以下步骤：

步骤4.1，从邻接图集中任取两个邻接图，计算两个邻接图对应主截面节点的几何相似性，若两个邻接图主截面节点的几何相似性大于阈值，则将两个邻接图归为一类，得到粗分类集合gs＝{gs1,gs2,…gsh}，h表示粗分类的聚类数目；

步骤4.2，针对每一类gsi∈gs，任取两个邻接图，以超图分割的方法对两个邻接图进行匹配，找到两个邻接图中对称节点，并删除非对称节点。

其中，步骤4.1中节点的几何相似性计算公式如下：

式中，v^a和分别为两个邻接图对应的主截面节点，表示对应的主截面节点描述子f1，f2，f3，f4和f5的相似度大小，f1为节点的外接矩形面积l*w，f2为节点的外接矩形长宽比：w/l，f3为节点的点簇面积与节点外接矩形面积的比值，f4为节点的点簇法向量，f5为节点中心点到外接矩形的距离比该距离与外接矩形长、宽之和：d/(l+w+d)，其中d表示主截面中心点沿主截面法向量到外接矩形的距离，w表示主截面节点的外接矩形的宽度，l表示主截面节点的外接矩形的长度；为主截面节点v^a的节点描述子，为主截面节点的节点描述子。

其中，步骤4.2具体为：

针对每一类gsi∈gs，任取两个邻接图g(e,v)与gt(et,vt)，g(e,v)的节点数v少于gt(et,vt)的节点数，将两个邻接图g(e,v)与gt(et,vt)的主截面通过一个边连接，形成超图sg，从超图sg对称轴的两方分别对应任选多对节点进行节点相似度计算，找到超图sg中相似度最大的对应节点，若相似度最大值大于阈值，则该对应节点为同一类型物体，保留两个邻接图中为同一类型物体的对应节点，删除两个邻接图中的其余节点；

节点相似度score计算公式如下：

ws+wc＝1

其中，g为从超图sg对称轴的两方分别对应选取的节点对数；为对应节点的几何相似度，为对应节点的连接相似度，v和vt分别为两个邻接图中的节点，e和et为两个邻接图中的边；

连接相似度计算公式如下；

计算方法如下：

其中，c1、c2和c3为边描述子，c1为边的长度，c2为边的朝向与主截面法向量夹角，c3为边的朝向与另一连接面法向量夹角，

相对于现有技术，本发明的有益效果为：

本发明提供了一种基于激光测量点云的室内物体提取方法，该方法不需要额外的训练样本集，也可以取得良好的分割分类效果。本发明可以快速高效的实现点云的室内物体分割与分类，可以更好的处理室内复杂环境下的点云数据。

附图说明

图1为本发明实施例的方法流程图；

图2为本发明实施例的输入的原始点云；

图3为本发明实施例中室内物体邻接图示意图；

图4为本发明实施例中两物体邻接图组成超图的示意图；

图5为本发明实施例中输入点云的房间一分类结果图；

图6为本发明实施例中输入点云的房间二分类结果图。

具体实施方法

为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。

通过点云进行室内场景理解一直以来都是研究的热点与难点。在此背景下，本发明提供了一种基于室内三维点云的不使用训练样本直接进行室内物体的分割与分类的方法。

见图1，本发明实施例提供的一种基于激光测量点云的室内物体提取方法，包括以下步骤：

步骤1，利用区域生长算法对室内点云进行类平面分割，得到分割后的类平面点簇；具体为：

输入室内点云p＝{p1,p2,…pn}，输入的原始点云如图2所示，n表示输入点云的点数，对输入点云，利用区域生长算法对点云进行类平面分割，并保存分割后的类平面点簇pc＝{pc1,pc2,…pcm}，m表示分割后的点簇数量，区域生长算法见下述文件[1]使用平滑约束对点云进行分割。

[1]rabbani,t.；heuvel,f.a.v.d.；vosselman,g.segmentationofpointcloudsusingsmoothnessconstraint.int.arch.photogramm.remotesens.spat.inf.sci.2006,36,248–253

步骤2，对每一个类平面点簇做邻域分析，并连接相邻类平面点簇，生成邻接图集；

顺序选取一个类平面点簇pci，再选取另外一个类平面点簇pcj，判断两个类平面点簇的距离是否小于阈值ε，若是则连接相邻类平面点簇，否则重新选取类平面点簇进行判断，直至遍历完成所有的类平面点簇，生成邻接图集；其中，阈值ε范围在0.05-0.5m，默认值为0.2m；

步骤3，提取邻接图的节点属性与边属性，节点属性主要由主截面标识(若该节点在邻接图中，连接节点最多，则该节点为主截面，主截面一般是桌子的上表面、椅子的座面等)与节点描述子组成，边属性由主边标识(连接主截面与普通节点的边)与边描述子组成；包括以下子步骤：

步骤3.1，根据观察，一个物体中的主截面(如桌子的上表面，椅子的座面，灯管的管状面)均占据相对较大的面积，能覆盖大多物体结构，且在物体中大多结构与其相连接。因此，对每个邻接图，提取邻接图中连接节点最多的节点，并将该节点定义为主截面，与主截面相连的定义为主边，生成的室内物体邻接图如图3所示。

步骤3.2，提取邻接图的节点属性，节点属性包含三个部分：主截面标识(0为非主截面，1为主截面)，类平面中心点center以及主截面节点描述子。主截面节点描述子包含属性如表1所示。表中d表示主截面中心点沿节点的法向量到外接矩形的距离，w主截面节点的表示外接矩形的宽度，l表示主截面节点的外接矩形的长度。

表1节点描述子

步骤3.3，提取邻接图的边属性，边属性包含主边标识与边描述子。主边标识为连接主截面中心点与其他类平面中心点的边，边描述子包含属性如表2所示。

表2边描述子

步骤4，根据提取的每个邻接图的节点属性与边属性，计算任意两个邻接图的相似性，根据相似性以超图分割的方法进行图分割和图匹配；

包括以下子步骤：

步骤4.1，从邻接图集中任取两个邻接图，计算两个邻接图对应主截面节点的几何相似性，若两个邻接图主截面节点的几何相似度大于阈值γ，则将两个邻接图归为一类，得到粗分类集合gs＝{gs1,gs2,…gsh}，h表示粗分类的聚类数目；其中，阈值γ选定范围在0.1-0.3，默认值为0.1；

节点的几何相似性计算公式如下：

式中，v^a和分别为两个邻接图对应的主截面节点，表示对应的主截面节点描述子f1，f2，f3，f4和f5的相似度大小；

步骤4.2，针对每一类gsi∈gs，任取两个邻接图g(e,v)与gt(et,vt)，g(e,v)的节点数|v|少于gt(et,vt)的节点数，将两个邻接图g(e,v)与gt(et,vt)的主截面通过一个边连接，形成超图sg，如图4所示；从超图sg对称轴的左右两方分别对应任选多对节点进行节点相似度计算，找到超图sg的相似度最大值，以及使相似度最大的对应(多对)节点对，若最大相似度大于阈值w，则该对应节点为同一类型物体，保留两个邻接图中为同一类型物体的对应节点，删除两个邻接图中的其余节点；其中，阈值w选定范围在0.1-0.3，默认值为0.1；

节点相似度score计算公式如下：

ws+wc＝1

其中，g为从超图sg对称轴的左右两方分别对应选取的节点对数；为对应节点的几何相似度，为对应节点的连接相似度，v和vt分别为两个邻接图中的节点，e和et为两个邻接图中的边；

连接相似度计算公式如下；

计算方法如下：

其中，c1、c2和c3为边描述子，c1为边的长度，c2为边的朝向与主截面法向量夹角，c3为边的朝向与另一连接面法向量夹角，

步骤5，最终，对步骤4处理后的邻接图，沿主截面的法向量方向上下延伸，做主截面外接长方体，将包围盒内未分类的点云簇，加入主截面所在类别，完成室内物体分类，最终的室内物体的分类结果如图5，图6所示。

以上所述仅为本发明中的一个实施例，并不用于限制本发明。凡在本发明的精神与原则之内，所做的任何修改，改进等，均应包含在本发明的保护范围之内。