CN108648277A - 一种激光雷达点云数据的快速重建方法 - Google Patents
一种激光雷达点云数据的快速重建方法 Download PDFInfo
- Publication number
- CN108648277A CN108648277A CN201810376549.9A CN201810376549A CN108648277A CN 108648277 A CN108648277 A CN 108648277A CN 201810376549 A CN201810376549 A CN 201810376549A CN 108648277 A CN108648277 A CN 108648277A
- Authority
- CN
- China
- Prior art keywords
- point cloud
- point
- cloud data
- laser radar
- line
- 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
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/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Image Generation (AREA)
Abstract
本发明公开了一种激光雷达点云数据的快速重建方法,首先获取激光雷达发送的距离数据和二维振镜发送的角度数据,转换为三维点云坐标,并进行格式化操作;然后对点云数据进行聚类标记和拟合,以剔除离群点;接着遍历相邻扫描线寻找基准线,并以两基准线为一组进行三角网格化;最后结合三角网格和设定的物体模型光照参数进行上色渲染,完成物体重建。本发明避免了传统方法的空洞效应,并且具有速度快、精度高和鲁棒性好等突出优势。
Description
技术领域
本发明涉及一种点云数据处理方法,特别涉及一种激光雷达点云数据的快速重建方法。
背景技术
三维物体重建是指通过传感器获得空间物体的二维图像,进而将其恢复成物体的三维结构的过程。获取二维图像的传感器可以选择激光雷达,其具有扫描精确度高、速度快的优势。但是通过激光雷达扫描获得的被测对象的点云数量庞大,且分布密度不均匀,这给三维重建带来了很多的挑战。
现有技术通过对点云数据进行三角化来重建物体的面,但是大部分已知的三角化方式都存在局限性。专利201710216985.5使用贪婪投影三角化算法对三维点云数据进行三角化网格处理,首先选取一个样本三角形作为初始曲面,然后不断选择新点,根据新点来构造新的三角形,扩张曲面边界,形成一个完整的三角网格曲面,最后根据投影点云的连接关系确定原始点云间的拓扑关系。此方法可以有效的找出散乱的数据点云之间蕴含的拓扑学关系,但是由于每个三角形面都要单独生成,所以需要不断计算比对,消耗大量时间,速度较慢。专利201610101813.9以扫描线对点云数据按照点云矩阵的方式保存,对点云矩阵的m×n数据进行划分,得到(m-1)×(n-1)个小的点云矩阵,划分后的每个点云矩阵包含四个点云数据,采用连接对角线方式对其划分,完成三角网格构建。此方法有效的利用了点云矩阵之间的联系,操作简单,重建速度快,但是无法识别点云矩阵中出现的激光空洞点进行正确的三角化,此外此方法对于点云数据的质量要求很高。
发明内容
本发明的目的在于提供一种激光雷达点云数据的快速重建方法,克服现有方法进行二维振镜激光雷达点云数据重建时遇到的效率低、重建效果差、无法解决空洞等问题。
实现本发明目的的技术解决方案为:一种激光雷达点云数据的快速重建方法,包括如下步骤:
步骤1、点云数据格式化:获取激光雷达发送的距离数据和二维振镜发送的角度数据,将其转换为三维点云坐标,并对点云数据进行格式化操作;
步骤2、数据聚类处理:根据数据所在的扫描线位置,进行点云数据的聚类标记;
步骤3、离群点剔除与滤波:拟合聚类标记后的点云数据,根据拟合曲线进行离群点剔除和平滑处理;
步骤4、非空洞点并行基准线规划:遍历相邻扫描线寻找相邻两行同时存在非空洞点的列,将两点连线作为基准线;
步骤5、基准线紧邻点网格化:以两基准线为一组,计算基准线上点与周围散乱点之间的距离,根据紧邻原则,选择距离最短的点云数据进行三角网格化;
步骤6、网格曲面渲染:结合三角网格和设定的物体模型光照参数进行上色渲染,完成物体重建。
本发明与现有技术相比,其显著优点为:1)本发明通过遍历相邻扫描线寻找出上下相邻两行同时存在非空洞点的列,并将两点连线作为基准线,以两基准线为一组实现点云数据分组,这种分组构建方式有效避免了传统方法面临的空洞点问题,使重建后的物体表面更加平滑、整齐,物体边缘更加清晰;2)本发明基于二维振镜激光雷达蛇形扫描以及电机驱动振镜角度的特性,按照扫描线位置聚类的方式,确定了同一条扫描线上的点云数据处于同一平面,并依此对点云数据进行聚类标记,在聚类标记的基础上可对同一分组的数据进行剔除离群点与平滑数据等相关处理,改善了网格化重建后整个曲面的平整性;3)本发明选取并行的基准线进行规划,遍历相邻扫描线寻找出上下相邻两行同时存在非空洞点的列,存储列号的数组重复使用,且在对各个分组进行三角化的时候根据紧邻原则,选择距离最短的点云数据进行三角网格化,简化了三角网格数据判别算法,大幅度提高了重建效率,且能保证点云网格面的平滑。
附图说明
图1为本发明的方法步骤示意图;
图2为本发明基于的二维振镜激光雷达工作原理示意图;
图3为本发明的基于并行基准线的分组示意图;
图4为本发明的分组三角化处理空洞点示意图;
图5为本发明的可能出现于首尾的散乱数据点网格化示意图。
具体实施方式
下面结合附图和具体实施例,进一步说明本发明方案。
如图1所示,激光雷达点云数据的快速重建方法,包括如下步骤:
步骤1、点云数据格式化:获取激光雷达发送的距离数据和二维振镜发送的角度数据,将其转换为三维点云坐标,并对点云数据进行格式化操作。
本发明实施例中,使用串口连接或网口连接的方式接收激光雷达和二维振镜发送的数据,利用如下矩阵关系将其转换为三维点云坐标:
式中,表示点云的三维坐标,L表示激光雷达发送的距离数据,θx和θy分别表示二维振镜发送的水平和竖直方向偏转角。
格式化后的数据结构中的点云数据类建成如下情况:
转化成三维点云坐标后,按照扫描线行列关系存储点云数据参数,包括点云数据的三维坐标、探测点横向纵向的偏转角和探测点位于原扫描线的行列。
步骤2、数据聚类处理:根据数据所在的扫描线位置进行点云数据的聚类标记。基于二维振镜激光雷达的蛇形扫描特性,可知同一条扫描线上的点云处于同一平面上,根据扫描线位置可以将同一条扫描线上点云数据的三维坐标转化成扫描线平面上的二维坐标。
步骤3、离群点剔除与滤波:拟合聚类标记后的点云数据,根据拟合曲线进行离群点剔除和平滑处理。
本发明实施例采用四阶最小二乘法对点云数据进行拟合,采用空间均值滤波对聚类的数据进行平滑处理。四阶最小二乘法具体如下:假设所在扫描平面上扫描线的四阶拟合多项式为:
y=a0+a1x+a2x2+a3x3+a4x4
每点到这条曲线的距离平方之和,即偏差平方和为:
对等式右边求ai偏导,得到:
求解方程组得到各系数,可以确定扫描线,进而计算每个点相对于曲线的偏离值,结合动态阈值确定是否为离群点,其中阈值并非固定值,阈值大小为该点所在扫描线上距此点最近的八个点的偏离值的平均值。
步骤4、非空洞点并行基准线规划:遍历相邻扫描线寻找相邻两行同时存在非空洞点的列,并将两点连线作为基准线。
本发明实施例中相邻扫描线的遍历过程为:如图3所示,从第一行扫描线开始,依次遍历相邻两条扫描线,记录下所有有效数据点的列号,同时用一个数组记录上下相邻两行同时存在非空洞点的列,像这样上下两条扫描线同一列都存在的两个点称作基准点,此两点连成的线称为基准线,以两基准线为一组实现点云数据分组。
步骤5、基准线紧邻点网格化:以两基准线为一组,计算基准线上点与周围散乱点之间的距离,根据紧邻原则,选择距离最短的点云数据进行三角网格化。
本发明实施例中对各个并行基准线组进行三角化的具体过程为:如图4所示,构造每组第一个三角形时先选取第一条基准线上两点,第三个点选取第i条或第i+1条扫描线上下一个紧邻的有效数据点,然后去除和紧邻点同一条扫描线上靠前的那个点,以剩下两个点为基础继续寻找下一个紧邻的有效数据点,直到第二条基准线上的两点也被构造进三角形为止。
如图5所示,若第一条基准线之前或者位于最后一条基准线之后存在散乱点,对散乱点单独进行网格化的方法与之前相同,区别为终止条件改为全部点都参与构造成三角形即停止。
步骤6、网格曲面渲染:结合三角网格和设定的物体模型光照参数进行上色渲染,完成物体重建,具体先计算每个三角网格面的法向量,然后据此计算物体表面法向量与观察视角向量的夹角,接着结合此夹角与设定的物体模型环境光、镜面反射光、漫反射光计算当前的光照条件,最后进行上色渲染,完成物体重建。
Claims (8)
1.一种激光雷达点云数据的快速重建方法,其特征在于,包括如下步骤:
步骤1、点云数据格式化:获取激光雷达发送的距离数据和二维振镜发送的角度数据,将其转换为三维点云坐标,并对点云数据进行格式化操作;
步骤2、数据聚类处理:根据数据所在的扫描线位置,进行点云数据的聚类标记;
步骤3、离群点剔除与滤波:拟合聚类标记后的点云数据,根据拟合曲线进行离群点剔除和平滑处理;
步骤4、非空洞点并行基准线规划:遍历相邻扫描线寻找相邻两行同时存在非空洞点的列,将两点连线作为基准线;
步骤5、基准线紧邻点网格化:以两基准线为一组,计算基准线上点与周围散乱点之间的距离,根据紧邻原则,选择距离最短的点云数据进行三角网格化;
步骤6、网格曲面渲染:结合三角网格和设定的物体模型光照参数进行上色渲染,完成物体重建。
2.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤1中三维点云坐标与距离、角度的矩阵关系为:
式中,表示点云的三维坐标,L表示激光雷达发送的距离数据,θx和θy分别表示二维振镜发送的水平和竖直方向偏转角。
3.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤1的格式化操具体为:按照扫描线行列关系存储点云数据参数,包括点云数据的三维坐标、探测点横向纵向的偏转角和探测点位于原扫描线的行列。
4.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤3采用四阶最小二乘法对聚类标记后的点云数据进行拟合,以剔除离群点,具体方法如下:
假设所在扫描平面上扫描线的四阶拟合多项式为:
y=a0+a1x+a2x2+a3x3+a4x4
每点到这条曲线的距离平方之和,即偏差平方和为:
对等式右边求ai偏导,得到:
求解方程组得到各系数即确定扫描线,进而计算每个点相对于曲线的偏离值,结合动态阈值确定是否为离群点,其中动态阈值大小为该点所在扫描线上距此点最近的八个点的偏离值的平均值。
5.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤3采用空间均值滤波对聚类的数据进行平滑处理。
6.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤4规划基准线的具体方法为:从第一行扫描线开始,依次遍历相邻两条扫描线,记录所有有效数据点的列号,筛选相邻两行同时存在非空洞点的列,选择其对应的两个点为基准点,此两点的连线为基准线。
7.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤5对基准线组进行三角化的具体过程为:构造每组第一个三角形时先选取第一条基准线上两点,第三个点选取第i条或第i+1条扫描线上下一个紧邻的有效数据点,然后去除和紧邻点同一条扫描线上靠前的那个点,以剩下两个点为基础继续寻找下一个紧邻的有效数据点,直到第二条基准线上的两点也被构造进三角形或者全部点都参与构造三角形为止。
8.根据权利要求书1所述的激光雷达点云数据的快速重建方法,其特征在于,步骤6首先计算每个三角网格面的法向量,然后据此计算物体表面法向量与观察视角向量的夹角,接着结合此夹角与设定的物体模型环境光、镜面反射光、漫反射光计算当前的光照条件,最后进行上色渲染,完成物体重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376549.9A CN108648277B (zh) | 2018-04-25 | 2018-04-25 | 一种激光雷达点云数据的快速重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376549.9A CN108648277B (zh) | 2018-04-25 | 2018-04-25 | 一种激光雷达点云数据的快速重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108648277A true CN108648277A (zh) | 2018-10-12 |
CN108648277B CN108648277B (zh) | 2021-11-09 |
Family
ID=63747658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810376549.9A Active CN108648277B (zh) | 2018-04-25 | 2018-04-25 | 一种激光雷达点云数据的快速重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108648277B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109345592A (zh) * | 2018-11-01 | 2019-02-15 | 中国矿业大学(北京) | 基于探地雷达的地下空洞三维坐标提取算法 |
CN109584294A (zh) * | 2018-11-23 | 2019-04-05 | 武汉中海庭数据技术有限公司 | 一种基于激光点云的路面点云提取方法和装置 |
CN111912346A (zh) * | 2020-06-30 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 一种适用于飞机表面机器人钻铆系统的窝孔在线检测方法 |
CN112506959A (zh) * | 2020-11-30 | 2021-03-16 | 智慧航海(青岛)科技有限公司 | 智能船舶数据库检索的数据调度方法、装置和检索系统 |
CN113291847A (zh) * | 2021-03-31 | 2021-08-24 | 湖南千盟工业智能系统股份有限公司 | 一种基于三维成像的散料智能堆取料方法 |
CN114119710A (zh) * | 2021-11-23 | 2022-03-01 | 燕山大学 | 一种计算敞车车厢剩余冻煤体积的方法及系统 |
CN117687027A (zh) * | 2023-12-08 | 2024-03-12 | 中国科学院空天信息创新研究院 | 一种面向机载阵列干涉sar三维点云多径抑制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170053438A1 (en) * | 2014-06-13 | 2017-02-23 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Method and system for reconstructing a three-dimensional model of point clouds |
CN107123162A (zh) * | 2016-02-24 | 2017-09-01 | 中国科学院沈阳自动化研究所 | 基于二维激光传感器的三维环境表面三角网格构建方法 |
-
2018
- 2018-04-25 CN CN201810376549.9A patent/CN108648277B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170053438A1 (en) * | 2014-06-13 | 2017-02-23 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Method and system for reconstructing a three-dimensional model of point clouds |
CN107123162A (zh) * | 2016-02-24 | 2017-09-01 | 中国科学院沈阳自动化研究所 | 基于二维激光传感器的三维环境表面三角网格构建方法 |
Non-Patent Citations (3)
Title |
---|
ANGELO L D等: "A new mesh-growing algorithm for fast surface", 《COMPUTER-AIDED DESIGN》 * |
FLOATER M S等: "Meshless parameterization and surface reconstruction", 《COMPUTER AIDED GEOMETRIC DESIGN》 * |
高莉: "改进的Delaunay三角剖分算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109345592A (zh) * | 2018-11-01 | 2019-02-15 | 中国矿业大学(北京) | 基于探地雷达的地下空洞三维坐标提取算法 |
CN109345592B (zh) * | 2018-11-01 | 2021-09-21 | 中国矿业大学(北京) | 基于探地雷达的地下空洞三维坐标提取算法 |
CN109584294A (zh) * | 2018-11-23 | 2019-04-05 | 武汉中海庭数据技术有限公司 | 一种基于激光点云的路面点云提取方法和装置 |
CN109584294B (zh) * | 2018-11-23 | 2020-08-28 | 武汉中海庭数据技术有限公司 | 一种基于激光点云的路面点云提取方法和装置 |
CN111912346A (zh) * | 2020-06-30 | 2020-11-10 | 成都飞机工业(集团)有限责任公司 | 一种适用于飞机表面机器人钻铆系统的窝孔在线检测方法 |
CN111912346B (zh) * | 2020-06-30 | 2021-12-10 | 成都飞机工业(集团)有限责任公司 | 一种适用于飞机表面机器人钻铆系统的窝孔在线检测方法 |
CN112506959A (zh) * | 2020-11-30 | 2021-03-16 | 智慧航海(青岛)科技有限公司 | 智能船舶数据库检索的数据调度方法、装置和检索系统 |
CN112506959B (zh) * | 2020-11-30 | 2022-11-08 | 智慧航海(青岛)科技有限公司 | 智能船舶数据库检索的数据调度方法、装置和检索系统 |
CN113291847A (zh) * | 2021-03-31 | 2021-08-24 | 湖南千盟工业智能系统股份有限公司 | 一种基于三维成像的散料智能堆取料方法 |
CN114119710A (zh) * | 2021-11-23 | 2022-03-01 | 燕山大学 | 一种计算敞车车厢剩余冻煤体积的方法及系统 |
CN114119710B (zh) * | 2021-11-23 | 2024-05-07 | 燕山大学 | 一种计算敞车车厢剩余冻煤体积的方法及系统 |
CN117687027A (zh) * | 2023-12-08 | 2024-03-12 | 中国科学院空天信息创新研究院 | 一种面向机载阵列干涉sar三维点云多径抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108648277B (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108648277A (zh) | 一种激光雷达点云数据的快速重建方法 | |
CN107025685B (zh) | 拓扑感知下的机载建筑屋顶点云建模方法 | |
CN105205866B (zh) | 基于密集点云的城市三维模型快速构建方法 | |
WO2020134082A1 (zh) | 一种路径规划方法、装置和移动设备 | |
CN110230998B (zh) | 基于线激光和双目相机的快速精密三维测量方法和装置 | |
CN108803621B (zh) | 一种钢轨焊缝修磨机器人路径规划方法 | |
CN104966317B (zh) | 一种基于矿体轮廓线的三维自动建模方法 | |
Zhu et al. | Seed point selection method for triangle constrained image matching propagation | |
CN102682103A (zh) | 一种面向海量激光雷达点云模型的三维空间索引方法 | |
CN111581776B (zh) | 一种基于几何重建模型的等几何分析方法 | |
CN105303616A (zh) | 基于单张照片的浮雕建模方法 | |
CN113868476B (zh) | 一种基于局部密度的八叉树点云预处理方法 | |
CN105046735A (zh) | 一种基于基元分布的三维表面纹理合成方法 | |
CN107146279A (zh) | 一种基于共生相关函数的多孔介质三维建模方法 | |
CN109544672A (zh) | 一种三维建筑模型纹理映射方法及装置 | |
CN115661374B (zh) | 一种基于空间划分和模型体素化的快速检索方法 | |
CN113421338A (zh) | 一种将二维矢量进行三维渲染的方法 | |
CN107247926A (zh) | 一种人体检测方法及装置 | |
CN114898043A (zh) | 一种激光点云数据瓦片构建方法 | |
CN110415362B (zh) | 重叠点云融合方法 | |
CN112465991B (zh) | 一种隧道点云去噪和生成可视化模型的方法 | |
CN110176053A (zh) | 一种大规模实景三维整体匀色方法 | |
CN110738726B (zh) | 一种基于八叉树的机器人视觉引导三维物体重建方法 | |
CN116659460A (zh) | 一种公路横断面激光点云切片快速生成方法 | |
CN111028349A (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 |