CN116681653A - 一种三维点云提取方法及提取系统 - Google Patents
一种三维点云提取方法及提取系统 Download PDFInfo
- Publication number
- CN116681653A CN116681653A CN202310546415.8A CN202310546415A CN116681653A CN 116681653 A CN116681653 A CN 116681653A CN 202310546415 A CN202310546415 A CN 202310546415A CN 116681653 A CN116681653 A CN 116681653A
- Authority
- CN
- China
- Prior art keywords
- point cloud
- dimensional
- cloud data
- tunnel
- particle
- 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.)
- Pending
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000012706 support-vector machine Methods 0.000 claims abstract description 47
- 238000001914 filtration Methods 0.000 claims abstract description 34
- 238000012360 testing method Methods 0.000 claims abstract description 31
- 238000012549 training Methods 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 239000002245 particle Substances 0.000 claims description 75
- 239000013598 vector Substances 0.000 claims description 42
- 230000008569 process Effects 0.000 claims description 22
- 238000000926 separation method Methods 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 14
- 239000007788 liquid Substances 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 9
- 239000004744 fabric Substances 0.000 description 8
- 238000010276 construction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 125000003275 alpha amino acid group Chemical group 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000009423 ventilation Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- 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/10024—Color image
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- 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/20—Special algorithmic details
- G06T2207/20081—Training; Learning
Abstract
本申请提供了一种三维点云提取方法及提取系统,提取方法包括以下步骤:获取初始点云数据;根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理;基于CSF算法对下采样后的部分初始点云数据进行地面滤波;对滤波后的点云数据构建三维空间特征和二维平面特征,并将三维空间特征和二维平面特征划分为训练集和测试集,利用训练集创建支持向量机,采用支持向量机提取隧道点云。本申请具有准确性高、鲁棒性好以及处理速度快等优势,还可应用在铁路装备检测领域,能够准确剔除铁路罐车罐体点云中梯子、液位计、标靶等附件点云,保留罐壁数据,提升后续容积计算准确性,为铁路罐车智能检定检测奠定基础。
Description
技术领域
本申请属于隧道工程测量及铁路装备检测领域,具体涉及一种基于多特征的三维点云提取方法及提取系统。
背景技术
在隧道施工过程中,需要经常对不同时期的隧道进行断面质量检测。采用全站仪或激光断面仪的传统检测方法存在效率低、手段复杂、检测不全面等问题。随着三维扫描技术的发展,利用扫描仪获得隧道内海量三维点云数据,借助高效的点云处理算法,可以快速、准确、轻松地完成隧道断面质量检测任务。
如图1所示,三维激光扫描仪获取的初始点云数据包括隧道点云和非隧道点云,其中,非隧道点云种包括地面、车辆和人员等。因此采取高效的点云滤波算法剔除非隧道点云是提取干净、高质量隧道点云的关键。常用的点云滤波方法包括:(1)基于人工选择的滤波方法,该方法通过人机交互剔除非目标点云,不适合噪声点分布散乱的情况;(2)基于距离阈值的滤波方法,该方法通过设定具体阈值划分目标点云和非目标点云;因此该方法通用性差,且容易出现误判;(3)基于统计的滤波方法,该方法通过统计每个点与K近邻范围内各点之间的距离剔除离群点,但该方法无法剔除距离目标点云较近的噪声点;(4)基于随机采样一致性(Random Sample Consensus,RANSAC)算法的滤波方法通过对不同平面进行拟合实现点云分割,该方法对参数依赖性较强,且针对特定物体平面。
发明内容
为至少在一定程度上克服相关技术中存在的问题,本申请提供了一种三维点云提取方法及提取系统。
根据本申请实施例的第一方面,本申请提供了一种三维点云提取方法,其包括以下步骤:
获取初始点云数据;
根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理;
基于CSF算法对下采样后的部分初始点云数据进行地面滤波;
对滤波后的点云数据构建三维空间特征和二维平面特征,并将三维空间特征和二维平面特征划分为训练集和测试集,利用训练集创建支持向量机,采用支持向量机提取隧道点云。
上述三维点云提取方法中,所述获取初始点云数据的方式为:
使用全站式扫描仪对隧道进行扫描,直接获取绝对坐标点云数据;
或者,使用地面式三维激光扫描仪与标靶结合的方式对隧道进行扫描,首先间接获取相对坐标点云数据,然后通过坐标配准转换为绝对坐标点云数据。
上述三维点云提取方法中,所述基于CSF算法对下采样后的部分初始点云数据进行地面滤波的过程为:
对下采样后的部分初始点云数据进行过滤,滤除离群点,倒置点云;
初始布栅格化,根据实际点云密度确定栅格大小,初始位置位于点云纵坐标的最大值上;
将点云与栅格粒子投影至同一水平面,计算当前粒子能抵达的最低位置;
计算每个粒子在仅考虑重力时外力影响下移动的距离;
计算每个粒子在内力影响下移动的距离;
不断计算每个粒子在仅考虑重力时外力影响下移动的距离以及在内力影响下移动的距离,直到所有粒子的最大高度变化小于预设阈值或者迭代次数到达指定次数,计算停止;
计算栅格化粒子与原始点云的距离;
判断栅格化粒子与原始点云的距离是否小于预设的距离阈值,如果是,则判定当前点云为地面点云;否则,判定当前点云为非地面点云;
判断是否遍历完所有点云,如果是,则结束;否则,返回重新计算栅格化粒子与原始点云的距离。
进一步地,所述粒子在仅考虑重力时外力影响下移动的距离为:
式中,X(t+Δt)表示粒子在外力影响下移动的距离,X(t)表示粒子在t时刻的位置,X(t-Δt)表示粒子在t-Δt时刻的位置,Δt表示时间间隔,G表示重力常数,m表示粒子质量,通常为1。
进一步地,所述粒子在内力影响下移动的距离为:
式中,表示粒子的位移向量,当粒子为可移动粒子时,b=1;否则,b=0;/>表示当前粒子,/>表示/>的相邻粒子,/>表示垂直方向的归一化向量,
进一步地,所述采用支持向量机提取隧道点云的过程为:
构建隧道点云特征向量;
根据构建的隧道点云特征向量,将滤波后的点云数据的三维空间特征和二维平面特征划分为训练集和测试集;
利用训练集创建支持向量机,获取最优参数下的支持向量机模型;
利用训练好的支持向量机模型,将测试集作为输入,对隧道点云进行二分类判别;如果输出结果为1,则判定当前输入为隧道点云,保留当前输入的测试集中的点云数据;如果输出结果为-1,则表示当前输入为非隧道点云,丢弃当前输入的测试集中的点云数据。
更进一步地,所述构建隧道点云特征向量的过程为:
在点云数据上构造三维KD-树,以将离散的三维点云转换为有关联的拓扑结构;
在KD-树中寻找每个点的邻近点云簇,点云簇的个数根据实际的点云密度进行确定;
计算邻近点云簇内所有点云的协方差特征值,构建三维空间特征和二维平面特征。
更进一步地,所述三维空间特征包括发散性指数Sλ、平面性指数Pλ、线性指数Lλ、各向同性Oλ、各向异性Aλ、特征熵Eλ、特征和∑λ、曲率变化Cλ、垂直度V、点云簇z轴极差Rz、点云簇z轴方差和点云簇密度D;
发散性指数Sλ为:
平面性指数Pλ为:
线性指数Lλ为:
各向同性Oλ为:
各向异性Aλ为:
特征熵Eλ为:
特征和∑λ为:∑λ=λ1+λ2+λ3;
曲率变化Cλ为:
垂直度V为:V=1-nz;
点云簇z轴极差Rz为:Rz=Max(z)-Min(z);
点云簇z轴方差为:/>
点云簇密度D为:
上述各式中,λ1,λ2,λ3表示协方差特征值,nz表示每个点法向量的z轴分量,k表示使用KD-树近邻搜索点云簇的k值,r表示搜索半径;
二维平面特征包括隧道垂直横截面内每个点与轨顶之间的距离dis(i′,j);
式中,表示实际轨顶位置/>的坐标,(x(i′,j),z(i′,j))表示第i′个横截面下第j个点。
更进一步地,所述利用训练集创建支持向量机,获取最优参数下的支持向量机模型的过程为:
定义分离超平面H(x),当H(x)=0时,表明样本点位于分离超平面上;当H(x)>0时,表明样本类别为y=1;当H(x)<0时,表明样本类别为y=-1;
H(x)的表达式为:
H(x)=wTx+b;
式中,w表示参数向量,b表示平移向量;
构建最优分离超平面的表达式为:
式中,i5=1,2,…,m;
对最优分离超平面的表达式进行转换并引入拉格朗日乘子,化解后最终优化目标函数为:
式中,α表示向量,i5和j均表示序号;
定义核函数和惩罚系数C>0,最终优化目标函数为:
式中,j=1,2,…,m;
利用SMO算法求解出向量α,并获得对应的参数向量w:
根据KKT条件找出所有的S个支持向量,计算每个支持向量的平均值b:
确定最优超平面,分类决策函数为:
f(x)=sign(wTx+b)。
根据本申请实施例的第二方面,本申请基于上述权利要求1~9任一项所述的三维点云提取方法还提供了一种三维点云提取系统,其包括扫描模块、下采样模块、滤波模块和提取模块;
所述扫描模块,用于获取初始点云数据;
所述下采样模块,用于根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理;
所述滤波模块,基于CSF算法对下采样后的部分初始点云数据进行地面滤波,以对隧道内的地面进行识别和剔除;
所述提取模块,用于根据滤波后的点云数据得到训练集和测试集,并利用训练集创建支持向量机,采用支持向量机提取隧道点云。
根据本申请的上述具体实施方式可知,至少具有以下有益效果:本申请提供的三维点云提取方法通过采用CSF算法对隧道内获取的初始点云进行地面滤波,随后将保留的离散三维点云转换为有关联的拓扑结构,构建多种三维空间特征与二维平面特征,以用于扩展点云数据,最后通过训练最优支持向量机模型快速、准确地提取出隧道点云,具有准确性高、鲁棒性好以及处理速度快等优势,本申请还可应用至铁路装备检测领域,可以准确剔除铁路罐车罐体点云中梯子、液位计、标靶等附件点云,保留罐壁数据,提升后续容积计算准确性,为铁路罐车智能检定检测奠定基础。
应了解的是,上述一般描述及以下具体实施方式仅为示例性及阐释性的,其并不能限制本申请所欲主张的范围。
附图说明
下面的所附附图是本申请的说明书的一部分,其示出了本申请的实施例,所附附图与说明书的描述一起用来说明本申请的原理。
图1为现有技术中隧道初始点云数据的示意图。
图2为本申请具体实施方式提供的一种三维点云提取方法的流程图。
图3为本申请具体实施方式提供的一种三维点云提取方法中获取的初始点云数据示意图。
图4为本申请具体实施方式提供的一种三维点云提取方法中基于CSF算法对下采样后的部分初始点云数据进行地面滤波的流程图。
图5为本申请具体实施方式提供的一种三维点云提取方法中采用支持向量机提取隧道点云的流程图。
图6为本申请具体实施方式提供的一种三维点云提取方法中线性可分支持向量机图解。
图7为本申请具体实施方式提供的一种三维点云提取方法中基于CSF算法对点云数据中的地面进行滤波的结果示意图;其中,图(a)表示初始点云,图(b)表示按里程切分后的点云,图(c)表示下采样后的点云,图(d)表示采用CSF算法剔除地面后的点云。
图8为本申请具体实施方式提供的一种三维点云提取方法中不同测试集隧道点云提取结果示意图。
图9为本申请具体实施方式提供的一种三维点云提取方法在铁路装备检测中的应用示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚明白,下面将以附图及详细叙述清楚说明本申请所揭示内容的精神,任何所属技术领域技术人员在了解本申请内容的实施例后,当可由本申请内容所教示的技术,加以改变及修饰,其并不脱离本申请内容的精神与范围。
本申请的示意性实施例及其说明用于解释本申请,但并不作为对本申请的限定。另外,在附图及实施方式中所使用相同或类似标号的元件/构件是用来代表相同或类似部分。
关于本文中所使用的“第一”、“第二”、…等,并非特别指称次序或顺位的意思,也非用以限定本申请,其仅为了区别以相同技术用语描述的元件或操作。
关于本文中所使用的“包含”、“包括”、“具有”、“含有”等等,均为开放性的用语,即意指包含但不限于。
关于本文中所使用的“及/或”,包括所述事物的任一或全部组合。
关于本文中的“多个”包括“两个”及“两个以上”;关于本文中的“多组”包括“两组”及“两组以上”。
某些用以描述本申请的用词将于下或在此说明书的别处讨论,以提供本领域技术人员在有关本申请的描述上额外的引导。
由于隧道现场施工环境较为复杂,除扫描到的隧道点云数据外,还包括以下点云数据:(1)地面;(2)正在施工的台车、除渣车和喷浆车等;(3)现场的施工人员;(4)连通洞内的通风管、排水管;(5)隧道爆破后残存的水汽、尘土等;因此需要采用快速、准确的滤波方式将以上噪声点云剔除,保留隧道点云。
三维扫描仪在隧道内获取的初始点云包含着大量数据,其中只有隧道点云才是后续隧道质量检测使用的正确数据,因此有必要对隧道三维点云数据进行提取。
如图2所示,本申请提供的三维点云提取方法包括以下步骤:
S1、获取初始点云数据;
为保证获取的点云数据为真实世界的大地坐标,可以采用以下两种方式获取点云数据:
一种方式是使用全站式扫描仪对隧道进行扫描,直接获取绝对坐标点云数据。采用这种方式,人员操作较为简单,但是扫描时间较长。
另一种方式是使用地面式三维激光扫描仪与标靶结合的方式对隧道进行扫描,首先间接获取相对坐标点云数据,然后通过坐标配准转换为绝对坐标点云数据。采用这种方式,人员操作较为复杂,但是扫描时间较短。
在实际操作中,可以根据现场的既有设备和工序流程,选择合适的方式,以获取初始点云数据。获取的初始点云数据如图3所示。
需要说明的是,真实世界的大地坐标是指CGCS2000国家大地坐标系下的三维坐标,实际工程计算中要按照大地坐标进行计算。
S2、根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理。
需要说明的是,里程范围是根据现场需求确定的,例如,里程范围可以是100m~200m。
对待计算的部分初始点云数据进行下采样,过滤掉冗余点云,保留局部范围内代表性点云数据,能够降低点云数据的数量,减少计算量,提升后续点云数据处理速度,进一步保证后续隧道三维点云质量检测的准确性。
S3、基于CSF(Cloth Simulation Filter,布模拟滤波)算法对下采样后的部分初始点云数据进行地面滤波,以对隧道内特征较为连贯的地面进行识别和剔除。
CSF算法是假设将一块布铺至原始倒置点云上,通过分析布上各个节点与相应点云之间的关系,最终迭代出布的形态,该布的形态即为地面形态。根据该形态即可划分出地面与非地面。
通过分析三维空间内各粒子受到的作用力(例如,重力和碰撞力组成的外力、内部牵引的内力)情况,定义如下基本公式:
式(1)中,Fext(X,t)表示粒子在时刻t位置X时受到的外力,Fint(X,t)表示粒子在时刻t位置X时受到的内力,m表示粒子质量(通常为1)。
假设粒子相对于地面只有垂直方向的运动,初始状态时布铺在倒置的点云上方。随着受外力影响,布的粒子会下沉,出现与点云相交的情况。判断哪些布粒子与点云相交,通过判断某处粒子的垂直高度和当前点云的高度是否都低于地面点云的高度,如果是,则判定此处粒子为不可移动粒子,否则判定此处粒子为可移动粒子。
对可移动粒子考虑内力情况,从而最终拟合出布的形状。
具体地,如图4所示,基于CSF算法对下采样后的部分初始点云数据进行地面滤波的过程为:
S31、对下采样后的部分初始点云数据进行过滤,滤除离群点,倒置点云。
S32、初始布栅格化,根据实际点云密度确定栅格大小(Grid Resolution,GR),初始位置位于点云纵坐标的最大值上。
S33、将点云与栅格粒子投影至同一水平面,寻找每个粒子的最近点(Corresponding Point,CP),将该点投影前的高度值定义为相交高度值(IntersectionHeight Value,IHV),IHV表示当前粒子可抵达的最低位置。
S34、计算每个粒子在外力影响(仅考虑重力)下移动的距离,将位置改变后的高度与IHV比较,如果位置改变后的高度小于或等于IHV,则将该粒子高度设置为等于IHV高度,并将该粒子标记为不可移动粒子。
粒子在外力影响下移动的距离用如下公式进行计算:
式(2)中,X(t+Δt)表示粒子在外力影响下移动的距离,X(t)表示粒子在t时刻的位置,X(t-Δt)表示粒子在t-Δt时刻的位置,Δt表示时间间隔,G表示重力常数。
S35、计算每个粒子在内力影响下移动的距离,相邻两个粒子间存在三种情况:1)均受力为可移动粒子,则向相反方向移动相同的距离;2)二者受力不同,一个移动,另一个不移动;3)均不受力为不可移动粒子,原地不动。粒子在内力影响下移动的距离用如下公式进行计算:
式(3)中,表示粒子的位移向量,当粒子为可移动粒子时,b=1;否则,b=0。/>表示当前粒子,/>表示/>的相邻粒子,/>表示垂直方向的归一化向量,即/>
通过定义刚性因子(Rigidness,RI)来表示粒子的移动次数,常规的刚性因子为1、2、3,刚性因子越大,布的表面越硬,表面越趋向平坦。
S36、按照步骤S34和步骤S35不断迭代,直到所有粒子的最大高度变化小于预设阈值或者迭代次数到达指定次数,计算停止。
S37、计算栅格化粒子与原始点云的距离。
S38、判断栅格化粒子与原始点云的距离是否小于预设的距离阈值,如果是,则判定当前点云为地面点云;否则,判定当前点云为非地面点云。
S39、判断是否遍历完所有点云,如果是,则结束;否则,返回至步骤S37,重新计算栅格化粒子与原始点云的距离。
在隧道内获取的点云数据,正常情况下地面都处于较平状态,相对于其他非地面点云更具有连通性和相似性,因此采用CSF的方式可以有效地分割地面点云与非地面点云,保留非地面点云,为后续提取隧道点云奠定基础。
S4、对滤波后的点云数据构建三维空间特征和二维平面特征,并将三维空间特征和二维平面特征划分为训练集和测试集,利用训练集创建支持向量机,采用支持向量机提取隧道点云。
三维点云一般包括空间(x,y,z)坐标、RGB颜色、反射率等信息,利用点云数据可以生成特征数据,机器学习算法可以对特征数据进行学习训练,生成稳定模型,最终快速、准确识别目标点云。
SVM(Support Vector Machine,支持向量机)是机器学习算法中用于完成分类工作的经典算法,通过输入需要提取隧道的点云数据,构建特征向量,生成模型最优参数,实现对输入点云是否为隧道点云的二分判别,达到提取隧道点云的目的,为后续隧道三维点云断面质量检测任务提供高质量数据。
如图5所示,采用支持向量机提取隧道点云的过程为:
S41、构建隧道点云特征向量,其具体过程为:
首先,在点云数据上构造三维KD-树,目的是将离散的三维点云转换为有关联的拓扑结构。
其次,在KD-树中寻找每个点的邻近点云簇,点云簇的个数根据实际的点云密度进行确定。
最后,计算邻近点云簇内所有点云的协方差特征值,构建三维空间特征和二维平面特征,以供后续支持向量机输入使用。
假设p(x,y,z)∈U3为三维点云空间中的一个点云,NU(p)表示点云p(x,y,z)在KD-树近邻搜索中搜索的最近点云簇,n表示点云簇的个数,则该点云簇形成的矩阵为:
式(4)中,表示第i2个点云簇的坐标,n表示点云簇的个数,i2=1,2,…,n。
为点云簇的中心点,构建与矩阵A的行数和列数均相等的矩阵B:
点云簇的协方差矩阵为:
根据点云簇的协方差矩阵计算出的协方差特征值为λ1,λ2,λ3,利用协方差特征值λ1,λ2,λ3构建如下三维空间特征:
发散性指数Sλ:
平面性指数Pλ:
线性指数Lλ:
各向同性Oλ:
各向异性Aλ:
特征熵Eλ:
特征和∑λ:
∑λ=λ1+λ2+λ3 (13)
曲率变化Cλ:
垂直度V:
V=1-nz (15)
式(15)中,nz表示每个点法向量的z轴分量。
点云簇z轴极差Rz:
Rz=Max(z)-Min(z) (16)
点云簇z轴方差
点云簇密度D:
式(18)中,k表示使用KD-树近邻搜索点云簇的k值,r表示搜索半径。
除此之外,将点云按照平行z轴方向进行切片,每个切片即为隧道实际的横截面。为确保每个横截面形状完整,因此可以将切片厚度设置为10cm。根据设计横截面及设计轨顶位置,确定每个实际横截面Si′={p1,p2,…,pn′}的实际轨顶位置其中,p1,p2,…,pn′表示实际横截面Si′上的各个点,n′表示点数。将实际横截面投影至xoz平面,计算第i′个横截面下第j点与第i′个横截面的轨顶距离dis(i′,j),该距离即为二维平面特征,计算公式如下:
式(19)中,表示实际轨顶位置/>的坐标,(x(i′,j),z(i′,j))表示第i′个横截面下第j个点。
S42、根据构建的隧道点云特征向量,将滤波后的点云数据的三维空间特征和二维平面特征划分为训练集和测试集。
其中,训练集用于训练SVM模型,确定最优参数,因此需要对样本集中的每个样本进行数据标注。
确定mtrain个输入样本/>
其中,表示第i4个样本的n维特征向量,i4=1,2,…,mtrain。
根据步骤S41构建的点云特征向量,此处n=13,表示第i4个样本的类别,-1表示非隧道点云,1表示隧道点云。
测试集作为SVM最优参数模型的输入,将mtest个测试样本输入到训练好的SVM模型中进行隧道点云判别,如果判定为隧道点云,则将该测试样本保留;否则,将该测试样本丢弃。
S43、利用训练集创建支持向量机,获取最优参数下的支持向量机模型;
支持向量机是通过解决凸优化问题,寻找最合适的支持向量,构建最优分离超平面,从而进行二分类决策。在遇到线性可分的样本时,可以通过引入软间隔降低异常样本点带来的影响,提升最优分离超平面的鲁棒性和泛化性;在遇到线性不可分的样本时,通过提升样本特征空间维度构建最优分离超平面,引入合适的核函数以解决高维度的计算灾难问题。
将步骤S42得到的训练集作为输入样本创建支持向量机,获取最优参数下的支持向量机模型的具体过程为:
S431、如图6所示,定义分离超平面H(x),当H(x)=0时,表明样本点位于分离超平面上;当H(x)>0时,表明样本类别为y=1;当H(x)<0时,表明样本类别为y=-1。
H(x)的表达式为:
H(x)=wTx+b (20)
式(20)中,w表示参数向量,b表示平移向量。
S432、构建最优分离超平面,即寻找图4中示出的间隔最大时的超平面,即:
式(21)中,i5=1,2,…,m。
S433、对式(21)进行转换并引入拉格朗日乘子,化解后最终优化目标函数为:
式(22)中,α表示向量,i5和j均表示序号。
S434、为解决线性不可分样本,定义核函数和惩罚系数C>0,最终优化目标函数为:
式(23)中,j=1,2,…,m。
S435、利用SMO(Sequential minimal optimization,序列最小优化算法)算法求解出向量α,则可获得对应的参数向量w:
S436、根据KKT条件找出所有的S个支持向量,计算每个支持向量的平均值b:
S437、确定最优超平面,分类决策函数为:
f(x)=sign(wTx+b) (26)
S44、利用训练好的支持向量机模型,将测试集作为输入,对隧道点云进行二分类判别;如果输出结果为1,则判定当前输入为隧道点云,保留当前输入的测试集中的点云数据;如果输出结果为-1,则表示当前输入为非隧道点云,丢弃当前输入的测试集中的点云数据。
下面采用一个具体的实施例对本申请提供的隧道三维点云提取方法进行说明。
采用Leica RTC360地面式三维激光扫描仪分别在宣绩铁路路口隧道、金建铁路上徐隧道及西十铁路西岭隧道采集点云数据。
基于CSF算法对点云数据中的地面进行滤波:
如图7(a)所示,初始点云长度为153.9m。在扫描过程中距离扫描仪越近,扫描的点云密度越高,获得的点云数据越好,因此在计算前需要将扫描仪附近的点云按照里程分割出来。
如图7(b)所示,本实施例分割出长度为21m的部分点云进行计算,此时点云数量为23215278。为了加速后续计算速度,对当前点云进行下采样处理。如图7(c)所示,下采样后点云数量变为196132。最后,采用CSF算法对地面进行剔除,保留剩余点云。如图7(d)所示,可以看出此时隧道点云与非隧道点云连通性较差,离散性较强。
基于SVM对隧道点云进行提取:
构建三维空间特征和二维平面特征,对每组特征归一化处理,防止特征值过大或过小产生的影响。
将70%的特征数据划分为训练集,30%的特征数据划分为测试集,对训练集中的隧道点云和非隧道点云进行标注。采用高斯核函数并确定其他最优参数,将五折交叉验证的平均值作为当前训练模型的评分结果,定义准确率为其中N表示预测正确的点云数,M表示总点云数。最终在不同测试集上的预测准确率如表1所示。
表1支持向量机预测结果
由表1可知,该支持向量机模型在不同测试集上的预测准确率为99%左右,准确率较高。这是因为由CSF算法准确剔除地面后,非隧道点云与隧道点云处于离散状态,且具有明显的特征不一致性,因此可以准确提取出如图8所示的隧道点云提取结果。
与隧道距离非常接近的少部分管道、线路、矮边墙点云会被误判为隧道点云,但这些数量占总点云数的比重非常小,对后续隧道三维点云断面质量检测影响可忽略不计。
表1中的三组测试数据均为形状相似的不同隧道,但可以使用相同的训练模型进行隧道点云提取,因此本申请具有较好的通用性。
本申请提供的三维点云提取方法首先采用CSF算法对隧道内获取的初始点云进行地面滤波,随后将保留的离散三维点云转换为有关联的拓扑结构,构建多种三维空间特征与二维平面特征以用于扩展点云数据,最后通过训练最优支持向量机模型快速、准确地提取出隧道点云,本申请具有以下优势:
(1)准确性高;
由于隧道点云与非隧道点云具有明显的特征不一致性,先采用CSF滤波将地面进行剔除,再构建点云特征空间,使用支持向量机对隧道点云进行提取,可以准确地提取出隧道点云,降低了后续对隧道点云进行检测时出现的误判情况。
(2)鲁棒性好;
由于隧道内施工情况较为复杂,经常会将正在施工的台车、除渣车、喷浆车、人员、通风管、排水管等点云数据一起获取到,但隧道点云相对于该类噪声点云特征相对稳定,且形状相似的隧道特征基本一致,因此模型可以迁移至相似隧道使用,具有较好的鲁棒性。
(3)速度较快;
与现有技术相比,本申请提供的三维点云提取方法的处理速度具有显著提升,能够满足现场施工规范,具有较好的市场前景。
(4)具备较强迁移性;
本申请提供的三维点云提取方法可成功迁移至铁路装备检测领域,尤其在铁路罐车三维激光扫描点云数据中进行罐体提取具备较好的识别能力,可以剔除罐体内部的梯子、液位计、标靶等附件点云,准确提取罐壁,为后续计算罐体容积,铁路罐车智能检定检测奠定良好的基础。
基于本申请实施例提供的三维点云提取方法,本申请实施例还提供了一种三维点云提取系统,其包括扫描模块、下采样模块、滤波模块和提取模块。
其中,扫描模块,用于获取初始点云数据。
下采样模块,用于根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理。
滤波模块,基于CSF算法对下采样后的部分初始点云数据进行地面滤波,以对隧道内特征较为连贯的地面进行识别和剔除。
提取模块,用于根据滤波后的点云数据得到训练集和测试集,并利用训练集创建支持向量机,采用支持向量机提取隧道点云。
需要说明的是:上述实施例提供的三维点云提取系统仅以上述各程序模块的划分进行举例说明,实际应用中,可以根据需要而将上述处理分配由不同的程序模块完成,即将三维点云提取系统的内部结构划分成不同的程序模块,以完成以上描述的全部或者部分处理。另外,上述实施例提供的三维点云提取系统与三维点云提取方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
在示例性实施例中,本申请实施例还提供了一种存储介质,是计算机可读存储介质,例如,包括计算机程序的存储器,上述计算机程序可由处理器执行,以完成前述三维点云提取方法中的所述步骤。
上述的本申请实施例可在各种硬件、软件编码或两者组合中进行实施。例如,本申请的实施例也可为在数据信号处理器中执行上述方法的程序代码。本申请也可涉及计算机处理器、数字信号处理器、微处理器或现场可编程门阵列执行的多种功能。可根据本申请配置上述处理器执行特定任务,其通过执行定义了本申请揭示的特定方法的机器可读软件代码或固件代码来完成。可将软件代码或固件代码发展为不同的程序语言与不同的格式或形式。也可为不同的目标平台编译软件代码。然而,根据本申请执行任务的软件代码与其他类型配置代码的不同代码样式、类型与语言不脱离本申请的精神与范围。
以上所述仅为本申请示意性的具体实施方式,在不脱离本申请的构思和原则的前提下,任何本领域的技术人员所做出的等同变化与修改,均应属于本申请保护的范围。
Claims (10)
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λ、平面性指数Pλ、线性指数Lλ、各向同性Oλ、各向异性Aλ、特征熵Eλ、特征和∑λ、曲率变化Cλ、垂直度V、点云簇z轴极差Rz、点云簇z轴方差和点云簇密度D;
发散性指数Sλ为:
平面性指数Pλ为:
线性指数Lλ为:
各向同性Oλ为:
各向异性Aλ为:
特征熵Eλ为:
特征和∑λ为:∑λ=λ1+λ2+λ3;
曲率变化Cλ为:
垂直度V为:V=1-nz;
点云簇z轴极差Rz为:Rz=Max(z)-Min(z);
点云簇z轴方差为:/>
点云簇密度D为:
上述各式中,λ1,λ2,λ3表示协方差特征值,nz表示每个点法向量的z轴分量,k表示使用KD-树近邻搜索点云簇的k值,r表示搜索半径;
二维平面特征包括隧道垂直横截面内每个点与轨顶之间的距离dis(i′,j);
式中,表示实际轨顶位置/>的坐标,(x(i′,j),z(i′,j))表示第i′个横截面下第j个点。
9.根据权利要求7所述的三维点云提取方法,其特征在于,所述利用训练集创建支持向量机,获取最优参数下的支持向量机模型的过程为:
定义分离超平面H(x),当H(x)=0时,表明样本点位于分离超平面上;当H(x)>0时,表明样本类别为y=1;当H(x)<0时,表明样本类别为y=-1;
H(x)的表达式为:
H(x)=wTx+b;
式中,w表示参数向量,b表示平移向量;
构建最优分离超平面的表达式为:
式中,i5=1,2,…,m;
对最优分离超平面的表达式进行转换并引入拉格朗日乘子,化解后最终优化目标函数为:
式中,α表示向量,i5和j均表示序号;
定义核函数K(xi5,xj)和惩罚系数C>0,最终优化目标函数为:
式中,j=1,2,…,m;
利用SMO算法求解出向量α,并获得对应的参数向量w:
根据KKT条件找出所有的S个支持向量,计算每个支持向量的平均值b:
确定最优超平面,分类决策函数为:
f(x)=sign(wTx+b)。
10.一种基于如权利要求1~9任一项所述的三维点云提取方法的提取系统,其特征在于,包括扫描模块、下采样模块、滤波模块和提取模块;
所述扫描模块,用于获取初始点云数据;
所述下采样模块,用于根据里程范围在初始点云数据中切分出待计算的部分初始点云数据并对该部分初始点云数据进行下采样处理;
所述滤波模块,基于CSF算法对下采样后的部分初始点云数据进行地面滤波,以对隧道内的地面进行识别和剔除;
所述提取模块,用于根据滤波后的点云数据得到训练集和测试集,并利用训练集创建支持向量机,采用支持向量机提取隧道点云。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310546415.8A CN116681653A (zh) | 2023-05-16 | 2023-05-16 | 一种三维点云提取方法及提取系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310546415.8A CN116681653A (zh) | 2023-05-16 | 2023-05-16 | 一种三维点云提取方法及提取系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116681653A true CN116681653A (zh) | 2023-09-01 |
Family
ID=87786320
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310546415.8A Pending CN116681653A (zh) | 2023-05-16 | 2023-05-16 | 一种三维点云提取方法及提取系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116681653A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117291845A (zh) * | 2023-11-27 | 2023-12-26 | 成都理工大学 | 一种点云地面滤波方法、系统、电子设备及存储介质 |
-
2023
- 2023-05-16 CN CN202310546415.8A patent/CN116681653A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117291845A (zh) * | 2023-11-27 | 2023-12-26 | 成都理工大学 | 一种点云地面滤波方法、系统、电子设备及存储介质 |
CN117291845B (zh) * | 2023-11-27 | 2024-03-19 | 成都理工大学 | 一种点云地面滤波方法、系统、电子设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Deng et al. | Review on computer vision-based crack detection and quantification methodologies for civil structures | |
CN109658387B (zh) | 电力列车的受电弓碳滑板缺陷的检测方法 | |
Yan et al. | Automated extraction of structural elements in steel girder bridges from laser point clouds | |
Movaghati et al. | Road extraction from satellite images using particle filtering and extended Kalman filtering | |
US9619691B2 (en) | Multi-view 3D object recognition from a point cloud and change detection | |
Sánchez-Rodríguez et al. | Automated detection and decomposition of railway tunnels from Mobile Laser Scanning Datasets | |
Xu et al. | Reconstruction of scaffolds from a photogrammetric point cloud of construction sites using a novel 3D local feature descriptor | |
CN110889318B (zh) | 利用cnn的车道检测方法和装置 | |
Lee et al. | Semantic segmentation of bridge components based on hierarchical point cloud model | |
CN112184736B (zh) | 一种基于欧式聚类的多平面提取方法 | |
CN110879961B (zh) | 利用车道模型的车道检测方法和装置 | |
CN111582380A (zh) | 一种基于时空特征的船舶轨迹密度聚类方法及装置 | |
CN105809113A (zh) | 三维人脸识别方法及应用其的数据处理装置 | |
CN116681653A (zh) | 一种三维点云提取方法及提取系统 | |
Liang et al. | An extraction and classification algorithm for concrete cracks based on machine vision | |
Özdemir et al. | Automatic extraction of trees by using multiple return properties of the lidar point cloud | |
König et al. | What's Cracking? A Review and Analysis of Deep Learning Methods for Structural Crack Segmentation, Detection and Quantification | |
CN113609895A (zh) | 基于改进Yolov3的获取道路交通信息方法 | |
Mirzaei et al. | End-to-end point cloud-based segmentation of building members for automating dimensional quality control | |
Xing et al. | An improved automatic pointwise semantic segmentation of a 3D urban scene from mobile terrestrial and airborne LiDAR point clouds: A machine learning approach | |
Zeybek | Inlier point preservation in outlier points removed from the ALS point cloud | |
Gehrung et al. | A fast voxel-based indicator for change detection using low resolution octrees | |
CN116452826A (zh) | 基于机器视觉的遮挡情况下煤矸石轮廓估计方法 | |
CN116052099A (zh) | 一种面向非结构化道路的小目标检测方法 | |
CN115100464A (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 |