基于点线特征的建筑物点云配准方法.pdf

第4 6卷第6期2 0 2 3年1 1月现 代 测 绘M o d e r nS u r v e y i n ga n dM a p p i n gV o l.4 6,N o v.,2 0 2 3 第一作者简介:尹盘飞,高级工程师,研究方向为测绘地理信息、航测遥感数据解译等。基于点线特征的建筑物点云配准方法尹盘飞,刘 霞,魏 军(三和数码测绘地理信息技术有限公司,甘肃 天水7 4 1 0 0 0)摘 要 针对建筑物表面具有明显的几何特征这一特点,提出基于点线特征的建筑物点云配准方法。首先基于点云配准基本理论,利用提取的建筑物点特征使用对偶四元数实现建筑物不同测站点云的粗配准,获取初始配准参数以及配准后的建筑物点云数据,然后将粗配准获取的参数作为待求精确参数值的初始值,利用建筑物中的线特征将共线方程作为精配准的数学模型,最后通过平差迭代获取参数的精确值,实现不同测站建筑物点云的高精度配准。实验结果表明,获取的配准后同名特征距离中误差为2.11 0-3m,证明了将点线特征相结合可以有效改善建筑物点云配准质量,提高建筑物点云配准的精度,对建筑物三维重建具有重要意义。关键词 建筑物点云;对偶四元数;共线方程;点云配准中图分类号:P 2 3 7;T P 3 9 1.4 1 文献标志码:A 文章编号:1 6 7 2-4 0 9 7(2 0 2 3)0 6-0 0 2 1-0 40 引 言随着三维扫描技术的快速发展,三维重建技术也得到了广泛的应用,主要在历史珍贵文物数字化保护领域1、医学领域2、自动驾驶或无人机飞行领域3等,三维重建技术为人类的生活提供了极大的便利,推动了各个行业的快速发展。由于三维扫描仪的视角有限,每次仅能扫描到物体的部分点云,需要将获取的多个视角的点云通过旋转、平移的步骤,将各站点云配准到同一个基准坐标系下,从而获得完整的信息,因此配准的质量和效率对三维重建具有重要意义,是三维重建过程中的关键步骤。点云配 准 一 般 包 括 粗 配 准 和 精 配 准 两 个 过程4,粗配准主要是将未知位置姿态的基准站点云和待配准站点云转移拼接到基准站点云所在的坐标系下,一定程度上减小基准站点云和待配准站点云在空间位置上的旋转与平移错位,为精配准提供良好的参数初值。粗配准常用的方法主要有4 P C S算法、快速点特征直方图、基于统计概率学的正态分布变换(N D T)算法等。精配准是在粗配准的基础上通过多次迭代优化将基准站点云和待配准站点云进行精确的配准,获取更加精确的配准参数值,提高点云数据配准的精度,使两片点云实现完好的拼接。其中被广泛应用的是由B e s l和M k c y a提出的最近点迭代算法(I t e r a t i v eC l o s e s tP o i n t,I C P),该方法对初值精度要求较高,在迭代过程中容易陷入局部最优解。点云配准在配准过程中根据选择的配准基元不同,将配准方法分为基于点特征基元、线特征基元、面特征基元,由于点特征基元在目标点云中大量存在,因此使用较为广泛,当建筑物部分点云被遮挡时,在基准站点云和待配准站点云中无法提取到精确的同名特征点,所以在点特征进行配准时,点云的配准质量会受到一定程度的影响。然而同名线特征和同名面特征受获取的点云质量影响较小,可以有效的减小点云配准误差。针对上述问题,本文提出利用提取的点特征,通过对偶四元数构建粗配准模型,解算配准参数,为精配准提供初值;然后从建筑物中提取线特征,并将共线方程作为精配准的数学模型,通过平差迭代解算精确配准参数。本文方法在粗配准过程中可以同时解算出旋转矩阵和平移参数,有效降低了平移矩阵和旋转矩阵之间的耦合误差,提高了点云的粗配准精度;在精配准过程中使用同名线特征,避免了选取精确的同名点,而且相对于仅使用单一特征,使用点线特征相结合可以获取更好的点云配准精度。1 点云基本理论1.1 点云配准基本理论不同测站点云配准的实质是将点云转换到同一个坐标系下,获取不同测站之间的转换关系,从而获取一个完整的物体。传统的点云转换模型为:X1Y1Z1=RX2Y2Z2+T(1)式中,X1、Y1、Z1为基准站坐标系下的点云坐标,X2、Y2、Z2为待配准站点云的坐标,R为旋转矩阵,T为平移矩阵。使用提取的点特征利用公式(1)实现不同测站点云拼接时,需要从不同测站点云中分别提取3个特征点。在求解点云配准参数时,为了同时解算出配准参数,减小分别解算配准参数的误差,使用对偶四元数解算配准参数旋转矩阵R以及平移矩阵T。1.2 对偶四元数对偶四元数是由C l i f f o r d提出的一种数学运算工具,由对偶数和四元数两部分组成,可以同时表示刚体的旋转和平移运动,对偶四元数的表达形式如式(2)所示:p=k+r=k0 kx ky kzT+r0 rx ry rzT(2)式中,k和r分别是对偶四元数p的实部和对偶部,为对偶单位。k0和kx、ky、kz为k的实部和虚部,r0和rx、ry、rz为r的实部和虚部。由于对偶四元数的实部和对偶部都属于四元数的对偶数,因此对偶数和四元数的运算法则也适用于对偶四元数。同时,对偶四元数的几何意义可以表示为2个三维集围绕着一个轴做刚体运动,先沿着n方向做平移运动,然后平移距离为d,平移到p点,然后在p点沿着n的方向旋转角度,如图1所示。图1 对偶四元数描述的刚体运动令q1=k1 k2 k3T,基于对偶四元数的运算法则以及性质得到q1与四元数矩阵U(k)和V(k)的关系如式(3)所示。U(k)=k0I+Sq1()q1-qT1k0V(k)=k0I-Sq1()q1-qT1k0(3)式中,S为反对称矩阵。根据式(2)可得,三维点集刚体运动的旋转矩阵和平移矩阵 可 以 用 对 偶 四 元 数表 示,如式(4)所示。R00T1=V(k)TU(k),T=V(k)Tr(4)通过计算,三维点云的配准参数旋转矩阵和平移矩阵可以用对偶四元数表示,如式(5)所示。R=a1a2a3b1b2b3c1c2c3=k20+k21-k22-k232k1k2-2k0k32k1k3+2k0k22k1k2+2k0k3k20-k21+k22-k232k2k3-2k0k12k1k3-2k0k22k2k3+2k0k1k20-k21-k22+k23T=X Y Z=2k0r1-2k1r0+2k2r3-2k3r22k0r2-2k1r3-2k2r0+2k3r12k0r3+2k1r2-2k2r1-2k3r0(5)式中,a1、a2、a3、b1、b2、b3、c1、c2、c3为旋转角、构成的方向余弦;(X,Y,Z)为平移向量。通过式(2)(3)可以得到对偶四元数描述的基于点特征的建筑物点云配准数学模型,如式(6)所示。Qi=P(k)TG(k)Qi+P(k)Ts(6)式中,Qi为基准站坐标系下点云的位置;Qi为待配准站坐标系下的点云的位置;P(k)TG(k)Qi为旋转矩阵;P(k)Ts为平移向量。将式(6)进行线性化,并从基准站点云和待配准站点云中提取3对同名特征点,得到的误差方程如式(7)所示。V=AX+L(7)式中,V=vxvyvz,A=A1 1A1 2A1 3A1 4A1 5A1 6A1 7A1 8A1 9A2 1A2 2A2 3A2 4A2 5A2 6A2 7A2 8A2 9A3 1A3 2A3 3A3 4A3 5A3 6A3 7A3 8A3 9,X=d k0d k1d k2d k3d r0d r1d r2d r3T,L=L1L2L3T通过平差解算出参数k0、k1、k2、k3、r0、r1、r2、r3,并将求解的参数值代入式(4),即可求出点云粗配准参数:旋转矩阵R和平移矩阵T。1.3 基于线特征的点云配准数学模型从基准站点云和待配准站点云分别提取同名直线,但是提取的同名直线中对应的端点不一定是同名特征点。因此本文选用共线方程作为基础方程来描述同名直线间的三维相似变换,并将此作为约束条件,如式(8)所示。22现 代 测 绘 第4 6卷XTYTZT+S RXaYaZa=XcYcZc+1Xd-XcYd-YcZd-ZcXTYTZT+S RXbYbZb=XcYcZc+2Xd-XcYd-YcZd-Zc(8)式中,XT、YT、ZT为待配准站点云到基准站点云平移参数,R为待配准站点云到基准站点云的旋转矩阵,S、1、2为待配准站点云到基准站点云的缩放因子。每一对同名直线都可以列出两个方程,如式(8)所示,将式(8)的两个公式相减,并获得到点云配准的数学模型,如式(9)所示,并利用四元数方法解算配准参数。Xd-XcYd-YcZd-Zc=v RXbYbZb-XaYaZa(9)2 实验分析为了验证本文方法 的 有效 性,采 用WHU-T L S基准数据集中一座建筑物点云数据,选取的建筑物的两站点云数据如图2所示,点云信息如表1所示。图2 配准前两测站点云的基本位置表1 实验数据基本信息激光点云个数基准站2 4 24 0 0待配准站1 2 75 1 5 在图2中,灰色的点云是基准站点云,黑色的点云是待配准站点云。由图2可知,在配准前两测站点云数据在空间位置上存在较大的差异。本文首先对两站点云进行了点云噪声处理,为了证明本文方法基于点线特征的建筑物点云配准效果,先从处理后的点云中均匀选取4对同名点进行点云粗配准,然后再在粗配准的基础上从点云中选取3对同名线特征进行精配准。4对同名点的坐标信息如表2所示。表2 两站中提取的同名特征点坐标同名点基准站待配准站X/mY/mZ/mX/mY/mZ/m1-3 6.5 3-1 2 9.7 6 6 2.1 9-1 3 9.3 7-1 0 4.2 56 2.5 02-4 5.0 9-9 0.8 96 7.7 4-1 1 7.2 0-7 1.6 3 6 8.0 43-1 0.7 6-1 2 6.3 9 4 7.1 1-1 1 9.5 3 7-1 2 0.9 34 7.3 44-3 5.3 4-1 3 1.8 3 4 0.4 8-1 4 0.0 2-1 0 7.0 14 1.0 2 使用提取的4对同名特征点通过对偶四元数将待配准测站点云配准到基准站点云所在的坐标系下,获取的粗配准效果如图3所示。由图3可知,经过粗配准后,很大程度上减小了待配准站点云与基准站点云在空间位置上的误差,使两站点云拼接到了一起,但可以发现,两站点云在有些地方融合得不是很好,仍然存在误差。图3 粗配准效果图在粗配准之后,将粗配准获取的配准参数值作为精配准的初始参数值。然后利用从两侧站中提取的4对同名线特征,基于共线条件方程实现了待配准测站和基准站的精配准,获取的配准参数如表3所示,精配准效果如图4所示。表3 获取的配准参数值旋转矩阵平移矩阵缩放因子0.6 8 80.7 2 70.0 0 1-2 0.1 0 3-0.7 2 70.6 8 8-0.0 0 2-4 1.8 7 70.9 8 7-0.0 0 20.0 0 11.0 0 10.3 6 7 由图4可知,待配准点云和基准站点云很好的融合在了一起,边界处也较为贴切,证明了本文方法的有效性。32第6期尹盘飞等:基于点线特征的建筑物点云配准方法图4 精配准效果图为了进一步分析基于点、基于点线的点云配准精度,本文分别选用点特征和点线特征进行建筑的点云配准,选用同名特征配准后的平均距离(距离为配准后的待配准站点云坐标与基准站点云坐标之差)中误差作为精度的评价指标。两种方法的配准精度如表4所示。表4 两种配准方法获取的同名特征配准中误差方法距离中误差/m基于点特征4.51 0-3基于点线特征2.11 0-3 由表4可知,采用点线特征相结合的建筑点云配准方法获取的同名特征点云的中误差比仅使用点特征的配准方法的精度高,证明了本文方法的可行性与有效性。3 结 语在建筑物的点云数据中,存在大量的几何特征,本文针对建筑存在的这一特性,提出基于点线特征结合的建筑物点云配准方法,该方法将选取不同的特征进行了建筑物不同测站的粗配准和精配准。实验结果表明,该方法获取的配准后同名特征的距离中误差为2.11 0-3m,使用多种特征相结合的配准方法比仅使用一种特征进行配准获取更好的配准精度,一定程度上也可以提高点云配准的质量。本文此次仅选用了建筑物中的点线特征,在后续研究中,可以加入面特征,从而进一步提高建筑物点云配准的精度。参考文献1 夏国芳,胡春梅,范亮.一种面向造像类文物的真三维模型精细重建方法J.敦 煌研究,2 0 1 8(3):1 3 1-1 4 0.2 赵夫群,周明全.颅骨点云模型的优化配准J.光学精密工程,2 0 1 7,2 5(7):1 9 2 7-1 9 3 3.3 刘晓文,付莉娜,徐工.基于3 D-S I F T与S I C P的激光点云与影像配准方法J.北京测绘,2 0 2 2,3 6(5):5 5 7-5 6 2.4 郭瑞科,王立,朱飞虎,等.空间非合作目标的多视角点云配准算法研究J.中国空间科学技术,2 0 1 6,3 6(5):3 2-3 9.B u i l d i n gP o i n tC l o u dR e g i s t r a t i o nM e t h o dB a s e do nP o i n tL i n eF e a t u r eY I NP a n-f e i,L I UX i a,W E I J u n(S a n h eD i g i t a lS u r v e y i n ga n dM a p p i n gG e o g r a p h i c I n f o r m a t i o nT e c h n o l o g yC o.L t d.,T i a n s h u iG a n s u7 4 1 0 0 0,C h i n a)A b s t r a c t Ap o i n t c l o u d r e g i s t r a t i o nm e t h o db a s e do np o i n t l i n e f e a t u r e s i sp r o p o s e d f o rb u i l d i n gs u r f a c ew i t ho b v i o u sg e o m e t r i cf e a t u r e s.F i r s t l y,b a s e do nt h eb a s i c t h e o r yo fp o i n t c l o u dr e g i s t r a t i o n,t h ee x t r a c t e db u i l d i n gp o i n t f e a t u r e sa r eu s e dt or e a l i z et h ec o a r s er e g i s t r a t i o no fd i f f e r e n tb u i l d i n gs i t ec l o u d sb yd u a lq u a t e r n i o n,a n dt h ei n i t i a lr e g i s t r a t i o np a r a m e t e r sa n dt h er e g i s t e r e db u i l d i n gp o i n t c l o u dd a t aa r eo b t a i n e d.T h e n,t h ep a r a m e t e r so b t a i n e db yc o a r s er e g i s t r a t i o na r eu s e da st h e i n i t i a lv a l u eo f t h ea c c u r a t ep a r a m e t e rv a l u e st ob ec a l c u l a t e d,a n dt h ec o l l i n e a re q u a t i o ni su s e da st h em a t h e m a t i c a lm o d e lo ff i n er e g i s t r a t i o nb yu s i n gt h el i n ec h a r a c t e r i s t i c si nt h eb u i l d i n g.F i n a l l y,t h ea c c u r a t ev a l u e so fp a r a m e t e r sa r eo b t a i n e db ya d j u s t m e n t i t e r a t i o n,a n dt h eh i g h-p r e c i s i o nr e g i s t r a t i o no fp o i n tc l o u d so fb u i l d i n g sa td i f f e r e n ts t a t i o n si sr e a l i z e d.T h ee x p e r i m e n t a l r e s u l t ss h o wt h a tt h ee r r o ro ft h es a m en a m ef e a t u r ed i s t a n c eo b t a i n e db yt h i sm e t h o di s2.11 0-3m,w h i c hp r o v e s t h a t t h ec o m b i n a t i o no f p o i n t a n d l i n e f e a t u r e s c a ne f f e c t i v e l y i m p r o v e t h eq u a l i t yo f b u i l d i n gp o i n t c l o u dr e g i s t r a t i o na n di m p r o v e t h ea c c u r a c yo fb u i l d i n gp o i n t c l o u dr e g i s t r a t i o n,w h i c hi so fg r e a ts i g n i f i c a n c et ot h r e e-d i m e n s i o n a l r e c o n s t r u c t i o no fb u i l d i n g s.K e yw o r d s p o i n t c l o u do fb u i l d i n g s;d u a l q u a t e r n i o n;c o l l i n e a re q u a t i o n;p o i n t c l o u dr e g i s t r a t i o n42现 代 测 绘 第4 6卷