基于平面多边形匹配的3D岩体点云配准方法技术

技术编号:21005006 阅读:47 留言:0更新日期:2019-04-30 21:45
本发明专利技术涉及一种基于平面多边形匹配的3D岩体点云配准方法。本发明专利技术主要包含3个阶段(1)基于体素的多次区域生长的平面检测(2)基于PCA的多边形匹配(3)基于匹配点对的点云配准。首先,初始岩体点云被划分为体素网格,在每个体素内部利用区域生长算法提取可靠的生长单元;基于生长单元进行二次生长得到初始平面集合;对于初始平面集合对每个平面周围邻域内的剩余点云进行生长,提取岩体点云中的平面集合;然后提取平面点云的多边形边界,利用PCA将多边形投影到二维平面并进行多边形匹配;利用匹配关系计算初始变换矩阵,最后利用ICP算法得到精确的配准结果,本发明专利技术适用于大场景岩石激光扫描点云的配准过程,可以高效、准确的实现岩体点云配准。

【技术实现步骤摘要】
基于平面多边形匹配的3D岩体点云配准方法
本专利技术属于计算机视觉和三维点云数据处理的交叉领域,涉及点云处理、岩体三维重建等技术,特别涉及3D岩体点云配准方法。
技术介绍
激光扫描仪以高分辨率、高准确度的获取目标物体的三维坐标信息并广泛的应用于三维重建,三维目标识别,虚拟现实等领域。但是在获取大场景目标时,往往一次扫描难以获取目标的全貌,需要进行多站扫描,点云配准的目标是将激光扫描仪获取的多站数据统一到同一个坐标系下。点云配准作为点云处理的初始环节,是目前研究的热点问题。目前在三维点云配准领域,最近邻迭代(ICP)是最常用的方法,其通过迭代的使误差最小化来得到变换矩阵,ICP算法可以获得非常精确的配准效果;不必对处理的点集进行分割和特征提取;在较好的初值情况下,可以得到很好的算法收敛性。但是算法性能跟点云初始位置相关性较大,容易陷入局部最优。为了解决这个问题,初始配准主要是在没有任何关于待配准点云之间相对位置关系信息的基础上计算两个点云之间刚体变换的粗略估计,该变换只是点云数据匹配的一个近似匹配,点云经过初始配准会得到一个相近的位置,然后在通过精细配准得到最终的配准结果。目前较多的利用点特征直方图(FPFH)、SIFT特征等实现点云的初始配准。在岩体工程领域,岩体往往是大场景环境,因此需要利用点云配准过程将激光扫描仪获取的多站扫描数据统一到同一坐标系下,这为后续的岩体三维建模以及数值计算过程提供零基础。而且岩体表面不规则、粗糙程度以及尺度不一,传统的特征提取方法难以在岩体点云中提取有效的特征用于配准过程。因此对岩体点云的配准的研究具有重要的意义。专利技术内容本专利技术针对岩体表面呈现出粗糙、形状不规则等特性,提出了一种基于平面多边形匹配的3D岩体点云配准方法。由于岩体表面主要有形状不规则的平面组成,平面特征能够很好适用于岩体点云的配准过程,并具有较高的精度。本专利技术提供一种基于平面多边形匹配的3D岩体点云配准方法,主要包括以下步骤:步骤1、基于体素的多次区域生长的平面检测:根据岩体点云的包围盒尺寸将点云快速体素化,然后通过基于体素的区域生长算法检测岩体点云中平面。此方法分别应用于源点云与目标点云中;步骤2、基于主成分分析的多边形匹配:将提取后的多边形利用主成分分析方法投影到同一二维坐标系下,利用平面的多边形轮廓实现多边形匹配;步骤3、建构匹配点对,实现点云配准:利用多边形的匹配关系,在三维空间中构建匹配点对实现点云的初始配准;最后利用最近邻迭代实现点云的精细配准。本专利技术所述的一种基于平面多边形匹配的3D岩体点云配准方法,作为优选方式,步骤1进一步包括以下步骤:步骤1.1、根据给定体素长度阈值,将点云快速的体素化;步骤1.2、在每个体素内部使用基于点的区域生长算法提取可靠的生长单元;步骤1.3、利用步骤1.2提取的生长单元,实现基于体素的区域生长过程,得到平面的初始集合;步骤1.4、对于步骤1.3提取的初始平面集合,以每个初始平面为基础,对该平面周围邻域部分进行一次生长过程,得到最终的平面集合。本专利技术所述的一种基于平面多边形匹配的3D岩体点云配准方法,作为优选方式,步骤2进一步包括以下步骤:步骤2.1、根据给定的参数利用凹包算法,对岩体点云中提取的平面计算其多边形边界;步骤2.2、从源点云多边形选取一个多边形si,包含n个点pi(1≤i≤n),将坐标中心平移到期质心则每个点的新坐标为pi=pi-Os;步骤2.3、利用主成分分析(PCA)将该多边形投影到二维坐标系;首先构建协方差矩阵Σ,λi和vi是该协方差矩阵对应的特征值和特征向量(λ1≤λ2≤λ3)。构建变换矩阵T={v1,v2},则多边形每个点pi所对应的二维空间坐标为pi_2d=pi*T;Σ·vj=λj·vj,j∈{1,2,3}步骤2.4、在目标点云中顺序选取多边形tj,将选取的多边形利用步骤2.2-2.3投影到同一个二维坐标系,计算多边形相交围成的面积若其占比Sr超过阈值则认为两个多边形是相匹配的;步骤2.5、重复在源点云中选取多边形,直到所有源点云多边形被处理,保存所有匹配关系。本专利技术所述的一种基于平面多边形匹配的3D岩体点云配准方法,作为优选方式,步骤3进一步包括以下步骤:步骤3.1、顺序从步骤2选取对应匹配关系,找到多边形所对应的三维平面点云,分别计算两个平面点云的质心Cs和Ct,构建匹配点对;步骤3.2、由于点云配准以刚性变换为主,包含旋转矩阵R和平移矩阵t。假设总共有K个对应点对。则利用解最优化问题,使配准误差ε最小,求得初始变换矩阵;步骤3.3、将经过初始变换的两个点云放入最近邻迭代(ICP)算法中得到最终的配准结果。本专利技术针对岩体点云特性的分析,利用基于体素的多次区域生长算法,快速的提取岩体点云中的平面并具有较高的准确性;利用PCA将平面投影到同一二维坐标系下,利用多边形匹配寻找对应点对实现点云的初始配准,最后利用ICP算法实现点云的精细配准,该算法提升了算法在应对不规则以及粗糙程度不一的岩体环境,也提升了算法的准确性。附图说明图1是基于平面多边形匹配的3D岩体点云配准方法流程图。图2是提取生长单元示意图。图3是平面提取结果示意图。图4是多边形匹配结果示意图。图5是最终配准结果示意图。具体实施方式本专利技术主要基于计算机视觉和点云处理技术,提出一种基于平面多边形匹配的3D岩体点云配准方法。本专利技术利用平面提取算法提取了岩体中的平面特征,利用多边形匹配寻找点对匹配关系,提高了配准的效率和准确性,使该算法能够在岩体等复杂场景下获得准确的点云配准结果。为了使本专利技术的目的、技术方案及优点更加清楚,以下结合附图和实施例,对本专利技术作进一步详细说明。本专利技术提供的方法能够用计算机软件技术实现流程,整体技术流程图参见图1,包括以下步骤:步骤1、为了实现快速的点云平面提取过程,首先需要设定体素的边长,实例的具体实施过程说明如下:步骤1.1、根据给定的体素边长voxel_size(本实例中取值0.8),将点云体素化;步骤1.2、选取体素内部区域生长算法准则:1邻域点法线与种子点集法线的夹角小于阈值θ1(本实例中取值10°),2邻域点与种子点集之间的距离小于阈值d1(本实例取值0.08m);步骤1.3、对于每个体素内部,如果点数大于一定的阈值Nr(本实例中取值30)则利用基于点的区域生长算法在平面内部提取平面作为生长单元。首先选取与体素中心点距离最近的点作为种子点,加入生长队列Q1;如果生长队列不为空,则在Q1中取一个点在其邻域向外生长,满足生长规则的点将会被加入种子点集,并加入生长队列Q1;直到生长队列为空则生长结束。该过程持续数次(本实例中取值10),选取最大的平面作为生长单元。提取后的生长单元结果如图2所示。步骤1.4、从步骤1.3得到的生长单元选取未被处理的、包含最大点数的平面作为种子平面,加入到生长队列Q2中,并标记为已处理;步骤1.5、建立区域生长准则。在本实例中选取准则:1平面与种子平面的夹角小于阈值θ2(本实例中取值15°),2平面与种子平面之间的距离小于阈值d2(本实例取值0.12m);步骤1.6、若生长队列不为空,则从生长队列Q2中取出一个平面。判断该平面的未被处理的邻域体素是否满足生长准则,如果满足准则,该邻域平面合并到种子点集中,同时本文档来自技高网...

【技术保护点】
1.一种基于平面多边形匹配的3D岩体点云配准方法,其特征在于主要包括以下步骤:步骤1、基于体素的多次区域生长的平面检测:根据岩体点云的包围盒尺寸将点云快速体素化,然后通过基于体素的区域生长算法检测岩体点云中平面。此方法分别应用于源点云与目标点云中;步骤2、基于主成分分析的多边形匹配:将提取后的多边形利用主成分分析方法投影到同一二维坐标系下,利用平面的多边形轮廓实现多边形匹配;步骤3、建构匹配点对,实现点云配准:利用多边形的匹配关系,在三维空间中构建匹配点对实现点云的初始配准;最后利用最近邻迭代实现点云的精细配准。

【技术特征摘要】
1.一种基于平面多边形匹配的3D岩体点云配准方法,其特征在于主要包括以下步骤:步骤1、基于体素的多次区域生长的平面检测:根据岩体点云的包围盒尺寸将点云快速体素化,然后通过基于体素的区域生长算法检测岩体点云中平面。此方法分别应用于源点云与目标点云中;步骤2、基于主成分分析的多边形匹配:将提取后的多边形利用主成分分析方法投影到同一二维坐标系下,利用平面的多边形轮廓实现多边形匹配;步骤3、建构匹配点对,实现点云配准:利用多边形的匹配关系,在三维空间中构建匹配点对实现点云的初始配准;最后利用最近邻迭代实现点云的精细配准。2.根据权利要求1所述的一种基于平面多边形匹配的3D岩体点云配准方法,其特征在于:所述步骤1进一步包括以下步骤:步骤1.1、根据给定体素长度阈值,将点云快速的体素化;步骤1.2、在每个体素内部使用基于点的区域生长算法提取可靠的生长单元;步骤1.3、利用步骤1.2提取的生长单元,实现基于体素的区域生长过程,得到平面的初始集合;步骤1.4、对于步骤1.3提取的初始平面集合,以每个初始平面为基础,对该平面周围邻域部分进行一次生长过程,得到最终的平面集合。3.根据权利要求1所述的一种基于平面多边形匹配的3D岩体点云配准方法,其特征在于:所述步骤2进一步包括以下步骤:步骤2.1、根据给定的参数利用凹包算法,对岩体点云中提取的平面计算其多边形边...

【专利技术属性】
技术研发人员:胡亮肖俊王颖
申请(专利权)人:中国科学院大学
类型:发明
国别省市:北京,11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1