CN115170950A - 基于多特征约束的室外场景建筑物提取方法 - Google Patents

基于多特征约束的室外场景建筑物提取方法 Download PDF

Info

Publication number
CN115170950A
CN115170950A CN202210717889.XA CN202210717889A CN115170950A CN 115170950 A CN115170950 A CN 115170950A CN 202210717889 A CN202210717889 A CN 202210717889A CN 115170950 A CN115170950 A CN 115170950A
Authority
CN
China
Prior art keywords
point
points
point cloud
building
ground
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210717889.XA
Other languages
English (en)
Inventor
宁小娟
孙泽乾
王兰兰
王琳
王映辉
金海燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to CN202210717889.XA priority Critical patent/CN115170950A/zh
Publication of CN115170950A publication Critical patent/CN115170950A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/176—Urban or other man-made structures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/273—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion removing elements interfering with the pattern to be recognised
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/36—Applying a local operator, i.e. means to operate on image points situated in the vicinity of a given point; Non-linear local filtering operations, e.g. median filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/762—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using clustering, e.g. of similar faces in social networks
    • G06V10/7625—Hierarchical techniques, i.e. dividing or merging patterns to obtain a tree-like representation; Dendograms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/762—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using clustering, e.g. of similar faces in social networks
    • G06V10/763—Non-hierarchical techniques, e.g. based on statistics of modelling distributions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/60—Type of objects
    • G06V20/64—Three-dimensional objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Evolutionary Computation (AREA)
  • Databases & Information Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Probability & Statistics with Applications (AREA)
  • Nonlinear Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了基于多特征约束的室外场景建筑物提取方法，首先使用体素网格滤波方法对大规模点云场景进行下采样，对于下采样后的场景用布料滤波算法去除地面点。然后，结合基于密度聚类算法和KD树对非地面点进行聚类，将非地面点聚为不同的点簇。最后，结合点云高度、体积、方向、维度以及颜色特征将建筑物点云与非建筑物点云进行分离，最终从室外场景中提取出建筑物点云。本发明有效提高了使用传统方法提取精度低的问题，对室外场景中大部分建筑物都能成功提取，为大规模场景中的建筑物提取提供了相应的处理思路。

Description

基于多特征约束的室外场景建筑物提取方法
技术领域
本发明属于点云场景分割技术领域，涉及基于多特征约束的室外场景建筑物提取方法。
背景技术
数字化城市的最终目标是将提取到的建筑物信息转化并描述成计算机能够理解的数字化模型。从点云三维重建的规模来看，当前对于小规模室内场景的重建和小型物体的重构已经取得了很多研究成果，但对于大规模室外场景中的物体重建仍然存在许多亟需解决的问题。
建筑物提取作为建筑物重建的前提和基础，是指将给定的数据集划分为不重叠的同质区域，并从这些区域中识别出建筑物的任务，使得建筑物在重建过程中不会受到其它因素的影响。建筑物提取需要充分考虑建筑物对象区别于其它类别物体的特征，从而将建筑物与非建筑物进行分离。长期以来，大规模复杂室外场景中包含着大量不同类别的地物，如植被、车辆、与建筑物高度类似的树木、在空间距离上非常接近建筑物甚至与建筑物连在一起的树木等，导致建筑物的自动提取一直很难实现，传统的提取方法存在提取精度差和效率低的问题，这将进一步影响建筑物重建模型的质量。
发明内容
本发明的目的是提供基于多特征约束的室外场景建筑物提取方法，有效地提高了传统方法存在的提取精度低的问题。
本发明所采用的技术方案是，基于多特征约束的室外场景建筑物提取方法，具体按照如下步骤实施：
步骤1，数据预处理：输入原始点云数据，使用体素网格滤波方法对原始场景进行下采样得到预处理后的点云数据；
步骤2，地面点滤除：使用布料滤波算法将场景中的地面点进行滤除；
步骤3，非地面点聚类：构建KD(K-Dimensional)树，使用索引方法查找邻域点，通过基于密度聚类算法实现非地面点聚类；
步骤4，建筑物点提取：从点云的高度特征、点云体积特征、方向特征、维度特征和颜色特征五个方面考虑，区分建筑物和非建筑物，完成建筑物点提取。
本发明的特点还在于，
步骤1具体为：
步骤1.1，计算每个体素网格的中心点
在原始点云数据上绘制体素网格，假设当前输入点为pi(xi,yi,zi)，体素网格边长Lgrid，则体素网格的中心点坐标W(Xcen,Ycen,Zcen)为，
Figure BDA0003704655240000021
根据体素网格的边长和体素网格的中心点坐标得到每个体素网格在空间中的位置；
步骤1.2，计算每个体素网格的重心点
假设网格的重心点为G(X,Y,Z)，则对于含有n个点的体素网格，每个体素网格的重心点表示为：
Figure BDA0003704655240000031
步骤1.3，使每个体素网格的重心点替换对应体素网格内的所有点，过滤掉重心点以外的剩余点，得到预处理后的点云数据。
步骤2具体为：
步骤2.1，将预处理后的点云数据进行反转，设置初始布料网格边长，构造初始布料格网，假设粒子只在垂直方向上移动；
步骤2.2，将格网粒子与点云投影到同一水平面，找到格网中每个粒子的邻近点并记录对应点的高程值HIV；
步骤2.3，对每一个可移动的粒子，计算其位置受到重力影响后产生的移动距离，计算公式如式(3)；
Figure BDA0003704655240000032
式(3)中，X(t)表示粒子在t时刻的位置，Δt表示时间步长，G表示万有引力常数，m表示粒子质量；
步骤2.4，将移动后粒子所在的位置与其对应近邻点的高程值HIV进行比较：如果粒子的所在高度小于等于邻近点的高程值HIV，则把粒子移动到HIV高度处并设置为不可移动点；
步骤2.5，对于每个粒子，计算受到内部驱动因素影响所产生的移动距离；如果两个相连接的粒子都为可移动粒子，则将两粒子移动到两粒子的平均高程值处；若一个为不可移动点，另一个为可移动粒子，则移动可移动粒子；若两个粒子都为不可移动点，则不移动；
粒子移动产生的距离由式(4)计算得到；
Figure BDA0003704655240000041
式(4)中，
Figure BDA0003704655240000042
表示粒子的位移向量，b表示为可移动粒子时，该值为1否则为0，
Figure BDA0003704655240000043
表示可移动粒子的位置，
Figure BDA0003704655240000044
表示当前粒子邻近粒子的位置，
Figure BDA0003704655240000045
表示垂直方向上的单位向量；
步骤2.6，迭代步骤2.3至2.5，直到所有粒子的最大高程变化小于阈值Mhv或者迭代次数到达预设的值，则停止模拟过程；
步骤2.7：计算点云与粒子之间的高度差异，如果点云中的点与粒子之间的距离小于高度阈值，则认为点云中的点是地面点，反之则认为是非地面点。
步骤3具体为，
步骤3.1，构建KD树；
步骤3.2，从点云中选取一个未被处理的点p0，未被处理的点指未被标记已处理的点，根据索引查找p0ε邻域内的点，如果邻域内至少包含Tmin个点，证明p0为核心点，此时，将该点标记为已处理并创建新簇C，将p0的ε邻域点加入到C中，否则将该点标记为已处理并重新选取新的点，直到选取到核心点；
步骤3.3，遍历簇C中未被处理的点，若遍历过程中存在点在给定的阈值下也是核心点，则将该点的ε邻域内的所有点加入到簇C中，并将该点标记为已处理；
步骤3.4，重复步骤3.3直到没有新的点加入到簇C；
步骤3.5，重复步骤3.2至3.4，直至点云中的点都被处理完毕，即可完成聚类。
步骤4具体为：
步骤4.1，使用高程值剔除场景中的低矮植被：高程值小于2m的点簇为低矮地物，需将其剔除，而高于2m的点簇为建筑物和树木等较高地物，将这些点簇保留进行下一步计算，高程值计算如式(5)；
H＝Zmax-Zground (5)
式(5)中，H表示高程值，Zmax表示点云在垂直于地面方向上的最大值，Zground表示点云地面点Z坐标，用地面点Z坐标的平均值代替；
步骤4.2，使用最小凸包体积剔除场景中面积较小且高度高于低矮植被但不高于建筑物的地物：求解聚类后点簇的最小凸包，进一步计算最小凸包的体积来区分建筑物与其他地物，体积阈值Vb根据不同场景中的建筑物体积而定；
步骤4.3，使用方向特征剔除场景中的高大树木：构造并分析点云法向量与Z轴方向向量
Figure BDA0003704655240000051
夹角的余弦值来确定当前点簇所属类别；
步骤4.4，使用维度特征剔除场景中表现为线状结构和球面结构的地物，根据主成分分析方法对点云进行分析；
步骤4.5，使用颜色特征剔除场景中密度较稀疏的树木以及和建筑物连接在一起的树木。
步骤4.3具体为：
步骤4.3.1，首先将点云法向量与Z轴方向向量夹角的余弦值划分为多个区间，区间大小在[0,1]范围内，设置组距μ＝0.1，然后构造夹角余弦值的分布直方图；
步骤4.3.2，计算点云法向量与Z轴方向向量夹角的余弦值分布中每个区间的概率以及所有概率的标准差，统计概率大于标准差的区间个数Pnum和概率小于标准差的区间个数Nnum，如果Pnum>Nnum，证明该点簇余弦值分布较为离散为非建筑物，将其剔除。
步骤4.4具体为：
步骤4.4.1，计算每个点簇的协方差矩阵的特征值γ1、γ2、γ3(γ1≥γ2≥γ3)；
步骤4.4.2，令
Figure BDA0003704655240000061
使用式(6)计算空间中的点分别属于三类维度的概率；
Figure BDA0003704655240000062
式(6)中，
Figure BDA0003704655240000063
表示线性度，即点云属于一维特征的概率，
Figure BDA0003704655240000064
表示平面度，即点云属于二维特征的概率，
Figure BDA0003704655240000065
表示散乱度，即点云属于三维特征的概率；
步骤4.4.3，根据不同的概率值大小获得当前点的维度特征，进一步得到该点所属的类别，如果一个点簇中面状点占比最高，则认定属于建筑物，反之如果散乱点或者线状点占比高于面状点，则认定属于非建筑物，由此提取出场景中的所有建筑物点云。
步骤4.5具体为：
步骤4.5.1，将过绿指数EXG(Excess Green)引入三维点云场景中作为剔除稀疏树木的依据，其计算公式如式(7)；
EXG＝2g-r-b(7)
式(7)中，r＝R/(R+G+B)、g＝G/(R+G+B)、b＝B/(R+G+B)，R、G、B分别表示点云的颜色值；
步骤4.5.2，如果EXG指数大于阈值E，则判断该点r邻域内所有点的EXG指数并统计其个数n，如果n大于点数阈值N，则该点为树木点，否则为建筑物点，将阈值E设置为0.08，阈值N自适应的设置为邻域内点云个数的一半，保留建筑物点完成建筑物点提取。
本发明的有益效果是：
本发明基于多特征约束的室外场景建筑物提取方法，将下采样后的场景通过布料滤波算法将地面点滤除，其适用性更强且能适用于多种具有复杂特征的地形，在对非地面点进行聚类后选取有较强区分能力的特征进行建筑物提取，可以有效剔除低矮且面积较小的地物，对密度稀疏以及和建筑物相连接的树木也有很好的剔除效果，有效提高了使用传统方法提取精度低的问题，并且对室外场景中大部分建筑物都能成功提取，为大规模场景中的建筑物提取提供了相应的处理思路。
附图说明
图1是基于多特征约束的室外场景建筑物提取方法框图；
图2是基于多特征约束的室外场景建筑物提取方法通过布料滤波算法滤除地面点的结果；
图3是基于多特征约束的室外场景建筑物提取方法结合密度聚类算法和KD树进行非地面点聚类的结果；
图4是基于多特征约束的室外场景建筑物提取方法最小凸包构建的结果；
图5是基于多特征约束的室外场景建筑物提取方法是不同点簇夹角余弦值分布直方图；
图6是基于多特征约束的室外场景建筑物提取方法不同点簇中的点类别分布情况直方图；
图7是基于多特征约束的室外场景建筑物提取方法对建筑物点提取的结果。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
实施例
本实施例提供基于多特征约束的室外场景建筑物提取方法，如图1所示，具体按照如下步骤实施：
步骤1，数据预处理：输入原始点云数据，使用体素网格滤波方法对原始场景进行下采样得到预处理后的点云数据；
步骤1.1，计算每个体素网格的中心点
在原始点云数据上绘制体素网格，假设当前输入点为pi(xi,yi,zi)，体素网格边长Lgrid，则体素网格的中心点坐标W(Xcen,Ycen,Zcen)为，
Figure BDA0003704655240000081
根据体素网格的边长和体素网格的中心点坐标得到每个体素网格在空间中的位置；
步骤1.2，计算每个体素网格的重心点
假设网格的重心点为G(X,Y,Z)，则对于含有n个点的体素网格，每个体素网格的重心点表示为：
Figure BDA0003704655240000082
步骤1.3，使每个体素网格的重心点替换对应体素网格内的所有点，过滤掉重心点以外的剩余点，得到预处理后的点云数据。
步骤2，地面点滤除：使用布料滤波算法将场景中的地面点进行滤除；
步骤2.1，将预处理后的点云数据进行反转，设置初始布料网格边长，构造初始布料格网，假设粒子只在垂直方向上移动；
步骤2.2，将格网粒子与点云投影到同一水平面，找到格网中每个粒子的邻近点并记录对应点的高程值HIV；
步骤2.3，对每一个可移动的粒子，计算其位置受到重力影响后产生的移动距离，计算公式如式(3)；
Figure BDA0003704655240000091
式(3)中，X(t)表示粒子在t时刻的位置，Δt表示时间步长，G表示万有引力常数，m表示粒子质量；
步骤2.4，将移动后粒子所在的位置与其对应近邻点的高程值HIV进行比较：如果粒子的所在高度小于等于邻近点的高程值HIV，则把粒子移动到HIV高度处并设置为不可移动点；
步骤2.5，对于每个粒子，计算受到内部驱动因素影响所产生的移动距离；如果两个相连接的粒子都为可移动粒子，则将两粒子移动到两粒子的平均高程值处；若一个为不可移动点，另一个为可移动粒子，则移动可移动粒子；若两个粒子都为不可移动点，则不移动；
粒子移动产生的距离由式(4)计算得到；
Figure BDA0003704655240000092
式(4)中，
Figure BDA0003704655240000093
表示粒子的位移向量，b表示为可移动粒子时，该值为1否则为0，
Figure BDA0003704655240000094
表示可移动粒子的位置，
Figure BDA0003704655240000095
表示当前粒子邻近粒子的位置，
Figure BDA0003704655240000096
表示垂直方向上的单位向量；
步骤2.6，迭代步骤2.3至2.5，直到所有粒子的最大高程变化小于阈值Mhv或者迭代次数到达预设的值，则停止模拟过程；
步骤2.7：计算点云与粒子之间的高度差异，如果点云中的点与粒子之间的距离小于高度阈值，则认为点云中的点是地面点，反之则认为是非地面点，布料滤波的结果如图2所示。
步骤3，非地面点聚类：构建KD(K-Dimensional)树，使用索引方法查找邻域点，通过基于密度聚类算法实现非地面点聚类；
步骤3.1，构建KD树；
步骤3.2，从点云中选取一个未被处理的点p0，未被处理的点指未被标记已处理的点，根据索引查找p0ε邻域内的点，如果邻域内至少包含Tmin个点，证明p0为核心点，此时，将该点标记为已处理并创建新簇C，将p0的ε邻域点加入到C中，否则将该点标记为已处理并重新选取新的点，直到选取到核心点；
步骤3.3，遍历簇C中未被处理的点，若遍历过程中存在点在给定的阈值下也是核心点，则将该点的ε邻域内的所有点加入到簇C中，并将该点标记为已处理；
步骤3.4，重复步骤3.3直到没有新的点加入到簇C；
步骤3.5，重复步骤3.2至3.4，直至点云中的点都被处理完毕，即可完成聚类，聚类的结果如图3所示。
步骤4，建筑物点提取：从点云的高度特征、点云体积特征、方向特征、维度特征和颜色特征五个方面考虑，区分建筑物和非建筑物，完成建筑物点提取；
步骤4.1，使用高程值剔除场景中的低矮植被：高程值小于2m的点簇为低矮地物，需将其剔除，而高于2m的点簇为建筑物和树木等较高地物，将这些点簇保留进行下一步计算，高程值计算如式(5)；
H＝Zmax-Zground (5)
式(5)中，H表示高程值，Zmax表示点云在垂直于地面方向上的最大值，Zground表示点云地面点Z坐标，用地面点Z坐标的平均值代替；
步骤4.2，使用最小凸包体积剔除场景中面积较小且高度高于低矮植被但不高于建筑物的地物：求解聚类后点簇的最小凸包，最小凸包的构建结果如图4所示，进一步计算最小凸包的体积来区分建筑物与其他地物，体积阈值Vb根据不同场景中的建筑物体积而定；
步骤4.3，使用方向特征剔除场景中的高大树木：构造并分析点云法向量与Z轴方向向量
Figure BDA0003704655240000111
夹角的余弦值来确定当前点簇所属类别；
步骤4.3具体为：
步骤4.3.1，首先将点云法向量与Z轴方向向量夹角的余弦值划分为多个区间，区间大小在[0,1]范围内，设置组距μ＝0.1，然后构造夹角余弦值的分布直方图如图5所示；
步骤4.3.2，计算点云法向量与Z轴方向向量夹角的余弦值分布中每个区间的概率以及所有概率的标准差，统计概率大于标准差的区间个数Pnum和概率小于标准差的区间个数Nnum，如果Pnum>Nnum，证明该点簇余弦值分布较为离散为非建筑物，将其剔除。
步骤4.4，使用维度特征剔除场景中表现为线状结构和球面结构的地物，根据主成分分析方法对点云进行分析；
步骤4.4具体为：
步骤4.4.1，计算每个点簇的协方差矩阵的特征值γ1、γ2、γ3(γ1≥γ2≥γ3)；
步骤4.4.2，令
Figure BDA0003704655240000121
使用式(6)计算空间中的点分别属于三类维度的概率；
Figure BDA0003704655240000122
式(6)中，
Figure BDA0003704655240000123
表示线性度，即点云属于一维特征的概率，
Figure BDA0003704655240000124
表示平面度，即点云属于二维特征的概率，
Figure BDA0003704655240000125
表示散乱度，即点云属于三维特征的概率；
步骤4.4.3，根据不同的概率值大小获得当前点的维度特征，进一步得到该点所属的类别，不同点簇中的点类别分布情况如图6所示，点簇1和点簇2中面状点占比最高，而其它三个点簇中的散乱点占比较高，因此认为点簇1和点簇2属于建筑物，而点簇3、点簇4和点簇5为非建筑物，由此提取出场景中的所有建筑物点云。
步骤4.5，使用颜色特征剔除场景中密度较稀疏的树木以及和建筑物连接在一起的树木；
步骤4.5具体为：
步骤4.5.1，将过绿指数EXG(Excess Green)引入三维点云场景中作为剔除稀疏树木的依据，其计算公式如式(7)；
EXG＝2g-r-b(7)
式(7)中，r＝R/(R+G+B)、g＝G/(R+G+B)、b＝B/(R+G+B)，R、G、B分别表示点云的颜色值；
步骤4.5.2，如果EXG指数大于阈值E，则判断该点r邻域内所有点的EXG指数并统计其个数n，如果n大于点数阈值N，则该点为树木点，否则为建筑物点，将阈值E设置为0.08，阈值N自适应的设置为邻域内点云个数的一半，保留建筑物点完成建筑物点提取。
通过上述方式可知，本发明基于多特征约束的室外场景建筑物提取方法，首先使用体素网格滤波方法对大规模点云场景进行下采样，对于下采样后的场景用布料滤波算法去除地面点。然后，结合基于密度聚类算法和KD树对非地面点进行聚类，将非地面点聚为不同的点簇。最后，结合点云高度、体积、方向、维度以及颜色特征将建筑物点云与非建筑物点云进行分离，最终从室外场景中提取出建筑物点云。本发明有效提高了使用传统方法提取精度低的问题，技术路线完整，并且对室外场景中大部分建筑物都能成功提取，为大规模场景中的建筑物提取提供了相应的处理思路。

Claims (8)

1.基于多特征约束的室外场景建筑物提取方法，其特征在于，具体按照如下步骤实施：
步骤1，数据预处理：输入原始点云数据，使用体素网格滤波方法对原始场景进行下采样得到预处理后的点云数据；
步骤2，地面点滤除：使用布料滤波算法将场景中的地面点进行滤除；
步骤3，非地面点聚类：构建KD树，使用索引方法查找邻域点，通过基于密度聚类算法实现非地面点聚类；
步骤4，建筑物点提取：从点云的高度特征、点云体积特征、方向特征、维度特征和颜色特征五个方面考虑，区分建筑物和非建筑物，完成建筑物点提取。
2.根据权利要求1所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤1具体为：
步骤1.1，计算每个体素网格的中心点
在原始点云数据上绘制体素网格，假设当前输入点为pi(xi,yi,zi)，体素网格边长Lgrid，则体素网格的中心点坐标W(Xcen,Ycen,Zcen)为，
Figure FDA0003704655230000011
根据体素网格的边长和体素网格的中心点坐标得到每个体素网格在空间中的位置；
步骤1.2，计算每个体素网格的重心点
假设网格的重心点为G(X,Y,Z)，则对于含有n个点的体素网格，每个体素网格的重心点表示为：
Figure FDA0003704655230000021
步骤1.3，使每个体素网格的重心点替换对应体素网格内的所有点，过滤掉重心点以外的剩余点，得到预处理后的点云数据。
3.根据权利要求2所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤2具体为：
步骤2.1，将预处理后的点云数据进行反转，设置初始布料网格边长，构造初始布料格网，假设粒子只在垂直方向上移动；
步骤2.2，将格网粒子与点云投影到同一水平面，找到格网中每个粒子的邻近点并记录对应点的高程值HIV；
步骤2.3，对每一个可移动的粒子，计算其位置受到重力影响后产生的移动距离，计算公式如式(3)；
Figure FDA0003704655230000022
式(3)中，X(t)表示粒子在t时刻的位置，Δt表示时间步长，G表示万有引力常数，m表示粒子质量；
步骤2.4，将移动后粒子所在的位置与其对应近邻点的高程值HIV进行比较：如果粒子的所在高度小于等于邻近点的高程值HIV，则把粒子移动到HIV高度处并设置为不可移动点；
步骤2.5，对于每个粒子，计算受到内部驱动因素影响所产生的移动距离；如果两个相连接的粒子都为可移动粒子，则将两粒子移动到两粒子的平均高程值处；若一个为不可移动点，另一个为可移动粒子，则移动可移动粒子；若两个粒子都为不可移动点，则不移动；
粒子移动产生的距离由式(4)计算得到；
Figure FDA0003704655230000031
式(4)中，
Figure FDA0003704655230000032
表示粒子的位移向量，b表示为可移动粒子时，该值为1否则为0，
Figure FDA0003704655230000033
表示可移动粒子的位置，
Figure FDA0003704655230000034
表示当前粒子邻近粒子的位置，
Figure FDA0003704655230000035
表示垂直方向上的单位向量；
步骤2.6，迭代步骤2.3至2.5，直到所有粒子的最大高程变化小于阈值Mhv或者迭代次数到达预设的值，则停止模拟过程；
步骤2.7：计算点云与粒子之间的高度差异，如果点云中的点与粒子之间的距离小于高度阈值，则认为点云中的点是地面点，反之则认为是非地面点。
4.根据权利要求3所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤3具体为，
步骤3.1，构建KD树；
步骤3.2，从点云中选取一个未被处理的点p0，未被处理的点指未被标记已处理的点，根据索引查找p0ε邻域内的点，如果邻域内至少包含Tmin个点，证明p0为核心点，此时，将该点标记为已处理并创建新簇C，将p0的ε邻域点加入到C中，否则将该点标记为已处理并重新选取新的点，直到选取到核心点；
步骤3.3，遍历簇C中未被处理的点，若遍历过程中存在点在给定的阈值下也是核心点，则将该点的ε邻域内的所有点加入到簇C中，并将该点标记为已处理；
步骤3.4，重复步骤3.3直到没有新的点加入到簇C；
步骤3.5，重复步骤3.2至3.4，直至点云中的点都被处理完毕，即可完成聚类。
5.根据权利要求4所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤4具体为：
步骤4.1，使用高程值剔除场景中的低矮植被：高程值小于2m的点簇为低矮地物，需将其剔除，而高于2m的点簇为较高地物，较高地物包括建筑物和树木，将这些较高地物的点簇保留进行下一步计算，高程值计算如式(5)；
H＝Zmax-Zground (5)
式(5)中，H表示高程值，Zmax表示点云在垂直于地面方向上的最大值，Zground表示点云地面点Z坐标，用地面点Z坐标的平均值代替；
步骤4.2，使用最小凸包体积剔除场景中面积较小且高度高于低矮植被但不高于建筑物的地物：求解聚类后点簇的最小凸包，进一步计算最小凸包的体积来区分建筑物与其他地物，体积阈值Vb根据不同场景中的建筑物体积而定；
步骤4.3，使用方向特征剔除场景中的高大树木：构造并分析点云法向量与Z轴方向向量
Figure FDA0003704655230000041
夹角的余弦值来确定当前点簇所属类别；
步骤4.4，使用维度特征剔除场景中表现为线状结构和球面结构的地物，根据主成分分析方法对点云进行分析；
步骤4.5，使用颜色特征剔除场景中密度较稀疏的树木以及和建筑物连接在一起的树木。
6.根据权利要求5所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤4.3具体为：
步骤4.3.1，首先将点云法向量与Z轴方向向量夹角的余弦值划分为多个区间，区间大小在[0,1]范围内，设置组距μ＝0.1，然后构造夹角余弦值的分布直方图；
步骤4.3.2，计算点云法向量与Z轴方向向量夹角的余弦值分布中每个区间的概率以及所有概率的标准差，统计概率大于标准差的区间个数Pnum和概率小于标准差的区间个数Nnum，如果Pnum>Nnum，证明该点簇余弦值分布较为离散为非建筑物，将其剔除。
7.根据权利要求6所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤4.4具体为：
步骤4.4.1，计算每个点簇的协方差矩阵的特征值γ1、γ2、γ3，γ1≥γ2≥γ3；
步骤4.4.2，令
Figure FDA0003704655230000051
使用式(6)计算空间中的点分别属于三类维度的概率；
Figure FDA0003704655230000052
式(6)中，
Figure FDA0003704655230000053
表示线性度，即点云属于一维特征的概率，
Figure FDA0003704655230000054
表示平面度，即点云属于二维特征的概率，
Figure FDA0003704655230000055
表示散乱度，即点云属于三维特征的概率；
步骤4.4.3，根据不同的概率值大小获得当前点的维度特征，进一步得到该点所属的类别，如果一个点簇中面状点占比最高，则认定属于建筑物，反之如果散乱点或者线状点占比高于面状点，则认定属于非建筑物，由此提取出场景中的所有建筑物点云。
8.根据权利要求7所述的基于多特征约束的室外场景建筑物提取方法，其特征在于，所述步骤4.5具体为：
步骤4.5.1，将过绿指数EXG引入三维点云场景中作为剔除稀疏树木的依据，其计算公式如式(7)；
EXG＝2g-r-b (7)
式(7)中，r＝R/(R+G+B)、g＝G/(R+G+B)、b＝B/(R+G+B)，R、G、B分别表示点云的颜色值；
步骤4.5.2，如果EXG指数大于阈值E，则判断该点r邻域内所有点的EXG指数并统计其个数n，如果n大于点数阈值N，则该点为树木点，否则为建筑物点，将阈值E设置为0.08，阈值N自适应的设置为邻域内点云个数的一半，保留建筑物点完成建筑物点提取。
CN202210717889.XA 2022-06-21 2022-06-21 基于多特征约束的室外场景建筑物提取方法 Pending CN115170950A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210717889.XA CN115170950A (zh) 2022-06-21 2022-06-21 基于多特征约束的室外场景建筑物提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210717889.XA CN115170950A (zh) 2022-06-21 2022-06-21 基于多特征约束的室外场景建筑物提取方法

Publications (1)

Publication Number Publication Date
CN115170950A true CN115170950A (zh) 2022-10-11

Family

ID=83487118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210717889.XA Pending CN115170950A (zh) 2022-06-21 2022-06-21 基于多特征约束的室外场景建筑物提取方法

Country Status (1)

Country Link
CN (1) CN115170950A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116071661A (zh) * 2023-04-06 2023-05-05 厦门大学 一种基于激光点云的城市道路场景语义分割方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116071661A (zh) * 2023-04-06 2023-05-05 厦门大学 一种基于激光点云的城市道路场景语义分割方法

Similar Documents

Publication Publication Date Title
CN107292276B (zh) 一种车载点云聚类方法及系统
CN112070769B (zh) 一种基于dbscan的分层点云分割方法
CN106022381B (zh) 基于车载激光扫描点云的路灯杆自动提取方法
CN110992341A (zh) 一种基于分割的机载LiDAR点云建筑物提取方法
CN111191628B (zh) 基于决策树与特征优化的遥感影像震害建筑物识别方法
CN112132840B (zh) 一种车载行道树点云分类与特征信息提取方法
CN112529044B (zh) 基于车载LiDAR的铁路接触网提取分类的方法
CN115205690B (zh) 基于mls点云数据的行道树单体化提取方法及装置
CN114926699A (zh) 一种室内三维点云语义分类方法、装置、介质及终端
CN114241217B (zh) 一种基于圆柱特征的树干点云高效提取方法
CN113255677B (zh) 一种岩体结构面及产状信息快速提取方法、设备及介质
CN113920360A (zh) 一种道路点云杆状物提取与多尺度识别方法
CN108765478A (zh) 一种建筑点云中分离单体建筑的密度聚类算法
CN104268600A (zh) 一种基于Minkowski距离的矿物浮选泡沫图像纹理分析及工况识别方法
CN110348478B (zh) 一种基于形状分类与组合的室外点云场景中树木提取方法
CN113657216A (zh) 基于形状特征的点云场景树木树冠及木质点分离方法
CN116258857A (zh) 一种面向室外树木激光点云分割与提取方法
CN114608476B (zh) 一种复杂岩体三维点云结构面智能分析提取方法
CN115170950A (zh) 基于多特征约束的室外场景建筑物提取方法
CN114463338B (zh) 一种基于图割与后处理的建筑物激光脚点自动提取方法
CN112200248B (zh) 一种基于dbscan聚类的城市道路环境下的点云语义分割方法、系统及存储介质
CN114170149A (zh) 一种基于激光点云的道路几何信息提取方法
CN110458111B (zh) 基于LightGBM的车载激光点云电力线的快速提取方法
CN116824379A (zh) 一种基于多维特征的激光点云建筑物轮廓渐进优化方法
CN116385659A (zh) 一种点云建筑物建模方法、系统、存储介质及电子设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination