一种三维点云提取方法及提取系统技术方案

本申请提供了一种三维点云提取方法及提取系统，提取方法包括以下步骤：获取初始点云数据；根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理；基于CSF算法对下采样后的部分初始点云数据进行地面滤波；对滤波后的点云数据构建三维空间特征和二维平面特征，并将三维空间特征和二维平面特征划分为训练集和测试集，利用训练集创建支持向量机，采用支持向量机提取隧道点云。本申请具有准确性高、鲁棒性好以及处理速度快等优势，还可应用在铁路装备检测领域，能够准确剔除铁路罐车罐体点云中梯子、液位计、标靶等附件点云，保留罐壁数据，提升后续容积计算准确性，为铁路罐车智能检定检测奠定基础。检定检测奠定基础。检定检测奠定基础。

【技术实现步骤摘要】
一种三维点云提取方法及提取系统


[0001]本申请属于隧道工程测量及铁路装备检测领域，具体涉及一种基于多特征的三维点云提取方法及提取系统。

技术介绍

[0002]在隧道施工过程中，需要经常对不同时期的隧道进行断面质量检测。采用全站仪或激光断面仪的传统检测方法存在效率低、手段复杂、检测不全面等问题。随着三维扫描技术的发展，利用扫描仪获得隧道内海量三维点云数据，借助高效的点云处理算法，可以快速、准确、轻松地完成隧道断面质量检测任务。
[0003]如图1所示，三维激光扫描仪获取的初始点云数据包括隧道点云和非隧道点云，其中，非隧道点云种包括地面、车辆和人员等。因此采取高效的点云滤波算法剔除非隧道点云是提取干净、高质量隧道点云的关键。常用的点云滤波方法包括：(1)基于人工选择的滤波方法，该方法通过人机交互剔除非目标点云，不适合噪声点分布散乱的情况；(2)基于距离阈值的滤波方法，该方法通过设定具体阈值划分目标点云和非目标点云；因此该方法通用性差，且容易出现误判；(3)基于统计的滤波方法，该方法通过统计每个点与K近邻范围内各点之间的距离剔除离群点，但该方法无法剔除距离目标点云较近的噪声点；(4)基于随机采样一致性(Random Sample Consensus，RANSAC)算法的滤波方法通过对不同平面进行拟合实现点云分割，该方法对参数依赖性较强，且针对特定物体平面。

技术实现思路

[0004]为至少在一定程度上克服相关技术中存在的问题，本申请提供了一种三维点云提取方法及提取系统。
[0005]根据本申请实施例的第一方面，本申请提供了一种三维点云提取方法，其包括以下步骤：
[0006]获取初始点云数据；
[0007]根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理；
[0008]基于CSF算法对下采样后的部分初始点云数据进行地面滤波；
[0009]对滤波后的点云数据构建三维空间特征和二维平面特征，并将三维空间特征和二维平面特征划分为训练集和测试集，利用训练集创建支持向量机，采用支持向量机提取隧道点云。
[0010]上述三维点云提取方法中，所述获取初始点云数据的方式为：
[0011]使用全站式扫描仪对隧道进行扫描，直接获取绝对坐标点云数据；
[0012]或者，使用地面式三维激光扫描仪与标靶结合的方式对隧道进行扫描，首先间接获取相对坐标点云数据，然后通过坐标配准转换为绝对坐标点云数据。
[0013]上述三维点云提取方法中，所述基于CSF算法对下采样后的部分初始点云数据进
行地面滤波的过程为：
[0014]对下采样后的部分初始点云数据进行过滤，滤除离群点，倒置点云；
[0015]初始布栅格化，根据实际点云密度确定栅格大小，初始位置位于点云纵坐标的最大值上；
[0016]将点云与栅格粒子投影至同一水平面，计算当前粒子能抵达的最低位置；
[0017]计算每个粒子在仅考虑重力时外力影响下移动的距离；
[0018]计算每个粒子在内力影响下移动的距离；
[0019]不断计算每个粒子在仅考虑重力时外力影响下移动的距离以及在内力影响下移动的距离，直到所有粒子的最大高度变化小于预设阈值或者迭代次数到达指定次数，计算停止；
[0020]计算栅格化粒子与原始点云的距离；
[0021]判断栅格化粒子与原始点云的距离是否小于预设的距离阈值，如果是，则判定当前点云为地面点云；否则，判定当前点云为非地面点云；
[0022]判断是否遍历完所有点云，如果是，则结束；否则，返回重新计算栅格化粒子与原始点云的距离。
[0023]进一步地，所述粒子在仅考虑重力时外力影响下移动的距离为：
[0024][0025]式中，X(t+Δt)表示粒子在外力影响下移动的距离，X(t)表示粒子在t时刻的位置，X(t
‑
Δt)表示粒子在t
‑
Δt时刻的位置，Δt表示时间间隔，G表示重力常数，m表示粒子质量，通常为1。
[0026]进一步地，所述粒子在内力影响下移动的距离为：
[0027][0028]式中，表示粒子的位移向量，当粒子为可移动粒子时，b＝1；否则，b＝0；表示当前粒子，表示的相邻粒子，表示垂直方向的归一化向量，
[0029][0030]进一步地，所述采用支持向量机提取隧道点云的过程为：
[0031]构建隧道点云特征向量；
[0032]根据构建的隧道点云特征向量，将滤波后的点云数据的三维空间特征和二维平面特征划分为训练集和测试集；
[0033]利用训练集创建支持向量机，获取最优参数下的支持向量机模型；
[0034]利用训练好的支持向量机模型，将测试集作为输入，对隧道点云进行二分类判别；如果输出结果为1，则判定当前输入为隧道点云，保留当前输入的测试集中的点云数据；如果输出结果为
‑
1，则表示当前输入为非隧道点云，丢弃当前输入的测试集中的点云数据。
[0035]更进一步地，所述构建隧道点云特征向量的过程为：
[0036]在点云数据上构造三维KD
‑
树，以将离散的三维点云转换为有关联的拓扑结构；
[0037]在KD
‑
树中寻找每个点的邻近点云簇，点云簇的个数根据实际的点云密度进行确定；
[0038]计算邻近点云簇内所有点云的协方差特征值，构建三维空间特征和二维平面特征。
[0039]更进一步地，所述三维空间特征包括发散性指数S
λ
、平面性指数P
λ
、线性指数L
λ
、各向同性O
λ
、各向异性A
λ
、特征熵E
λ
、特征和∑λ、曲率变化C
λ
、垂直度V、点云簇z轴极差R
z
、点云簇z轴方差和点云簇密度D；
[0040]发散性指数S
λ
为：
[0041]平面性指数P
λ
为：
[0042]线性指数L
λ
为：
[0043]各向同性O
λ
为：
[0044]各向异性A
λ
为：
[0045]特征熵E
λ
为：
[0046]特征和∑λ为：∑λ＝λ1+λ2+λ3；
[0047]曲率变化C
λ
为：
[0048]垂直度V为：V＝1
‑
n
z
；
[0049]点云簇z轴极差R
z
为：R
z
＝Max(z)
‑
Min(z)；
[0050]点云簇z轴方差为：
[0051]点云簇密度D为：
[0052]上述各式中，λ1,λ2,λ3表示协方差特征值，n
z
表示每个点法向量本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种三维点云提取方法，其特征在于，包括以下步骤：获取初始点云数据；根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理；基于CSF算法对下采样后的部分初始点云数据进行地面滤波；对滤波后的点云数据构建三维空间特征和二维平面特征，并将三维空间特征和二维平面特征划分为训练集和测试集，利用训练集创建支持向量机，采用支持向量机提取隧道点云。2.根据权利要求1所述的三维点云提取方法，其特征在于，所述获取初始点云数据的方式为：使用全站式扫描仪对隧道进行扫描，直接获取绝对坐标点云数据；或者，使用地面式三维激光扫描仪与标靶结合的方式对隧道进行扫描，首先间接获取相对坐标点云数据，然后通过坐标配准转换为绝对坐标点云数据。3.根据权利要求1所述的三维点云提取方法，其特征在于，所述基于CSF算法对下采样后的部分初始点云数据进行地面滤波的过程为：对下采样后的部分初始点云数据进行过滤，滤除离群点，倒置点云；初始布栅格化，根据实际点云密度确定栅格大小，初始位置位于点云纵坐标的最大值上；将点云与栅格粒子投影至同一水平面，计算当前粒子能抵达的最低位置；计算每个粒子在仅考虑重力时外力影响下移动的距离；计算每个粒子在内力影响下移动的距离；不断计算每个粒子在仅考虑重力时外力影响下移动的距离以及在内力影响下移动的距离，直到所有粒子的最大高度变化小于预设阈值或者迭代次数到达指定次数，计算停止；计算栅格化粒子与原始点云的距离；判断栅格化粒子与原始点云的距离是否小于预设的距离阈值，如果是，则判定当前点云为地面点云；否则，判定当前点云为非地面点云；判断是否遍历完所有点云，如果是，则结束；否则，返回重新计算栅格化粒子与原始点云的距离。4.根据权利要求3所述的三维点云提取方法，其特征在于，所述粒子在仅考虑重力时外力影响下移动的距离为：式中，X(t+Δt)表示粒子在外力影响下移动的距离，X(t)表示粒子在t时刻的位置，X(t
‑
Δt)表示粒子在t
‑
Δt时刻的位置，Δt表示时间间隔，G表示重力常数，m表示粒子质量，通常为1。5.根据权利要求3所述的三维点云提取方法，其特征在于，所述粒子在内力影响下移动的距离为：
式中，表示粒子的位移向量，当粒子为可移动粒子时，b＝1；否则，b＝0；表示当前粒子，表示的相邻粒子，表示垂直方向的归一化向量，6.根据权利要求3所述的三维点云提取方法，其特征在于，所述采用支持向量机提取隧道点云的过程为：构建隧道点云特征向量；根据构建的隧道点云特征向量，将滤波后的点云数据的三维空间特征和二维平面特征划分为训练集和测试集；利用训练集创建支持向量机，获取最优参数下的支持向量机模型；利用训练好的支持向量机模型，将测试集作为输入，对隧道点云进行二分类判别；如果输出结果为1，则判定当前输入为隧道点云，保留当前输入的测试集中的点云数据；如果输出结果为
‑
1，则表示当前输入为非隧道点云，丢弃当前输入的测试集中的点云数据。7.根据权利要求6所述的三维点云提取方法，其特征在于，所述构建隧道点云特征向量的过程为：在点云数据上构造三维KD
‑
树，以将离散的三维点云转换为有关联的拓扑结构；在KD
‑
树中寻找每个点的邻近点云簇，点云簇的个数根据实际的点云密度进行确定；计算邻近点云簇内所有点云的协方差特征值，构建三维空间特征和二维平面特征。8.根据权利要求7所述的三维点云提取方法，其特征在于，所述三维空间特征包括发散性指数S
λ
、平面性指...

【专利技术属性】
技术研发人员：刘云锴，庞庆，张志鹏，邵学君，武慧杰，郭有为，李祉萱，
申请(专利权)人：中国铁道科学研究院集团有限公司，
类型：发明
国别省市：