一种基于三维点云的建筑物平面分割方法技术

本发明专利技术为一种基于三维点云的建筑物平面分割方法，其将区域生长算法与基于边界点的距离算法相结合，首先，利用区域生长算法对建筑物平面进行粗分割，并从原始建筑物点云中去除粗提取的平面点得到剩余点云数据。然后，利用在切面上的点与其相邻点之间形成的向量之间夹角的最大值来提取粗平面的边界点。最后，根据剩余点到提取的粗平面边界点的距离阈值，得到两平面相交的角点，将粗提取的平面点与对应的角点相结合，得到最优的平面分割。在不引起过度分割的前提下，本发明专利技术对同一平面上的不同目标能进行精确分割。本发明专利技术还可以精确地提取出区域生长算法无法分离的两个平面相交的角点。点。点。

【技术实现步骤摘要】
一种基于三维点云的建筑物平面分割方法


[0001]本专利技术属于三维点云分割领域，涉及一种基于三维点云数据的建筑物平面分割方法。

技术介绍

[0002]随着我国城镇化的快速发展、城市人口的高速增长以及智慧城市建设的迅猛需求， 三维建筑物模型已成为“数字城市”的重要组成部分，在城市规划、环境模拟、文化遗产保护 等领域得到广泛的应用。建筑物平面包含丰富和直观的细节信息，是建筑物三维建模不可或 缺的组成部分。建筑物平面点云分割是实现数字城市三维模型重建的前提和基础，如何高效 率、高精度地从点云数据中提取建筑物平面信息，成为建筑物三维重建的重要任务之一。建 筑物平面数据的主要通过三种方法获取：航空摄影测量、深度相机和激光雷达扫描仪，这些 方法可以提供建筑物平面的稠密点云数据。通过航空摄影测量可以获得建筑物平面特征的详 细纹理图案。深度相机可以记录物体到相机的距离。激光雷达扫描仪不仅可以获取在不同位 置物体的高密度、高精度几何结构信息，还可以提供高质量点云数据的三维建筑物模型。
[0003]传统分割算法根据点云的一些基本特性如几何特征、颜色信息等设计特征描述子， 然后基于这些特征描述子建立特征滤波器进行分割。传统的方法主要包括：1)霍夫变换 (Hough Transform)：Hough变换主要用于检测具有给定特征的几何形状(如直线、平面和 圆)。2)随机采样一致性算法(RANSAC)：RANSAC算法通过估计数学模型参数来检测平 面，进行点云数据平面分割。3)区域增长算法：区域增长算法是检测平面、圆和其他参数 化形状的经典方法。该算法的基本思想是首先在分割区域中随机找出一个种子像素，将其作 为初始点，然后将邻近区域中与种子像素具有相似特征的像素聚集在初始点区域中得到一个 新像素，将得到的新像素作为新的种子像素进行生长，直到附近邻域中没有种子像素特征相 似的像素为止。虽然，区域生长算法可以分割建筑平面。但会造成两平面相交点的缺失，而 RANSAC算法会引起平面的过度分割。Hough变换算法由于每个点需要参与多次投票，计 算量大、占用内存多且计算复杂度和存储量随参数和数据量的增加呈指数增长。因此，大多 数研究将区域生长算法与其他方法相结合来检测建筑平面，虽然这些方法可以提取建筑平面 点云数据，但操作复杂。

技术实现思路

[0004]为了从三维点云数据中能精确对建筑物平面进行分割，本专利技术将区域生长算法与基 于边界点的距离算法相结合，提出了一种基于三维点云的建筑物平面分割方法。
[0005]本专利技术所采用的技术方案是：一种基于三维点云的建筑物平面分割方法，其特征在 于，包括以下步骤：
[0006]步骤1：利用区域生长算法对建筑平面进行粗分割，并从原始建筑物点中去除粗提 取平面点得到剩余点数据。
[0007]步骤2：利用在切面上的点与其相邻点之间形成的向量之间夹角的最大值来提取粗 平面的边界点。
[0008]步骤3：根据剩余点到提取的粗平面边界点的距离阈值，得到两平面相交的角点， 将粗提取的平面点与对应的角点相结合，得到最优的平面分割。
[0009]对于建筑物平面分割算法，由于Hough变换的计算量大、占用内存多，RANSAC 算法容易产生过度分割，而区域生长算法容易缺失两平面的相交点。本专利技术将区域生长算法 与基于边界点的距离算法相结合，提出一种建筑物平面点云分割算法。与RANSAC算法相 比，该方法不仅可以在不产生过度分割的情况下对同一平面上的不同目标进行分割，还可以 获取区域增长算法无法提取两个平面相交的角点。
附图说明
[0010]图1为本专利技术实施例的流程图；
[0011]图2为本专利技术实施例的边界点提取图，(a)原始数据、(b)边界点提取结果；
[0012]图3为本专利技术实施例中cottage的原始点云；
[0013]图4为本专利技术实施例中pantry的原始点云；
[0014]图5为本专利技术实施例中cottage的平面分割结果图；
[0015]图6为本专利技术实施例中pantry的平面分割结果图。
具体实施方式
[0015]下面将结合附图对本专利技术作进一步的详细描述，所描述的实例仅旨在便于对本专利技术 的理解，并不用于限定本专利技术。
[0016]请见图1，本专利技术提供了一种基于三维点云的建筑物平面分割方法，本专利技术的具体 步骤如下：
[0017]1.平面的粗提取
[0018]区域生长是一种基于空间邻近性和同质性的分割算法。本专利技术采用区域生长算法对 建筑平面进行粗提取。法向量和曲率是曲面的基本特征，是进行曲面特征识别的重要基础。 利用区域生长算法进行平面的粗提取，需要估计采样点的法向量和曲率来获得平面分割。本 专利技术利用由相邻点组成的曲面估计点的法向量和曲率，并利用主成分分析或由相邻点生成的 协方差矩阵的特征向量和特征值计算点的法向量。p
i
为平面上任意一点，则对应的协方差矩阵 C为:式中，n为p
i
的相邻点个数，p
s
为相邻点的 中心。p
i
的法向量为：其中λ
j
为协方差矩阵的第j个特征值，为第j 个特征向量。法向量是对应于最小特征值的特征向量。由于特征值λ0<λ1<λ2代表一个点随三 个主要方向的变化程度，故p
i
的曲率σ可近似为：其中λ0,λ1,λ2为协方差矩 阵的特征值。
[0019]区域生长分割算法的具体流程：1)计算所有点的曲率，并将最小曲率点添加到种子点集。 2)计算当前种子点及其相邻点的法向量。如果当前种子点的相邻点不仅满足该点的法向量与当 前种子点法向量的夹角值，而且其曲率值小于设置的曲率值，则该相邻点到该种子点序列。3) 删除当前的种子点，将相邻的点视为新的种子点，继续区域的生长。4)这
个过程将继续下去， 直到种子点序列为空为止。
[0020]2.边界点提取
[0021]粗提取平面的边界点是根据切平面上点与其相邻点所形成的向量夹角的最大值来提取 的。在得到所有点的法向量后，利用点与其邻点在切平面上形成的向量的最大夹角准则来判断其 是否位于平面的边界上。需要提取适当数量的边界点来计算剩余点到粗提取平面边界点的距离。 如果角度阈值太小，很多内部点将被视为边界点。如果角度阈值太大，可能会丢失一些边界点(如 图2)。
[0022]粗提取平面的边界点提取过程：1)计算由一个点及其邻近点构成的曲面的法向量。切平 面是由法向量计算出来的。2)将这些点投影到切平面上。3)计算这些相邻点在切平面上构成的 向量之间的夹角。4)若相邻向量之间的最大夹角大于设定的阈值，则视为边界点；否则，视为 内点。
[0023]3.最优平面获取
[0024]获得边界点后，利用剩余点到粗提取平面边界点的距离阈值提取两个平面相交的共面点。 剩余点从建筑的原始点中去除粗提取平面点后得到。剩余点用P(x,y,z)表示。粗提取平面的边 界点表示为Q(l,m,n)。剩余点到本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于三维点云数据的建筑物平面分割方法，其特征在于如下步骤：步骤1：利用区域生长算法对建筑物平面进行粗分割，并从原始建筑物点中去除粗提取的平面点得到剩余点数据。步骤2：利用在切面上的点与其相邻点之间形成的向量之间的夹角的最大值来提取粗平面的边界点。步骤3：根据剩余点到提取的粗平面边界点的距离阈值，得到两平面相交的角点，将粗提取的平面点与对应的角点相结合，得到最优的平面分割。2.根据权利要求1所述的基于三维点云数据的建筑物平面分割方法，其特征在于，步骤1进一步包括：采用区域生长算法对建筑平面进行粗提取。该算法的思想是选择曲率最小的点作为初始种子点。该种子点的相邻点满足该点的法向量与当前种子点的法向量之间的夹角值。此外，它们的曲率值小于设定的曲率值。这些被视为新的种子点，继续区域生长，直到种子点序列为空。邻接点被认为属于同一个平面。具体流程如下:(1)计算所有点的曲率，并将最小曲率点添加到种子点集。(2)计算当前种子点及其相邻点的法向量。如果当前种子点的相邻点不仅满足该点的法向量与当前种子点法向量的夹角值，而且其曲率值小于设置的曲率值，则该相邻点到该种子点序列。(3)删除当前...

【专利技术属性】
技术研发人员：周贵云，苏中花，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：