基于层次化闭合环约束的多视点云全局优化配准方法技术

本发明专利技术公开了一种基于层次化闭合环约束的多视点云全局优化配准方法，属于计算机视觉和激光扫描数据处理的交叉领域。本发明专利技术针对手不同扫描视角产生的配准基元邻域信息的空间差异，根据配准基元的特征一致性和几何一致性对点云进行连通性分析，对多视点云构建连通图，将成对点云之间的相似变换统一到基准坐标系。基于层次化闭合环约束，通过最小二乘联合平差的误差分析与修正，得到全局优化精配准结果。果。果。

【技术实现步骤摘要】
基于层次化闭合环约束的多视点云全局优化配准方法


[0001]本专利技术属于计算机视觉和激光扫描数据处理的交叉领域，尤其涉及点云数据预处理，多视点云配准，针对激光点云测量数据处理自动化研究。

技术介绍

[0002]激光扫描测量技术是获取三维目标空间信息的重要手段，利用激光测距的原理创建被测目标几何表面的点云，通过记录大量的密集点的空间坐标、色彩纹理和反射强度等信息，可快速重建被测目标的三维实体及线、面、体等各种模型基元数据。随着三维点云获取技术的发展与推广，用户可以快捷地获取场景目标的多分辨率、多尺度、多时态点云数据。其中，手持式激光扫描仪具有良好的便携性、较高的几何精度(亚毫米级)，可快速地完成准确、高分辨率的三维扫描，在精密零部件测量、文化遗产数字化保护等应用中发挥着重要作用，为多尺度、多场景目标的三维空间信息获取提供了强有力的技术支撑。
[0003]三维点云记录了目标的全方位几何、属性信息，其对细节的描述也十分精确，但点云的高密度、海量、空间无序性以及场景中三维目标的数据不完整性，目标间的重叠性、遮挡性、相似性等现象给手持激光点云要素提取和增强带来了巨大的挑战。因此需要对多视点云数据进行拼接与配准，从而获取完整的目标点云，兼顾场景的宏观信息和细节特征。对于存在不同程度噪声干扰、大范围对称区域的大尺度点云，迫切需要一种鲁棒、高效、普适的算法来获得全局最优的配准结果，为后续点云分类分割、目标提取、曲面重建等处理过程提供便利。

技术实现思路

[0004]本专利技术针对不同扫描视角产生的配准基元邻域信息的空间差异，根据连通性对多视点云构建连通图，将成对点云之间的相似变换统一到基准坐标系。基于层次化闭合环约束，通过最小二乘联合平差的误差分析与修正，得到全局优化精配准结果。
[0005]本专利技术针对不同扫描视角下捕获的手持式激光扫描点云，提出一种基于层次化闭合环约束，对已完成特征提取和成对配准的点云对，通过最小二乘联合平差的误差分析与修正，进行多视点云全局优化配准方法。本方法分为4个关键步骤：全局连通性分析、层次化遍历闭合环、闭合环误差剔除和累积空间变换。
[0006]本专利技术通过以下技术流程解决上述问题：
[0007]一种基于层次化闭合环约束的多视点云全局优化配准方法，包括以下步骤：
[0008]步骤1，选定基准点云坐标系。首先，读取多视点云文件、同名特征点对信息以及闭合环粗差探测后的成对配准结果。然后对输入的多视点云，采用同名特征点对数衡量点云之间的配准关系强弱，将相应同名特征点数量定义为全局点云连通图中相应节点之间的边缘权重。继而，判断配准关系是否足够稳健，当同名特征点对数少于给定最小点对集合数，则认为点云之间的配准关系弱并给予拒绝，否则接纳配准关系。最后，计算连通图各节点的边权值之和，并按照降序排列，以边权值之和最大的节点作为根节点，即边权值之和最大节
点代表的点云为基准点云。
[0009]步骤2，层次化遍历闭合环。首先，对步骤1确定的根节点各边按权值从大至小遍历。若边不构成闭合环，则根据输入的成对配准结果与根节点进行配准；若边构成闭合环，则重新调整局部配准中的空间变换误差，进行闭合环误差重分配。
[0010]步骤3，闭合环误差重分配；构成闭合环的点云P
i
、P
j
和P
k
两两之间分别存在同名特征点对，由特征一致性约束和几何一致性检验得P
i
和P
j
之间的同名特征点对BSC
i
和BSC
j
，满足条件方程(1)：
[0011]BSC
i
＝R
i
→
j
BSC
j
+T
i
→
j
ꢀꢀ
(1)
[0012]同理，P
j
和P
k
之间的同名特征点对BSC
j
和BSC
k
，P
i
和P
k
之间的同名特征点对BSC
′
i
和BSC
′
k
，分别满足条件方程(2)和(3)：
[0013]BSC
j
＝R
j
→
k
BSC
k
+T
j
→
k
ꢀꢀ
(2)
[0014][0015]其中BSC
′
j
为P
k
中BSC
′
k
在P
j
中对应的虚拟同名特征点；
[0016]对以上三组方程整理并构造法方程(4)，采用联合平差法迭代计算出两组旋转矩阵R
i
→
j
、R
j
→
k
和平移矩阵T
i
→
j
、T
j
→
k
，其中观测值为同名特征点对，未知数由旋转角Δω、Δκ和平移参数ΔX、ΔY、ΔZ构成；
[0017][0018]其中，V表示观测方程未知数求解未知数的改正值，L
′
、L
″
表示观测值与近似值之差，A
ij
、B
jk
、A
′
ik
、B
′
ik
分别表示同名特征点对BSC
i
和BSC
j
、BSC
j
和BSC
k
、BSC
′
i
和BSC
′
k
、BSC
′
j
和BSC
′
k
对应关系的系数矩阵，N
ij
表示P
i
和P
j
之间的同名特征点对数量，N
jk
表示P
j
和P
k
之间的同名特征点对数量，N
ik
表示P
k
和P
i
之间的同名特征点对数量；
[0019]步骤3中法方程(4)的求解如下；
[0020]当V
T
PV最小时，求得未知数的解，此处P为单位阵，记B为求解参数建立的误差方程式的系数，则有X＝(B
T
PB)
‑1B
T
L；
[0021]其中，B对应法方程(4)中同名特征点对对应关系的系数矩阵，L对应法方程(4)中观测值与近似值之差。
[0022]步骤4，点云融合并更新连通图。将完成与根节点配准的点云融合为新的跟节点，更新连通图，并进行下一轮配准。重复上述步骤1
‑
3，直至所有点云都完成配准，通过累积空间变换，可以将每个点云配准到全局坐标系中，从而获得最终的多视点云配准结果。<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于层次化闭合环约束的多视点云全局优化配准方法，其特征在于，包括以下步骤：步骤1，选定基准点云坐标系；首先，读取多视点云文件、基于特征一致性和几何一致性提取的同名特征点对信息以及闭合环粗差探测后的成对配准结果；然后对输入的多视点云，采用同名特征点对数衡量点云之间的配准关系强弱，将相应同名特征点数量定义为全局点云连通图中相应节点之间的边缘权重；继而，判断配准关系是否足够稳健，当同名特征点对数少于给定最小点对集合数，则认为点云之间的配准关系弱并给予拒绝，否则接纳配准关系；最后，计算连通图各节点的边权值之和，并按照降序排列，以边权值之和最大的节点作为根节点，即边权值之和最大节点代表的点云为基准点云；步骤2，层次化遍历闭合环；首先，对步骤1确定的根节点各边按权值从大至小遍历；若边不构成闭合环，则根据输入的成对配准结果与根节点进行配准；若边构成闭合环，则重新调整局部配准中的空间变换误差，进行闭合环误差重分配；步骤3，闭合环误差重分配；构成闭合环的点云P
i
、P
j
和P
k
两两之间分别存在同名特征点对，由特征一致性约束和几何一致性检验得P
i
和P
j
之间的同名特征点对BSC
i
和BSC
j
，满足条件方程(1)：BSC
i
＝R
i
→
j
BSC
j
+T
i
→
j
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)同理，P
j
和P
k
之间的同名特征点对BSC
j
和BSC
k
，P
i
和P
k
之间的同名特征点对BSC
′
i
和BSC
′
k
，分别满足条件方程(2)和(3)：BSC
j
＝R
j
→
k
BSC
k
+T
j
→
k
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中BSC
′
j
为P
k
中BSC
′
k
在P
j
中对应的虚拟同名特征点；对以上三组方程整理并构造法方程(4)，采用联合平差法迭代计算出两组旋转矩阵R
i
→
j
、R
j
→
k
和平移矩阵T
i
→
j
、T
j
→
k
，其中观测值为同名特征点对，未知数由旋转角Δω...

【专利技术属性】
技术研发人员：张卉冉，董震，杨必胜，
申请(专利权)人：武汉大学，
类型：发明
国别省市：