
摘要:针对传统基于特征的粗配准效率低、误匹配较多的不足,提出一种基于特征空间匹配的配准方法。利用简化的PointNet模型实现特征空间的提取,以优化的点云PPF信息作为输入,根据提取的特征空间向量计算欧氏距离以筛选匹配点,通过RANSAC剔除误匹配点对完成粗配准,利用ICP实现精配准。实验结果表明,本文算法相比FPFH和SHOT算法与ICP结合可有效提升配准效率,且配准结果的均方根误差较小。
受测量条件限制,激光扫描技术往往需要架设多站点进行扫描,并对获取的点云数据进行配准以消除不同站点的相对位置变化。目前点云配准算法主要分粗配准和精配准2个阶段。迭代最近点算法(iterativeclosestpoint,ICP)[1]为目前应用最广泛且十分有效的一种精配准方法,但计算成本大、易陷入局部最小值。高效且鲁棒性优良的粗配准方法可为ICP提供较好的初始值[2,3],但传统的基于特征的粗配准方法效率低、误匹配点对较多。本文基于PointNet模型[4]提取点云深度特征的优点,将一种改进的点对特征(pointpairfeature,PPF)信息[5]作为输入,通过模型直接提取点云局部信息的特征空间而无需训练,提取的特征空间经实验验证具有较优的局部表达能力,然后与传统算法相结合依次完成粗配准与精配准。
1、基于特征空间匹配的配准算法
基于特征空间匹配的配准算法的主要步骤为:首先计算源点云与目标点云的内部形态描述子(intrinsicshapesignatures,ISS)关键点,求取关键点k邻域内优化的PPF信息,并用改进的PointNet模型进行特征提取;然后根据特征空间向量的欧氏距离完成匹配点对的筛选,通过随机采样一致性(randomsampleconsensus,RANSAC)剔除误匹配点对;最后利用k-d树加速ICP精配准。
1.1ISS特征点提取
首先由ISS算法计算源点云A与目标点云B(以下称两站点云)的关键点,该步骤阈值的选取较为关键,阈值过大会导致特征点数目多而不利于计算,阈值过小会使得局部信息缺失。阈值选取依赖于经验值,主要原则为调整阈值大小及邻域半径得到适合于算法的关键点数量。本文将关键点数控制在10000左右,阈值ε1和ε2采用点云库(pointcloudlibrary,PCL)官方文档的参数0.975,提取的关键点直接进入下一步计算,特征点提取步骤如下:
1)为点云中每个点pi建立局部坐标系,设定每个点的搜索半径r。
2)以搜索半径r建立k-d树,计算其范围内所有点的权值:
wij=1/|pi−pj|,|pi−pj|
式中,pj为pi半径r内的邻域点坐标。
3)计算所有点pi的协方差矩阵:
cov(pi)=∑|pi−pj|
4)计算各点pi协方差矩阵的特征值{λ1ii1,λ2ii2,λ3ii3},由大到小进行排列。
5)设置阈值ε1与ε2,满足式(3)则认为该点为ISS特征点:
λ2i/λ1i≤ε1,λ3i/λ2i≤ε2 (3)λi2/λi1≤ε1,λi3/λi2≤ε2 (3)
1.2改进的PointNet模型
PointNet为一种新型的处理点云数据的深度学习模型,可直接在点云数据上应用。该模型在点云分类和任务分割中具有重要作用,但所获取的点云特征未充分考虑邻域点的信息,因此难以直接应用于点云配准。图1为PointNet主要架构,关键流程包括:输入为点云数据的坐标集合;通过网络得到转换矩阵,将原始点云与转换矩阵相乘进行点云对齐;再通过多层感知器(multilayerperceptron,MLP)提取点云数据特征,对特征进行对齐并升维;执行最大池化操作来获取各组点云数据的全局特征;最后针对点云分类或任务分割,分别对特征进行MLP处理。
图1PointNet简略架构图
PointNet处理点云数据具有可提取深度信息的优点,但其在点云配准上仍具有一定挑战:PointNet对于姿态错位较为敏感,旋转不变性的前提是点云相对于规范坐标框架错位较小;需要大量的训练数据来提高抗噪性;原始模型针对不同样本需要等量点数等。
针对PointNet模型无法获取纯粹局部信息的缺陷,将模型输入的点云数据改变为优化的PPF信息,点云PPF信息计算公式为:
F12=(∥d∥2,∠(n1,d),∠(n2,d),∠(n1,n2)) (4)F12=(∥d∥2,∠(n1,d),∠(n2,d),∠(n1,n2)) (4)
式中,d为两点间的坐标差向量,n1和n2分别为两点的法向量,‖d‖2为欧氏距离,角度计算公式为:
∠(n1,n2)=arccos(n1⋅n2/(|n1|×|n2|)) (5)∠(n1,n2)=arccos(n1⋅n2/(|n1|×|n2|)) (5)
由于PPF信息中d含坐标信息,而匹配点对依赖的局部信息应与坐标无关,故只保留PPF中欧氏距离与法向量夹角,优化的PPF信息如下:
F12=(∥d∥2,∠(n1,n2)) (6)F12=(∥d∥2,∠(n1,n2)) (6)
将PPF信息由两点扩展至k个点即可获得局部范围的特征表达,文中模型的输入即为关键点在其k邻域内的优化PPF信息,如第r个关键点的优化PPF信息为:
Fr=(∥d1r∥2,∠(n1,nr),∥d2r∥2,∠(n2,nr)⋯∥dkr∥2,(∠(nk,nr)) (7)Fr=(∥d1r∥2,∠(n1,nr),∥d2r∥2,∠(n2,nr)⋯∥dkr∥2,(∠(nk,nr)) (7)
PointNet模型首先采用训练网络学习从输入的点云坐标中得到转换矩阵,将原始点云与转换矩阵相乘进行对齐,以保证整块点云空间变换的旋转不变性,但无法在配准中提取局部特征,因此将模型中该处的训练网络去除以提高运行效率。PointNet模型中最后获得的特征为1024维,研究结果表明,当维度取128和256时可获得与1024维相当甚至更优的效果[6],故本文将提取特征由1024维改为256维,可提高后续计算的效率。改进后的模型架构如图2所示。
图2改进的PointNet架构图
根据模型得到的特征空间向量计算两站点云间关键点的特征空间向量欧氏距离,设置阈值后可初步筛选出匹配点对。将初步筛选的匹配点对集合记为M(An,Bn),然后通过RANSAC作进一步筛选,求取的内点即为最终的粗匹配点对。根据上述点对计算旋转和平移矩阵,作用于原始两站点云即可得到粗配准结果。RANSAC过程的计算时间与匹配点对集合M的数量正相关,因此需要综合考虑效率与准确度2个因素,本文将点对集合M控制在原始点数的1/30左右,RANSAC迭代上限设置为10000次。
1.3k-d树加速的ICP精配准
粗配准的结果主要是为后续的精配准作准备,因此将本文算法与传统算法粗配准进行对比后,仍需对精配准的结果作进一步比较。考虑到效率因素,选取k-d树加速的ICP算法进行精配准,通过精配准的实验结果反映不同粗配准算法提供的初始化条件的优劣。ICP算法通过最小二乘原理实现最优匹配进行精配准,主要流程为寻找源点云与目标点云间的旋转和平移矩阵,重复寻找点集和计算刚体变换,当2次迭代的目标函数值差异小于设定的阈值或迭代达到一定次数时停止迭代,目标函数为:
f(Rk,Tk)=1N∑i=1N∥a¯¯i−Rkb¯k−1i−Tk∥2 (8)f(Rk,Τk)=1Ν∑i=1Ν∥a¯i-Rkb¯ik-1-Τk∥2 (8)
式中,Rk为旋转矩阵,Tk为平移矩阵,a¯¯a¯i和b¯b¯i分别为源点云和目标点云坐标与各自质心差值后的坐标。
传统ICP算法在确定对应点对,即在目标点云中搜索源点云相应最近点时,每次迭代都需遍历一次目标点云,导致效率低下,因此可采用k-d树对ICP算法进行加速。在迭代之前,预先对两站点云建立k-d树,进而一次性计算源点云中各点在目标点云中相应的最近点。
2、实验与分析
2.1实验数据及环境
实验数据分为3组:第1组为Stanford3DScanningRepository下载的Lucy点云模型,用来测试算法在不同噪声水平及参数下的粗配准效果;第2组为实测的汉武大帝雕像数据,数据由Z+F5016三维激光扫描仪获取,命名为Han,目的是对实测轮廓性较好的对象进行配准实验;第3组为Robotic3DScanRepository下载的实测数据,由RieglVZ-400扫描仪获取,命名为Street,用来测试大规模室外场景下该算法的有效性。
实验笔记本电脑处理器为AMDA10-5745M、内存8GB、Windows10操作系统,算法基于Python3.6、Tensorflow2.0,FPFH和SHOT算法在VisualStudio2013、PCL点云库1.8中实现。
2.2Lucy数据粗配准实验
Lucy数据实验仅测试粗配准效果,将从2个方面对该组数据的配准结果进行评价:1)模型在不同噪声水平下的配准均方根误差(rootmeansquareerror,RMSE);2)与传统配准方法的效率和精度的比较。均方根误差通过查询源点云各点至目标点云中最近点获得,计算公式为:
RMSE=(1Np∑i=1Np∥Rpi+T−qi∥22)12 (9)RΜSE=(1Νp∑i=1Νp∥Rpi+Τ-qi∥22)12 (9)
式中,Np为源点云点数,pi为源点云,qi为目标点云对应pi中的最近点集。
Lucy模型原始数据点云个数为1300万,通过曲率压缩后为20万,再通过ISS检测出9166个关键点,在此数据基础上随机删除1/100、1/50并进行刚性变换,添加递增的高斯白噪声以分成多组数据,得到的待配准点云个数分别为9075和8983。由图3(a)可知,本文算法在随机删除1/100和1/50点的情况下,添加1~20dB强度噪声时均能完成配准,且精度均在2.5mm以内,噪声增强及缺失点数增多均会使配准误差变大。采用本文算法选择不同邻域点数计算法向量夹角及欧氏距离时,最终配准误差会产生约0.1mm的差异,当点数为9075和8983时,选取20个邻域点的配准效果较好(图3(b))。从表1和表2可以看出,在Lucy数据中添加不同强度噪声时,本文算法平均配准时间比FPFH减少61.1s,提高73.2%,比SHOT算法减少53.3s,提高70.4%;平均配准精度比FPFH算法提高89.9%,比SHOT算法提高89.9%。图4(a)为原始点云,4(b)、4(c)、4(d)为本文算法、FPFH算法、SHOT算法对Lucy模型添加1dB噪声并删除1/50点的配准效果。从图中可以看出,本文算法配准后的点云重叠度优于FPFH、SHOT算法。
图3多组实验配准误差对比图
表1不同噪声配准误差对比
表2不同噪声配准时间对比
图4Lucy配准图
2.3Han数据精配准实验
Han两个站点曲率压缩和去噪后均约为20万点云,经ISS特征检测得到关键点数分别为11249和11890,站点间重叠度约为70%,计算PPF信息时选择20个邻域点。从图5(a)和5(b)可以看出,本文算法粗配准相比于FPFH算法能够提供更好的初始状态,以满足后续精配准的需求;图5(c)和5(d)为FPFH算法与本文算法在ICP精配准后的效果,两者在视觉上无明显区别。由表3可知,本文算法配准过程相比于FPFH与ICP算法,运行时间减少100.5s,提高68.6%,RMSE减小3.5882mm,提高28.1%,具有一定的优势。
图5Han配准图
表3Han配准时间及误差对比
2.4Street数据精配准实验
Street两个站点曲率压缩和去噪后均约为30万点云,经ISS特征检测得到关键点数分别为9005和8512,站点间重叠度约为80%,计算PPF信息时选择20个邻域点。从图6(b)可以看出,ISS关键点的分布位置在边缘处较为集中,能够较好地表达场景特征。从图6(c)和6(d)可以看出,本文算法的粗配准效果优于FPFH算法,能保持较好的初始配准姿态。图6(e)和6(f)为2种方法的精配准效果,两者差异不明显,需采用表4中的量化结果进行比较。由表4可知,本文算法相比于FPFH算法,运行时间减少74.4s,提高65.6%,RMSE减少1.2429mm,提高7.8%。
图6Street配准图
表4Street配准时间及误差对比
3、结语
本文对点云配准算法进行深入研究,针对基于特征的配准方法存在误匹配点对多且计算成本高的问题,提出一种改进的PointNet模型代替传统的特征提取和匹配方法,相对于传统方法可大幅提高处理效率,同时由于提取过程中信息利用率较大,可在一定程度上减少误匹配问题。
研究结果表明,本文算法在粗配准过程与FPFH、SHOT算法相比,可减少配准时间,同时提高匹配点对的正确率,为后续的精配准提供更好的初始状态,有效减少精配准匹配失败的情况。实验数据表明,本文算法相比于FPFH与ICP算法,运行时间平均减少67.1%,均方根误差平均减小17.95%,能够高速且有效地完成物体配准。
该算法的不足之处在于整个配准流程涉及多种算法的融合,其中参数和邻域点的选择较为关键,需要一定的人工干预,因此对算法的自适应性作出改进将是下一步的研究重点。
参考文献:
[3]张崇军,许烨璋,郑善喜,等.改进权重的迭代最近点算法在点云配准中的应用[J].大地测量与地球动力学,2019,39(4):417-420.
[6]白静,司庆龙,秦飞巍.轻量级实时点云分类网络LightPointNet[J].计算机辅助设计与图形学学报,2019,31(4):612-621.
陈强,岳东杰,陈健.基于特征空间匹配的激光雷达点云配准算法[J].大地测量与地球动力学,2020,40(12):1303-1307.
分享:
近年来,各种分类或聚类算法应用于地球物理资料处理和反演,主要包括模糊C均值(FCM)聚类[1]、支持向量机(SVM)[2-3]和神经网络(NN)[4-5]等。其中,SVM和NN算法均属于监督学习分类方法,需要足够数量的样本维持结果的稳定性;而FCM聚类算法是一种无监督学习方法,具有样本量需求少的优势[6-7],更适用于地球物理方法联合反演。
2024-08-05咸潮是发生在滨海河口地区的一种水文现象,它是由于太阳和月球对地表海水的引力作用,海水倒灌进入内陆河道,导致水体变咸。咸潮危害巨大,例如影响水厂供水安全及部分水生生物的正常生存。因此,加强对咸潮的监测,并据此及时制定应对措施,对当地防灾减灾和城市供水保障等具有重要意义。
2024-01-02在新时期背景下,对国土空间规划功能进行准确定位,全面落实各项规划方案,不仅可以保护周围生态环境,而且能够实现统筹治理目标。结合现阶段我国国土空间规划功能定位现状可知,环境修复难度不断增加,各项管控措施落实效果较差,严重影响生态环境质量。鉴于此,本文重点分析新形势下国土空间规划功能定位和实施途径。
2021-06-01日益严重的环境污染、自然资源无节制开发,不仅开始制约人类社会的发展,同时威胁人类的生存和繁衍,如何改变这种情况,成为当今人类亟待解决的问题。基于卫星遥感的国土空间监测的发展与应用可以帮助人类全面调查全球自然资源和生态环境,确保人类更合理的开发、分配、利用自然资源,使实现可持续发展成为可能。
2021-05-04文章围绕着国土空间规划中对地理信息大数据的有效应用展开充分讨论,首先简要介绍了地理信息大数据的概述,然后表明了国土空间规划的特点,最后着重探讨了地理信息大数据在国土空间规划中的具体应用,包括提供合理数据信息、优化国土空间规划、强化信息的数据源、明确空间具体位置,以此完善地理信息的内容,并形成重要的数据,促进国土空间规划的良好发展。
2021-04-27在当前这一社会与经济快速发展背景下,地震预警已经成为了人民群众关心的问题,因为地震为社会与经济带来的威胁远超于其他地质灾害,再加上地震的前期预兆不明显,所以就需要利用合理的方式进行监测预警,进而有效提高防控效率,降低灾害为人民群众生命财产带来的威胁。那么在当前这一社会背景下,信息技术成为了各行业广泛应用的技术之一。
2021-02-19受测量条件限制,激光扫描技术往往需要架设多站点进行扫描,并对获取的点云数据进行配准以消除不同站点的相对位置变化。目前点云配准算法主要分粗配准和精配准2个阶段。迭代最近点算法(iterativeclosestpoint,ICP)[1]为目前应用最广泛且十分有效的一种精配准方法,但计算成本大、易陷入局部最小值。
2020-12-05高铁列车可当做一种全新的震源类型,其确定的长度和荷载,固定的行驶路线以及近乎匀速运动的特性,使之具有可重复震源的特征,另外,地震仪记录的高铁震源记录,具有宽频带分立谱的特征,为高铁高架桥附近地表结构和物性高精度探测提供了条件.因此,近年来,以高铁列车震源为基础的高铁地震学已成为了一个新的研究热点.
2020-12-04第三极不同地区(如季风区和西风影响区)和不同冰川表面(如老雪区、新雪区、裸冰区)雪冰中高分辨率的碳质组分的观测与分析势在必行,特别是消融期开始后碳质组分在短时间尺度内的浓度、黑碳粒径以及有机碳组成变化对于模型的改进至关重要,可为我们全面认识碳质组分对青藏高原不同地区冰川消融影响的机制提供科学依据。
2020-10-24喜马拉雅地区的冰川是南亚地区许多大江大河的源头,如印度河、恒河和布拉马普特河(上游中国境内河段称雅鲁藏布江)。喜马拉雅高大的山体阻挡来自印度洋和阿拉伯海的湿润水汽,南坡温暖湿润,北坡寒冷干旱,剧烈的海拔梯度差异使得在短距离空间内水热急速变化,为各种动植物生存繁衍提供了场所,成为全球生物多样性最丰富的地区之一。
2020-10-24人气:3316
人气:2758
人气:2463
人气:2148
人气:1951
我要评论
期刊名称:大地测量与地球动力学
期刊人气:1938
主管单位:中国地震局
主办单位:中国地震局地震研究所,地壳运动监测工程研究中心
出版地方:湖北
专业分类:科学
国际刊号:1671-5942
国内刊号:42-1655/P
邮发代号:38-194
创刊时间:1981年
发行周期:月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:0.435
影响因子:1.417
影响因子:0.700
影响因子:0.925
影响因子:1.760
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!