CN110009726B - 一种根据平面基元间的结构关系对点云提取平面的方法 - Google Patents
一种根据平面基元间的结构关系对点云提取平面的方法 Download PDFInfo
- Publication number
- CN110009726B CN110009726B CN201910175983.5A CN201910175983A CN110009726B CN 110009726 B CN110009726 B CN 110009726B CN 201910175983 A CN201910175983 A CN 201910175983A CN 110009726 B CN110009726 B CN 110009726B
- Authority
- CN
- China
- Prior art keywords
- plane
- planes
- point cloud
- psi
- extracting
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种根据平面基元间的结构关系对点云提取平面的方法,该方法首先根据邻域协方差矩阵求出各点的法向量与曲率;其次采用区域生长法从点云中提取平面基元;然后从上述提取的平面中抽象出一系列通用的结构知识,包括共面、平行、正交和对称;最后检测平面,如果发现两个平面有正交或对称的关系,则可被认定为正确的平面,并以这些正确的面片为种子面片合并其他的面片。本发明利用平面间提取的结构知识,大大提高了密集匹配点云平面提取的可靠性。
Description
技术领域
本发明属于无人机倾斜影像建模的技术领域,具体涉及一种根据平面基元间的结构关系对点云提取平面的方法。
背景技术
无人机倾斜摄影测量技术是近几年新兴的一种高测量技术,通过多台传感器从不同的角度进行数据的采集,快速高效获取丰富的数据信息,真实地反映地面的客观情况,将用户引入符合人眼视觉的真实直观世界。倾斜摄影测量具有周期短、效率高,成本低、操作简单,机动性、灵活性,可形成更加直观的3D模型等其他测量方式所不具备的优点。
随着无人机倾斜摄影的广泛应用,如何通过无人机获取的点云数据自动提取精细化的面片模型并重构其相互之间复杂的拓扑关系一直是倾斜摄影测量建筑物三维自动重建的挑战性问题之一。建筑物三维模型种类多样形态各异、城市环境错综复杂,受点云噪声、缺漏等因素影响,平面基元提取精度较差,易产生缺漏和拓扑结构错误。模型驱动方法受上述因素影响较小,因为该方法将平面提取隐式的转换成平面拟合,拟合精度仅受数据精度和模型种类影响,而不会导致拓扑结构错误。然而数据驱动方法受限于点云数据精度和分布密度,会导致空间上不连续的共面平面结合成一个面,难以精确提取平面集合。
发明内容
为了解决上述问题,本发明的目的在于提供一种根据平面基元间的结构关系对点云提取平面的方法,基于平面间提取的结构知识进行平面配准,大大提高了密集匹配点云平面提取的可靠性。
为了实现上述的目的,本发明采用了以下的技术方案:
步骤一,计算点云对应的法向量与曲率:
对于点云P,根据领域协方差矩阵求出各点的法向量与曲率;
步骤二,从点云中提取平面基元:
对于步骤一的处理结果,通过区域生长法提取平面基元;
步骤三,提取上述平面基元间的结构关系,如果两个平面有正交或对称的关系,则可被认定为正确的平面,并以这些正确的平面为种子平面合并其他的平面:
根据平面法向量间关系从基元平面集合ψ中抽象出平面间的结构关系,包括共面、平行、正交和对称,得到平面合集A和平面合集B,其中平面合集A中的平面之间为正交或对称关系,平面B中的合集为平行或共面关系;
按区域生长法,以A中的平面为种子平面,将B中的平面合并到A中。
作为优选方案,所述步骤二的具体步骤包括:
点云P中的一个点的K领域为其中K为点的K个邻近点;根据K邻域计算协方差矩阵C,该矩阵有三个特征值λ0,λ1,λ2,其中λ0>λ1>λ2,最小的特征值λ2对应的特征向量即为该点的法向,曲率为依次求出N个点对应的法向集合为曲率集合为{k1,k2,...,kN}。
作为优选方案,步骤四中所述的根据平面法向量间关系从基元平面集合ψ中抽象出平面间的结构关系的具体步骤包括:
本发明的根据平面基元间的结构关系对点云提取平面的方法,该方法围绕区域增长法实现对平面提取,主要在初始种子点的选取上进行改进,通过对源点云数据进行初步处理使其可以用于提取平面间通用的结构知识的计算,以具有正交或对称关系的平面作为种子面片合并其他面片。本方法利用平面间提取的结构知识,大大提高了密集匹配点云平面提取的可靠性。
上述方法可以通过计算机程序实现,因此,本发明相应提出一种存储介质,存储介质上存储有计算机程序,程序被处理器执行实现如上任一项所述的根据平面基元间的结构关系对点云提取平面的方法。
本发明还公开了一种建筑物精细化三维自动重建方法,该方法包括如上任一项所述的根据平面基元间的结构关系对点云提取平面的方法。
建筑物精细化三维自动重建方法一般包括以下步骤:1、建筑物三维自动化建模点云的获取;2、点云数据的处理,包括点云数据的配准拼接、去噪简化;3、建筑物平面的提取及拓扑重构。
其中,建筑物三维自动化建模点云的获取方式有多种,一般来说需要通过三维激光扫描设备获取建筑物平面的三维数据,设置三维激光扫描设备的扫描参数以及网络参数,使其能够获取目标建筑物的扫描数据并将扫描数据传输到数据处理装置进行数据处理,于实施例中,三维激光扫描设备可以通过无人机搭载。[李清泉.三维空间数据的实时获取、建模与可视化[M].武汉大学出版社,2003]详细介绍了三维空间数据获取的相关技术,本领域技术人员借助该文献即可实现本方法所述的建筑物三维自动化建模点云的获取。具体实施例中,数据处理装置一般采用电脑设备。
由于在获取建筑物三维自动化建模点云时,为了获取建筑物完整的三维数据,一般从多站点对目标建筑物进行扫描,因此,对于初始获取的点云数据需要对其进行处理提取目标建筑物的完整数据。一般来说,对于该步骤中对于点云数据的处理包括点云配准拼接、去燥简化等,上述点云配准拼接、去燥简化等均具有已公开的具体实现方法。例如[蒋荣华.地面三维激光扫描点云配准研究综述[J].科技创新与生产力,2016(12)]一文详细介绍了点云配准相关的内容。
建筑物平面的提取是三维重建的一个重要环节,该环节的质量与效率对三维重建而言至关重要。本申请中的建筑物三维自动重建方法,所述建筑物平面的提取采用上述的根据平面基元间的结构关系对点云提取平面的方法,首先通过邻域协方差求出各点的法向量与曲率,保证点云噪声较大时依然能求出较可靠的法向量与曲率;区域生长法保证了提取平面内的点共面的准确性;最后利用平面间的共面、正交与对称关系,可确定正确无误的平面,再将剩余的平面合并到正确的平面上得出最终的平面集合,这样可大大提高了密集匹配点云平面提取的可靠性。
具体实施方式
下面对本专利的优选实施方案作进一步详细的说明。
本发明的实施例公开一种根据平面基元间的结构关系对点云提取平面的方法,用于从建筑物三维点云中提取平面,包括以下步骤:
步骤一,根据邻域协方差矩阵计算点云对应的法向量与曲率:
对于通过三维激光扫描设备获取的点云数据P,点云其中N表示顶点个数, R3指三维空间;为点云P中的一个点,其K邻域为其中K为点的K个邻近点。根据K邻域计算协方差矩阵C,该矩阵有三个特征值,假设为λ0,λ1,λ2,其中λ0>λ1>λ2,则最小的特征值λ2对应的特征向量即为该点的法向,曲率为用此方法依次求出N个点对应的法向集合为曲率集合为{k1,k2,...,kN}。
步骤二,从点云中提取平面基元:
这里采用区域生长法提取平面基元,将曲率为0的顶点作为种子顶点,如果两个种子顶点间的法向相同且两顶点是彼此的最相邻点,则这两个种子顶点可合并,假设合并后的种子顶点为剩余顶点集合为reminder=P-seed。生长规则为若某剩余顶点与某种子顶点相似,则将从reminder中取出放入到某种子顶点对应的顶点集合Pq中,相似性采用顶点间的欧式距离及曲率衡量,本实施例中,两个顶点的欧式距离差及曲率差的加权平均越小则这两个顶点越相似。终止条件为即直至不存在满足生长规则的顶点则生长停止。最终得到平面基元集合ψ={ψ1,ψ2,...,ψN},其中平面ψi对应的顶点记作且
步骤三:提取上述平面基元间的结构关系
从基元平面集合ψ中抽象出平面间的结构关系,包括共面、平行、正交和对称。如果发现两个面正交或对称,则可认定这两个面为正确的平面,将该平面放入种子平面集合A中。这里的判断依据是平面法向量间的关系,如面ψ1和ψ2对应的法向量分别为有以下几种关系:
检测完ψ中的平面后得出种子平面A与剩余平面B,其中B=ψ-A。以A中的平面为种子平面,将B中的平面合并到A中,即如果A中的平面A1与B中的平面B1共面且两平面所包含的顶点在空间位置上邻近,则将面B1从B中取出合并到A中的平面A1里得到新的平面A′1。经过合并处理后得到最终的面片集合为A′={A′1,A′2,...,A′N}。
上述的根据平面基元间的结构关系对点云提取平面的方法,首先通过邻域协方差求出各点的法向量与曲率,保证点云噪声较大时依然能求出较可靠的法向量与曲率;区域生长法保证了提取平面内的点共面的准确性;最后利用平面间的共面、正交与对称关系,可确定正确无误的平面,再将剩余的平面合并到正确的平面上得出最终的平面集合,这样可大大提高了密集匹配点云平面提取的可靠性。
上述方法可使用软件编程实现,相应的,本发明还请求保护一种存储介质,该存储介质上存储有计算机程序,程序被处理器执行实现如上所述的根据平面基元间的结构关系对点云提取平面的方法。
建筑物精细化三维自动化重建方法一般包括以下步骤:1、建筑物三维自动化建模点云的获取;2、点云数据的处理,包括点云数据的配准拼接、去噪简化;3、建筑物平面的提取及拓扑重构。其中所述的建筑物平面的提取采用上述的根据平面基元间的结构关系对点云提取平面的方法,利用平面间提取的结构知识,大大提高了密集匹配点云平面提取的可靠性。
Claims (5)
步骤一,计算点云对应的法向量与曲率:
对于点云P,根据领域协方差矩阵求出各点的法向量与曲率;
步骤二,从点云中提取平面基元:
对于步骤一的处理结果,通过区域生长法提取平面基元;
步骤三,提取上述平面基元间的结构关系,如果两个平面有正交或对称的关系,则可被认定为正确的平面,并以这些正确的平面为种子平面合并其他的平面:
根据平面法向量间关系从基元平面集合ψ中抽象出平面间的结构关系,包括共面、平行、正交和对称,得到平面合集A和平面合集B,其中平面合集A中的平面之间为正交或对称关系,平面B中的合集为平行或共面关系;
按区域生长法,以A中的平面为种子平面,将B中的平面合并到A中。
4.一种建筑物精细化三维自动重建方法,包括以下步骤:
建筑物三维自动化建模点云的获取;
点云数据处理,包括,配准拼接和去噪简化;
建筑物平面的提取即拓扑重建;
其特征在于,所述建筑物平面的提取采用如权利要求1-3中任一项所述的根据平面基元间的结构关系对点云提取平面的方法。
5.一种存储介质,其特征在于,存储介质上存储有计算机程序,程序被处理器执行实现如权利要求1-3中任一项所述的根据平面基元间的结构关系对点云提取平面的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910175983.5A CN110009726B (zh) | 2019-03-08 | 2019-03-08 | 一种根据平面基元间的结构关系对点云提取平面的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910175983.5A CN110009726B (zh) | 2019-03-08 | 2019-03-08 | 一种根据平面基元间的结构关系对点云提取平面的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110009726A CN110009726A (zh) | 2019-07-12 |
CN110009726B true CN110009726B (zh) | 2022-09-30 |
Family
ID=67166668
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910175983.5A Active CN110009726B (zh) | 2019-03-08 | 2019-03-08 | 一种根据平面基元间的结构关系对点云提取平面的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110009726B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112489207B (zh) * | 2021-02-07 | 2021-07-13 | 深圳大学 | 一种空间约束的密集匹配点云平面基元提取方法 |
CN114821013B (zh) * | 2022-07-01 | 2022-10-18 | 深圳大学 | 基于点云数据的基元检测方法、装置和计算机设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104395932A (zh) * | 2012-06-29 | 2015-03-04 | 三菱电机株式会社 | 对准数据的方法 |
CN104463872A (zh) * | 2014-12-10 | 2015-03-25 | 武汉大学 | 基于车载LiDAR点云数据的分类方法 |
CN105631929A (zh) * | 2014-11-28 | 2016-06-01 | 富泰华工业(深圳)有限公司 | 点云精简方法及系统 |
CN105869149A (zh) * | 2016-03-24 | 2016-08-17 | 大连理工大学 | 基于主向量分析的断骨截面分割及断骨模型配准方法 |
CN106780509A (zh) * | 2016-12-01 | 2017-05-31 | 山东交通学院 | 融合多维特征的建筑物点云层次聚类分割方法 |
CN107223268A (zh) * | 2015-12-30 | 2017-09-29 | 中国科学院深圳先进技术研究院 | 一种三维点云模型重建方法及装置 |
CN107845099A (zh) * | 2017-10-18 | 2018-03-27 | 安徽佳通乘用子午线轮胎有限公司 | 一种面向轮胎点云的数据分割方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5465128B2 (ja) * | 2010-08-11 | 2014-04-09 | 株式会社トプコン | 点群位置データ処理装置、点群位置データ処理システム、点群位置データ処理方法、および点群位置データ処理プログラム |
JP5963353B2 (ja) * | 2012-08-09 | 2016-08-03 | 株式会社トプコン | 光学データ処理装置、光学データ処理システム、光学データ処理方法、および光学データ処理用プログラム |
-
2019
- 2019-03-08 CN CN201910175983.5A patent/CN110009726B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104395932A (zh) * | 2012-06-29 | 2015-03-04 | 三菱电机株式会社 | 对准数据的方法 |
CN105631929A (zh) * | 2014-11-28 | 2016-06-01 | 富泰华工业(深圳)有限公司 | 点云精简方法及系统 |
CN104463872A (zh) * | 2014-12-10 | 2015-03-25 | 武汉大学 | 基于车载LiDAR点云数据的分类方法 |
CN107223268A (zh) * | 2015-12-30 | 2017-09-29 | 中国科学院深圳先进技术研究院 | 一种三维点云模型重建方法及装置 |
CN105869149A (zh) * | 2016-03-24 | 2016-08-17 | 大连理工大学 | 基于主向量分析的断骨截面分割及断骨模型配准方法 |
CN106780509A (zh) * | 2016-12-01 | 2017-05-31 | 山东交通学院 | 融合多维特征的建筑物点云层次聚类分割方法 |
CN107845099A (zh) * | 2017-10-18 | 2018-03-27 | 安徽佳通乘用子午线轮胎有限公司 | 一种面向轮胎点云的数据分割方法 |
Non-Patent Citations (1)
Title |
---|
基于平面提取的点云数据分割算法;李宝顺等;《计算机应用于软件》;20140731;第31卷(第7期);第145-148页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110009726A (zh) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102154470B1 (ko) | 다중 특징 검색 및 변형에 기반한 3차원 인체 헤어스타일 생성 방법 | |
CN112002014B (zh) | 面向精细结构的三维人脸重建方法、系统、装置 | |
CN106780619B (zh) | 一种基于Kinect深度相机的人体尺寸测量方法 | |
Brenner et al. | Coarse orientation of terrestrial laser scans in urban environments | |
CN110599528A (zh) | 一种基于神经网络的无监督三维医学图像配准方法及系统 | |
CN110796694A (zh) | 一种基于KinectV2的果实三维点云实时获取方法 | |
CN108154531B (zh) | 一种计算体表损伤区域面积的方法和装置 | |
Sinko et al. | 3D registration of the point cloud data using ICP algorithm in medical image analysis | |
CN110009726B (zh) | 一种根据平面基元间的结构关系对点云提取平面的方法 | |
CN110009745B (zh) | 根据平面基元和模型驱动对点云提取平面的方法 | |
CN107481274A (zh) | 一种三维作物点云的鲁棒性重构方法 | |
Gonzalez-Aguilera et al. | New tools for rock art modelling: automated sensor integration in Pindal Cave | |
CN112329726B (zh) | 人脸识别方法及装置 | |
CN113689535A (zh) | 基于无人机影像的建筑物模型的生成方法和装置 | |
Piras et al. | Current options for visualization of local deformation in modern shape analysis applied to paleobiological case studies | |
Zhou et al. | Point cloud registration for agriculture and forestry crops based on calibration balls using Kinect V2 | |
CN112712590B (zh) | 一种动物点云生成方法及系统 | |
Zvietcovich et al. | A novel method for estimating the complete 3D shape of pottery with axial symmetry from single potsherds based on principal component analysis | |
CN112581511A (zh) | 基于近似直立扫描点云快速配准的三维重建方法及系统 | |
CN116797733A (zh) | 实时三维物体动态重建方法 | |
CN114399547B (zh) | 一种基于多帧的单目slam鲁棒初始化方法 | |
CN110232695A (zh) | 基于混合模态图像的左心室图像分割方法及系统 | |
CN115861532A (zh) | 基于深度学习的植被地物模型重建方法与系统 | |
WO2014106747A1 (en) | Methods and apparatus for image processing | |
CN108921908B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |