基于结构关联的三维点云配准研究

基于结构关联的三维点云配准研究

ID:10310770

大小:8.77 MB

页数:72页

时间:2018-07-05

上传者:文档小小白
基于结构关联的三维点云配准研究_第1页
基于结构关联的三维点云配准研究_第2页
基于结构关联的三维点云配准研究_第3页
基于结构关联的三维点云配准研究_第4页
基于结构关联的三维点云配准研究_第5页
资源描述:

《基于结构关联的三维点云配准研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

学校代码:10004密级:公开又4乂攀BEIJINGJIAOTONGUNIVERSITY硕士学位论文基于结构关联的三维点云配准研究作者姓名张哲学科专业计算机科学与技术指导教师许宏丽教授培养院系计算机与信息技术学院f二零一八年四月 硕士学位论文基于结构关联的三维点云配准研究Registrationof3DointcloudsbasedonStructuralAssociationp作者:张哲导师:许宏丽北京交通大学2018年4月 学位论文版权使用授权书本学位论文作者完全了解北京交通大学有关保留、使用学位论文的规定。特授权北京交通大学可以将学位论文的全部或部分内容编入有关数据库进行检索,提供阅览服务,并采用影印、缩印或扫描等复制手段保存、汇编以供查阅和借阅。同意学校向国家有关部门或机构送交论文的复印件和磁盘。学校可以为存在馆际合作关系的兄弟高校用户提供文献传递服务和交换服务。(保密的学位论文在解密后适用本授权说明)导师签名:学位论文作者签名:m?曰期年3月曰签字曰期?M年芩月方曰签字豸V :学校代码:00410密级公开北京交通大学硕士学位论文基于结构关联的三维点云配准研究Reistrationof3DPointCloudsBasedonStructuralAssociationg?15120470作者姓名.张哲学号::导师姓名许宏丽职称:教授学位类别:工学学位级别:硕士学科专业:计算机科学与技术研究方向:计算机视觉北京交通大学2018年4月i 致谢,至今仍记得刚上研,两年多的研宄生求学生活即将结束岁月如梭,转眼间、自己的担心,感觉不知道如何下手同学宄生确定研究课题时。但是在各位老师一一的关心和帮助下,我慢慢走近课题研宄内容,我才能步,是在他们的鼓励下步完成研宄生毕业设计。值此毕业论文完成之际,我谨向所有关心、爱护、帮助我的人们表示最诚挚的感谢与最美好的祝愿。本论文是在导师许宏丽教授的悉心指导下完成的。两年多的时间里,严谨的治学态度,工作作,许宏丽老师渊博的专业知识精益求精的风,诲人不倦的高尚师德,朴实无华、平易近人的人格魅力对我影响深远。许老师不仅授我以文,,,而且教我做人虽历时三载却赋予我终生受益无穷之道。本一论文从论文的选题、修改直至定稿,倾注了老师大量的时间和精力,每步都是在导师的指导下完成的。至此论文完成之际,我向我的导师许宏丽教授表示深切的谢意与祝福。在本文研究过程中,,我还要感谢实验室的同学包括侯琳琳、刘静静、白宏运、高佳月等,他们在研宄过程中对我的帮助和支持非常多。另外也感谢我挚爱的爸爸妈妈,他们的理解和支持是我学习上最大的动力。正是他们的支持,才使我能够在学校专心完成学业。本文在国家自然科学基金No.6147202961473031、科技部国家重点研发计划(,)(No.2016YFB1200100)、北京市自然科学基金No.4152042和中央高校基本科研业()务费专项资金@〇.2016见20052016圯^1019,2016阳1^016的资助下完成。,)V 北京交通大学硕士学位论文摘要摘要随着三维点云处理技术的快速发展,三维点云处理技术被广泛应用于各项研究领域中。在实际工程中由于检测物体易受物体遮挡、环境、表面布局和获取方式等因素的影响,需要多次从不同角度获取点云。因此,要将不同视角的点云完一一成配准,个整体形成,选取种好的配准方法尤为重要。本文综合考虑点云配准流程中每个步骤对配准结果精度和效率的影响,从不同角度出发,保证配准精度,提高配准效率。本文主要工作如下:1本文对点云特征描述子做了研究,。针对点云数据量大,特征描述能力弱()特征计算效率低的问题,本文提出了基于法向量分布的特征描述子和基于统计特征直方图的特征描述子。其中,基于法向量分布的特征描述子计算效率快,能快速提取点云中的特征点;基于统计特征直方图的特征描述子综合考虑点云法向量分布统计特征、深度统计特征和密度统计特征,能够全面描述点云的结构关联特性,,。实验表明使用基于法向量分布的特征描述子提取出的特征点具有较高的鲁棒性;使用基于统计特征直方图的特征描述子,能够准确描述探寻点所在区域内点的结构关联关系。2本文对基于特征描述子的配准算法做了研究,提出了基于统计特征直方图()的点云初始配准算法和基于ICP的精确配准的两步配准算法,。在初始配准中针对初始点云数据过大问题,提出了可变边长的体素网格算法,在保证点云质量的基础上,,精简点云。针对精简点云中存在的特征相近的点影响配准精度的问题云的关键点集本文提出了基于法向量分布特征的关键点提取算法,提取出每个点合,。针对提取关键点后细节的丢失问题使用统计特征直方图来增强点云的描述能力。本文结合RANSAC算法的基本思想对两视角的关键点完成初始配准。对初始配准的结果,使用ICP算法实现最终的精确配准。实验表明,本算法与ICP算法相比,能够对任意初始位置的点云完成正确配准。与使用FPFH特征的配准算法相比,,提高配准速度本算法能够在相近的配准精度下。3实现了点云三维重构系统的设计与开发。系统运用本文提出的两步配准算()法,完成高精度配准elauna算法完成最终的三维重构系统界面分别;并使用Dy;对配准结果和重构结果可视化。关键词:点云配准;;;多视角点云精简关键点;特征描述子iii 北京交通大学硕士学位论文ABSTRACTABSTRACTttoWiththeraiddevelomentof3Doincloudrocessinechnolo3Dintppppggy,pcloudprocessingtechnologyhasbeenwidelyusedinvariousresearchfields.Duetothedetectedobjectsinfluencedbyobjectsenvironment,surfacelaoutandacuisition,yqmethodsinactualroectitisnecessartoobtainointcloudsfromdiferentanlespj,ypgmanytimesThereforeinordertocomletethereistrationofointcloudswith,pgpdiferentersectivesiisarticularlimortanochooseaoodreistrationmethod.pp,typttggpIntheaerweconsidertheinfluenceofeachsteintheointcloudreistrationpp,ppgprocessontheaccuracyandefficiencyoftheregistrationresults.Besides,thereistrationtimeisimprovedonthebasisofensurinthereistrationaccuracyfromgggdifferenterspectives.Themainworkofthispaerisasfollows:pp-Firstltheointcloudstructurerelatedfeaturedescritorhasbeenstudied.Aiminy,ppgattheproblemoflargeamountofdatainpointcloud,weakfeaturedescriptionabilityandlowfeaturecomutinefficienctheaerresentsafeaturedescritorbasedonpgy,ppppnormalvectordistributionandafeaturedescritorbasedonstatisticalfeaturehistogram.pAmongthemtheformerhashihcomutationaleficiencandcanuicklextractthe,gpyqyfeatureointsintheointcloud.Thelattercomrehensivelyconsidersthestatisticalpppcharacteristicsofointcloudnormalvectordistribution,depthstatisticalfeaturesandpdensitystatisticalfeaturesandcancomprehensivelydescribethestructurecorrelation,characteristicsofointclouds.Exerimentsshowthatthefeatureointsextractedbasedppponthenormalvectordistributionhavehighrobustness.Usingthefeaturedescriptorsbasedonthehistogramofstatisticalfeaturescanaccuratelydescribethestructuralrelationshisbetweenheointsinthesearcharea.ptpture-tSecondlythereistrationalorithmbasedonstrucrelatedfeaturedescriors,ggphasbeenstudied,andputforwardthealgorithmthatapointcloudinitialregistrationbasedonstatisticalfeaturehistogramSFandanaccurateregistrationbasedonICRIn()theinitialreistrationaiminattheroblemoftoolareinitialointclouddata,g,gpgpavariable-edevoxelraorihmroosedtoensureheualiofoincloudandggidlgtispptqtytpsmtecoudminattheroemthatteointsnearthefeaureointsinhiplifyhl.Aigpblhptptestreamlinedointaffectthereistrationaccuracyakeointextractionalorithmbasedg,ypgpiv 北京交通大学硕士学位论文ABSTRACTonthenormalvectordistributionfeatureisroosedinthisaer.Atthesametimethepppp,statisticalfeaturehistoramisusedtoenhancethedescritionabilityoftheointcloudgpptosolvetheroblemoflossofdetailaftertheextractionofthekeyoints.WeborrowpptheideaofRANSACalgorithmtocompletetheinitialregistrationofthekeyointsofpthetwoviewsorheiniialreistrationresulthefinalreisationisdoneusinthe.Fttgttr,ggICPalgorithm.ResultsshowthatcomparedwithICPalgorithm,thisalgorithmcancorrectlyregisterthepointcloudinanyinitialposition.Inaddition,theproposedalgorithmcanimprovetheregistrationspeedwithsimilarregistrationaccuracyncomarsottereaonFHreipinohgistrtialgorithmusingtheFPfeatu.,BasedonVisualC++6.0andQT5.6.3developmentplatformthe3D,econmmcan-rstructionsysteisdesignedandimlemented.Thesysteusethetwosteppretrtonttotettgisaialorithmroosedinheaerfinishhointcloudofhecddaagppppppmat-consforhighreisireitrationthentheDelaunaalorithmisusedtorealizethepg,yg-andthreedimensionalreconstructionandtheresultsofthekeointsextraction,ypreconstructionaredisplayedonthesystemae.pg-KEYWORDS:PointCloudReistrationMultiviewPointCloudSimlificationKeg;;p;yPointsFeatureDescritor;pV 北京交通大学硕士学位论文^目录觀iiiABSTRACTiv1雜3U课题背景及研宄意义31_2三维点云配准的研宄现状61.3论文主要内容及组织结构81.3.1论文的主要内容81.3.2论文的组织结构92基于结构关联的点云特征描述子研宄102.1相关工作102.2基于法向量分布的特征描述子162.3基于统计特征直方图的特征描述子182.4实验结果与分析212.4.1基于法向量分布的特征描述子实验结果与分析212.42基于统计特征直方图的特征描述子实验结果与分析23.2.5本章小结263基于特征描述子的配准算法273.1相关工作273.2算法流程303.3基于SF的点云初始配准算法313.3.1基于可变边长的体素网格点云精简算法313.3.2基于NV的关键点提取算法353.3.3基于SF的初始配准算法363_4基于ICP的精确配准算法383.5实验结果及分析:393.5.1散乱点云数据配准实验393.5.2自获取深度点云数据配准实验443.6本章小结471 北京交通大学硕士学位论文4三维重构演示系统494.1基于Delaunay的三维重建494.2三维重构系统软件总体设计504.2.1三维重构系统软件搭建环境504.2三维重构系统软件设计思路1.254.3三维重构系统软件操作流程524.4本章小结545工作总结与展望555.1工作总结55一5.2进步的研究建议56参考文献57作者简历及攻读硕士学位期间取得的研究成果62独创性声明63学位论文数据集642 北京交通大学硕士学位论文^1绪论本章阐述了所研宄课题的研宄背景和研究意义,并对三维点云配准研宄领域内的国内外研究现状进行综述,介绍本文的主要内容及组织结构。1.1课题背景及研究意乂在逆向工程中,通过测量设备得到的物体表面点的数据集合称之为点云。点一组数据的集合可以存储每个点的三维方向、位置坐标值、法向量等云是,不仅一、的几何信息,还可以将被测物体的表面的诸如颜色纹理及透明度等属性并存储,。现有的点云获取的方式呈现多样化趋势针对不同应用可以选择相对应的点123[][][]云获取方案-。,如图11所示三维点云数据获取技术II1接触式非接触式fi?謇罾手持式三维标测量光学声学电磁学*z可见光不可见光£IIII干涉测量图像分析测距法X射线红外线-图11点云数据获取方式分类F-i11Theointclouddataacuisiedclaifiigitonmthosscatonpq在实际应用中,常采用扫描设备获取点云,不同扫描设备获取的点云千差万别,根据能否确定邻近点位置,点云数据可以分为无序点云和有序点云。无序点云可以根据点云数据之间是否存在拓扑关系或者组织关系分为多类,如散乱点云、扫描点云、矩阵点云和多边形点云等。其中,散乱点云是指整个点-2云数据是散乱、无序的,点与点之间没有明显结构关系。如图1(a)所示;扫3 北京交通大学硕士学位论文绪论描线点云是指点云排列分布整体呈扫描线形态,通常需要高精度的扫描设备来确-保曲面点云获取质量和获取效率。如图12(b)所示矩阵点云相较于前两种,;一个坐标点整体上分布接近于矩形阵列,矩阵中每个顶点对应点云数据中的,通常是经过网格化插值计算获得-)。如图12(c所示;多边形点云不同于前几种存一,多边形点云中的所有点存储在同个平面上储方式,点之间通过连接形成多个-)。平面多边形,点云结构如图12(d所示??..*#??*???..KrTyK广1爾a)散乱点云b)扫描线点云c)矩阵点云d)多边形点云3图-2[]1无序点云的分类3F-[】i12ThecifiiiPglasscatonofdsorderedointClouda)ScatteredPointCloudb)ScanlinePointCloudc)MatrixPointCloud(d)PolygonPointCloud有序点云指点云通过矩阵进行组织,数据由杂乱无序进而被划分为行与列,如从立体相机或TOF相机等设备获得的点云数据。该类点云的优点为每个点云能够快速找到相邻点云的位置,最近邻节点查找效率高。此外,点云按其点的几何分布稠密度可以分为稀疏型点云和密集型点云。通常稀疏型点云在几十到几千点数据量之间,结构信息较弱,在不增加其他约束的情况下,很难处理;密集型点云数据量大,典型的密集型点云为三维激光设备捕获的点云,该类点云蕴含丰富的三维点特征,具有超强的描述能力。三维点云处理技术己成为计算机图形学中的研宂热点,其优越性主要体现在如下几个方面。首先,三维点云数据能够精确记录物体表面结构信息,直接存储空间点的多维信息,描述能力大大高于平面点云。其次,,能够全视角展现物体信息对捕获到的一个视角的点云数据,通过对点云数据的平移、旋转和缩放操作,使人们可以全方面的查看物体结构。与二维图形相比,通过空间变换避免了图像中物体的遮挡现象。三维点云处理可以获取局部的目标线、面、体、空间等三维数据,最终建立目标体的三维数字模型给人们带来视觉上的真实感^,可以构造虚,体验三维物体的视觉效应拟物体形象。在三维点云处理技术中,高精度快速点云配准方法是其研究热点和难点,在45[][]计算机视觉、逆向工程等领域,均有着十分重要的意义。在获取物体点云数字化模型过程中,需要在物体完整点云数据的基础上进行4 北京交通大学硕士学位论文绪论,重建但在实际应用领域中,由于测量设备的差异和被测物体大小等因素的影响,一无法通过次测量就获得物体的完整点云数据,需要从多个视角对被测物。因此一体或空间场景进行捕获,以得到多个在相邻视角之间有定的重叠区域的局部点云数据,利用点云配准技术对上述的不同点云数据处理,进而得到物体的准确且完整的点云数据。然而,从多个视角获得的点云所在的坐标系不同,不能将得到的多个视角的点云简单地叠加在一起一。点云配准过程就是力求求解个最佳的旋转和平移变换矩阵一,将多个从不同视角下捕获的局部点云数据整合到个相同的一坐标系下,通过对不同点云之间重叠部分的匹配,近似配准出满足定精度要求的完整点云的过程,三维物体坐标信息的获取已经变。随着三维扫描技术的发展得更快速更精确,而三维点云配准算法成为三维配准与点云重构技术中的关键环,精确配准后的数据才能构成整体模型节,应用于更多领域中。因此,点云配准是后续工作的基础,,点云配准的精度决定了后续重建的精度以及测量的准确性对物体三维模型的精度具有直接的影响。在制造业中,利用配准算法将生产的实物模型与标准三维模型比对,以逆向工程为主,,检测工件是否达标,改善了质量检测效率和精确度。在建筑领域中可以借助大型扫描设备的多个视角获取建筑物的局部区域,通过配准技术将局部6[]。区域构建完整的建筑物模型,完成建筑物的建模点云配准从建,用于破损检测78口[][]筑港和路桥隧道等工程检测,到地形地质研宄与测量等应用,使工程检测过。程更简捷、更迅速以及更准确,综合效益得到显著提高此外,三维配准技术在91()11[][][]目标识别、医学、图像匹配等领域中也广泛得到应用,具有无与伦比的优越性。12三维点云配准也促进了虚拟现实[]、电子商务无人驾驶等诸多方面的发展。例如,激光雷达产品可以自动对车速进行控制,,作为主动巡航控制系统的核心其技术是否成熟直接决定了未来无人驾驶汽车能否走出理论研宄,被应用到曰常生活中。谷歌对无人驾驶汽车的研制中,己经使用了激光雷达探测器来探测车况,将多视角下获取的雷达点云数据精确配准,对汽车周围环境实现三百六十度无死角感应,动态形成运动汽车周围环境的实时监测,执行内部指令实现汽车平稳安全行驶。综上所述,点云配准技术具有重要的研究意义和丰富的应用价值。目前,现有配准算法难以满足现实中人们的应用要求,在效率和精度两个层面上都有很大?的改进空间,现有的配准算法都i各自的针对性和局限性,所以导致其中。此外一多数方法在实现应用中受到了极大的限制。因此,如何获得种适用性广泛且精确度高的配准算法成为研宄的主要方向。5 北京交通大学硕士学位论文绪论1.2三维点云配准的研究现状现阶段点云配准方法得到学术前沿的普遍关注,按配准原理可分为迭代计算最小误差配准、基于数学模型变换配准和基于特征寻求对应关系三类配准方法。一第类是迭代最小误差法。其以迭代最近点算法为基础,在不同阶段有很多不同的改进算法。迭代最近点算法(IterativeClosestPoints,ICP)是经典的配准算14[]等在法,由说31992年首次提出。ICP算法的简单易懂,且精度较高,但由于一需要多次迭代,该算法的效率较低,且对点云数据具有较高的限制条件。其,一待配准的两个点云之间要有定重叠区域,否则会存在较多错误的对应点对,影响配准结果,不能相距;其二,待配准的两个点云之间要有相对合理的初始位置,否则可能导致算法的收敛方向不准确太远,致使配准失败。上述两个条件不只保证了该算法计算的有效性,同时也极大地限制了该算法的应用范围。因此,多位学者对原始ICP算法的各个阶段进行研宄,提出了很多改进方法。15[]乙11〇1^等在原始ICP算法的基础上,提出欧氏距离阈值去噪和点云方向矢量夹角阈值两种改进方法,提高配准精确度和速度。然而这种方法存在不可回避的问题,Zhang的配准环境太过理想化,在实际操作扫描过程中很难做到两个点集完一,同时,达到预期的配准计算速度对于硬件配置来说也是个很大的挑战全相同。1617[][]Chen等提出了HT-ICP算法,有效的去除了错误点对,提升配准精度。Yang等将BnB算法与ICP算法结合,通过重新划分旋转空间和平移空间来加快对应点的查找速度。然而,此种ICP改进算法在两片点云初始位置相差太大的情况下,很容易陷入局部最优解,所以在实际应用中常在二次配准时使用。二一第类是基于数学模型变换法。基于数学模型变换法主要包括随机采样致性(RandomSampleConsensus,RANSAC)算法和基于统计计算刚体变换算法。一^一,致性算法最早是由Fchier其中随机采样is提出,是种从样本中准确拟19[el合出数学模型的方法。Chii首先将RANSAC思想应用于点云配准中,由于RANSAC算法对存在大量外点的数据集也同样适用,因此这种配准方法对重叠区2Q域较小的两个点云的配准具有较高的鲁棒性[]。Papazov提出了基于RANSAC的一四点法,在每次采样中,随机选取同平面上的四个点,这四个点的对角线的交点将两条对角线划分为四条线段,依据在刚体变换下这些线段的长度和比例不变特性来减少搜索对应点的时间,降低算法的时间复杂度,在噪声较多的点云配准'21[】中也能使用。Meng提出了基于RANSAC思想的采样球配准方法,该方法通过一建立采样球模型,提高对应点对的查找效率,,进步降低算法的时间复杂。综上,而且对噪声的鲁棒性较高基于RANSAC的配准方法不必计算复杂的几何特征,但确定对应点对的过程相对比较复杂,当点云的特征非常明显时,适当加。通常6 北京交通大学硕士学位论文绪论入一些几何特征约束来减少搜索对应关系的时间复杂度。这类算法通常用于大型复杂点云数据的初始配准过程。基于统计计算刚体变换算法原理是将空间离散化,定义配准误差评价标准,22使用标准最优化技术确定点云变换【]提出的法向量分。典型的算法是\1叫1111%011布变换(NDT)算法,该算法的缺点是在初始参数误差较大的情况下,会使配准23[]陷入局部极小值-。胡修祥等将NARF特征与3DNDT特征结合,实现对三维场24[]-NDT算法景的快速高精度配准3D,实现了对百。骆林等通过融合RGB影像与万级大型复杂溶洞点云数据的快速高精度配准。这类算法目前常用于大型复杂点云数据的二次配准。第三类是基于特征寻求对应关系方法。基于特征寻求对应关系的配准方法根据源点云与目标点云中对应点特征的相似性,寻求两幅点云中特征相同的点与点的对应关系,根据对应关系计算从源点云到目标点云的坐标变换。这类方法的核心是选取一个可区分性强一、不变性强的特征描述子。点的几何特征般包括点的法向量、曲率、积分不变量、矩不变量、结构关联描述子、旋转图像及特征直方图等等。上述的几何特征均为刚体变换下的不变量。其中最常见的是利用法向量和曲率结合来确定两个点云之间的对应关系。25基于此e[],学术界对点云特征描述子进行了更深入的研宂。Ba等提出利用点26[]的曲率变化并增加对应点对间的法向量夹角约束来匹配点对,朱延娟和徐金婷27[]也采用该方式得到点云对应关系,不能保,但由于曲率和法向量是低维描述子证对应关系的准确性,算法的精度较低。梅元刚等利用点的旋转图像来搜索确定对应关系,首先利用低维描述子确定候选对应点对,再根据旋转图像查找正确的对应关系。该方法利用低维描述子来减少高维描述子对应点的搜索区域,提高了对应关系的准确性和计算效率[M。Rusu等同样利用特征直方图查找对应关系完成点云配准,由于特征直方图具有不受点云位置、方向及点云密度的影响的特点,所以该算法具有较好的鲁棒性一。基于几何特征的配准算法般适用于具有明显特3()[]征的点云,11,否则使用特征确定对应点易出错导致配准错误。了〇^〇11等提出了-旋转图像描述子(sinimae)该描述子是使用最为广泛的三维形状配准描述子pg,;332tnt]Rusu等分别提出了点特征直方图描述子(PFH)和快速点特征直方图描述子33[](),FPFH。其中FPFH和Tombari等提出的直方图方向特征描述子(SHOT)计34算速度快Guo[l,特征描述维度低;等提出的旋转投影统计局部描述子(R〇PS),一35步提高了三维点云配准的精度[]进;\31^等提出了局部特征统计直方图36][(LFSH),该描述子对于场景点云的配准,可以获得较高的配准精度。陆军等37一提出了FPFH的改进算法[],定程度上降低了配准时间。0^11等提出了区域方38一向包围盒描述子t],使用全局特征实现两片点云的精确配准。陈茂霖等提出了7 北京交通大学硕士学位论文种无附加信息的点云自动配准方法,根据地面激光点云PFH描述子各维度的正态分布统计特性来简化点云,,并使用香农熵选择不确定性提炼配准点集替身配准精度。基于特征寻求对应关系配准方法结果可信度高。基于特征寻求对应关系配准方法最为充分的描述每个点云的区域特征,具有更好的可区分性。因此,本课题选用基于特征寻求对应关系的方法来实现多视点云的配准。其中,针对特征提取,本课题使用探寻点与邻近点的空间结构关联信息,构建点云的基本结构约束,利用结构特征构建特征描述子,实现点云的配准研究。1.3论文主要内容及组织结构1.3.1论文的主要内容针对在多视角下捕获点云数据的配准问题,分别对配准的各个阶段进行了深入的研宄,在分析各个步骤实现方法的基础上,主要利用探寻点与邻近点的空间结构关联特性构建描述子,实现高精度的关键点提取和特征描述,从而达到快速、准确的对不同视角的点云数据配准的目的。本文主要内容包括:通过点云空间结构关联关系,构建点云的基本结构约束,利用结构特征构建综合特征描述子,从而准确的描述探寻点云与邻接点云的位置关系。本文提出了基于法向量分布的特征描述子,并将法向量分布统计特征与点云深度统计特征和点云密度统计特征相结合,形成基于统计特征直方图的特征描述子,相对充分的描述点云结构关系,通过实验,验证本文所提出的描述子具有较好的描述能力。本文提出的两步配准法,,在点云初始配准过程中的点云预处理阶段利用可变边长的体素网格对点云精简,考虑到可变边长体素网格下采样对点云模型形状一有限的敏感程度和点云特征描述子的高维计算问题,通过再次精简点云数据进步提升效率,定义点云。本文提出基于法向量分布特征描述子的关键点提取算法的邻接区域和法向量分布特征计算模型,使用该模型提取关键点。本文提出基于统计特征直方图描述子的初始配准算法,对提取的关键点计算统计特征直方图描一述子,使用采样致性配准算法去除错误匹配点,准确求解变换矩阵,完成初始。。配准阶段在点云精确配准过程中,使用ICP算法求取最佳变换本文使用VisualC++6.0和QT5.6.3开发平台,PCL、FLANN、BOOST、VTK等第三方依赖库,,实现三维重构系统该系统可实现点云关键点提取结果和重构结果的可视化。8 北京交通大学硕士学位论文1.3.2论文的组织结构本文共分为五个章节,各章节的内容如下:一第:本文绪论,简要介绍了课题的研究背景及研究意义章,点云配准国内外研究现状以本文的主要内容和组织结构。一第二章:本文提出种基于法向量分布的特征描述子,在原始点云中,使用点云法向量分布特征描述点云与邻近点的结构关系,较好的描述点云基本空间结构。综合考虑法向量分布统计特征、密度统计特征、深度统计特征三种结构分量一变化对配准误差的影响,从,构建基于统计特征直方图的特征描述子而进步更全面地对点云空间结构进行描述。对提出的以上两种特征描述子的描述性进行实验分析。第三章:在第二章的基础上,将提出的结构特征描述子运用于点云配准中。在点云预处理阶段,提出使用可变边长的体素网格点云精简算法,并对精简的点云使用法向量分布特征描述子来提取关键点,弥补可变边长体素网格方法对点云形状的不够敏感的缺陷,保留关键点。对提取的关键点使用基于统计特征直方图一的特征描述子完成后序配准,步提升配准效率进。通过实验证明算法的有效性。第四章:本章使用VisualC++6.0和QT5.6.3开发平台,利用PCL及FLANN等第三方依赖库,实现完整的三维重构演示系统的设计和开发。实现了本文第三章提出的配准算法,并使用Delaunary三角化实现最终的三维重构,系统界面分别对配准前点云的三维结构图、点云的NV关键点三维结构图、点云的配准结果三维结构图和基于Delaunary重构的重构结果图实现了可视化。一第五章:本章对全文总结,说明主要研宄成果,提出下步的工作及展望。9 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究2基于结构关联的点云特征描述子研究本章主要研宄点云的特征描述子,以现有点云的低维特征为基础,参考点云的,高维特征定义方式构建基于空间结构关联关系的结构约束,利用结构特征构建综合特征描述子,从而准确的描述探寻点云与邻接点云的位置关系,提高每个点的描述能力,对后续点云配准具有重要意义。本章的主要工作如下:1、基于法向量分布的特征描述子的定义和意义。2、基于统计特征直方图的特征描述子的定义和意义。3、通过实验,初步验证基于法向量分布的特征描述子和基于统计特征直方图的特征描述子两种描述子的描述能力。2.1相关工作在不同视角下采集的两片点云配准过程中,如何实现重合部分的自动识别是一394()[][],,实现难度也最大。常见的方法是标志法决定配准成果的关键因素之同时,一即在物体表面附加层具有特殊纹理或者颜色信息的标志,通过识别特殊标识来实现配准一步多片点云之间的配准问题就简化成为计算参数使多个点云中。进,标志的重合的问题,,对。但标志法的局限明显首先,被测物体必须允许贴标志,珍惜文物,或者流态表面等都是明显不适用的由于附加标志必然会遮挡被测物一一体的部分表面特征,配准完成后需要通过算法补充缺失部分的信息类常。另41[],当计算机展示多个待配准点云时见方法是通过人工指定,通过人工分别向两片点云的图片输入重合范围,或者重合点,计算机再使用具体算法从选定的范围进行计算,对配准过程中的计算效率有着明,这种算法由于需要过多的人工干预显的制约。上述两种方法虽然都可以有效的缩小两片点云直接的匹配范围,但是“”都没有从根本上解决如何完全自动处理的问题,利用特征描述子建立。综上曲面间的对应关系,在曲面识别与重建过程中有非常广泛的应用。点云的特征描述子是指探寻点与邻域的结构关联特征,结构关联特性常使用相邻规则来约束,经典描述子的3D结构关联描述子具有辨别力强、结构简单、低敏感度等特点。一点云的3D特征描述与计算是点云信息处理中至关重要的步,特征描述的好,,,配准坏对后续的步骤包括曲面识别分割,重采样等产生很大影响。特征描01 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究述子选择的优劣,直接决定着配准最终结果的好坏。从几何特征描述范围来讲,,包括全局特征和局部特征,曲较复杂的包括由,简单的局部法线率等几何特征多个维度特征组合形成的特征描述等。424344[][][]点云的几何特征主要包括低维特征描述子,例如:法向量、曲率、线4531性规划[]及积分不变量等[]。高维特征描述子包括形状内容描述子、旋转图像、pfh32[]等及fpfh。低维特征描述子只需要很少的几个参数值就可以表示邻近的几何特征,虽然估算过程简单高效,却难以完整地表征点邻域内的几何信息,容易受到噪声干扰,点云中仍会存在许多低维特征描述子相似的点。高维特征描述子包含丰富的点的邻域结构关联关系,例如,FPFH特征描述子使用33维的直方图来描述点的邻域内的特征信息,这些特征不仅与点的三维坐标有关,同时也与法向量,但估算过程中时间复杂度高。下面概括说明几种常用的特征描述子有关。31m[]1点特征直方图(PointFeatureHistoras,PFH)()gPFH通过参数化探寻点与邻域点之间的空间差异,构建多维直方图。直方图中存储探寻点与区域内邻近点法线方向之间的两两相互关系,表示物体表面变化一情况-。PFH特征的数学模型如图21所示,表示的是任意个探寻点/V在其局部区域内,与邻近点间的结构关系,并以红色点A为球心。Pn5翁〇-312[]图1PFH数学模型31-FePFma】ig21ThHmatheticalmodel定义固定的局部坐标轴,计算P,和乃两点之间的法向方向偏差,其中N,表-示八的法向量,,凡表示八的法向量,如图22所示11 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究=WUxV皆|wt-31图22局部坐标系的定义[]31-1Teefna[Finecacoordig22hditionofthlolitesstemy-==x-P/—=图22中wxva,w?vw/wuvw坐标中,将法向间的s,[(乃s)(,||f^)]一-关系定义为个三元组,如式(21)所示:a=v奶-—(21)6=arctanw7w7g,gr)(其中,d是两点P,和P的欧氏距离。,之间-计算的特征示意图如图23所示。供'.........,愈二、__顯一氣31图2-3点云PFH特征示意图[】31F2-H[1reig3TheointcloudPFfeatusketchmappm32[]2快速点特征直方图(FastPointFeatureHistograms,FPFH)()一FPFH特征描述子,对于任意个探寻点IV计算这个点和它的邻域点之间的一个关系元组,使用,定义为简单点特征直方图SPFH。对重新确定邻近点的邻域12 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究邻域的SPFH值来构造最终的FPFH直方图。FPFH与PFH类似,保留大部分的描-述能力,使算法复杂度为FPFH特征的数学模型如图24所示。<''許PJ':\kJPk2^7l_^;VPlO/”11?…-23-4【]图2FPFH数学模型32F-mamod[1ig24TheointcloudFPFHmatheticalelp对于探寻点IV在其邻域半径内,,首先计算其与邻近点间点对的SPFH特征2-4如图中红线所示。其次,对邻近点P&的所有SPFH值和探寻点匕的SPFH值-4中黑线所示的FPFH特征统计分析,如图2,从而得到最终。其中,重要的关联-信息被重复计算,如图24中加粗黑线表示。其中,FPFH计算公式可以通过公式2-2)计算出FPFH特征的计算结果示意图如图2-5(,所示。11k—=-FPFH/SPFH+xS-PFHP(^)(^)X(,)(22)^/=i….klI__32-]图25点云FPFH特征示意图[32-F25PoudFPFHeae[】igintclofturesktchmap46[](3)视点特征直方图描述子(ViewpointFeatureHistogram,VFH)13 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究VFH源于FPFH,具有FPFH的识别能力,并且在计算中加入视点变量信息,增加位姿判断。它的数学模型和计算。从而可应用于位姿估计和目标识别等问题-6-()示意图如图2(a)和26b所示。p、fv'QnV',pP',二^j^nw=u*v46图2-6乂?11示意图[】a)VFH特征示意图b)VFH计算示意图46F-^1i26VFHschematicdiaramgga)VFHcharacteristicschematicdiagramb)VFHcalculationschematicdiagram一VFH特征包含两个部分,计算中心点,,第部分扩展FPFH描述子。首先。设置视点,得到中心点的视点方向其次,计算中心点的视点方向与点云中所有点的FPFH。第二部分包含了视点方向信息,视点特征分量可描述点云的方向特性。-VFH特征的示意图如图27所示。VFH.*|展的FPFH元素|扩|视点元素I|\\.iH!1/^!!vhi\I!尸"\/\I?I/U丨/丨r.Bins462-7[】图VFH特征示意图46-H[]Fig27VFfeaturesketchmap47[](4)NARF特征14 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究NARF特征基于深度图像,需要从其中提取NARF特征点位置的NARF特征。描述子。主要由两个部分组成,关键点提取和关键点信息描述,首先,对点云提取边缘并进行打分通过结果来判断该点作为边缘点的概率。边缘和点仍的位置关系分为上、下。选择A和其右侧邻点求相、左和右四个方向对点间距离,并把该值和5比较,来判断边缘是否在该点右边。若距离大于表一示仍右侧点和外不在同个平面上。为了增加对噪声的适应能力,以A右侧若干(-点的中点来替代开右侧邻点23)所示。接下来依据此信息对该点进行,如公式-)打分,如公式(24。=—-FP(23)right^jx+i,ym=p<1fs)--—=max0-S1(24)r,ighl^rightVy一一其次个重要的参考依据,但并不唯。,在提取关键点时,边缘仅仅作为对于某个物体来说,除了边缘点之外,还应该有若干能表达其某些特征的点来做关键点,所以在设计关键点提取算法时,需要从边缘和曲面结构两个方面考虑,,并且关键点要能重复。同时关键点最好落在比较稳定的区域方便提取法线。一。某点的曲率与点云所在的曲面中,某点的曲率是最重要的结构描述因素之该点处曲面变化的剧烈程度成正比。在二维深度图像中,先找到开点及其周边与可以得到一:之距离小于25的点并进行PCA主成分分析,组值主方向V以及曲一=率值A。其中主方向V必然是个三维向量。对于边缘点而言,取权重wl,则V=,取权重w为边缘方向。对于其他点而言则方向为V在平面p上的--5)至(28)是。投影,平面垂直于仍与原点的连线将权重与方向带入公式(2。某点为特征点的可能性,其结果可得到若干特征点。最后进行极大值抑制-—s-四^ ̄0-)min1max1(25々(ja'|aVVJJn,^=- ̄nw(2'6)/l()?\^J'---ax/?/?lcosa(27)n(;)7)(^j)()J(I=-pI(28){)\{p)^2{p)因此,计算NARF关键点的NARF特征的步骤如下:a)计算关键点在所在patch的对齐方向15 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究b将星形图案叠加到atch中,每个光束对应于描述在该光束下的像素变化)p一c从描述子中提取个不变的方向)d)旋转描述子到这个值,直到旋转不变一一e法向量和唯的方向起定义了该点局部区域的一个6D0F坐标帧。)-按照该步骤得到的NARF特征如图28所示。纖顯b(c)()47图2-8NARF特征示意图【]4?F-RFur[1ig28TheNAfeatesketchmap目前的描述子虽然考虑了点云局部间的结构关联关系,但定义结构关联约束(中,多采用复杂的局部参考轴Localreferenceaxis,LRA)作为关联基准,计算一复杂且有进步简化的空间,。因此,以上述特征描述子为基础本文中提出了基于法向量分布的特征描述子,并从点云法向量、深度以及密度三个层面的结构关联关系出发,分别构建法向量分布统计特征直方图、深度统计特征直方图、密度统计特征直方图,兼顾以上三种特征直方图,构建统计特征直方图。2.2基于法向量分布的特征描述子点云法向量是描述点云局部区域几何表面形状的重要特征,能够准确的表现一出区域姿态,的法向量特征虽,是描述点云区域几何形态的重要数学特征而单然估算过程简单高效,却只是单点特征,难以表征点的邻域内几何信息,容易受到噪声的干扰。因此,,本文基于法向量,构建了法向量分布特征描述子并将其定义为NV(NormalVector)特征,NV特征能够有效加强邻域点法向量的空间结构关联关系,加强特征描述。16 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究根据邻域的概念,,计算点云中探寻点f在邻域的法向量。因此对邻域的定义尤为重要。由于点云局部区域的邻近点分布可以体现探寻点表面特征和非特征一区域,若探寻点与周围邻近点距离比较小,反之,,般为特征区域若邻近点分布稀疏,探寻点与邻近点距离比较大,通常为模型中平滑区域,此时,若仍考虑较远的邻近点,很容易引入噪声点或影响探寻点局部特征正确性的点。因此,本文将探寻点f的邻近点定义为,以£为中心,半径为搜索空间的之个与点欧氏距离最近的点,即为点云在探索半径内的邻近点个数£t,设f在;?=内点云个数为《,邻近点个数为A。若只考虑及半径内的的点,即,<=,灸否则选取半径内距离^点最近的t个点作为邻近点,即弋。1.潘???????:2-9V图N特征示意图F2-Thig9eNVfeaturesketchmap-那么,将NV特征定义为探寻点与近邻点法向量夹角的均值描述,如图29所示,其中,箭头表示其近邻点的法向量方向示意图。可以看出£点和,左面平缓一变化区域内点云法向量方向接近致,而右侧曲面区域内点云法向量方向变化较大。--设/V为点云中f点的法向量,计算公式如式(29)、(210)么,那,定义f点的法相均值特征为《v+-。,计算公式如式(211),= ̄ ̄C^P)(-ZLtkPi\PiP,29))JC-^A-enn123-PP,j{,,}(210)lJjU其中,/表示乃点的第y个特征向量,冗表示严点邻近点集合的质心,表示点的第个特征点,表示协方差矩阵中的第个特征值,为个邻近点中的y'_/之17 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究/第/个点的索引,^为f点在局部区域内的V个邻近点的法向量,为向量;C的丨卜|模长,arcco为JC的反余弦值。由于特征是对邻近点的点云法向夹角均值描述,只要空间f点与其邻近点之间的空间结构关系不变,空间的缩放和旋转,对£点与其邻近点之间的法向夹角特征是不会改变的,,因此特征具有尺度不变性和旋转不变性。v.特征来描述该点所在区域特性;使用《,若m特征值较大。则说明尸点所在,,;区域的几何结构特征变化较大,在其邻近附近的点曲率变化较大,存在突出的几何特征。相反,则说明尸点所在区域几何结构特征变化平缓,该区域内点云的特,.征接近配准过程中容易相互干扰,可特征,。因此根据每个点的《v表征该点所,在区域的变化程度。并将每个点的特征统称为点云的NV特征。2.3基于统计特征直方图的特征描述子基于法向量分布的特征描述子能够快速表示点云,但对于精度要求极其高的配一准或三维重构应用中,,其只能描述单个点的特征信息描述能力有限,若进步描述结构特征更复杂的空间区域,,则需要增加更多的统计特征。因此本文综合考虑点云中每个探寻点在各自邻域中,与邻近点之间的结构关联关系,从点云法向量分布统计特征、点云的深度统计特征以及点云的密度统计特征三个方面来定义完整的统计特征描述子,并将其定义为SF(Statisticalfeatures)特征。其中,点云法向量分布统计特征主要体现单个点的特征信息,而在探寻点邻域内的邻近点集的空间特性信息,,主要通过深度统计特征和密度统计特征来体现,具体来说深度统计特征主要体现探寻点邻域结构的凹凸性,而密度统计特征主要体现探寻点邻域结构的稠密度。在现有结构关联关系的定义中,基于局部参考轴(LocalReferenceAxis,LRA)一定义空间点云特征的方法被广泛使用,但它的精准度是建立在LRA唯并且稳定,若LRA不精确,的基础上,后续严苛的定义约束也不会得到好的描述性能并且LRA的定义复杂,基于LRA的描述子计算效率低下。因此,本文将SF特征设计为不涉及复杂的LRA描述,从而减少计算的复杂度。本文使用多种二维结构关联关系同时来约束点云间的三维关系,通过探寻点的NV统计特征、深度统计特征和密度统计特征来构建特征描述子。为T加快计算效率,本文使用直方图累加的;-方法将其组合,构建统计特征直方图。如公式(212)所示:s=thdtit-fcomcaihNV,e,hdens(212)t{p[ytj)^))){18 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究?六?/2其中,甙为第Z点的祆特征,表示z点的NV统计特征直方图,碑吻表()示/点的深度统计特征直方图,表示点的密度统计特征直方图,cowcaf〇表示直方图的连接函数。探寻点NV特征如2.2节定义,将探寻点领域内所有点的NV特征统计到直方图中,使用直方图存储探寻点与邻近点间的法向量关联关系,构成点云的NV统计特征。NV统计特征描述了曲面变化程度,能有效区分不同弯曲程度的曲面和平面。点云的深度统计特征定义为探寻点的邻近点与其在法向量投影方向的新点集-间的距离关联关系。如图210所示。LLthjpf图2-10点云深度特征F-i210Theolethfregintcouddeatupp=L,假设输入点云为P,其中"表示输入点云的总个数。对于探灼,,丹4寻点;?,,,r为半径的球,,将邻域定义为以;^点为球心可以表示为■==?12,,叹《表不/?的法向量。将作为p/,,表示点的法向量{\j//?探寻点的局部坐标系的正方向,过仍点的法向量,做与球表面垂直的切平面I作为局部投影平面,,将外在半径r内的所有邻近点投影到切平面Z上得到新的'"=……12点云集乂,?和间的欧式距离作为深度信息办7,,,将|_/4;/A;办办-办的计算公式如式(213)所示。将所有邻近节点的深度关联关系统计到直方图中,构成p点的深度统计特征。=--depthr2-gp/p(13)ji()-其中,w.8表不%与尸的点成关系。,^/A)(/19 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究/卜|||fv,?髮^z(?fx:A???'???‘:、????夕W**?*-K%??f::!:vy:〇yvl;^y))2-图1〗凹凸点特征示意图F-vig211Thecharacteristicschematicofconcaeandconvexointsp根据深度统计特征,能够很好的鉴别探寻点所在区域结构的凹凸性。对于不同区域的两个探寻点〗和,区域结构凹凸的局部深度的统计值存在强烈差异:丨点所_/-在区域为凸区域,投影距离为私,,其邻近点投影距离均大于等于&如图211中的+4^4+&;/点所在区域为凹区域,投影距离为其邻近点投影----,。,距离均小于等于 ̄,如图211中的欠心'&上\^1因此深度统计特征能够较好的定义该点所在曲面的区域结构凹凸性。点云的密度统计特征定义为探寻点^在投影面Z中的投影点的分布情况。如图2-r12所示,以半径为最大边缘,将等分成多组等间距的长度〇。那么,Z区域可以划分为以p的投影点为圆心,倍数为半径的同心圆。每个相邻区域内点云数-量可根据只和W的欧氏距离计算得到,将其定义为p。计算如式(214)所示:/= ̄--2'P/(14)PPPPiii^t\^{))^其中代表;?点的距离,%为仍点的法向量,/?的取值范围为〇r。密度只与//[,]||I特征可以区分点云分布区域的稠密度和分布不均匀度,从而区分不同的点云区域。20 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究a\-I!、乂.\、\jH'?\\:—?....--图212点云密度特征F-1inldifeaeig22potcoudensttury在给出反应探寻点局部关联关系的三种结构特征后,通过统计得到三个子直方图。子图包含Nu>^和>^3的分量个数可根据输入点云的几何结构信息,设置具体数值,确定每个特征对统计直方图的约束比重。2.4实验结果与分析、为验证本章提出的NV特征和SF特征的描述能力,分别对NV特征SF特征中的各分量特征进行实验验证。2.4.1基于法向量分布的特征描述子实验结果与分析48D[]利用斯坦福Bunny、ragon以及Happy数据集来验证NV特征的描述性能。对不同视角下的点云数据bunl和bun2、dragl和drag2、hapl和hap2分别提取特。征点,点云数目如表1所示表1点云数量比较Table1pointcloudnumbercomparison数量bunlbun2dragldrag2haplhap2原始数据400974025641841434677805675583NV关键点数据9484954197089731144041467821 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究一根据输入点云集合先验结构信息.设置个初始阈值与其比较,若,对《v,<q,表示所在区域起伏变化平缓,特征不明显,不作为特征点考虑;反之,若则表示该点与邻近点的法向量夹角均值较大,特征鲜明,将其作为提取的特征变化剧烈的点集。阈值A直接决定着特征点提取结果。综合考虑点云先验结构信息和NV特征提取后的点云数量对点云整体的描述能力,设置数据集Bunny-的&等于0.08,得到Bunny的NV特征点示意图如图213(b)所示:mfmmmm丨nu图2-13原始点云与NV特征点示意图a)bunl原始点云b)bunlNV特征点c)hal原始点云d)happylNV特征点ppy-racFiure213thesketchmaoftheoriinalointcloudandNVchateristicointsgpgppa)Bunloriinalpointcloudb)bunlNVcharacteristicpoint;gc)happyloriginalpointcloudd)happylNVcharacteristicpointNV.同理,将Happy和Dragon数据集的阈值&分别设为0.08和0.06,得到的2---4特征点如表1所示,可视化结果如图13(d)、214(b)、21(e)所示:22 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究圖國國rll画I國图2-14特征点选择结果49[]a)drag1原始点云b)dra1NV特征点c)dral特征点gg49)26)[](2]^(2(1(1原始点云1特征点0^8特征点%%Fure-214thechraceristiconectionresulsigatpitselta)drag1theoriginalpointcloudb)draglNVcharacteristicpoints49[Jcdralcharacteristicointsddra2oriinalointcloud)gp)ggp49NV[]edra2characteristicointsf)dra2characteristicoints)gpgp在三组实验数据集中,drag数据集的点云分布变化更加剧烈,更难提取点云49an[1a中的特征点。使用Wg等提出的特征提取算法分别对drgl、dmg2提取特征-4e)--点,结果如图21(、214(f)所示。与NV特征提取结果比较,如图214(b)、-4e21(),由于使用到了不同的点云数据集,需要重新考虑距离、曲率、夹角等条一件对判别特征点的约束,,。因此如何设置多个参数的复杂组合成为了个难题另外,多个影响因素的计算降低了特征点选择效率。相比而言,本文对NV特征点的提取方法简单,效率高,能够提取出点云模型中曲面剧烈变换的点云。2.4.2基于统计特征直方图的特征描述子实验结果与分析利用斯坦福Bunny数据集来验证统计特征描述子中三个分量的描述能力。实验对bunny的局部点云数据集中的代表点A、B和C展开说明,A、B和C点的位一一-置如图215所示。A点为凹曲面上的点;B为平缓曲面上的点C点为凸曲;面上的一点。23 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究’.‘f丨_图2-15局部点云图2-Fi15Thelocalointcloudgp一一一NV统计首先个特征。对比A点和B点的,分析统计特征中的第特征NV统计特征,NV特征反映夹角变化程度,这里将NV的范围设置为(0,tt/4),并将其划分为,若局部区域夹角大于20度,。如对于12个区间其空间变化显著-Bunn巴连接区域,A点的NV统计特征直方图如图216(a)所示y尾,大多数法向量间的夹角都分布在(68)区间中,法向量间角度变化较大,而对于平缓变化,的周围B点-(04),其NV统计特征如图216(b)所示,大多数点夹角都分布在,一,法向量朝向几乎致,,区间中。因此可根据NV统计特征较好的描述探寻点的点法向量特性。160130&140—-——^———110:1j90一-I一墨100醫70>I80I25060I4〇30I11201〇■壓謹I§||'10123456789101112123456789101112binbin-图216NV统计直方图a)A点b)B点.???-Finoni216NVstatistcahsoramaAotbBitgilitg)p)p一一分析统计特征中的第二个特征深度统计特征。由于NV统计特征只描述了点所在区域曲面变化程度,并不具备曲面结构凹凸结构的区分性,因此,比较A、24 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究C两点的深度统计特征,将深度特征范围设置为(02r3),并将其划分为13个区,/,间,对局部区域凹曲面上的点A其深度统计特,其空间变化显著。对于Bunny-征如图217(a)所示,可以看出,该点所在区域的邻域点多分布在直方图bin<7-1,平均深度值较小,7的区域中;而对于凸曲面中的C点其深度统计特征如图2(b)所示,该点所在区域的邻域点多分布在直方图中bin多7的区域中,平均深度值较大。220200漏_170150需圓120胃1〇00?:MI_|I2。…HD■■IIIi祖爆|羅ss123456789101112131245678910111213_3〇3binbin图2-直方图17深度统计a)A点b)C点F-rAig217ThedethstatisticalhistoamaointbCointpg)p)p一一最后,分析统计特征中的第三个特征密度统计直方图。对于点云密度分一布不均匀的点云,,凹凸性并不能唯的描述该点因此,考虑描述探寻点的密度特征,对比A点和B点的密度,设置密度范围为(0,r),并将r均匀的划分为8-个b,in计算得到的直方图分别如图218所示,其不同的直方图中邻域点的数量存在较大差异,并且分布可以用于区分点云的曲面分布特征。具体分析,由于密度特征反应的是点云分布区域的稠密度和分布不均匀度,所以考虑计算密度的期望==-值和方差值。对于图218(a),五denA24.6250,AienA6.4457,可以判断A点为稠-==密度高且分布不均匀度高的点;对于图218(b),五denB16.1250,DdenB4.9355,可以判断B点为稠密度和分布不均匀度均适中的点。???25 北京交通大学硕士学位论文基于结构关联的点云特征描述子研究4040303〇11BIlium.iiiim1234567812345678binbin-图218密度统计直方图a)A点b)B点F-1TheaaAiBig28densitstatisticlhistoramontbointyg)p)p综合每个描述子的特征,由NV统计特征、深度统计特征和密度统计特征构一建的SF特征描述子,能够较充分的描述探寻点所在区域的结构特性,从而唯的区分探寻点。2.5本章小结本章首先介绍了四种常用的特征描述子,即PFH特征、FPFH特征、VFH特征和NARF特征,这些描述子多采用LRA作为结构关联坐标系,其精准度建立在一LRA唯并且稳定的基础上。但LRA常存在定义复杂和计算不准确的问题,导致特征计算不准确,计算效率低下。鉴于此,本文基于现有描述子的构造方法,,首先提出了NV特征描述子该描述子使用法向量间的关联关系来描述点云中探寻点的点特征信息。其次,考虑单个点的特征信息对于精度要求极高的配准应用中描述能力的有限性,本文提出了统计特征描述子,该描述子在二维结构中使用局部区域内探寻点与邻近点的法向量间关联关系、深度间关联关系以及密度间的关联关系,约束点云空间三维结构关系一体的SF,形成将法向量统计特征、深度统计特征和密度统计特征整合为特征描述子,能够较充分的描述探寻点所在区域的结构特征,从而加强探寻点区分性。在本章实验中NV特征、SF特征中的各分量特征分别进行实验,初步证,对明了两个特征描述子的描述能力。26 北京交通大学硕士学位论文基于特征描述子的配准算法3基于特征描述子的配准算法一在上章中初步验证了NV特征与SF特征的描述能力。在配准算法中,基于特征寻求对应关系的配准方法最为充分的描述每个点云的区域特征,具有更好的可区分性一一。因此,本章将上章中的两种特征应用于点云配准中,进步说明NV特征与SF特征的描述能力。在本章中,提出了基于特征描述子的两步配准算法,两步分别为基于SF特征的点云初始配准算法和基于ICP的点云精确配准算法。在基于SF特征的点云初始配准算法中,提出了基于可变变长的体素网格下采样算法,能够近似等比例的化-简点云,进;提出了基于NV特征的关键点筛选算法,保留点云的关键点步降一低点云数据量,准确;提出基于SF特征的采样致性配准算法,去除错误匹配点求解变换矩阵,完成初始配准。在点云精确配准过程中,使用ICP算法,对初始配准结果完成精确配准。3.1相关工作目标全表面三维重构需要解决多视角点云的配准问题。现阶段点云配准研究已经取得了一定成果,目前配准过程的改进方法主要集中在如下的几个方面:1点集预处理()在扫描过程中,物体点的数量级非常高被扫描进视线,背景也会,被存到点云数据里。因此,对获得的点云数据,首先需要对大量的点云数据降级,减少点,提高配准效率。其次云数据量,降低运算负荷,需要对物体背景点、异常点和5G[]噪声点进行去除处理,,避免非物体上的点影响配准精度。在预处理过程中刘斌等将八叉树与ICP算法结合,提出OTOAICP算法,该算法使用八叉树求出重叠的51[]部分区域,根据部分点云使用最近点迭代算法完成配准。靳洁使用小波变化,,实现点云数据滤波重新定义划分空间。(2)点云的选择配准过程的前提是探寻点位于两块待配准点云的重叠区域,若位于非重叠区域则肯定无法得到准确的对应点,现有国内外对点云配准方法的改进多是为了提高点云对应点选择的效率和精度。52一[]针对效率,&1^1〇1^%&2等使用迭代最小误差法,从定程度上提高了ICP53一[]算法效率,但是并没有提高配准精度;Shini等提出了种三维匹配方法,使用j27 北京交通大学硕士学位论文基于特征描述子的配准算法奇异值分解方法,再不需要迭代的情况下,就能确定点云匹配的最佳旋转矩阵。54一[]Ya等利用CTNC(c-o-nexnglosettt)算法,提取局部点集,定程度上提高了ICP55算法效率[]-D-tree方,但是匹配精度还有待于提升。徐万鑫等使用分割最长轴K法,加快曲率特征配准,提髙ICP配准速度。梁群仙等使用改进的扫描线剖分算法,,完成点云快速重建根据扫描线之间点的相对位置关系。57-[]针对精度,J〇hnson等提出了旋转图像描述子(sinimage)该描述子是目p,一三维人脸识别的应用中效果显著前使用最为广泛的三维结构配准描述子之,在。58〇[]11〇等人提出了旋转投影统计局部描述子(RoPS),该描述子具有较强的描述能一59[]力能够保证点云的刚性变化,进步提升了三维点云配准的精度。3611丨11等将《:?算法与点云增强技术相结合,对原始点云数据增加新的点集来提,在配准过程中升配准性能,,减少测量错误完成了对工业现场中多个传感器获取点云数据的配6()准[]。郑德华等针对具有线性特征的扫描对象,提出了线性特征描述子,不需要61[]求取主方向,,使用最少数量的共轭线对数实现初始的点云配准。秦绪佳等提出了基于法向量直方图的配准算法,并根据三维不变矩特征,将描述子与ICP的%结合使用,能获得较好的配准结构。Men等针对彩色点云中特殊的颜色特性,将一P算法归化的点云数据与加权RGB色调值结合,加,提出了四维IC强点云描述tM特征。Richard等将点云模型表面细分为拟线性表面,轮廓细分为大量的拟线性,再用ICP算法对各部分分别进行处理曲线,采用分治的思想以达到更好的拼接效果。(3)最近点对查找一个环节最近点对查找是整个配准过程中耗时较长的,计算的对应点示意图果3-如图1所示。最简单的方式是使用欧氏距离作为最近点计算的度量,除了欧氏距离约束条件,常用多维特征来描述点云,对于多维特征,将曲面上探寻点及邻域一,的结构特征用个特定向量描述,邻域则采用A邻域查询和半径搜索两种方法通过匹配向量值来建立不同曲面中点的对应关系。图3-1对应点计算结果F-i31Theresultofthecorresondinointcalculationgpgp28 北京交通大学硕士学位论文基于特征描述子的配准算法一64[]八叉树(Octree)是种可以描述三维空间的树状结构。树中每个节点表一-。32示个立方体,每个节点如图。若立方体不空又可以被分为八个子立方体所示,右侧为对应的空间数据结构。,左侧为视觉上立方体的划分结果64-t图32八叉树数据结构]64F-2Thercu[]i3datastutreofoctreeg一-A::2A:维树(KDTree)同样是种可以描述A2维空间的树状结构,是基于坐()65[]。e,标轴的二叉分隔法,常用于多维点云数据的空间划分和检索与Octre相比-Tree划分点云数据,使用KD,三维空间划分相对均匀更适用于初始密度不等的点云。它的缺点为整个数据结构都是二分法,深层次的划分会缩小数据点数。因-,不是对所有类型的点云Tree划分的搜索速度都是最优。此,KD66[]在点云最近点查找中。Xin等提出动态分层配准,空间划分方法还有很多g算法,使用分层思想得到的配准,引入动态距离矩阵记录层次搜索距离完成配准67一000[]85mm。孙军误差达到.华等采用种分层块状全局搜索到邻近局部搜索的改进ICP算法,提高了ICP的速度及消除了点云缺失对配准的影响。目前,空间划分策略仍然是配准中的研究热点。4去除错误对应点对()对于计算得到的对应点对,其中含有大量不符合对应关系的错误匹配点,可以适用增加约束方法,去除错误对应点,优化对应点集合。常使用的约束方法包括'.。.点间距离约束、高斯曲率约束等.点间距离约束是从获取的初始点对中选取任意两个点对^,6和,如果(,)勺次)(正确匹配点对,根据刚体变换距离不变性可知成实际情,况下,所以这里取近似对应点对只,在两片离散点云中很难找到完美的对应点29 北京交通大学硕士学位论文基于特征描述子的配准算法需满足近似等于,,可选择适当的阈值4,计算每个初始点对从而得到符合距离约束的点对;高斯曲率约束是在微分中一点的高斯曲率值反应该点所在曲面的,使用曲面上变化程度,从而提取出符合条件的对应点集合。3.2算法流程本文提出了基于特征描述子的两步配准算法,两步分别为基于SF的点云初始配准算法和基于。ICP的精确配准算法在基于SF的点云初始配准过程中,首先使用基于可变边长的体素网格点云。精简算法,使,精简点云集合中的数据并保留重要的点云特征信息对精简结果用2.2节提出的NV描述子提取精简结果中的NV关键点。NV关键点集合能较一些细节信息好的表示整个点云的真实结构特征,但可能丢失。因此,使用2.3节。提出的SF特征计算方法,计算每个关键点的SF特征SF特征通过参数化关键点与邻近点的空间结构关联关系一,形成个多维的统计特征直方图。使用一RANSAC算法的基本思想计算SF特征对应点,进步计算出初始变换矩阵,完成基于SF特征的初始配准。在基于ICP的精确配准中,对初始配准结果,使用ICP算法快速迭代,根据一距离约束,计算出变换矩阵,进步提高最终配准精度。-算法流程框图如图33所示。30 北京交通大学硕士学位论文基于特征描述子的配准算法获取点云P获取点云g^—-—:“:'i丨rHf'云c精简点云p.,精简点.-;1!!^'提取NV关键点/V、提取NV关键点S■*_}...?:.:1:!\,,:!计算"的统计特征直W涵兩涵逞:||方图(SF)方图(SF)..;:L丨————"一、‘:'Ii1?-i厲I,sf初始—i_ir—?±-¥」Sii色—■ESF丨10的精■■dHzi^Ji!:d图3-3配准流程图F-ig33Thereistrationflowchartg3.3基于SF的点云初始配准算法基于SF的点云初始配准过程又可以分解为精简点云、提取NV关键点、计算关键点的统计特征直方图(SF)、计算变换矩阵等步骤。在保证原算法有配准结果的基础上,减少稠密点云数据对点云配准的不利影响。其中,每个步骤的计算都会改变初始配准精度和效率一,影响最终的配准结果。下面逐介绍每个步骤的意义和方法。3.3.1基于可变边长的体素网格点云精简算法,由于获取点云方式不同,导致点云数据存在差异常会出现点云分布不均匀。此外,影响点云质量。如果直接计算探寻点处的,环境中也会存在稀疏的异常点特征描述,引入的异常点会增加特征计算时间,特征描述性差,并且点云数量级31 北京交通大学硕士学位论文基于特征描述子的配准算法大。为此,必须对,对点云的后期处理产生不利影响,大幅度降低点云配准时间获取的点云数据在进行配准前进行过滤和优化处理。针对点云化简的问题,使,原有的化简算法多不考虑点云的初始分布得点云。初始稠密的部分化简后变得稀疏,从而使点云分布情况失真戀68图3-4可变边长的体素网格示意图[】68F-r['ig34Voxelschematicofvaiableedeg68[]本文基于体素网格滤波器实现点云下采样,提出可变边长的体素网格点云精简算法,,,在保证点云分布密度的基础上修剪点云使点云近似于等密度简化。三维体素网格将所有点云包围在子体素中-,如图34所示。本文将可变边长体素网格定义为以Z为边长的子立方体网格,在每个网格内,将点云数量记为《;若网格内点云数量《小于某个阈值a,则将其删除;否则,判断当前网格的点云数量^是否满足划分点阈值6,若数量n在阈值a与阈值6之间,一个顶点为源点,对当前网格以立方体的,以当前网格各边长的二倍为扩充长度,若不满足阈值6的数量要求,继续扩充,重新统计扩充网格内点的数量;否则计算本扩充网格内所有点云数据的重心坐标,选取与重心位置最接近的坐标来代替网格区域内部其他点云数据3-5。可变边长体素网格的精简流程图如图所示。32 北京交通大学硕士学位论文基于特征描述子的配准算法入点计算包围盒计算栅格边长LX计算栅格内点云数量《要区IYl^边长以2倍增长一1|计算栅格重心,作为代表点i—图3-5可变边长体素网格的精简流程图F-EVig35TheFlowchartSimlifieddeoxelReductionpg小立方网格的边长Z的确定。在基于可变变长的体素化网格下采样点云简化,立方体网格边长Z的设置十分重要,算法中。若Z过大子网格中划分的点云数,,子网格设置过小量多降低搜索速度;若Z过小,网络中存在大量空节点。由于小立方体网格的变长与点云的密度成反比。若点云密度小,表示空间点云数量少,L需要取大些来提高搜索效率,Z;若点云密度大,表示空间内点云数量多需要取小些来降低搜索时间。设小立方体网格边长为Z,单位小网格中点云数据个数为《,点的总数为iV。-可由公式(31)计算得出。L=a^[s[gNn=-/(31)yV=LLLxyz式中:《为权重,可以修改网格边长;s为比重系数;g为小网格中点云个数;设--。,、网格的体积为F那么,可以得到点云的体素网格边长Z如公式(32)(33)所示。33 北京交通大学硕士学位论文基于特征描述子的配准算法=JfWi(3-2)IVN=^-^+^4(maxmin)=-*L+;l-y(33)y(ma,yrain)_zz+A(raaxmin)式中:4为点云X轴上的最大值;为点云Y轴上的最大值;I为点云Z轴上&,的最大范围。l23/x/x/根据网格边长z,将点云数据划分为个子立方体,将点云数据映射到子立方体中,根据网格内数量与阈值关系,调整网格边长,实现高质量的点云精简。具体操作如下。一+网格编号为式-)对于任点P,,其所属的(34,x—^17Pirain_y ̄?yi=--^?(3-4)fpLL_2 ̄3Z7Pi_ima?,三元组为p点的网格编码表示,w表示x向下取整。设置阈值A其中,代表每个子立方体内最少点云数量;阈值^代表需要满足的最低上限6;统计第灸个子立方体中所有点云个数若6<-/?,计算点云集合的重心,如式(35),"V-<=-Y(35)ckncn??丨,咖、:个的X、Y、Z坐标;;fet、辦、其中、外分别表示第A网格中点云数据zcA:、Y、Z坐标。若《<a,将网格内的所有点lt分别代表第个网格重心点a的X去除。若a<?<6,表明当前区域内的点云较少,更新第A:个网格的边长,将其设一置为原始边长的2倍,即将其与周围三个网格合并,组成个新网格。重新统计3-5网格内点云数量,若大于阈值6,使用公式()计算重心,否则继续扩充,直,至满足精简条件,计算网格内点云重心坐标。遍历所有网格,动态更新网格边长34 北京交通大学硕士学位论文基于特征描述子的配准算法选取最接近重心的点云数据代替该子立方体网格内部其他点云数据,实现点云下采样。一基于可变边长的体素网格精简算法,可以定程度上去除点云集合中的数据并保证重要的点云特征信息,完成点云数据精简。优化原始点云数据。3.3.2基于NV的关键点提取算法对于使用可变边长的体素网格点云精简结果,若原始点云数据为稠密点云,仍会保留大量的点云,在计算特征时,计算每个点的特征会严重影响点云配准的整体效率并且稠密点云中包含大量特征相近的点云数据,;影响不同视角下的点云对应点匹配。因此,考虑对精简后的点云提取关键点。4769RF特征[][]现有的关键点特征提取方法有NA、Harris特征等,但多数方法7()68][][。只适用于处理深度图像对散乱点云的关键点检测,马麵溟等提出先计算曲,利用梯度法搜索曲面上的高斯曲率极值点,来提取点云的关键点面片模型;王49[]等将曲率和密度信息融合丽辉,实现对稀疏点云的准确关键点检测;这些方法,但模型存在较多参数描述性强,不同的参数值组合获得的关键点差别较大,若将其应用到配准中,不能保证配准算法的稳健性。?v的大小直接反映了该点所在真实物体表面局部区域起伏情况,因此可以利,用nv.值来提取特征鲜明的点。,一本文首先计算每个点的nv特征,其次,根据输入点云集合P和Q的先验结,构信息,对设置阈值若表示所在区域起伏变化平缓,不作为特征点考虑,,;反之若则表示该点与邻近点的法向量夹角均值较大将其作7374[][]为提取的关键点。^可根据配准误差变化动态调整。配准误差定义有多种,本文使用RMSE作为配准误差,即采用均方根误差作为配准变换后两片点云间的-对应点间的距离误差,可由公式(36)得到。''RMSE= ̄-x--Yx+y^+zz(3-6)jHIlllll=ii其中,;,《代表对应点的总对数和xz分别代表配准后两片点集中的对(,>,)应点的三维坐标。首先为^设置初始估计值,使用值选择出关键点集合和进行初始匹,计算匹配的均方根误差RMSE配s,并将其作为配准误差。若此时配准误差不满足约束条件公式-(37),减少^值,计算相应^值的关键点及RMSEa,直到得到最优关键点集合&和35 北京交通大学硕士学位论文基于特征描述子的配准算法-RMSERMSE<dRMSE<-V()_,37£,ElQmeX\\其中,(3_8)itk^%7娜a▽式中,为设定的配准误差收敛值,是允许的最大配准误差,为阈值&下点;RMSE和K中的对应点集合,,云关键点集合K,e*卜、^1p9;)为具体的对应点对阈值q下的配准误差。选择关键点的算法伪代码如下:-算法31(关键点选择算法)1:P=L=k,R,E输入a,,,p?u,>q,L,g,p2Q9xl2{}{j2计算点云法向量3构建-reeKDT4do初始化K,Kpg-for/Itonu计算点/与邻近点的法向夹角2-8计算《v?特征按式(),n■v>£ifi,更新关键点集合Kp,endifendfor基于关键点K,的初始配准^,得到对应点对-=-按式(37)计算RMSE&,qA0.01wh-<dileabsRMSERMSEC+002e+Qm(ii)5输出K=LK=kLk:,,,,Dk,U^},ak?a2j分别对待配准点云和参考点云选择关键点,并将获取的关键点称为NV关键点。3.3.3基于SF的初始配准算法对获取两视角下的点云数据P和Q,使用基于可变边长的体素网格算法精简点云,,分别为点云集I;对精简点云计算NV特征并基于NV特征选择关键点V36 北京交通大学硕士学位论文基于特征描述子的配准算法’一和Qa。NV关键点集合能较好的表征整个点云的真实结构特征,但可能丢失些细节信息。因此,使用基于统计特征直方图的特征描述子来增强点云的描述能力,统计特征描述子主要有四个关键的参数值,分别为三个子描述特征的直方图bin的尺寸N!、N2、N3和搜索半径r。其中价、N2和N3的大小决定了统计特征直方图的描述能力和计算效率,根据模型的先验知识和每个子描述子贡献特性,选取合适的,为每个分量做出不同的贡献。半径r决定了探寻点邻域范围冲、沁和处值,并影响计算效率。若r设置较大,会给探寻点邻域带来额外的样本;反之,r设置较小会使探寻点的邻域包含较少的点云,造成描述子可区分度低。由于输入三维模型的规模通常是未知的,因此很难直接求出r值的数量级,因此,本文采用模型大小的百分比来赋值。模型的大小被定义为模型包围盒对角线的长度。SFANV统计特征深度统计特征密度统计特征llBlMilyflllk,,biinlbnl+bin2bnl+bin2+bin3i图3-6SF直方图F-ig36SFhistogram使用统计特征直方图通过参数化查询点云与邻近点的空间差异一,形成个多维的统计特征直方图,从而较好的捕捉点云局部区域变化情况,弥补点云数量较少的缺陷-6,构造的特征直方图如图3所示。通过将NV统计特征、深度统计特征和密度统计特征相结合,加强点云的特征描述能力,保证配准精度。接下来结合采样一致性算法(RANSAC)的基本思想来完成点云的初始配准RANSAC算法的思想是对候选对应关系中大量采样,力求求出最优的采样结果。具体采样步骤如下:一(1)设置配准中对应点约束的最小误差为dmin,从其中个点云数据集乂中选取s个点作为样本点。一2对每个样本点,在另个数据集A中找到与样本点SF特征相似的点),保(存这些所有求出的点对集合一,从这些对应点中选择出些代表的对应点。37 北京交通大学硕士学位论文基于特征描述子的配准算法(3)使用(2)得到的对应点计算采样点的旋转矩阵R和平移矢量T,并计算点云的配准误差作为配准评价转换质量。重复以上三步直至得到最佳迭代结果或满足迭代总次数,使用-levenbergMarquardt算法做非线性局部优化。经过上述过程,得到最优旋转和平移矩阵,将点云4变换到点云A坐标系。3.4基于ICP的精确配准算法针对初始配准后的结果,,使用距离作为特征提高最终配准结果的正确性。对初始配准结果使用ICP算法进行二次精确配准。在不显著增加运算时间的同时,一进步提高配准精度3-7。配准流程如图所示。/初始配准点云数据A/初始配准点云数据b//士计算最近点^计算坐标变换和误差X初始配准点云数据A配准到初始配准点云数据B否差是否收^>结束V,/图3-7ICP算法精确配准流程图F-7Thwi3eICPalorithmforaccuratereistratonrgggiflochat38 北京交通大学硕士学位论文基于特征描述子的配准算法基于ICP的精确配准算法描述如下:-2算法3(基于ICP的精确配准算法) ̄ ̄1设得到初始配准点云数据集A和B2-TreetA使用KD寻找第次迭代中,集合在集合B中的欧氏距离最小集合1^Y=CABA=A---39,i(i)(o)()3计算坐标变换向量和误差。求取对应点的数据庞大,采用最小二乘法求解最优的旋转平移得到矩阵R和T=〇〇、Aa与B的平均误差也必)4点集坐标变换,将计算的旋转平移矩阵作用于初始配准的点云数据At,得到Ai+1一5定义误差收,若必.算法收敛敛条件成,跳转到6;否则跳转到26将最优旋转平移矩阵作用于初始配准的点云数据A7结束3.5头验结果及分析3.5.1散乱点云数据配准实验为验证本方法对散乱点云数据的有效性,选用斯坦福Bunny、Dragon以及Happy数据集来说明本方法的配准过程。将不同视角下的点云数据biml和bun2、dragl和drag2、hapl和hap2分别作为待配准点云和参考点云,可视化结果如图3-8所示,点云数目如表2所示。圓■遲图3-8点云可视化a)bunl和bun2b)dral和dra2c)hal和2ggphap-8Fig3PointcloudvisualizationaBunlandbun2bdraglanddra2cHalandha2))g)pp表2点云数量比较Table2ointcloudnumbercomarisonpp数据集合bunlbun2draldra2halha2ggpp原始点云个数400974025641841434677805675583下采样点云个数1111411607115231215492598664NY关键点个数48474133370847314404467839 北京交通大学硕士学位论文基于特征描述子的配准算法首先,使用可变边长的体素网格下采样后的数据集bunl和bun2来说明NV特征中对阈值的设置直接决定着提取关键点的个数,影响初始配准的&。阈值A时间和配准误差。设置不同的q值,得到q与初始配准时间和配准误差的变化曲-。,线如图39所示其中,实线代表初始配准过程的总时间变化曲线点划线代表初始配准结果误差变化曲线,,;可以看出随着巧值的减小剩余的NV关键点数量。逐渐增加,总的初始配准时间持而配,计算剩余点的特征计算时间变长续上升准误差的变化表现为下降趋势,当A值小于0.11时,配准误差曲线变化趋势趋于平稳。—?I91!60000-IeZZj\E:\\■U\--S0000\Z肩\/'?4SO00-\y■I\§*^4〇〇〇〇?-\4Q^i\Z--350CX)势\-*sco■*res.S\.5-25000^\'—?*'*?■**—?—-■.?.#?20000111^1?l1SQ〇D40〇81301201101009〇0.1.阈值图3-9阈值&与初始配准时间和配准误差变化曲线-nevewretrirrFiiialritrtiontimechacurithisatoneroi39Thethresholdandntesagggg-综合考虑NV特征对配准时间和配准精度的影响,选取图39中配准误差曲线的近似收敛处为设置的q值,即q等于0.11。同理,将Happy和Dragon数据集的。0,提取出所有视角下的点云关键点阈值q分别设为.08和0.06计算关键点的SF特征,首先,确定斯坦福Bunny数据集SF特征中N!、N2和N3的最优取值。其中,NV统计特征配准得到的均方根误差随bin数量变化的-RMSE来。统计特征直曲线如图310所示,使用均方根误差评价点云配准的结果方图中各分量可根据点云数据的不同动态选取。40 北京交通大学硕士学位论文基于特征描述子的配准算法6-.00E05误舅aso〇E-osA4-.00E05I\\::100-0.E5O.OOE+OO48011314617181920212223242526272829303112356791121151Bins-RM图310NV特征bin数量与SE变化曲线3-ThandRMchaFi10eNVcharacteristicbinuantitySEnecurvegqg,RMSE值,着bin个数增加n总个数,随其中,纵轴为可以看出,横轴为bi=RMSE值逐渐较小2时RMSE。因此,设置N1的值为12;,当bin1,值趋于平稳-点云的深度统计特征配准得到的均方根误差随bin数量变化的曲线如图311=所示RMSE值随bin数目的变大逐渐变小,当bin8时,配准误差趋;可以看出,。于平稳。因此,将N2设置为87-^.00E05差-6.00E05I5-1.00E05||4-.00E05\h||^3.00E^)5V\//■E-05\-1.00E050+.00E0071819202123227230311234567891011121314151612242526829Bins3-RMSE变化曲线图11点云的深度特征bin数量与Fig3-llTheointclouddepthcharacteristicbinquantityandRMSEchangecurvep41 北京交通大学硕士学位论文基于特征描述子的配准算法同样,点云的密度统计特征随配准得到的均方根误差随bin数量变化曲线-312所示n值,即N3等3。如图;选取RMSE值最小的bin数目作为最优bi于14.00M4误差3-.50E04A3-00EO4I20E-04.52-00H04+O.OOEOO34567901112131415161181920222324267893011281721522223Bins-b图312点云的密度特征in数量与RMSE变化曲线F-nRM1ThitcloudiracteristicbiuantitandSEchanecurveig32eondenstchaqygpy一8,分别确定N,N13,1等于12,N2等于3等于将三种特征组合在起构建一一3-出个33维的统计特征直方图,点云其中点的试特征如图13所示。20—.—’个T:数18NV特征深度密度丨丨|13579111313315171921232527293Bins图3-直方图13统计特征F-careramig313statistilfeatuhistog其中,横坐标bin取值的1到12区域代表NV统计特征分布,12到20表示点。云深度统计特征分布,20到33表示点密度统计特征分布从而,根据点云法向量42 北京交通大学硕士学位论文基于特征描述子的配准算法分布特统计征与点云深度统计特征以及邻域点的密度统计特征三个方面,得到完整的统计特征描述子。一使用统计特征直方图参数化查询点云与邻近点的空间差异,形成个多维的统计特征直方图,能较好的捕捉点云局部区域变化情况,弥补提取关键点后点云数量减少的缺陷。通过将NV统计特征与点云密度统计特征和点云深度统计特征相结合,加强点云的特征描述能力,保证配准精度。同理,得到HaDraon的最优b3ppy和ginl、bin2、bin的值分别为10、6、14和10、10、12。将本文所提算法应用于点云Bunny、Dragon和Happy数据集,初始配准结果-5_5-?xx分别如图314(a)(c)所示.20518l〇,初始配准误差分别为3.27605l〇、2_52-?和.06316xl〇。使用ICP算法完成精确配准后,配准结果如图314(d)(f)-5x_5-6X所示,最终配准误差分别为18438x.0。.4l〇、l.〇666l〇和9728951國圓■i】■棚【瞿3-图14点云配准结果a)Bunny初始配准b)Dragon初始配准c)Happy初始配准d)Bunny精确配准e)Dragon精确配准f)Happy精确配准F-igure314TheointcloudregistrationresultspaBunnyinitialreistrationbDraoninitialreistrationcHainitialreistrationgg)ppyg;))gd)Bunnypreciseregistratione)Dragonpreciseregistrationf)Happypreciseregistration.一为了进步验证本文方法的有效性,对本文提出的基于SF的初始配准与直接使用PFH、FPFH和SHOT的初始配准对比。如表3所示,相比于直接使用PFH、43 北京交通大学硕士学位论文基于特征描述子的配准算法FPFH和SHOT和RANSAC算法,本文优化了初始点云数据,克服了特征相近的点对求取对应点的干扰,并使用多种特征来增强关键点的描述能力,在迭代次数相同的情况下,,基于SF的初始配准能够保证配准精度提高配准效率。表3配准结果比较Table3Comarisonofreistrationresultspg ̄数据集Eiml时间/mS初始配准误x_5点云数量点云数量次数差/i〇mNV-S-FRANSACoursbunl、bun2(7296())4347413310033.27605FPFH-RANSACbunl、bun2(>111141160710011831010.0778PFH-RANSACbunl、bun2(>11114116071004681432.76851SHOT-RANSACbunl、bun2()1111411607100.5711611472482NV-SF-RANSACoursl、d2()drara(gg>37034731100120192220518.FPFH-RANSACldra、dra2(gg)52312154100519484.38611214PFH-RANSACdral.dra2(gg)1152312154100.8935719196571SHOT-RANSACdragl、dra2()523121541004831.79357g111201NV-SF-RANSACoursl、ha2()(happ>440446781001389302.06316FPFH-RANSAChal、ha2(pp)925986641003763892.33662PFH-RANSAChapl、ha2()925986641001.8368p14302781SHOT-RANSAClha、ha2(pp>9259866410034251591.53188aon和Ha点云数据集,对下釆样Bunny,直接使用ICP算法配准得、Drgppy--55到的配准误差分别为1.62231xl〇,7.72603x10'1.28465xl〇。可以看出,本文算法配准更加精确。综上,本文算法与基于其他特征描述子的算法相比,配准效率显著提升;与直接ICP算法配准相比,配准结果更加精确。3.5.2自获取深度点云数据配准实验7576[][],kineCt20,为验证本算法的普适性本节使用.获取两片真实点云深度图像一-a步验证本算法的有效性,进。图315()为获取的两片真实稠密点云数据数量分别为679779、512678。44 北京交通大学硕士学位论文基于特征描述子的配准算法國圓國3-15图kinect2.0获取点云数据a)多视角原始点云b)视角1下NV关键点c)视角2下NV关键点F-5Theig31ointclouddataobtainedusinkinect2.0pgamultivieworiinalointcloudbNVkeointofviewonecNVkeointofviewtwo)gp)yp)yp=,012NV关键点的数目为678、经下采样后选取NV特征的阈值A.,得到93-ARF和Harr65486。如图315(b)和(c)所示,而使用Nis方法检测点云的关键-点,结果分别如图316(a)和(b)所示。由于NARF和Harris检测的关键点数量较少,少量的关键点并不能表征点云空间的整体结构特性。因此,很难使用该类方法,准确计算出两视图下的点云对应关系,从而保证配准算法的正确性。Wfgl-图316关键点检测结果a)NARF关键点b)Harris关键点F-ig316ThekeointdetectionresultsaNARFkeointbHarriskeointyp)yp)yp首先,经滤波去除背景和噪声,使用基于可变边长的体素网格滤波器简化点云,得到下采样后的点云数目分别为119756、93214。使用ICP算法配准,由于两视角下的点云数据初始位置下的重叠区域相互影响,造成错误的配准结果。由于初始位置相差较大,并且两片点云数据中有重叠X,ICP算法,2区域配准误差陷入局部极小值,导致配准错误配准误差为.2290710A-配准结果如图317(a)所示。45 北京交通大学硕士学位论文基于特征描述子的配准算法3-)图17配准结果a直接使用ICPb)本章算法F-ranIrecig317ThereisttioresultsauseCPditlyb)thisarticleg)使用本文提出的两步配准算法完成配准-,配准结果如图317(b)所示。最终><弋对比图-的配准误差为1.117561〇318(a)、(b)可以看出,在点云初始位置相差较大并且点云数据中有重叠区域时,直接使用ICP算法会导致配准错误。-同理,对kinect2.0捕获的chair点云做配准实验,图318为两个视角下的chair数据可视化图片。3-图18两视角自获取chair点云-Fi318Twoogintsofviewfromtheacuisitionofchairointcloudspqp配准参数和配准结果数值如表4所示,由表可以看出,使用ICP算法直接作x为约束,在初始位置相差较大的点云中,得到的配准误差为37680.45411'由于初始位置相差较大,直接使用ICP配准,距离约束过于简单,配准结果错误,如-图219(a)所示。使用本章提出的两步配准算法,多考虑了点云的深度特征以及邻域点的密度对点云描述,配准误差为1518x10.10'配准结果如图319(b)所示。46 北京交通大学硕士学位论文基于特征描述子的配准算法rz!c^3图3-I)19配准结果a)直接使用CPb本章算法-ThFig319ereistrationresultsofusindistanceauseICPdirectlbthisarticlegg)y)表4点云配准结果比较Table4Comparisonofregistrationresults ̄ifii迭代次数时间/ms配准误差x_5点云数量点云数量/lOmICP8205389766631008831034.547681(错误)Ours8205389766631001801921.10518综上,,本章提出的配准方法综合考虑了点云的空间结构关联特征,能有效的区分探寻点,避免配准错误,保证配准精度。有力的验证了本算法的普适性。3.6本章小结本章提出了基于统计特征直方图(SF)的点云初始配准和基于ICP的精确配准的两步配准算法。在基于SF的点云初始配准过程中,首先提出了基于可变边长的体素网格点云。,精简算法,能精简点云集合中的数据并保留重要的点云特征信息对精简结果使用2,.2节提出的NV描述子提取精简结果中的关键点减少特征相近点间的相互干扰,从而降低不同视角下的点云对应点匹配的错误率。NV关键点集合能较好的一表征整个点云的真实结构特征.3,但可能丢失些细节信息。因此,使用2节提出的SF特征计算方法,计算每个关键点的SF特征。SF特征通过参数化关键点与邻近点的空间结构关联关系一,形成个多维的统计特征直方图,并根据模型先验知识和每个描述子的贡献能力,为SF特征选取合适的bin值。SF特征能较好的捕捉关键点局部区域变化情况,弥补点云数量较少的缺陷。使用RANSAC算法的基本思想计算一SF特征对应点,进步计算出初始变换矩阵,完成基于SF特征的初始47 北京交通大学硕士学位论文基于特征描述子的配准算法配准。在基于ICP的精确配准中,对初始配准结果,使用ICP算法快速迭代,进一步提高最终配准精度。本章实验分别对散乱点云数据和自获取深度点云数据完成配准,实验结果表明,,与基于其他特征描述子的配准算法相比本章配准算法效率高;与直接使用ICP的配准算法相比,本章配准算法避免了配准错误,配准结果更加精确。48 北京交通大学硕士学位论文三维重构演示系统4三维重构演示系统三维重构技术在逆向工程、虚拟现实、工业测量等领域应用广泛,在实现第三章的点云配准方法后,可得到完整物体的点云数据,基于此数据本章使用三维重构技术实现物体的可视化,进行后续应用的研宄。。由点云数据重构出三维模型的方法可以划分为两类,插值法和曲面逼近法插值法是在原始点云的基础上插入新的数据得到重构物体;曲面逼近法则是使用曲面面片或复杂形状来逼近原始点云,形成三维形体。目前主要的曲面重建方法有77N56[][]。参数曲面重建、隐函数曲面重建^nDelaunay曲面重建三类方法其中,基eauna,。于Dly方法灵活,对曲面边界更为敏感因此被广泛使用4.1基于Delaunay的三维重建是一在三维点云重构处理中,Delaunay种三角剖分算法,能够将散乱的点云数据集划分为较为均匀的三角形网格,使点的几何变为面的结合,三角剖分后的结-。果如图41所示、4-1Vorono图i图F4-1Voronoidiaramigg'56[]本节简要介绍£)11111116等提出的Delauna三角化算法。其基本思想是对(^说y_£>D/,/于空间中的点云集合S,当前第个点的Delaunay四面体为(〇向()四面体“”中加入第汗1个点,控球准则,新加入点对原有四面体分裂并使分裂结果满足79[?]£/+1直至点云中所有点加入其中,得到(),重复这个过程。49 北京交通大学硕士学位论文三维重构演示系统算法主要使用三种数据结构,、四面体数据结构和三角片数据结构点数据结构。点数据结构是对所有数据点X、Y与Z轴坐标的存储;在四面体数据结构中,定义四个指向顶点的指针和四个指向邻接四面体的指针,设四面体中的四个顶点分别为v。,^',四个三角面片分别为■?VVVSVVV'yVVV'y0(0,,,,,,,,,W2。l2)i(03i)0(230)i()在三角片数据结构中,定义指向所在四面体的指针以及为所属四面体面的编号。Delaunay三角化算法描述如下:4-1算法(基于Delaunay三角化重构算法)一1、读入点云数据中的所有点,选取数据中4个不共面的点初步构成个四面体2一,。、读取数据中下个点判断该点与当前四面体位置关系若在四面体外部,3跳转至;若在四面体内部,跳转4;若落在三角面片上,跳转至5;若该点属于某个边上,跳转至6;若与当前点重合,则忽略该点,重复执行2。3、。计算距离该点最近的四面体上的三角面片,跳转至74,,,7。、以该点为顶点将该点所在的四面体分裂形成新的4个四面体跳转至5、以该点为顶点,在所属三角片所属的四面体分裂,形成3个新的四面体;跳转至7。6,7。、以该点为顶点,与当前四面体非共线的顶点连接,分裂四面体跳转至“”7合。,,进8、检测不符控球准则的四面体若存在行调整;否则跳转至8,、若点云数据集中无剩余点结束。否则跳转至2。elauna三角化算法,即可得到三将第三章得到的点云数据配准结果,使用Dy。维重构形体,完成三维重构任务接下来。将本文所述第二、第三及本节所述理论内容,用软件实现其具体功能,搭建三维重构系统软件,验证本文所述算法的实用性。4.2三维重构系统软件总体设计4.2.1三维重构系统软件搭建环境本系统搭建环境为:Windows7操作系统、Intel(R)XeonCPU处理器、12GB内存和64位操作系统,基于VisualC++6.0软件开发平台,结合C++语§、PCL72[8G][]界面开发工具,点云库、QT及FLANN等第三方依赖库实现本文第三章、第。四章4.1节所述算法。下面介绍上述几种关键库50 北京交通大学硕士学位论文三维重构演示系统PCL是一个独立一家商、大规模、开放式的二维或三维点云处理平台。其由业公司提供资金支持的非营利组织,在BSD许可证的条款下发布的,为商业和研宄提供免费使用。其支持多种操作系统,实现对点云获取、特征提取、过滤、分害I、检索、识别、重构等多种点云处理算法,被广泛使用。J提供了图形用户界面程序开发框架,支持GUI程序开发,包括大量宏QT,易—。a++扩展,并且允许真正的组件编程常与VisulC起使用,实现可视化界面开发。一Qhull是个强大的程序,常用于计算二维、三维,以及四维以上维度点集的凸包、Delaunay三角网等,常用与应用科学、工程、统计和数学中。Matlab和Octave都使用它来提供计算几何,Mathematica使用它实现Delaunay构造网格。FLANN是一个快速计算近似最近邻搜索库。包含最适合于最近邻搜索的算法,以及一个根据数据集自动选择最佳算法和最佳参数的系统++。FLANN使用C写成。能够很容易地通过C,MTALAB和Python等绑定提供的库,用在很多环境中。Eigen是C++中的开源的、高层类库,是关于矩阵类的模板函数,有效支持矩阵、矢量、线性代数、多项式计算器,快速傅立叶变换、数值分析及其相关的算法,适。支持各种矩阵的标准数据操作,并支持扩展数据类型用范围广。4.2.2三维重构系统软件设计思路本软件命名为PCReconstration(PointCloudReconstration),界面可视化框架如图4-2所示。—「一=====—..====二=======二:==P|j]〒||||I|一二视角原始点云!第视角原始点云丨第i丨|可视化丨可视化丨!|I|!二==3\===:今控制台1最终重建结果可视化输入■iIi[I!||[|一!第视角关键点丨丨第二视角关键点可视化可视化丨丨IIIIII-.:1—1vI1—.._...一--1—函4:2涵碗F-i42Sstemdesginyg一二视角原始点云软件主要分为控制区域和显示区域。显示区域又分为第、一可视化模块、第、二视角关键点可视化模块、最终重构结果可视化模块。而对51 北京交通大学硕士学位论文三维重构演示系统一于第视角下的点云可视化使用红色,对第二视角下的点云可视化使用绿色显示,从而区分两视角的点云;输入的点云要求为PCD格式,点云显示中可以实现在3D空间中的点云平移、旋转和缩放等交互操作。控制台中的动作能控制相应可视化模块,实现求关键点、重构的功能。一系统中输入的点云数据格式均为PCD文件格式,每个PCD文件由文件头和数据组成。其中,文件头用于声明文件中存储点云数据的格式和特征。PCD文件头包含常见的字段解释如下:VERSION字段指定版本IZE;S字段指定每个维度的单位字节数;FIELDS字段指定字段名称和点云维度;TYPE字段指定类型;POINTS字段指定点云的总数;COUNT字段指定每个维度中元素个数,默认所有维度的数目为hDATA字段指定点云数据类型,支持ASCII和二进制。WIDTH字段指定点云数据集的宽度;HEIGHT字段指定点云数据集的高度。PCL库支持对PCD点云格式处理,,可以识别PCD格式定义通过与其他依赖库合作,实现点云的处理和显示。4.3三维重构系统软件操作流程使用斯坦福Bunny数据集和Dragon数据集对三维重构系统做简单实验。首先,一系统开发界面如图4-3所示。每个黑色区域为vtkWet个Qidg,己实现对点云的旋转、平移和缩放功能。“”“”—openviewone按钮用来选择第视角下的pcd文件,openviewtwocd-按钮用来选择第二个视角下的p文件,得到的初始可视化结果如界面43所示;'—-:*?ntrl:coo…HHHHHH—keonejvorowareoornttwoogtir>ykevtotcob>o*wvnonBVB:j4-3aBunnb图原始点云加载界面)y)DragonF4-TnBDri3heoriinalointcloudloadiniterfaceaunnbaonggpg)y)g“”“”选择出两片不同视角的点云后,可以通过keyone按钮和keytwo按钮,分别调用本文提出的基于法向量分布的特征描述子计算方法计算每个点的NV特52 北京交通大学硕士学位论文三维重构演示系统“”“”征,并使用该特征提取关键点,分别在keyone和keytwo窗口实现可视化,该可视化界面的显示结果如图4-4所示; ̄ ̄ ̄"一"- ̄**二二一nlnt!cd:tr0^;coraDncooImmIvi?Uyteivn??-^?w?MV,,,:■國n〇??rcotwneyoetwocofM&tfiok?〇tevt?ersbrtkwekyw?|图4-4关键点加载界面aBunnbDraon)y)gF-ni44keointresultinterfaceaBunnbDragyp)y)go“”reistraton,通过gi按钮,系统会对求得的两视角下的关键点点云集合使用本文提出的基于SF的点云初始配准算法,得到初始配准变换矩阵R和T,对原始点云变换后,得到新的初始点云,再使用基于ICP的精确配准算法进行精确配准,最终得到的配准结果在系统中显示在最右侧-5,如图4所示。“*r—R..>0!- ̄一!'X.?-■-?^tv>K办?'-?-r?'''?...,.W、■'、->-''?'■〇'、-y'y叫t(S^{.nrltlni:cotoni:conropoepae—HIE3mmBtm=H_國|^|-〇kO/rcwwtrotokewri?torf?cfwexjyoneketo*camtrt?yy-nnbD图45三维配准加载界面aBuragon)y)-DFiure45ReistrationresultinterfaceaBunnbraongg)y)g“”三维重构算法实现点击reconstruction按钮,系统使用本章介绍的Delaunay对配准得到的完整Bunny点云和Dragon点云的三维重构,并进行光顺处理,最终-得到的重构结果如图46所示。53 北京交通大学硕士学位论文三维重构演示系统 ̄ ̄? ̄ ̄-*'mm:,...r|ro^rrlanelotv^cstopc——iH^RMHorvI〇n?????V^w*^*tw>7Tl—?〇lJ|;|^_|11^1::_國0j,rw?ssr<nk?y?〇ei?ytwofew她《iankeyok?ytro?cwsta〇ij|’;图4-6三维重构加载界面aBunnybDraon))gF-BDiure46Reconstructionresultinterfaceaunnbrang)y)go界面每个QvtkWidget窗口中,均使点云显示在窗口中间,并实现了绕中心旋转和平移等操作。本系统可以对初始位置任意的两片不同点云实现高精度配准,并将配准结果和重构结构可视化显示。4.4本章小结本章首先介绍了基于Delaunay三角化的三维重构算法。其次介绍了演示系统的搭建环境,给,并介绍了其中几款重要的库文件。最后出本文三维重构演示系统的总体设计。本章所述三维重构演示系统实现了本文第三章提出的基于统计特征的两步配准算法,并使用本章所述的Delaunay三角化和光顺处理实现最终的三维重构,系统界面分别对配准前点云原始三维结构图、点云的NV关键点三维结构图、点云配准结果的三维结构图以及基于Delaunay重构的重构结果结构图实现了可视化。本章使用斯坦福B_y和Dragon数据集展示了演示系统实现点云重构的整体流程。本章所述演示系统实现了本文第三章提出的点云两步配准算法,可以完成对任意位置的特征点数据集的高精度配准和三维重构,并以界面形式可视化。54 北京交通大学硕士学位论文工作总结与展望5工作总结与展望一,本章对全文工作进行总结,论述本论文的主要研宄成果并提出了下步工作及展望。5.1工作总结,本文以不同视角的点云数据作为原始数据,结合点云空间结构关联关系构建基于法向量分布的特征描述子,(NV)和基于统计特征直方图的特征描述子(SF)并将两种描述子应用到点云配准过程中,提出了基于SF的点云初始配准和基于配准的两步配准算法,ICP的精确,完成多视角点云的快速、高精度配准。最后使用软件实现了三维重构演示系统。本文所做的工作如下:1、深入的研宄了点云配准方法,并对点云精简算法、点云特征描述方法、点云关键点提取算法,确定、点云匹配算法进行了深入分析后,整理自己的思路了自己的研究方向。一2些经典点云特征描述子进行总结分析,针对现有特征描述子使用、对现有LRA定义探寻点与邻近点间关联关系存在的问题,本文使用直方图来存储点云间,的结构关联关系,用多种二维结构关联关系来约束点云间的三维关系即将探寻点的NV统计特征,形成完整的基于统计、深度统计特征和密度统计特征相结合SF),证明了本文所提出的描述子具有较特征值直方图的特征描述子(。通过实验好的描述能力。、,3本文提出的两步配准法,在点云初始配准过程中的点云预处理阶段提出了基于可变边长的体素网格精简算法,在传统固定网格边长拆分方法的基础上,,对基本单元选择性扩充根据每个单元点云的数量,使用扩充后区域重心点代替区域内的所有点,。对精简后的点云使用NV特征提取关键点,避免特征相似点一云之间的干扰,加强关键点的描述性能,并使用采样致。计算关键点的SF特征,性配准算法去除错误匹配点,准确求解变换矩阵完成初始配准阶段。在点云精?一通过实,确配准过程中,使用ICP算法求取最佳变换进步提高最终配准精度。验验证算法的有效性。4、本文基于VisualC++6.0和QT5.6.3开发平台,PCL和FLANN、BOOST、VTK等第三方依赖库,实现了本文第三章提出的点云配准算法,并使用Delaunay算法实现对配准结果的三维重构,并以界面形式可视化。55 北京交通大学硕士学位论文工作总结与展望5一.2进步的研究建议本文对点云配准方法进行了深入的研宄,并将其应用在散乱点云和自获取点一些优化和改进的方面云的配准中,取得了不错的效果。但是还有存在。1、在继续研究基于结构关联的点云特征描述子的基础上,本文希望加入点云的形状先验,增加全局特征描述,加强点云的描述能力。通过形状先验对点云方位的描述,快速检测点云间的相对位置关系。2、在继续研宄基于结构关联特征描述子点云配准算法的基础上,可以在未来工作中,有针对性的研究更复杂的点云数据结构,如稀疏点云、带有大量噪声点,以支持更复杂结构点云的高难度配准,云的配准算法。并且考虑对计算的对应点一增加约束条件,从而进步提升配准精度。3、在继续研究点云配准算法的基础上,可以在未来工作中根据多个视角下的两两点云配准误差,,更正多视角配准后的完整点云配准结果提高多个视角下得到的完整点云配准精度。本论文错误或者不足之处在所难免,希望各位老师们和同学们指正批评。56 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果参考文献1v-udSehalACemeaDMakaeevaM.RealTimescaleinvariant3draneointclo,,[]ggpreii2012111220-229straton.6:.^,g][2]PomerleauF,ColasF,SiegwartR.AReviewofPointCloudRegistrationAlgorithmsforcsoundationsTrensinobotcs20-:MobileRoboti.F&dRi15411104.,,()[J][3]邱春丽.基于点云的曲面重建技术研究[D].北京交通大学,2014.“”-4G.GenStructuredliht3DsuaceimainatutorialAdvancesinOt.andPhoonicsrf:t,,[]ggggp-321281602011.,()()M''5S.GorthiandP.RastogiFrineroectiontechniues:Whitherweare?Ot.Laser.En.,gpqpg[]j481-142010330.,()D6..华中科技大[]吴庆华基于线结构光扫描的三维表面缺陷在线检测的理论与应用研究[]2013学.,-7oredossicKovakBtaliiiPPPovNe.Threedimensonalmeasurementsofbodesnmotion[],,,basedonmutile-laser-lanetrianulationevistaecnicaDeLaFacultadDel.RTppg[J]-61InenieriaUniversidadDelZulia2015382:53.g,,()8.何原荣.,郑渊茂,潘火平,等基于点云数据的复杂建筑体真三维建模与应用J遥感技术[][]-与应用2016316:10911099.,,()一LED和9张朝勇苏真伟乔丽等..种基于线激光的钢轨表面缺陷检测系统科学技术[],,,[J]2012123677-与工程:989880.,,()[10]王文标,马孜,李爱国,等.再制造工件缺损部位定位技术研究m.科技通报,2010,265-:652656.()11KBlaskSG.Reistrationof3DointclouddatabcreationoffiltereddensiMineart,[]gpyyimaes:USEP2272045B1P.2011.,g[]^-i-Diniii12S.ZhanRecentroressesonrealtme3shaemeasurementusdtalfrne[]g,pgpggg^-roectiontechniues/Ot.LaserEn.4821491582010.p,()()jqpg[13]彭仲涛.结构光测量点云数据配准研究[D].哈尔滨工程大学,2015.MeKaiono-DshaesIEEransactio14PaulBeslNeilD.Amethodforrestratif3.ETnson,[J][]JygpPiiIntience142-attern:AnalssandMachneell992123956.g,,y()15ZhonYZhanM.AutomaticReistrationTechnoloofPointCloudBasedonImroved,ggyp[]ggICPAlorithm.ControlEnineerinofChina2014.,g[J]ggusin--16ChenWuXWanMYetal.3Dshaemodelinaselfdeveloedhandheld3D,,,[]Jgpggp-ntclaserscannerandanefficientHTICPoiloudreistrationalorithm.Otics&Laserpgg!}]pTechno-lo201345l:414423.gy,,()17YLiH-IloballtimalluontoDICF-anCambellDetalPointtJ.GoCP:AGOSoti3Se,,,[]gpypReistration.[J].IEEETransactionsonPatternAnalysis&MachineIntelligence,g204-15:11.()18FischlerMABollesRC.Randomsamleconsensus:aaradimformodelfittinwith,[]ppggalicationstoimaeanalsisandautomatedcartorahM.ACM1981.gpy,ppgy[]57 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果9enHunenBRAN-BasedDARCES1ChCSYPCh:roachtost.SACANewAFa[],g,gJppAutomaticRegistrationofPartiallyOverlappingRangeImages[J].IEEETransactionsonPatternAnalsis&MachineIntellience20022111-:12291234.,,yg()[20]PapazovC,BurschkaD.AnEfficientRANSACfor3DObjectRecognitioninNoisyandccludedScenesomterison-ccv20106492-O:..CuViA,135148.[J]p,21Men-LIianPanRANSAC-basedstabDiSHPlelanefittinmethodofointW.A[]g,Jg,gpgpclouds[J].ScienceofSurveying&Mapping,2015.[22]MagnussonM,LilienthalA,DuckettT.ScanregistrationforautonomousminingvehiclesusinD-ournaoFieRotics201024-3.NDT.lfldbo10:803827g,,[J]J()23修祥RFD-NDT多1.3.205]胡,张良结合NA特征的改进型视点云配准[J]信号处理,,[3112-:16741679_()243D—NDT骆林.,杜宁显云等.融合点云RGB影像和算法溶洞的点云自动精确配准J[],张,[]2015-矿山测量4:6771.,()25Kwan-oerekehtiAmethodforautomatedreistrationHBaeDDranizedoint.Li.ofuno[]g,ggpcloudsournalof-PhotorammetrandRemoteSensin2008631:3654..,[J]Jgy()g26朱延娟周来水张丽艳.散乱点云数据配准算法.计算机辅助设计与图形学学报2006[],,,,184-:75481()4.[27]徐金亭,孙玉文,刘伟军复杂曲面加工检测中的精确定位方法[J].机械工程学报,2007,436-:175179.()[28]梅元刚,何玉庆.低维特征空间中基于旋转图像的三维环境模型配准方法[J].中国科学杂志社2014441-:108118.,,()[29]RusuRB,BlodowN,MartonZC,etal.AligningpointcloudviewsusingpersistentfeaturehIR-istorams.ProceedinsofOS2008:384391.,g[J]g30JohnsonAEHebertM.Usinsinimaesforefficientobectreconitionincluttered3D[],gpggjscenes[J].IEEETransactionsonPatternAnalysis&MachineIntellience,2002,g2-15:433449.()31RusuRBBlodowNMartonZCetal.Alininoinloudviewsusinersistentfeature,,,tc[]ggpgphistoramsC//Ieee/rsInternationalConferenceonIntellientRobotsandSstems.IEEE,g[]jgy20-13:33843391.32RusuRBBlodowNBeetzM.FastointfeaturehistoramsFPFHfor3DreistrationC,,//[]()[]pggIEEEInternationaitiP20091-1lConferenceonRobotcsandAutomaon.IEEEress:848853.,33SaltiSTombariFStefanoLD.SHOT:Uniuesinaturesofhistoramsforsurfaceand,,qgg[]riind-uteriiIitexturedesctonComVson&maeUerstandn20141258:251264ppgg,,()34GuoYSohelFABennamounMetal.RoPS:Alocalfeaturedescritorfor3Driidobects,,,[]pgjbasedonrotationalroectionstatisticsC//InternationalConferenceonCommunications,pj[]SinaiiiI-lProcessinandTherAlcatons.EEE2013:16.gg,pp,[35]YangJ,CaoZ,ZhangQ.Afastandrobustlocaldescriptorfor3Dpointcloudregistration[J],norma--IftionSciences2016s346347:163179.,,[36]陆军,彭仲涛,董东来,等.点云FPFH特征提取优化配准算法[J].新型工业化,20-1477581.()58 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果ndaneta-uarestocloud37CheWebl.Alicationofconstrainedtotalleastsoint,,[]ppqpnam-reistration.ournalofeodcs7Geodes&Gi31213141.,(,gJyy)38.陈茂霖.中,卢维欣,万幼川,等无附加信息的地面激光点云自动配准方法[J]国激光,[]20-164:212220.()王赫.39.视觉测量点云数据拼接方法及关键技术研究D.哈尔滨理工大学2011,[][]-40王力20:1...测1136244145,李广云,贺嘉,等点云拼接中标志自动匹配方法J绘科学,,[][]()41WanXZhaoZLCasAGetal.Aniterativeclosestointaroachforthereistration,,[]g,pppppgofvolumetrichumanretinaimaedataobtainedboticalcoherencetomorah.gypgpy[J]Mu20-ltimediaTools&Alications1:115.pp,6ZhIIreistin3-Dmodel42iianTUKeZYanC,etal.mrovedCPration[]qg,,gpgreconstruction.TransactionsoftheChinaWeldinInstitution2013.,[J]g43ElbazGAvrahamTFischerA.3DPointCloudReistrationforLocalizationUsinaDee,,[]ggp-iiernRecoon207NeuralNetworkni.AutoEncoderC//ComuterVsonandPattti.1[]pgi-Dimint-oudataBaon44GeBPenBTian.RestratiofThreeensionalPoClDsed,,no[]gQgTiinii20274-1CurvatureMa.ournalofanUnverst1346:180.,,p[J]Jjy()[45]ChenQjKoltunV.RobustNonrigidRegistrationbyConvexOptixnizationfC]//IEEEuteriicie-InternationalConferenceonComVson.IEEEComuterSot2015:20392047.ppy,iPiirarThree-Dimsiona46AldomaAartonZCTombariFetal.Tutoral:ontCloudLb:enlM,[],,yObectReconitionand6DOFPoseEstimation[J].Robotics&AutomationMagazineIEEE,jg201213-:8091.,9()47BastianStederRaduBodanRusuKurtKonolieetal.NARF:3DRaneImaeFeatures,,,[]ggggforbectReconition.1O200.fJ]jg48htts://www.cc.atech.edu/roects/laremodels/[]pgpjg_-49王丽辉.2.袁保宗.三维散乱点云模型的特征点检测信号处理2011276:93938,,[],[I]()50刘斌郭际明邓祥祥.测绘科学2016.结合八叉树和最近点迭代算法的点云配准,,[],,[J]-412:130132.()D.201351靳洁..基于小波分析的地面三维激光扫描点云数据的滤波方法研究长安大学,[][]52RusinkiewiczSLevoM.EfficientvariantsoftheICPalorithmC//International[],yg[]--ConferenceonDDiitalImainandModeli2001Pin.I1.3n.roceedsEEE2002:14552ggg,,gg[53]OomoriS,NishidaT,KuroiS.Pointcloudmatchinusinsingularvaluedecomposition^].gggrtificialLifeobotics212-A:154.&R2016149,,()[54]YangDE,MuWP,SunWK,etal.Coregistrationofterrestriallidarpointsbyadaptive-itutomascaleinvariantfeaturetransformationwhconstrainedeometr.Ationingy[J]onstrucon25-Cti201225:4958.,,()55徐万鑫许宏丽C//ChineseConferenceonSstem.改进的ICP算法在点云配准中的应用[],[]ySimulationTechnology&Application.2012.一2056梁群仙..计算机工程,13,,许宏丽种基于点云数据的快速曲面重构方法[J][]39-2:237240.()57ohnsonA.EHebertMUsinsinimaesforefficientobectreconitionincluttered3D.[]J,,,gpgjgscenes.IEEETransactionsonPatternAnalsisandMachineIntellience,2002,[J]yg2-15:433449.)(59 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果[58]GuoY,SohelFA,BennamounM,etal.RoPS:Alocalfeaturedescriptorfor3Drigidobjectsbasedonrotationalprojectionstatistics[C]//InternationalConferenceonCommunications,iandTher-SnalProcessinilicatonsEE20:Ai.IE1316.gg’pp,[59]SeninN,ColosimoBM,PacellaM,PointsetaugmentationthroughfittingforenhancedICPregistrationofointcloudsinmultisensorcoordinatemetrology.Roboticsandp[J]u--ComterInteratedManufacturin201329l:3952.pgg,,()[60]郑德华,岳东杰,岳建平.基于几何特征约束的建筑物点云配准算法田.测绘学报,2008,374464-:468.()61秦绪佳.ICP拼接.,徐菲,王建奇,等基于法向量直方图特征描述的点云小型微型计算[][J]2016373593-597机系统.,,()[62]MenH,GebreB,PochirajuK.Colorpointcloudregistrationwith4DICPalgorithmic]//IEEEIionalConferenceonRoboticsandAutomation201-nternat.IEEE1:15111516,.An-63BeareeRDieulotYRabateP.innovativesubdivisionICPreistrationmethodfor[],J,g-toolathcorrectionaliedtodeformedaircraftartsmachinin.Internationalournalppppg[J]Jof-AdvancedManufacturinTechnolo2010535:463471.g,,()gy64SchnabelRKleintree-basedont-cloudcomression//EurorahicsIEEEtcR.OciC/V[],pp[]gpg--ConferenceonPointBasedGrahics.Eurorahicsssocaton2006:.Aii111121pgp,65SchauerNiichterA.Collisiondetectionbetweenointcloudsusinanefficientk[]J,pg-dtreeiiAiIi2015293440-458mlementaton.dvancedEnneerinnformatcs:.p^,,]gg()[66]XingBB,ChengH,JingLW,etal.Studyonastitchingalgorithmoftheiterativeclosestontbasedondnamichierarchournaloftcalechnolocofticheski.OiTC/Oiipyy[J]JpgypZhuma-l201582l:2832.,,()[67]SunJH,XieP,LiuZ,etal.Automatic3Dpointcloudregistrationbasedonhierarchicalnnn-blocklobalsearch.Otics&PrecisioEieerin201321174.:1180g],,[Jpg()g68VetterMHofleBHollausMetal.VerticalVeetationStructureAnalsisandHdraulic,,,[]gyy-RouhnessDeterminationUsinDenseALSPointCloudDataaVoxelBasedAroach.ggppfJ]-iiISPRSInternationalArchivesofthePhotorammetrRemoteSensnandSatalgy,gpInformationSciences2011.,[69]DinaA.Hafiz,BayumyA.B.Youssef,WalaaM.Sheta,etal.Interestpointdetectionin3Dintc3DSobe-HisInolouddatausinlarroerator.ternationalournalofpgp[J]JPatternRecognition&ArtificialIntelligence,2015,29(07):150710022059009.70马骊溟,徐毅,李泽湘.基于高斯曲率极值点的散乱点云数据特征点提取.系统仿真学[][J]2008209-报:23412344.,,()[71]PankajDS,NidamanuriRR.Arobustestimationtechniquefor3DpointcloudreistrationImaenalsstereolo.Ai&S2016351.g[J]gygy,,()[72]朱德海.点云库PCL学习教程[M].北京航空航天大学出版社,2012.""Crauer-rchardtnknuhm73M.illM.d.BBuMoerC.MuetHeizeP.KstetandG.Notni[],,,,,,"<(Oofoicorresondencemei-iintheaccuracntthodsnthreedmensonalmeasurementyppMmsu-sstesifrieroectionOt.En.5260636010636012013nn.y,,,,ggpjpg()()74SalviMataboschCFofiDetal.Areviewofrecentraneimaereistrationmethods,,,[]Jggg-withaccuracevaluation.Imae&utin.VisionCom2007255:578596y,,[J]gpg()60 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果[75]AridT.IntroductiontoprogramminwithKinect:Understandinghand/arm/headgmotionandsokencommands[C]//SinalProcessinandCommunicationsAlicationspggppConferenceSIU20IEEE-1220th.2012:11.,,()76Zhan.AdacentScatterPointCloudReistrationTechnoloResearchBasedon[]gXjggy-ThreedimensionalCameraofTimeofFliht.urnalofMechanicalEnineerin2013J,,g[]Jogg12-49.:81336()[77]梁锡坤.B样条类曲线曲面理论及其应用研究[D].合肥工业大学,2003.78.D.2010莫堃华中科技大学.基于隐式函数的曲面重构方法及其应用,[][]79LawsonCLGenerationofaTrianularridWithAlicationtoContourlottn.GPi.[]gppg[J]MemoJetPropulsionLaboratory,1972.80鲍忠贵编程快速入门M.清华大学出版社2016.王涛,陈凌眼Q,[],j[]61 北京交通大学硕士学位论文作者简历及攻读硕士学位期间取得的研究成果作者简历及攻读硕士学位期间取得的研宄成果一、作者简历作者:张哲性别:女出生年月:1993年9月教育经历:-202011年9月15年7月北京工商大学计算机与信息工程学院软件工程获工学学士学位。-202015年9月18年4月北京交通大学计算机与信息技术学院计算机科学与技术攻读工学硕士学位。二、发表论文一1尹辉.[]张哲.种基于关键点选择的快速点云配准算法[J]激光与光电子学进,许宏丽,展,20175412):121002,(三、参与科研项目1国家自然科学基金No.6147202961473031)?(,[]2No.[]科技部国家重点研发计划.2016YFB1200100)(No.4152042?[3]北京市自然科学基金()4研业务费专项资金(N〇.2016BZ0052016JBM019,2016JBM016).[]中央高校基本科J,62 北京交通大学硕士学位论文独创性声明独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作和取得的研宄成果,除了文中特别加以标注和致谢之处外,论文中不包含其他人己经发表或撰写过的研宄成果,也一不包含为获得北京交通大学或其他教育机构的学位或证书而使用过的材料。与我同工作的同志对本研宄所做的任何贡献均已在论文中作了明确的说明并表示了谢意。:曰学位论文作者签名:签字日期七松年¥月《63 北京交通大学硕士学位论文学位论文数据集学位论文数据集表1:.1数据集页关键词*密级*论文资助中图分类号UDC1点云配准;多视公开TP39基于计算机视觉角点云;点云精的高速铁路轮轨简;关键点;特接触关系检测关征描述子键技术研宄(No:K16JB00240)学位授予单位名称*学位授予单位代学位类别*学位级别*^北京交通大学10004IL学硕士|_论文题名*并列题名论文语种*基于结构关联的维点云配准研宄I中文^"**5120470作者姓名学号1|张哲培养单位名称*培养单位代码*培养单位地址邮编|北京交通大学10004北京市海淀区西直100044门夕卜上园村3号^学科专业*研宂方向*学制*学位授予年*&机科学与技算机视觉U.5年2018,|计丨论文提交*日期2018.03*\*导师姓名午宏丽职称|教授|评阅人答辩委员会主席*答辩委员会成员常冬霞高博鲁凌云|电子版论文提交格式文本(V)图像()视频(()多媒体()其他())音频:推荐格式application/msword;aliti/dfppcaonp5版论文出版(发布)者电子版论文出版(发布)地权限声明||*仑文总页数3|6*33项,其中带为必填数据,为21项共。64

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
关闭