一种基于背包式三维激光点云数据的建筑物立面三维重建方法技术

本发明专利技术公开了一种基于背包式三维激光点云数据的建筑物立面三维重建方法，涉及地理信息技术领域。包括以下步骤：S1.建筑物点云数据获取；S2.建筑物立面点云数据自动提取；S3.单体建筑物自动分割；S4.建筑物立面几何位置边界获取；S5.建筑物立面三维重建。本发明专利技术采用基于体元投影密度的点云滤波算法在保留较为完整的建筑物目标的同时有效滤除地面和植被等非建筑物目标，然后利用图像全局搜索和剖面分析的方法实现单体建筑物的自动分割；针对单体建筑物点云利用RANSAC算法进行立面自动分割和冗余立面剔除获取建筑物立面几何位置边界；利用二维边界线约束原始点云数据并结合RANSAC算法进行立面三维边界直线拟合以获取建筑物立面几何框架模型。

A 3D reconstruction method of building facade based on knapsack 3D laser point cloud data

【技术实现步骤摘要】
一种基于背包式三维激光点云数据的建筑物立面三维重建方法
本专利技术涉及地理信息
，具体是涉及一种基于背包式三维激光点云数据的建筑物立面三维重建方法。
技术介绍
随着“智慧城市”等概念的逐渐深入，城市空间信息已经成为了当今社会中不可或缺的重要资源。建筑物作为城市最重要的组成部分，基于遥感数据的实现建筑物立面自动化三维重建的方法，将为“数字城市”提供重要的基础数据，同时也是目前研究领域的难点之一。目前的建模软件将建筑物视为刚体结构，通过在点云数据中量测建筑物部件的几何参数，生成类似的几何基元，拼凑出建筑物的模型。这种方法生成的模型过于理想化，与建筑物实际模型差别较大。且模型生成过程完全基于手动操作，效率极低。因此，移动式三维激光扫描系统在大规模城市场景的三维重建中具有越来越明显的优势。移动式三维激光扫描系统具有数据获取速度快、数据相对完整、精度高、点云数据密集和适应复杂场景的优点，但另一方面，其获取的数据量巨大，带有噪声并存在遮挡，这给点云数据的处理方法带来了巨大的挑战。公开号为CN107146280A的中国专利技术“一种基于切分的点云建筑物的重建方法”提供了一种基于切分的点云建筑物的重建方法，具体按照以下步骤实施：步骤1、利用区域增长以及基于距离的聚类算法对点云场景进行分割，并提取点云场景中构成建筑物的平面；步骤2、对步骤1提取的平面进行纵向切分，计算每个切片的长度，并依次比较每个切片与后续切片的高度差探测关键点，并依次连接关键点，完成平面的重建。步骤3、根据平面间的拓扑关系对建筑物初始重建模型进行优化；步骤4、利用墙面与窗户的拓扑关系，完成建筑物细节的重建。此专利技术解决了现有技术中存在的由于三维物体点云数据量大以及点云数据不完整，导致建筑重建细节不足的问题。公开号为CN108765568A的中国专利技术“一种基于激光雷达点云的多层次建筑物快速三维重建方法”提供了一种基于激光雷达点云的多层次建筑物快速三维重建方法，其步骤为：采集三维点云数据；点云滤波处理得到建筑物点云；使用优化的随机抽样一致性算法获得激光点云集；Delaunay三角剖分进行建筑物轮廓点的提取；冒泡排序对轮廓点进行排序，生成轮廓线；使用关键点提取算法，提取建筑物轮廓的关键点；连接关键点，并使用正交约束对轮廓线进行规则化处理；将建筑物轮廓线赋予点云中的高程信息，生成三维建筑物模型。此专利技术采用优化的算法能够提供更高的执行效率，并能有效的删除建筑物墙面点；使用Delaunay剖分方法精确的提取建筑物的轮廓线；使用轮廓线关键点提取算法较好的提取建筑物的关键点，并能抑制错误轮廓点的干扰，减少伪关键点的生成。随着卫星导航定位、惯性导航、激光扫描、近景摄像等多传感器集成技术与同步定位与制图(SLAM)技术的快速发展，背包式移动三维激光扫描系统应运而生。作为目前最热门的测绘新技术，设备通过人员背载进行数据扫描，人员可通过的地方都可进行数据采集，采集过程快速、便捷、低成本，为建筑物轮廓线的自动提取提供了全新的解决方案。背包式移动三维激光扫描系统作为一种新兴的测绘设备，其在国内外的应用还处于起步阶段，相关的应用研究相对较少，因此本专利技术旨在针对背包式三维激光点云数据的特点，提供一种行之有效的建筑物立面三维几何重建方法，实现建筑物立面的自动化高精度三维重建。
技术实现思路
本专利技术要解决的技术问题是提供一种实现建筑物立面的自动化高精度三维重建的基于背包式三维激光点云数据的建筑物立面三维重建方法。为解决上述问题，本专利技术的技术方案为：一种基于背包式三维激光点云数据的建筑物立面三维重建方法，包括以下步骤：S1.建筑物点云数据获取：通过背包式三维激光扫描系统获取建筑物点云数据；S2.建筑物立面点云数据自动提取：首先对建筑物点云数据进行体元化，然后对每个体元内的点云数据计算体元投影密度，设置阈值，提取体元投影密度大于阈值的该体元内的点云数据为建筑物立面点云数据；S3.单体建筑物自动分割：对建筑物立面点云数据进行栅格化，记录非空栅格坐标并构建k-d树，按照图像搜索的方法对点云数据进行分类，然后针对每个类别点云数据进行剖面分析，设置阈值，将剖面分析参数小于阈值的类别判别为单体建筑物点云数据；S4.建筑物立面几何位置边界获取：利用随机抽样一致性算法对单体建筑物点云数据进行立面分割，得到立面方程参数和立面点云数据，根据立面方程参数分析立面夹角和距离以合并冗余立面，然后利用主成分分析法分析立面点云数据的几何分布参数并拟合二维直线，通过直线拓扑求交进行直线规则化处理，获取建筑物立面几何位置边界；S5.建筑物立面三维重建:基于建筑物立面几何位置边界对步骤S1中获取的建筑物点云数据进行建筑物立面点云数据分割，利用随机抽样一致性算法拟合立面的参数方程，基于参数方程对立面点云数据进行垂直投影，然后利用凹凸包算法获取立面的三维边界数据，利用栅格化的方法获取立面的上边界数据，利用随机抽样一致性算法拟合立面上边界三维直线，对拟合的直线进行分段和排序，结合建筑物点云数据和直线求交对建筑物的上边界线进行规则化处理，立面上边界特征点向地面作垂直投影，投影点的平面坐标和上边界特征点的平面坐标相同，投影点的高程值为建筑物点云数据中的最小高程，最终连接上边界特征点和地面投影点获取建筑物立面三维几何框架模型。进一步地，所述步骤S2中的体元投影密度计算方法，包括以下步骤：S21.以分辨率dl对建筑物点云数据进行体元化和栅格化，计算体元投影密度，采用如下公式：Dp＝nvoxel/ngrid其中，nvoxel为非空体元数，ngrid为非空栅格数，Dp为体元投影密度；S22.设置阈值，将体元投影密度Dp小于阈值的点云数据去除。进一步地，所述步骤S3中图像搜索包括以下步骤：随机选取一个栅格作为起始搜索点，利用k-d树搜索其邻域半径内是否存在栅格坐标，若存在则将本次搜索到的邻域数据作为新的搜索点，若所有搜索点的邻域半径内都不存在未访问过的栅格坐标，则将已访问过的栅格坐标内部的点云数据作为同一类别，从剩余栅格中随机选择一个作为新的起始搜索点，重复上述过程，直至所有的数据都被访问。进一步地，所述步骤S3中剖面分析包括以下步骤：利用凸包算法对每个类别的栅格数据进行处理，获取对应的外边界，计算外边界栅格的面积与利用外边界栅格构造的多边形面积的比值作为剖面分析参数。进一步地，所述步骤S4中合并冗余立面包括以下步骤：S41.立面分割结果中每个立面的参数方程如下：Ax+By+Cz+D＝0立面水平投影直线与X轴夹角的计算公式如下：θ＝arctan(-A/B)S42.设置阈值，若两平面的夹角的差值和参数D的差值均小于对应的阈值，则判别这两个平面为同一个平面，并将二者对应的点云数据进行合并。进一步地，所述步骤S4中利用主成分分析法拟合直线包括以下步骤：根据立面点云数据XY坐标的第一主成分对应的单位特征向量e1＝[e11e12]来拟合直线的斜率k，计算公本文档来自技高网...

【技术保护点】
1.一种基于背包式三维激光点云数据的建筑物立面三维重建方法，其特征在于：包括以下步骤：/nS1.建筑物点云数据获取：通过背包式三维激光扫描系统获取建筑物点云数据；/nS2.建筑物立面点云数据自动提取：首先对建筑物点云数据进行体元化，然后对每个体元内的点云数据计算体元投影密度，设置阈值，提取体元投影密度大于阈值的该体元内的点云数据为建筑物立面点云数据；/nS3.单体建筑物自动分割：对建筑物立面点云数据进行栅格化，记录非空栅格坐标并构建k-d树，按照图像搜索的方法对点云数据进行分类，然后针对每个类别点云数据进行剖面分析，设置阈值，将剖面分析参数小于阈值的类别判别为单体建筑物点云数据；/nS4.建筑物立面几何位置边界获取：利用随机抽样一致性算法对单体建筑物点云数据进行立面分割，得到立面方程参数和立面点云数据，根据立面方程参数分析立面夹角和距离以合并冗余立面，然后利用主成分分析法分析立面点云数据的几何分布参数并拟合二维直线，通过直线拓扑求交进行直线规则化处理，获取建筑物立面几何位置边界；/nS5.建筑物立面三维重建:基于建筑物立面几何位置边界对步骤S1中获取的建筑物点云数据进行建筑物立面点云数据分割，利用随机抽样一致性算法拟合立面的参数方程，基于参数方程对立面点云数据进行垂直投影，然后利用凹凸包算法获取立面的三维边界数据，利用栅格化的方法获取立面的上边界数据，利用随机抽样一致性算法拟合立面上边界三维直线，对拟合的直线进行分段和排序，结合建筑物点云数据和直线求交对建筑物的上边界线进行规则化处理，立面上边界特征点向地面作垂直投影，投影点的平面坐标和上边界特征点的平面坐标相同，投影点的高程值为建筑物点云数据中的最小高程，最终连接上边界特征点和地面投影点获取建筑物立面三维几何框架模型。/n...

【技术特征摘要】
1.一种基于背包式三维激光点云数据的建筑物立面三维重建方法，其特征在于：包括以下步骤：
S1.建筑物点云数据获取：通过背包式三维激光扫描系统获取建筑物点云数据；
S2.建筑物立面点云数据自动提取：首先对建筑物点云数据进行体元化，然后对每个体元内的点云数据计算体元投影密度，设置阈值，提取体元投影密度大于阈值的该体元内的点云数据为建筑物立面点云数据；
S3.单体建筑物自动分割：对建筑物立面点云数据进行栅格化，记录非空栅格坐标并构建k-d树，按照图像搜索的方法对点云数据进行分类，然后针对每个类别点云数据进行剖面分析，设置阈值，将剖面分析参数小于阈值的类别判别为单体建筑物点云数据；
S4.建筑物立面几何位置边界获取：利用随机抽样一致性算法对单体建筑物点云数据进行立面分割，得到立面方程参数和立面点云数据，根据立面方程参数分析立面夹角和距离以合并冗余立面，然后利用主成分分析法分析立面点云数据的几何分布参数并拟合二维直线，通过直线拓扑求交进行直线规则化处理，获取建筑物立面几何位置边界；
S5.建筑物立面三维重建:基于建筑物立面几何位置边界对步骤S1中获取的建筑物点云数据进行建筑物立面点云数据分割，利用随机抽样一致性算法拟合立面的参数方程，基于参数方程对立面点云数据进行垂直投影，然后利用凹凸包算法获取立面的三维边界数据，利用栅格化的方法获取立面的上边界数据，利用随机抽样一致性算法拟合立面上边界三维直线，对拟合的直线进行分段和排序，结合建筑物点云数据和直线求交对建筑物的上边界线进行规则化处理，立面上边界特征点向地面作垂直投影，投影点的平面坐标和上边界特征点的平面坐标相同，投影点的高程值为建筑物点云数据中的最小高程，最终连接上边界特征点和地面投影点获取建筑物立面三维几何框架模型。


2.根据权利要求1所述的一种基于背包式三维激光点云数据的建筑物立面三维重建方法，其特征在于：所述步骤S2中的体元投影密度计算方法，包括以下步骤：
S21.以分辨率dl对建筑物点云数据进行体元化和栅格化，计算体元投影密度，采用如下公式：
Dp＝...

【专利技术属性】
技术研发人员：邢万里，蔡东健，周旺辉，岳顺，
申请(专利权)人：苏州工业园区测绘地理信息有限公司，
类型：发明
国别省市：江苏;32