CN111986322A - 一种基于结构分析的点云室内场景布局重建方法 - Google Patents
一种基于结构分析的点云室内场景布局重建方法 Download PDFInfo
- Publication number
- CN111986322A CN111986322A CN202010703972.2A CN202010703972A CN111986322A CN 111986322 A CN111986322 A CN 111986322A CN 202010703972 A CN202010703972 A CN 202010703972A CN 111986322 A CN111986322 A CN 111986322A
- Authority
- CN
- China
- Prior art keywords
- point
- points
- room
- indoor scene
- point cloud
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000012916 structural analysis Methods 0.000 title claims abstract description 28
- 239000013598 vector Substances 0.000 claims description 78
- 239000011159 matrix material Substances 0.000 claims description 37
- 238000001514 detection method Methods 0.000 claims description 20
- 230000011218 segmentation Effects 0.000 claims description 14
- 238000000354 decomposition reaction Methods 0.000 claims description 13
- 238000000513 principal component analysis Methods 0.000 claims description 12
- 238000002372 labelling Methods 0.000 claims description 7
- 230000003416 augmentation Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000010146 3D printing Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Computer Graphics (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于结构分析的点云室内场景布局重建方法,具体按照如下步骤实施:步骤1,对室内场景进行预处理操作;步骤2,对步骤1进行预处理的室内场景进行粗分割,提取房间的轮廓;步骤3,根据经步骤2提取的房间的轮廓对室内场景房间的布局进行划分,划分出房间布局;步骤4,分割经步骤1预处理的室内场景的不同面的墙壁并进行重建,然后在步骤3划分的房间布局进行室内场景的全局结构重建,重建完整的室内结构。本发明的一种基于结构分析的点云室内场景布局重建方法,解决了现有技术中存在的由于需要考虑场景中的大量数据特征,导致算法设计较为复杂且处理速度较慢的问题。
Description
技术领域
本发明属于计算机视觉和图像处理方法技术领域,涉及一种基于结构分析的点云室内场景布局重建方法。
背景技术
在计算机视觉与图像处理领域,对各种场景的三维建模一直都是一个具有重要研究意义和应用价值的课题。其在移动机器人自主导航、SLAM、虚拟现实和3D打印等诸多方面都发挥着十分重要的作用。现实世界的各种场景中,室内场景是与人们工作和生活最为密切相关的场所,是人类感知世界的最重要组成部分,随着机器人技术、计算机视觉以及人工智能等技术的飞速发展,室内场景的重建显得更为迫切。室内场景的三维建模可以应用于火灾救援、密室探险和考古挖掘等必不可少却能够对人类生命造成威胁的探测活动中。但基于三维激光点云数据的室内场景重建存在以下三个难点:①玻璃窗户、镜面等非漫反射物体造成三维点云数据采集的困难,容易干扰到对环境的理解;②室内场景杂物较多,环境相对复杂,存在遮挡的可能较大,这会导致信息不全使物体难以识别;③大范围采集场景数据产生大量噪声,点云碎片化严重,不利于场景分割与理解的展开。因此,对于室内场景的三维重建仍然是一个具有挑战性的研究。
目前较为常用的几种室内场景重建方法有基于数据驱动的方法、基于模型驱动的方法以及基于机器学习的方法。基于数据驱动的方法能够考虑到场景数据的全面性,通过分析建筑物点云的数据特征,唯一地确定建筑物形状,但是由于需要考虑场景中的大量数据特征,导致算法设计较为复杂,且处理速度较慢。基于模型驱动的方法对于建筑物重建模型时规则化程度较高,在数据采集过程中由局部遮挡引起的部分建筑数据缺失也能得到很好的弥补,处理方法较为简单高效,但是基于模型重建时并未完全依据数据,其结果可能会产生错误。基于机器学习的方法使用了先进的算法与技术,从大量训练集数据中学习并分析数据的各种特征,形成一系列场景重建的方法过程,此类方法对数据依赖程度较高,学习不同的数据集会获取对应的不同结果,数据集的选择越规范,重建结果更精确,但是该类方法一般所需数据量都较为庞大,且对设备性能要求较高,目前难以大量推广。
发明内容
本发明的目的是提供一种基于结构分析的点云室内场景布局重建方法,解决了现有技术中存在的由于需要考虑场景中的大量数据特征,导致算法设计较为复杂且处理速度较慢的问题。
本发明所采用的技术方案是,一种基于结构分析的点云室内场景布局重建方法,具体按照如下步骤实施:
步骤1,对室内场景进行预处理操作;
步骤2,对步骤1进行预处理的室内场景进行粗分割,提取房间的轮廓;
步骤3,根据经步骤2提取的房间的轮廓对室内场景房间的布局进行划分,划分出房间布局;
步骤4,分割经步骤1预处理的室内场景的不同面的墙壁并进行重建,然后在步骤3划分的房间布局进行室内场景的全局结构重建,重建完整的室内结构。
本发明的特征还在于,
步骤1具体为:
步骤1.1,通过PCA主成分分析法获取室内场景主平面,由主平面的法向量确定旋转方向,将原始室内场景点云数据进行场景旋转操作,使其处于自然竖直状态;
步骤1.2,去掉经步骤1.1旋转的室内场景点云数据的地板与天花板数据;
步骤1.3,对经步骤1.2去除地板与天花板数据的室内场景点云数据使用PCA主成分分析法拟合局部平面来确定法向量,根据法向量滤掉非墙壁数据点。
步骤1.1中将原始室内场景点云数据进行场景旋转操作具体按照如下步骤实施:
步骤1.1.1,室内场景点云数据为P={p1,p2,…,pn},其中,pi表示点云文件中第i个点的坐标信息(xi,yi,zi),将原始室内场景点云表示为n×3的矩阵A={XT,YT,ZT},其中,X={X1,X2,…,Xn},Y={Y1,Y2,…,Yn},Z={Z1,Z2,…,Zn},其中,Xi是第i个点的x轴坐标点,Yi是第i个点的y轴坐标点,Zi是第i个点的z轴坐标点,i=1,2,3…n;
步骤1.1.2,在矩阵A中每一位坐标值减去各自的平均值得到矩阵D;
步骤1.1.4,采用奇异值分解法求协方差矩阵的特征值与特征向量,对特征值按从大到小排序,选择前两个特征值v1和v2,对应的特征向量v1=(a1,b1,c1)和v2=(a2,b2,c2)即为所求得的场景主次趋势方向V;
V=v1×v2=(b1c2-c1b2,c1a2-a1c2,a1b2-b1a2) (3)
步骤1.1.5,通过旋转V,使得场景处于自然竖直状态,具体为:
V′为V在坐标平面XOY内的投影向量,与x轴之间的夹角为α,V″为向量V在坐标平面xoz内的投影向量,与z轴之间的夹角为β,将场景绕z坐标轴顺时针旋转角度α,确定第一个旋转矩阵为Rz,然后顺时针绕y轴旋转角度β,确定了第二个旋转矩阵为Ry,将点云场景转化为n×3的矩阵A后进行旋转,得到最终的旋转结果为Vres=(ARz)Ry=A(RzRy),其中,
步骤1.2具体按照下述步骤进行:
步骤1.2.1,统计经步骤1.1旋转的室内场景点云数据在z坐标方向上的分布构造直方图,其中,组距的设置是根据点的z坐标范围,将其均匀划分为200份,将其间隔作为组距μ,然后根据组距确定每组数据的坐标范围及处于该范围的点的数目完成直方图的构造;
步骤1.2.2,将直方图中z坐标的某个极值设为预期峰值peak;
步骤1.2.3,找出以peak为中心,半径为radius的区域中出现的全体点集M,将其归属于一个聚类C,其中,radius为μ的四倍;
步骤1.2.4,以peak为中心点,计算从peak到集合M中每个数据点构成的向量,叠加所有向量,可得到漂移向量shift,peak=peak+shift,即peak沿着shift的方向移动,移动距离是||shift||;
步骤1.2.5,重复步骤1.2.3、1.2.4,直到shift小于设定的阈值σ,即就是迭代到收敛或达到了最大迭代次数10次,记录当前的peak,迭代中所有出现在以peak为中心radius为半径圆内的点都归到一类C中;
步骤1.2.6,如果收敛时,当前簇C的peak与其它已经存在的簇C1中心的距离小于阈值σ,那么把C和C1合并,数据点出现次数也对应合并,其中,σ为μ的二倍,否则,把C作为新的聚类;
步骤1.2.7,找到天花板和地面的聚类点集后,将其去除。
步骤1.3具体按照如下步骤实施:
步骤1.3.1,对经步骤1.2去除地板与天花板数据的室内场景点云数据用KD-tree组织点云模型中的所有点;
步骤1.3.2,对于经步骤1.3.1组织的所有点中某个点p的法向量求解,首先利用KD-tree查询其k个近邻点;
步骤1.3.3,根据k个近邻点构造协方差矩阵;
步骤1.3.4,用奇异值分解算法求解步骤1.1.3中协方差矩阵最大的两特征向量的向量积,该向量积就是当前点的法向量;
步骤1.3.5,重复步骤1.3.2、1.3.3和1.3.4,求出所有点的法向量;
步骤1.3.6,遍历所有点,如果对于点pi,若|ni·nz|<Tangle,则删除点pi,其中,ni为点pi的法向量,nz=(0,0,1)是沿纵轴的单位向量,据此删除所有此类非墙壁点,其中,Tangle∈(0,2]。
步骤2具体为:
步骤2.1,对经步骤1预处理的室内场景点云数据投影到坐标平面xoy中,然后将x,y坐标不同的点作为数据投影的有效点,得到原始室内场景结构轮廓,即就是投影数据;
步骤2.2,应用基于改进的Ransac算法,对步骤2.1投影得到的原始室内场景结构轮廓进行直线检测,具体按照下述步骤进行:
步骤2.2.1,创建直线集合H,用于存放直线模型;
步骤2.2.2,构造投影数据的二维KD-tree;
步骤2.2.3,求解未测点p的k近邻,选取点p与其中一个邻近点求解直线模型参数,并将该邻近点标记为已检测;
步骤2.2.4,在最大误差范围内统计符合该直线模型的数据点数目,如果符合模型参数的样本点,即内点所占的比例大于设定的阈值0.8,则记录该直线,否则利用所有内点重新确定直线方程参数,然后进入下一次迭代;
步骤2.2.5,重复2.2.2~2.2.4,直到p的邻近点均已被检测;
步骤2.2.6,计算该直线模型的方向向量,在检测结果H中查询共线或接近共线的直线集合,加入其中,否则,无共线直线时新建类别H′,加入H中,然后标记该直线的内点;
步骤2.2.7,重复2.2.2~2.2.6,直到所有点已被检测,得到检测结果H={H1,H2,…,Hn};
步骤2.3,采用基于Mean shift方法将直线检测结果进行聚类,去除较短直线以及杂乱直线的影响,最后得到与墙壁一一对应的直线检测结果,从而得到室内房间轮廓。
步骤2.3具体为:
步骤2.3.1,选择经2.2步骤得到的直线分割结果簇Hi,在结果簇中随机选择一条作为起始中心直线,设定center=0,表示与该直线间的距离;
步骤2.3.2,找出与中心直线center距离小于radius的所有直线作为集合U,认为这些直线同属于一个聚类HCi,同时在该聚类中将该直线出现的次数加1;
步骤2.3.3,以center为基准,计算从center到集合U中每条直线的距离,为center其中一侧的直线距离设定为负值,将这些距离相加,得到漂移量distance;
步骤2.3.4,更新中心直线center=center+distance,即中心直线沿着distance的方向移动,移动距离是|distance|,根据center与直线结果簇Hi的直线方向确定新的中心直线,清空集合U,为下一次均值漂移做准备;
步骤2.3.5,重复上述2.3.2~2.3.4步骤,直到漂移量distance的值小于阈值μ/5,即为迭代到收敛,记录此时的center,将此迭代中经过的所有直线均归类为簇Hi;
步骤2.3.6,重复2.3.1~2.3.5直到归类完所有的结果簇H;
步骤2.3.7,如果簇Hi的center与其它已经存在的簇Hj的center之间距离小于墙壁厚度2μ,那么把Hj和Hi合并,直线出现次数也对应合并;
步骤2.3.8,分类:根据每个类,对每个点的访问频率,取访问频率最大的m个类,每个类别以其中心直线作为直线优化结果,得到的直线集合L={l1,l2,…,ln}作为潜在墙壁优化结果,即就是室内房间轮廓,其中,m值根据数据中房间数目预先设定。
步骤3具体为:
步骤3.1,由公式(5)确定扫描仪位置P,Min.Ptx为x坐标最小值,Min.Pty为y坐标最小值,Max.Pty为y坐标最小值,Min.Ptz为z坐标最小值,Max.Ptz为z坐标最大值;
P=(Min.Ptx,(Min.Pty+Max.Pty)/2,(Min.Ptz+Max.Ptz)/2) (5)
步骤3.2,求出步骤2房间粗分割结果中的每一个预期房间元素的右侧边界中心点;
步骤3.3,以点P为起始点,向右延伸穿过预期房间右侧边界中心点构造射线Ray;
步骤3.4,求解射线Ray与该房间的交点数目Number,若交点数目Number为偶数,将该房间标签化为1,表示为潜在房间,否则,将该房间标签化为0;
步骤3.5,将每个预期房间元素转化为无向图G中的一个顶点,对于能相互连接的预期房间添加边信息,其中边是具有权值的,边的权值通过与之连接的两个预期房间的中心点的欧氏距离定义,其中,房间中心点通过房间的顶点求均值得到,具体如公式(5)所示:
x,y,z为房间中心点的三维坐标,而xi,yi,zi是二维平面图中房间的顶点Pi的三维坐标;
步骤3.6,确定了每个房间的中心点之后,通过公式(7)求解房间中心点O和Q的欧氏距离,从而得到房间中心点及其欧氏距离的求解结果;
其中,xO,yO,zO是房间中心点O的三维坐标,xQ,yQ,zQ是房间中心点Q的三维坐标;
步骤3.7,向该无向图G中再添加了2个顶点,分别为源点“S”和汇点“T”,设置与源点相连的顶点权值为0,汇点相连的顶点权值为无穷大;
步骤3.8,建立无向图G的能量方程,如公式(8)所示;
E=∑i∈VDi(li)+Ri,j(li,lj) (8)
其中,Di为数据项,Ri,j为约束项,数据项由式(9)确定,约束项由式(10)确定;
其中,li为房间i的标签,系数c=0.6,Si为归一化后的面积,即预期房间面积除以整个预期房间的初步划分结果的面积;
Li,j为预期房间i、j的欧氏距离,该距离除以整个预期房间划分结果的最长边可做归一化处理,权重ωi,j为该预期房间占总房间初步划分结果总面积的比例;
步骤3.9,使用最大流/最小割算法对其进行能量最小化求解,找到从源点S到汇点T的一条增广路径,这条增广路径作为无向图G的最小割,最终将最小割映射到房间粗分割结果中,从而得到房间布局信息。
步骤4具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体按照下述步骤进行:
步骤4.1.1,为经步骤1处理后的室内场景点云数据构造KD_tree;
步骤4.1.2,在KD_tree中搜索当前点的K个邻近点;
步骤4.1.3,根据当前点以及其K个近邻点构造协方差矩阵;
步骤4.1.4,采用奇异值分解算法进行特征分解,最大的两特征向量的向量积即为当前点的法向量;
步骤4.1.5,重复上述4.1.2~4.1.4步骤,直到求出所有点的法向量;
步骤4.1.6,为点云数据中的每个点定义平面标记,初始标记为0,即FlatIndex=0,表示该点未被标记;
步骤4.1.7,遍历点云数据中的每个点,若该点的平面标记为0,为当前点设置墙壁类别为FlatIndex;
步骤4.1.8,将该邻近点入栈,在KD_tree中搜索该点的K个邻近点;
步骤4.1.9,判断该点与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.10,若栈不为空,取栈顶元素为搜索点,在KD_tree中搜索该点的K个邻近点,判断栈顶元素与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.11,重复上述步骤直到标记完所有的点;
步骤4.1.12,遍历点云数据中的每一个点,根据每个点的平面类别FlatIndex组织点云为不同墙壁,将点云数目小于阈值的墙壁删除得到结果Walls;
步骤4.2,对分割的墙壁进行栅格化,分析各墙壁是否含有门窗,对有门窗的墙壁,将各栅格区别分析并拼接,对没有门窗的墙壁,使用OBB包围盒的简单方式进行重建;
步骤4.3,对全局结构进行重建。
步骤4.2具体为:
步骤4.2.1,将墙壁划分为35×20的栅格,栅格尺寸的长度Gl和宽度Gw分别为:Gl=Wl/35,Gw=Ww/20,其中,Wl为墙壁长度,Ww为墙壁宽度;
步骤4.2.2,统计该栅格单元中点的数目Numg,根据式(11)求得单位面积的点密度Dg,并根据公式(12)计算分类阈值Dw:
其中,Dgi表示某栅格单元i的点密度;
步骤4.2.3,判断某个栅格单元i使得Dgi<Dw,则该墙壁为复杂墙壁,否则,该墙壁视为单一墙壁;
步骤4.2.4,对单一墙壁应用OBB包围盒进行墙壁重建;
步骤4.2.5,对复杂墙壁,将点密度小于阈值Dw的单元格看作门窗开口位置,将开口位置从墙壁中删除,将点密度大于阈值Dw的栅格单元拼接起来,将墙壁拟合为包含门窗细节的结果;
步骤4.3具体为:
步骤4.3.1,在步骤3获取到的室内场景房间布局信息的基础上,将房间布局的二维分割结果在竖直方向即z坐标轴上进行延伸,向z坐标轴负方向延伸至地面的上界坐标处,向z轴正方向延伸至天花板下界坐标处,将房间布局平面图中的每条直线看作墙壁的中心线,为其赋予合适的厚度2μ,μ为组距,即将z坐标平均分为200份,得到室内场景各个墙壁在三维空间中的粗略位置;
步骤4.3.2,将墙壁粗略位置与步骤4.2场景分割得到的墙壁分类重建结果相对应,并将墙壁的z坐标范围限定在地面与天花板位置之间,从而重建得到每个房间的不同墙壁结构,然后在墙壁两端分别补齐地面与天花板平面。
本发明的有益效果是:
本发明解决的现有技术中存在的由于需要考虑场景中的大量数据特征,导致算法设计较为复杂,且处理速度较慢以及需数据量都较为庞大,且对设备性能要求较高,目前难以大量推广的问题,本发明室内场景点云数据处理简单有序,布局重建工作的简单,在有效性和稳定性都有了很大的提高,并且拥有较好的鲁棒性,丰富了计算机图形学和视觉智能的方法体系。
附图说明
图1是本发明一种基于结构分析的点云室内场景布局重建方法步骤1中得到自然竖直室内场景数据;
图2是本发明一种基于结构分析的点云室内场景布局重建方法步骤1中去除室内场景的地板和天花板的数据;
图3是本发明一种基于结构分析的点云室内场景布局重建方法步骤1中滤掉室内场景非墙壁的数据;
图4是本发明一种基于结构分析的点云室内场景布局重建方法步骤2中将室内场景数据投影到坐标平面XOY中的数据;
图5是本发明一种基于结构分析的点云室内场景布局重建方法步骤2中对投影结果进行直线检测的结果;
图6是本发明一种基于结构分析的点云室内场景布局重建方法步骤2中对直线结果进行聚类的结果;
图7是本发明一种基于结构分析的点云室内场景布局重建方法步骤3中划分出房间布局的结果;
图8是本发明一种基于结构分析的点云室内场景布局重建方法步骤4中分割出不同面的墙壁的结果;
图9是本发明一种基于结构分析的点云室内场景布局重建方法步骤4中栅格化墙壁并重建墙壁的结果;
图10是本发明一种基于结构分析的点云室内场景布局重建方法步骤4中对全局结构进行重建的结果。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种基于结构分析的点云室内场景布局重建方法,具体按照如下步骤实施:
步骤1,对室内场景进行预处理操作;具体为:
步骤1.1,通过PCA主成分分析法获取室内场景主平面,由主平面的法向量确定旋转方向,将原始室内场景点云数据进行场景旋转操作,使其处于自然竖直状态;其中,将原始室内场景点云数据进行场景旋转操作具体按照如下步骤实施:
步骤1.1.1,室内场景点云数据为P={p1,p2,…,pn},其中,pi表示点云文件中第i个点的坐标信息(xi,yi,zi),将原始室内场景点云表示为n×3的矩阵A={XT,YT,ZT},其中,X={X1,X2,…,Xn},Y={Y1,Y2,…,Yn},Z={Z1,Z2,…,Zn},其中,Xi是第i个点的x轴坐标点,Yi是第i个点的y轴坐标点,Zi是第i个点的z轴坐标点,i=1,2,3…n;
步骤1.1.2,在矩阵A中每一位坐标值减去各自的平均值得到矩阵D;
步骤1.1.4,采用奇异值分解法求协方差矩阵的特征值与特征向量,对特征值按从大到小排序,选择前两个特征值v1和v2,对应的特征向量v1=(a1,b1,c1)和v2=(a2,b2,c2)即为所求得的场景主次趋势方向V;
V=v1×v2=(b1c2-c1b2,c1a2-a1c2,a1b2-b1a2) (3)
步骤1.1.5,通过旋转V,使得场景处于自然竖直状态,如图1所示,具体为:
V′为V在坐标平面XOY内的投影向量,与x轴之间的夹角为α,V″为向量V在坐标平面xoz内的投影向量,与z轴之间的夹角为β,将场景绕z坐标轴顺时针旋转角度α,确定第一个旋转矩阵为Rz,然后顺时针绕y轴旋转角度β,确定了第二个旋转矩阵为Ry,由于旋转矩阵数据量较小,而原始数据较大,为了简化计算,可根据矩阵运算的结合律(AB)C=A(BC),可以先将两旋转矩阵进行相乘,然后对原始数据进行旋转,由于原始数据是一系列包含三维坐标值的散乱点集,将点云场景转化为n×3的矩阵A后进行旋转,得到最终的旋转结果为Vres=(ARz)Ry=A(RzRy),其中,
步骤1.2,去掉经步骤1.1旋转的室内场景点云数据的地板与天花板数据;具体按照下述步骤进行:
步骤1.2.1,统计经步骤1.1旋转的室内场景点云数据在z坐标方向上的分布构造直方图,其中,组距的设置是根据点的z坐标范围,将其均匀划分为200份,将其间隔作为组距μ,然后根据组距确定每组数据的坐标范围及处于该范围的点的数目完成直方图的构造;
步骤1.2.2,将直方图中z坐标的某个极值设为预期峰值peak;
步骤1.2.3,找出以peak为中心,半径为radius的区域中出现的全体点集M,将其归属于一个聚类C,其中,radius为μ的四倍;
步骤1.2.4,以peak为中心点,计算从peak到集合M中每个数据点构成的向量,叠加所有向量,可得到漂移向量shift,peak=peak+shift,即peak沿着shift的方向移动,移动距离是||shift||;
步骤1.2.5,重复步骤1.2.3、1.2.4,直到shift小于设定的阈值σ,即就是迭代到收敛或达到了最大迭代次数10次,记录当前的peak,迭代中所有出现在以peak为中心radius为半径圆内的点都归到一类C中;
步骤1.2.6,如果收敛时,当前簇C的peak与其它已经存在的簇C1中心的距离小于阈值σ,那么把C和C1合并,数据点出现次数也对应合并,其中,σ为μ的二倍,否则,把C作为新的聚类;
步骤1.2.7,找到天花板和地面的聚类点集后,将其去除,如图2所示;
步骤1.3,对经步骤1.2去除地板与天花板数据的室内场景点云数据使用PCA主成分分析法拟合局部平面来确定法向量,根据法向量滤掉非墙壁数据点,具体按照如下步骤实施:
步骤1.3.1,对经步骤1.2去除地板与天花板数据的室内场景点云数据用KD-tree组织点云模型中的所有点;
步骤1.3.2,对于经步骤1.3.1组织的所有点中某个点p的法向量求解,首先利用KD-tree查询其k个近邻点;
步骤1.3.3,根据k个近邻点构造协方差矩阵;
步骤1.3.4,用奇异值分解算法求解步骤1.1.3中协方差矩阵最大的两特征向量的向量积,该向量积就是当前点的法向量;
步骤1.3.5,重复步骤1.3.2、1.3.3和1.3.4,求出所有点的法向量;
步骤1.3.6,遍历所有点,如果对于点pi,若|ni·nz|<Tangle,则删除点pi,其中,ni为点pi的法向量,nz=(0,0,1)是沿纵轴的单位向量,据此删除所有此类非墙壁点,如图3所示,其中,考虑到角度的误差,针对不同的数据为其设置角度阈值Tangle∈(0,2];
步骤2,对步骤1进行预处理的室内场景进行粗分割,提取房间的轮廓;具体为:
步骤2.1,对经步骤1预处理的室内场景点云数据投影到坐标平面xoy中,是通过舍弃第三个维度(z坐标)的方式将点云投影到地面,舍弃之后,需要判断其余两维坐标是否相同,只保留x,y坐标不同的点作为数据投影的有效点,得到原始室内场景结构轮廓,即就是投影数据,如图4所示;
步骤2.2,应用基于改进的Ransac算法,对步骤2.1投影得到的原始室内场景结构轮廓进行直线检测,具体按照下述步骤进行:
步骤2.2.1,创建直线集合H,用于存放直线模型;
步骤2.2.2,构造投影数据的二维KD-tree;
步骤2.2.3,求解未测点p的k近邻,选取点p与其中一个邻近点求解直线模型参数,并将该邻近点标记为已检测;
步骤2.2.4,在最大误差范围内统计符合该直线模型的数据点数目,如果符合模型参数的样本点,即内点所占的比例大于设定的阈值0.8,则记录该直线,否则利用所有内点重新确定直线方程参数,然后进入下一次迭代;
步骤2.2.5,重复2.2.2~2.2.4,直到p的邻近点均已被检测;
步骤2.2.6,计算该直线模型的方向向量,在检测结果H中查询共线或接近共线的直线集合,加入其中,否则,无共线直线时新建类别H′,加入H中,然后标记该直线的内点;
步骤2.2.7,重复2.2.2~2.2.6,直到所有点已被检测,得到检测结果H={H1,H2,…,Hn},如图5所示;
步骤2.3,采用基于Mean shift方法将直线检测结果进行聚类,去除较短直线以及杂乱直线的影响,最后得到与墙壁一一对应的直线检测结果,从而得到室内房间轮廓,具体为:
步骤2.3.1,选择经2.2步骤得到的直线分割结果簇Hi,在结果簇中随机选择一条作为起始中心直线,设定center=0,表示与该直线间的距离;
步骤2.3.2,找出与中心直线center距离小于radius的所有直线作为集合U,认为这些直线同属于一个聚类HCi,同时在该聚类中将该直线出现的次数加1;
步骤2.3.3,以center为基准,计算从center到集合U中每条直线的距离,为center其中一侧的直线距离设定为负值,将这些距离相加,得到漂移量distance;
步骤2.3.4,更新中心直线center=center+distance,即中心直线沿着distance的方向移动,移动距离是|distance|,根据center与直线结果簇Hi的直线方向确定新的中心直线,清空集合U,为下一次均值漂移做准备;
步骤2.3.5,重复上述2.3.2~2.3.4步骤,直到漂移量distance的值小于阈值μ/5,即为迭代到收敛,记录此时的center,将此迭代中经过的所有直线均归类为簇Hi;
步骤2.3.6,重复2.3.1~2.3.5直到归类完所有的结果簇H;
步骤2.3.7,如果簇Hi的center与其它已经存在的簇Hj的center之间距离小于墙壁厚度2μ,那么把Hj和Hi合并,直线出现次数也对应合并;
步骤2.3.8,分类:根据每个类,对每个点的访问频率,取访问频率最大的m个类,每个类别以其中心直线作为直线优化结果,得到的直线集合L={l1,l2,…,ln}作为潜在墙壁优化结果,即就是室内房间轮廓,其中,m值根据数据中房间数目预先设定,如图6所示;
步骤3,根据经步骤2提取的房间的轮廓对室内场景房间的布局进行划分,划分出房间布局;如图7所示,具体为:
步骤3.1,由公式(5)确定扫描仪位置P,Min.Ptx为x坐标最小值,Min.Pty为y坐标最小值,Max.Pty为y坐标最小值,Min.Ptz为z坐标最小值,Max.Ptz为z坐标最大值;
P=(Min.Ptx,(Min.Pty+Max.Pty)/2,(Min.Ptz+Max.Ptz)/2)(5)
步骤3.2,求出步骤2房间粗分割结果中的每一个预期房间元素的右侧边界中心点;
步骤3.3,以点P为起始点,向右延伸穿过预期房间右侧边界中心点构造射线Ray;
步骤3.4,求解射线Ray与该房间的交点数目Number,若交点数目Number为偶数,将该房间标签化为1,表示为潜在房间,否则,将该房间标签化为0;
步骤3.5,将每个预期房间元素转化为无向图G中的一个顶点,对于能相互连接的预期房间添加边信息,其中边是具有权值的,边的权值通过与之连接的两个预期房间的中心点的欧氏距离定义,其中,房间中心点通过房间的顶点求均值得到,具体如公式(5)所示:
x,y,z为房间中心点的三维坐标,而xi,yi,zi是二维平面图中房间的顶点Pi的三维坐标;
步骤3.6,确定了每个房间的中心点之后,通过公式(7)求解房间中心点O和Q的欧氏距离,从而得到房间中心点及其欧氏距离的求解结果;
其中,xO,yO,zO是房间中心点O的三维坐标,xQ,yQ,zQ是房间中心点Q的三维坐标;
步骤3.7,向该无向图G中再添加了2个顶点,分别为源点“S”和汇点“T”,设置与源点相连的顶点权值为0,汇点相连的顶点权值为无穷大;
步骤3.8,建立无向图G的能量方程,如公式(8)所示;
E=∑i∈VDi(li)+Ri,j(li,lj) (8)
其中,Di为数据项,Ri,j为约束项,数据项由式(9)确定,约束项由式(10)确定;
数据项综合考虑预期房间的不同尺寸与其标签化结果,为其设置不同权值。由于房间的面积一般不是特别大,所以对面积较大的预期房间赋予较大惩罚值,则用式(9)求解数据项,以增加数据本身对房间分割的有效约束,其中,根据先验知识:房间内部一般会采集到更多的点云数据,则为空单元格设置较大惩罚值,而使得非空预期房间更容易分类为房间,故当预期房间为空时乘以系数c=0.6,非空时乘以系数1-c,Si为归一化后的面积,即预期房间面积除以整个预期房间的初步划分结果的面积;
约束项用来降低最终模型复杂性,其根据相邻预期房间的欧氏距离按比例设置惩罚值,Li,j为预期房间i、j的欧氏距离,该距离除以整个预期房间划分结果的最长边可做归一化处理。为了结果更准确,可非空预期房间减小权重以降低权值,从而增加保留程度。约束项Ri,j的计算公式如式(10)所示,其中,权重为该预期房间占总房间初步划分结果总面积的比例;
Li,j为预期房间i、j的欧氏距离,该距离除以整个预期房间划分结果的最长边可做归一化处理,权重ωi,j为该预期房间占总房间初步划分结果总面积的比例;
步骤3.9,使用最大流/最小割算法对其进行能量最小化求解,找到从源点S到汇点T的一条增广路径,这条增广路径作为无向图G的最小割,最终将最小割映射到房间粗分割结果中,从而得到房间布局信息;
步骤4,分割经步骤1预处理的室内场景的不同面的墙壁并进行重建,然后在步骤3划分的房间布局进行室内场景的全局结构重建,重建完整的室内结构,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体按照下述步骤进行:
步骤4.1.1,为经步骤1处理后的室内场景点云数据构造KD_tree;
步骤4.1.2,在KD_tree中搜索当前点的K个邻近点;
步骤4.1.3,根据当前点以及其K个近邻点构造协方差矩阵;
步骤4.1.4,采用奇异值分解算法进行特征分解,最大的两特征向量的向量积即为当前点的法向量;
步骤4.1.5,重复上述4.1.2~4.1.4步骤,直到求出所有点的法向量;
步骤4.1.6,为点云数据中的每个点定义平面标记,初始标记为0,即FlatIndex=0,表示该点未被标记;
步骤4.1.7,遍历点云数据中的每个点,若该点的平面标记为0,为当前点设置墙壁类别为FlatIndex;
步骤4.1.8,将该邻近点入栈,在KD_tree中搜索该点的K个邻近点;
步骤4.1.9,判断该点与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.10,若栈不为空,取栈顶元素为搜索点,在KD_tree中搜索该点的K个邻近点,判断栈顶元素与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.11,重复上述步骤直到标记完所有的点;
步骤4.1.12,遍历点云数据中的每一个点,根据每个点的平面类别FlatIndex组织点云为不同墙壁,将点云数目小于阈值的墙壁删除得到结果Walls,如图8所示;
步骤4.2,对分割的墙壁进行栅格化,分析各墙壁是否含有门窗,对有门窗的墙壁,将各栅格区别分析并拼接,对没有门窗的墙壁,使用OBB包围盒的简单方式进行重建;具体为:
步骤4.2.1,,确定栅格单元的尺寸,由于常见的墙壁长宽比例有2:1和3:2,且门窗等开口占墙壁空间较多,因此,检测开口时,墙壁栅格化尺寸可以偏大一些,兼顾精确度、算法复杂度以及门窗在墙壁中的规则性,本文对常用的墙壁长宽比求均值,将墙壁划分为35×20的栅格,栅格尺寸的长度Gl和宽度Gw分别为:Gl=Wl/35,Gw=Ww/20,其中,Wl为墙壁长度,Ww为墙壁宽度;
步骤4.2.2,确定栅格单元属性值,即单位面积的点云密度,首先,在三维空间中,设置墙壁的厚度为2μ,由墙壁对应的平面图中线段确定其中心位置坐标,从而确定每个栅格单元的三维坐标范围,统计该栅格单元中点的数目Numg之后,根据栅格尺寸对点云密度进行归一化,根据式(11)求得单位面积的点密度作为栅格单元的属性值,分类阈值Dw由栅格单元点密度均值确定,如式(12)所示;统计该栅格单元中点的数目Numg,根据式(11)求得单位面积的点密度Dg,并根据公式(12)计算分类阈值Dw:
其中,Dgi表示某栅格单元i的点密度;
步骤4.2.3,判断某个栅格单元i使得Dgi<Dw,则该墙壁为复杂墙壁,否则,当所有单元格的单位面积点密度均不小于阈值时,将该墙壁视为单一墙壁;
步骤4.2.4,对单一墙壁应用OBB包围盒进行墙壁重建;具体为:
步骤4.2.4.1,若该类墙壁为不包含门窗的单一墙壁,首先构造点云数据的协方差矩阵;
步骤4.2.4.2,采用奇异值分解法(SVD)求出协方差矩阵的特征值与特征向量,对其进行排序并归一化特征向量,随后将三个特征向量施密特正交化之后便可作为OBB的三个主轴;
步骤4.2.4.3,将全体点云数据投影到三个主轴上估算出OBB的中心点以及长宽高信息,求出该OBB盒的八个顶点,完成OBB包围盒的求解;
步骤4.2.5,对复杂墙壁,将点密度小于阈值Dw的单元格看作门窗开口位置,将开口位置从墙壁中删除,将点密度大于阈值Dw的栅格单元拼接起来,将墙壁拟合为包含门窗细节的结果,如图9所示;
步骤4.3,对全局结构进行重建,具体为:
步骤4.3.1,在步骤3获取到的室内场景房间布局信息的基础上,将房间布局的二维分割结果在竖直方向即z坐标轴上进行延伸,向z坐标轴负方向延伸至地面的上界坐标处,向z轴正方向延伸至天花板下界坐标处,将房间布局平面图中的每条直线看作墙壁的中心线,为其赋予合适的厚度2μ,μ为组距,即将z坐标平均分为200份,得到室内场景各个墙壁在三维空间中的粗略位置;
步骤4.3.2,将墙壁粗略位置与步骤4.2场景分割得到的墙壁分类重建结果相对应,并将墙壁的z坐标范围限定在地面与天花板位置之间,从而重建得到每个房间的不同墙壁结构,然后在墙壁两端分别补齐地面与天花板平面,如图10所示。
本发明以仅包含顶点坐标信息的点云室内场景布局重建进行研究,首先进行了场景预处理,分析了场景的趋势从而旋转至自然竖直状态,根据墙壁投影线获取房间初步划分结果,据此,通过转化为无向图的方式求解出室内场景房间布局信息,根据区域增长法分割场景,并栅格化分类重建墙壁,最后根据场景布局平面图规范墙壁重建位置,补齐地面天花板而得到场景重建结果。本发明技术路线较为完整,在有效性和稳定性都有了很大的提高,并且拥有较好的鲁棒性,丰富了计算机图形学和视觉智能的方法体系。
Claims (10)
1.一种基于结构分析的点云室内场景布局重建方法,其特征在于,具体按照如下步骤实施:
步骤1,对室内场景进行预处理操作;
步骤2,对步骤1进行预处理的室内场景进行粗分割,提取房间的轮廓;
步骤3,根据经步骤2提取的房间的轮廓对室内场景房间的布局进行划分,划分出房间布局;
步骤4,分割经步骤1预处理的室内场景的不同面的墙壁并进行重建,然后在步骤3划分的房间布局进行室内场景的全局结构重建,重建完整的室内结构。
2.根据权利要求1所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤1具体为:
步骤1.1,通过PCA主成分分析法获取室内场景主平面,由主平面的法向量确定旋转方向,将原始室内场景点云数据进行场景旋转操作,使其处于自然竖直状态;
步骤1.2,去掉经步骤1.1旋转的室内场景点云数据的地板与天花板数据;
步骤1.3,对经步骤1.2去除地板与天花板数据的室内场景点云数据使用PCA主成分分析法拟合局部平面来确定法向量,根据法向量滤掉非墙壁数据点。
3.根据权利要求2所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤1.1中将原始室内场景点云数据进行场景旋转操作具体按照如下步骤实施:
步骤1.1.1,室内场景点云数据为P={p1,p2,…,pn},其中,pi表示点云文件中第i个点的坐标信息(xi,yi,zi),将原始室内场景点云表示为n×3的矩阵A={XT,YT,ZT},其中,X={X1,X2,…,Xn},Y={Y1,Y2,…,Yn},Z={Z1,Z2,…,Zn},其中,Xi是第i个点的x轴坐标点,Yi是第i个点的y轴坐标点,Zi是第i个点的z轴坐标点,i=1,2,3...n;
步骤1.1.2,在矩阵A中每一位坐标值减去各自的平均值得到矩阵D;
步骤1.1.4,采用奇异值分解法求协方差矩阵的特征值与特征向量,对特征值按从大到小排序,选择前两个特征值v1和v2,对应的特征向量v1=(a1,b1,c1)和v2=(a2,b2,c2)即为所求得的场景主次趋势方向V;
V=v1×v2=(b1c2-c1b2,c1a2-a1c2,a1b2-b1a2) (3)
步骤1.1.5,通过旋转V,使得场景处于自然竖直状态,具体为:
V′为V在坐标平面XOY内的投影向量,与x轴之间的夹角为α,V″为向量V在坐标平面xoz内的投影向量,与z轴之间的夹角为β,将场景绕z坐标轴顺时针旋转角度α,确定第一个旋转矩阵为Rz,然后顺时针绕y轴旋转角度β,确定了第二个旋转矩阵为Ry,将点云场景转化为n×3的矩阵A后进行旋转,得到最终的旋转结果为Vres=(ARz)Ry=A(RzRy),其中,
4.根据权利要求3所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤1.2具体按照下述步骤进行:
步骤1.2.1,统计经步骤1.1旋转的室内场景点云数据在z坐标方向上的分布构造直方图,其中,组距的设置是根据点的z坐标范围,将其均匀划分为200份,将其间隔作为组距μ,然后根据组距确定每组数据的坐标范围及处于该范围的点的数目完成直方图的构造;
步骤1.2.2,将直方图中z坐标的某个极值设为预期峰值peak;
步骤1.2.3,找出以peak为中心,半径为radius的区域中出现的全体点集M,将其归属于一个聚类C,其中,radius为μ的四倍;
步骤1.2.4,以peak为中心点,计算从peak到集合M中每个数据点构成的向量,叠加所有向量,可得到漂移向量shift,peak=peak+shift,即peak沿着shift的方向移动,移动距离是||shift||;
步骤1.2.5,重复步骤1.2.3、1.2.4,直到shift小于设定的阈值σ,即就是迭代到收敛或达到了最大迭代次数10次,记录当前的peak,迭代中所有出现在以peak为中心radius为半径圆内的点都归到一类C中;
步骤1.2.6,如果收敛时,当前簇C的peak与其它已经存在的簇C1中心的距离小于阈值σ,那么把C和C1合并,数据点出现次数也对应合并,其中,σ为μ的二倍,否则,把C作为新的聚类;
步骤1.2.7,找到天花板和地面的聚类点集后,将其去除。
5.根据权利要求4所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤1.3具体按照如下步骤实施:
步骤1.3.1,对经步骤1.2去除地板与天花板数据的室内场景点云数据用KD-tree组织点云模型中的所有点;
步骤1.3.2,对于经步骤1.3.1组织的所有点中某个点p的法向量求解,首先利用KD-tree查询其k个近邻点;
步骤1.3.3,根据k个近邻点构造协方差矩阵;
步骤1.3.4,用奇异值分解算法求解步骤1.1.3中协方差矩阵最大的两特征向量的向量积,该向量积就是当前点的法向量;
步骤1.3.5,重复步骤1.3.2、1.3.3和1.3.4,求出所有点的法向量;
步骤1.3.6,遍历所有点,如果对于点pi,若|ni·nz|<Tangle,则删除点pi,其中,ni为点pi的法向量,nz=(0,0,1)是沿纵轴的单位向量,据此删除所有此类非墙壁点,其中,Tangle∈(0,2]。
6.根据权利要求1一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤2具体为:
步骤2.1,对经步骤1预处理的室内场景点云数据投影到坐标平面xoy中,然后将x,y坐标不同的点作为数据投影的有效点,得到原始室内场景结构轮廓,即就是投影数据;
步骤2.2,应用基于改进的Ransac算法,对步骤2.1投影得到的原始室内场景结构轮廓进行直线检测,具体按照下述步骤进行:
步骤2.2.1,创建直线集合H,用于存放直线模型;
步骤2.2.2,构造投影数据的二维KD-tree;
步骤2.2.3,求解未测点p的k近邻,选取点p与其中一个邻近点求解直线模型参数,并将该邻近点标记为已检测;
步骤2.2.4,在最大误差范围内统计符合该直线模型的数据点数目,如果符合模型参数的样本点,即内点所占的比例大于设定的阈值0.8,则记录该直线,否则利用所有内点重新确定直线方程参数,然后进入下一次迭代;
步骤2.2.5,重复2.2.2~2.2.4,直到p的邻近点均已被检测;
步骤2.2.6,计算该直线模型的方向向量,在检测结果H中查询共线或接近共线的直线集合,加入其中,否则,无共线直线时新建类别H′,加入H中,然后标记该直线的内点;
步骤2.2.7,重复2.2.2~2.2.6,直到所有点已被检测,得到检测结果H={H1,H2,…,Hn};
步骤2.3,采用基于Mean shift方法将直线检测结果进行聚类,去除较短直线以及杂乱直线的影响,最后得到与墙壁一一对应的直线检测结果,从而得到室内房间轮廓。
7.根据权利要求6所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤2.3具体为:
步骤2.3.1,选择经2.2步骤得到的直线分割结果簇Hi,在结果簇中随机选择一条作为起始中心直线,设定center=0,表示与该直线间的距离;
步骤2.3.2,找出与中心直线center距离小于radius的所有直线作为集合U,认为这些直线同属于一个聚类HCi,同时在该聚类中将该直线出现的次数加1;
步骤2.3.3,以center为基准,计算从center到集合U中每条直线的距离,为center其中一侧的直线距离设定为负值,将这些距离相加,得到漂移量distance;
步骤2.3.4,更新中心直线center=center+distance,即中心直线沿着distance的方向移动,移动距离是|distance|,根据center与直线结果簇Hi的直线方向确定新的中心直线,清空集合U,为下一次均值漂移做准备;
步骤2.3.5,重复上述2.3.2~2.3.4步骤,直到漂移量distance的值小于阈值μ/5,即为迭代到收敛,记录此时的center,将此迭代中经过的所有直线均归类为簇Hi;
步骤2.3.6,重复2.3.1~2.3.5直到归类完所有的结果簇H;
步骤2.3.7,如果簇Hi的center与其它已经存在的簇Hj的center之间距离小于墙壁厚度2μ,那么把Hj和Hi合并,直线出现次数也对应合并;
步骤2.3.8,分类:根据每个类,对每个点的访问频率,取访问频率最大的m个类,每个类别以其中心直线作为直线优化结果,得到的直线集合L={l1,l2,…,ln}作为潜在墙壁优化结果,即就是室内房间轮廓,其中,m值根据数据中房间数目预先设定。
8.根据权利要求1所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤3具体为:
步骤3.1,由公式(5)确定扫描仪位置P,Min.Ptx为x坐标最小值,Min.Pty为y坐标最小值,Max.Pty为y坐标最小值,Min.Ptz为z坐标最小值,Max.Ptz为z坐标最大值;
P=(Min.Ptx,(Min.Pty+Max.Pty)/2,(Min.Ptz+Max.Ptz)/2) (5)
步骤3.2,求出步骤2房间粗分割结果中的每一个预期房间元素的右侧边界中心点;
步骤3.3,以点P为起始点,向右延伸穿过预期房间右侧边界中心点构造射线Ray;
步骤3.4,求解射线Ray与该房间的交点数目Number,若交点数目Number为偶数,将该房间标签化为1,表示为潜在房间,否则,将该房间标签化为0;
步骤3.5,将每个预期房间元素转化为无向图G中的一个顶点,对于能相互连接的预期房间添加边信息,其中边是具有权值的,边的权值通过与之连接的两个预期房间的中心点的欧氏距离定义,其中,房间中心点通过房间的顶点求均值得到,具体如公式(5)所示:
x,y,z为房间中心点的三维坐标,而xi,yi,zi是二维平面图中房间的顶点Pi的三维坐标;
步骤3.6,确定了每个房间的中心点之后,通过公式(7)求解房间中心点O和Q的欧氏距离,从而得到房间中心点及其欧氏距离的求解结果;
其中,xO,yO,zO是房间中心点O的三维坐标,xQ,yQ,zQ是房间中心点Q的三维坐标;
步骤3.7,向该无向图G中再添加了2个顶点,分别为源点“S”和汇点“T”,设置与源点相连的顶点权值为0,汇点相连的顶点权值为无穷大;
步骤3.8,建立无向图G的能量方程,如公式(8)所示;
E=∑i∈VDi(li)+Ri,j(li,lj) (8)
其中,Di为数据项,Ri,j为约束项,数据项由式(9)确定,约束项由式(10)确定;
其中,li为房间i的标签,系数c=0.6,Si为归一化后的面积,即预期房间面积除以整个预期房间的初步划分结果的面积;
Li,j为预期房间i、j的欧氏距离,该距离除以整个预期房间划分结果的最长边可做归一化处理,权重ωi,j为该预期房间占总房间初步划分结果总面积的比例;
步骤3.9,使用最大流/最小割算法对其进行能量最小化求解,找到从源点S到汇点T的一条增广路径,这条增广路径作为无向图G的最小割,最终将最小割映射到房间粗分割结果中,从而得到房间布局信息。
9.根据权利要求1所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤4具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体为:
步骤4.1,对经过步骤1处理后的数据,应用区域增长算法分割出不同面的墙壁,具体按照下述步骤进行:
步骤4.1.1,为经步骤1处理后的室内场景点云数据构造KD_tree;
步骤4.1.2,在KD_tree中搜索当前点的K个邻近点;
步骤4.1.3,根据当前点以及其K个近邻点构造协方差矩阵;
步骤4.1.4,采用奇异值分解算法进行特征分解,最大的两特征向量的向量积即为当前点的法向量;
步骤4.1.5,重复上述4.1.2~4.1.4步骤,直到求出所有点的法向量;
步骤4.1.6,为点云数据中的每个点定义平面标记,初始标记为0,即FlatIndex=0,表示该点未被标记;
步骤4.1.7,遍历点云数据中的每个点,若该点的平面标记为0,为当前点设置墙壁类别为FlatIndex;
步骤4.1.8,将该邻近点入栈,在KD_tree中搜索该点的K个邻近点;
步骤4.1.9,判断该点与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.10,若栈不为空,取栈顶元素为搜索点,在KD_tree中搜索该点的K个邻近点,判断栈顶元素与其K个邻近点法向量夹角的余弦值cos,若1-cos<0.1,则标记该邻近点墙壁类别为FlatIndex,将该邻近点入栈;
步骤4.1.11,重复上述步骤直到标记完所有的点;
步骤4.1.12,遍历点云数据中的每一个点,根据每个点的平面类别FlatIndex组织点云为不同墙壁,将点云数目小于阈值的墙壁删除得到结果Walls;
步骤4.2,对分割的墙壁进行栅格化,分析各墙壁是否含有门窗,对有门窗的墙壁,将各栅格区别分析并拼接,对没有门窗的墙壁,使用OBB包围盒的简单方式进行重建;
步骤4.3,对全局结构进行重建。
10.根据权利要求9所述的一种基于结构分析的点云室内场景布局重建方法,其特征在于,所述步骤4.2具体为:
步骤4.2.1,将墙壁划分为35×20的栅格,栅格尺寸的长度Gl和宽度Gw分别为:Gl=Wl/35,Gw=Ww/20,其中,Wl为墙壁长度,Ww为墙壁宽度;
步骤4.2.2,统计该栅格单元中点的数目Numg,根据式(11)求得单位面积的点密度Dg,并根据公式(12)计算分类阈值Dw:
其中,Dgi表示某栅格单元i的点密度;
步骤4.2.3,判断某个栅格单元i使得Dgi<Dw,则该墙壁为复杂墙壁,否则,该墙壁视为单一墙壁;
步骤4.2.4,对单一墙壁应用OBB包围盒进行墙壁重建;
步骤4.2.5,对复杂墙壁,将点密度小于阈值Dw的单元格看作门窗开口位置,将开口位置从墙壁中删除,将点密度大于阈值Dw的栅格单元拼接起来,将墙壁拟合为包含门窗细节的结果;
所述步骤4.3具体为:
步骤4.3.1,在步骤3获取到的室内场景房间布局信息的基础上,将房间布局的二维分割结果在竖直方向即z坐标轴上进行延伸,向z坐标轴负方向延伸至地面的上界坐标处,向z轴正方向延伸至天花板下界坐标处,将房间布局平面图中的每条直线看作墙壁的中心线,为其赋予合适的厚度2μ,μ为组距,即将z坐标平均分为200份,得到室内场景各个墙壁在三维空间中的粗略位置;
步骤4.3.2,将墙壁粗略位置与步骤4.2场景分割得到的墙壁分类重建结果相对应,并将墙壁的z坐标范围限定在地面与天花板位置之间,从而重建得到每个房间的不同墙壁结构,然后在墙壁两端分别补齐地面与天花板平面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010703972.2A CN111986322B (zh) | 2020-07-21 | 2020-07-21 | 一种基于结构分析的点云室内场景布局重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010703972.2A CN111986322B (zh) | 2020-07-21 | 2020-07-21 | 一种基于结构分析的点云室内场景布局重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111986322A true CN111986322A (zh) | 2020-11-24 |
CN111986322B CN111986322B (zh) | 2023-12-22 |
Family
ID=73437905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010703972.2A Active CN111986322B (zh) | 2020-07-21 | 2020-07-21 | 一种基于结构分析的点云室内场景布局重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111986322B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112380312A (zh) * | 2020-11-30 | 2021-02-19 | 重庆智行者信息科技有限公司 | 基于栅格检测的激光地图更新方法、终端及计算机设备 |
CN112833815A (zh) * | 2020-12-30 | 2021-05-25 | 广东工业大学 | 一种基于激光的房屋测量和墙面平整度检测方法 |
CN113112600A (zh) * | 2021-04-02 | 2021-07-13 | 浙江大学 | 基于结构的室内场景三维建模方法 |
CN113284228A (zh) * | 2021-05-27 | 2021-08-20 | 西安理工大学 | 一种基于点云的室内场景房间布局划分方法 |
CN113379898A (zh) * | 2021-06-17 | 2021-09-10 | 西安理工大学 | 一种基于语义分割的三维室内场景重建方法 |
CN113409448A (zh) * | 2021-06-17 | 2021-09-17 | 西安理工大学 | 一种室内点云场景的墙壁结构分析与重建方法 |
CN114596313A (zh) * | 2022-05-09 | 2022-06-07 | 深圳大学 | 一种基于室内点云的建筑物构件损伤检测方法及相关设备 |
JP2022142994A (ja) * | 2021-03-17 | 2022-10-03 | 株式会社構造計画研究所 | モデル生成装置、モデル生成方法、及びモデル生成プログラム |
CN118365746A (zh) * | 2024-04-29 | 2024-07-19 | 山东科技大学 | 三维点云语义引导的室内场景平面图生成方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107146280A (zh) * | 2017-05-09 | 2017-09-08 | 西安理工大学 | 一种基于切分的点云建筑物重建方法 |
CN108875133A (zh) * | 2017-05-09 | 2018-11-23 | 达索系统公司 | 确定建筑布局 |
CN109325998A (zh) * | 2018-10-08 | 2019-02-12 | 香港理工大学 | 一种基于点云数据的室内3d建模方法、系统及相关装置 |
CN110782524A (zh) * | 2019-10-25 | 2020-02-11 | 重庆邮电大学 | 基于全景图的室内三维重建方法 |
-
2020
- 2020-07-21 CN CN202010703972.2A patent/CN111986322B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107146280A (zh) * | 2017-05-09 | 2017-09-08 | 西安理工大学 | 一种基于切分的点云建筑物重建方法 |
CN108875133A (zh) * | 2017-05-09 | 2018-11-23 | 达索系统公司 | 确定建筑布局 |
CN109325998A (zh) * | 2018-10-08 | 2019-02-12 | 香港理工大学 | 一种基于点云数据的室内3d建模方法、系统及相关装置 |
US20200111251A1 (en) * | 2018-10-08 | 2020-04-09 | The Hong Kong Polytechnic University | 3d indoor modeling method, system and device based on point cloud data |
CN110782524A (zh) * | 2019-10-25 | 2020-02-11 | 重庆邮电大学 | 基于全景图的室内三维重建方法 |
Non-Patent Citations (1)
Title |
---|
张宏鑫;方雨桐;利明;: "结合视觉惯性模组的室内三维布局鲁棒重建方法", 计算机辅助设计与图形学学报, no. 02 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112380312A (zh) * | 2020-11-30 | 2021-02-19 | 重庆智行者信息科技有限公司 | 基于栅格检测的激光地图更新方法、终端及计算机设备 |
CN112380312B (zh) * | 2020-11-30 | 2022-08-05 | 北京智行者科技股份有限公司 | 基于栅格检测的激光地图更新方法、终端及计算机设备 |
CN112833815A (zh) * | 2020-12-30 | 2021-05-25 | 广东工业大学 | 一种基于激光的房屋测量和墙面平整度检测方法 |
JP2022142994A (ja) * | 2021-03-17 | 2022-10-03 | 株式会社構造計画研究所 | モデル生成装置、モデル生成方法、及びモデル生成プログラム |
JP7164652B2 (ja) | 2021-03-17 | 2022-11-01 | 株式会社構造計画研究所 | モデル生成装置、モデル生成方法、及びモデル生成プログラム |
CN113112600B (zh) * | 2021-04-02 | 2023-03-03 | 浙江大学 | 基于结构的室内场景三维建模方法 |
CN113112600A (zh) * | 2021-04-02 | 2021-07-13 | 浙江大学 | 基于结构的室内场景三维建模方法 |
CN113284228A (zh) * | 2021-05-27 | 2021-08-20 | 西安理工大学 | 一种基于点云的室内场景房间布局划分方法 |
CN113284228B (zh) * | 2021-05-27 | 2023-04-14 | 西安理工大学 | 一种基于点云的室内场景房间布局划分方法 |
CN113379898A (zh) * | 2021-06-17 | 2021-09-10 | 西安理工大学 | 一种基于语义分割的三维室内场景重建方法 |
CN113409448A (zh) * | 2021-06-17 | 2021-09-17 | 西安理工大学 | 一种室内点云场景的墙壁结构分析与重建方法 |
CN113379898B (zh) * | 2021-06-17 | 2022-11-11 | 西安理工大学 | 一种基于语义分割的三维室内场景重建方法 |
CN114596313A (zh) * | 2022-05-09 | 2022-06-07 | 深圳大学 | 一种基于室内点云的建筑物构件损伤检测方法及相关设备 |
CN114596313B (zh) * | 2022-05-09 | 2022-10-04 | 深圳大学 | 一种基于室内点云的建筑物构件损伤检测方法及相关设备 |
CN118365746A (zh) * | 2024-04-29 | 2024-07-19 | 山东科技大学 | 三维点云语义引导的室内场景平面图生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111986322B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111986322B (zh) | 一种基于结构分析的点云室内场景布局重建方法 | |
Xu et al. | Voxel-and graph-based point cloud segmentation of 3d scenes using perceptual grouping laws | |
Zhou et al. | 2.5 d dual contouring: A robust approach to creating building models from aerial lidar point clouds | |
CN109446691B (zh) | 基于激光点云与空气动力学的活立木抗风性能分析方法 | |
CN111696210A (zh) | 一种基于三维点云数据特征轻量化的点云重构方法及系统 | |
CN112288857A (zh) | 一种基于深度学习的机器人语义地图物体识别方法 | |
CN114266987B (zh) | 一种无人机高边坡危岩体智能识别方法 | |
CN113379898B (zh) | 一种基于语义分割的三维室内场景重建方法 | |
Zhou et al. | Decomposing polygon meshes by means of critical points | |
CN106127857A (zh) | 综合数据驱动与模型驱动的机载LiDAR数据建模方法 | |
CN110047036B (zh) | 基于极坐标格网的地面激光扫描数据建筑物立面提取方法 | |
CN113066162A (zh) | 一种用于电磁计算的城市环境快速建模方法 | |
CN109215112B (zh) | 一种单侧点云模型的标注方法 | |
CN114926699A (zh) | 一种室内三维点云语义分类方法、装置、介质及终端 | |
Han et al. | Urban scene LOD vectorized modeling from photogrammetry meshes | |
Hu et al. | Geometric feature enhanced line segment extraction from large-scale point clouds with hierarchical topological optimization | |
Yu et al. | Unsupervised ground filtering of airborne-based 3D meshes using a robust cloth simulation | |
CN114463338A (zh) | 一种基于图割与后处理的建筑物激光脚点自动提取方法 | |
CN117788735A (zh) | 基于栅格划分的动态点云去除方法 | |
Zhu et al. | Variational building modeling from urban MVS meshes | |
Sun et al. | Automated segmentation of LiDAR point clouds for building rooftop extraction | |
Kaushik et al. | Fast planar clustering and polygon extraction from noisy range images acquired in indoor environments | |
CN112200248A (zh) | 一种基于dbscan聚类的城市道路环境下的点云语义分割方法、系统及存储介质 | |
Sahebdivani et al. | Deep learning based classification of color point cloud for 3D reconstruction of interior elements of buildings | |
Zhu et al. | Feature line based building detection and reconstruction from oblique airborne imagery |
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 |