基于TerraScan软件的点云数据滤波处理研究

基于TerraScan软件的点云数据滤波处理研究TheResearchOffilteringofpointclouddataBasedonTerraScansoftware高志国(广州市城市规划勘测设计研究院，广东广州510060)Gaozhiguo(GuangzhouUrbanPlanningDesign&SurveyResearchInstitute,Guangzhou510060,China)摘要:点云滤波是地面三维激光扫描数据处理的重点和热点问题，也是点云数据处理的一项关键工作。阐述了目前几种主要的点云滤波方法的基础上，并着重探讨了基于TerraScan软件的点云数据处理方法，如误差点剔除、点云分类、模型重建等。经实验数据分析取得了良好的滤波效果，能有效滤除点云噪声，同时很好地保持了地形特征。关键词：三维激光扫描；数据处理；数据分类；点云滤波Abstract：PointCloudisathree-dimensionallaserscandataprocessinggroundfocusandhotspots，pointclouddataprocessingisakeywork．severalkeypointsofthecurrentcloudfilteringmethod，basedonfocusedonthesoftware-basedTerraScanpointclouddataprocessingmethods，suchastheerrorpointremoved，pointcloudclassification，modelreconstruction．Theexperimentaldataobtainedagoodfilteringeffect，caneffectivelyfilterthepointcloudnoise，whilemaintainingitsedgefeatures．Keywords：3Dlaserscanning；dataprocessing；Dataclassification；PointCloudFilter引言地面三维激光扫描技术（TerrestrialLaserScanning,TLS）是20世纪90年代发展起来的一种快速获取空间三维信息的新技术手段，并广泛应用于测绘技术行业，带动了测绘技术的全新发展，使得空间数据的获取从传统的单点数据采集向连续的、密集型自动数据获取技术转变，为测绘行业注入了新的活力，推动了测绘技术不断向前发展，同时也使专业的数据处理人员突破了传统模式，改进了数据处理的方法，并进一步探索新的数据处理理论和方法。然而相对三维激光扫描技术硬件的发展的成熟度而言，针对三维激光扫描技术所获取的实体的点云（PointCloud）数据的后处理等的发展还相对滞后，而且在测绘领域都还没有形成一套标准的数据处理方案，因此对数据处理的研究显得尤为重要，使数据处理更进一步地向自动化、智能化的方向发展[10]。1基于TerraScan模块的数据处理1.1Terrasolid系列产品介绍TerraSolid系列软件是第一套商业化LiDAR数据处理软件，基于Microstation开发的，运行于Micorstation系统之上，因此需要熟悉Microstation操作的用户才能很好的使用它。它包括：TerraMatch、TerraScan、TerraModeler、TerraPhoto、TerraSurvey、TerraPhotoViewer、TerraScanViewer、TerraPipe、TerraSlave、TerraPipeNet等模块。TerraSolid系列软件能够快速的载入lidar点云数据，在足够内存支持下（2G），载入39000000个点只需要40多秒。Terrasoild软件因受限于Microstation，浏览起来闪烁较为明显，速度不快。1.2TerraScan软件模块TerraScan是TerraSolid公司用来处理LiDAR点云数据的软件。他能够读取除了标准激光雷达数据格式*.las文件外，还可以读取*.xyz数据格式文件，因此为其他扫描仪提供了良好的数据接口。TerraScan的基本功能如下：·以xyz文本或类似于LAS和TerraScan的二进制文本读入原始的激光点云·三维方式浏览数据·自定义点类别·激光点自动手动分类·交互式判别三维目标（如：铁塔）·数字化地物·探测电力线·矢量化房屋·生成激光点的截面图·输出点分类TerraScan作为TerraSolid公司第一套成功、完善的商业化雷达点云数据处理软件的一个主要模块，具有能够完成大量的雷达点云数据的分类处理等工作，对该软件得操作简便易行。并且能根据用户对不同数据分类要求，完成数据读取、分类、处理、建模等相关工作。TerraScan软件还可以根据用户的不同需要对数据进行分类处理，结合该软件中的其他数据处理模块，（如TerraModel模块、TerraPhoto模块、TerraStreet模块等）完成对点云数据的二次分类，并通过人机交互等操作对地物进行建模等。1.3TerraScan的分类原理TerraScan是基于利用改进的不规则TIN网加密方法进行分类的。首先，通过迭代计算初始化一个稀疏不规则TIN网，事先设定TIN网阀值并通过计算将每次满足条件的点添加到TIN网中。然后，对不满足先前条件的点重新构建新的不规则TIN网，并重新计算新的阈值条件，对后续点进行同样的筛选判断。重复操作，直到不再有新点加入为止[7]。由于本次论文数据是采用徕卡地面三维激光扫描仪LeciaScanStation2获取的，其数据文件格式未知，而其工程文件格式为*.imp格式，因而不能直接读取到TerraScan软件进行后续的处理等操作，但可以通过数据格式之间的转换来完成相互的操作，TerraScan可读取的数据格式为*.dte文件、标准的激光雷达数据格式*.las和数据格式为*.XYZ的数据文件。因此可以将数据直接读入LeciaScanStation2扫描仪采集的点云数据进行相关操作，用徕卡本身的cyclone数据后处理软件能够提供的多种数据格式输出与转换接口，通过数据输出格式转换成TerraScan软件能够直接读取的*.XYZ数据格式读入软件见图,1示。图1数据格式转换与激光点云数据的读入1.4TerraScan的分类过程TerraScan工具箱为数据后处理提供了丰富的分类工具，使用户能够方便快捷地对点云数据进行分类工作。三维激光扫描仪接收返回的激光信号不仅包括地面反射，同时也包含多种误差信息等。因此需要将点云数据中各种类型的数据和误差等区分开来。TerraScan软件模块中提供了多种数据分类的类别。同时，用户还可以根据自己的实际或工程需要，添加新的所需的分类类别，每个激光信号分类类别都是唯一的，如图2所示。图2TerraScan点分类及其属性如图2所示，使用者可以根据实际需要，使用“Add”命令对分类的类别进行添加，使用“Edit”命令对上述的分类类别属性等信息进行修改等操作。(1)分离低点。所谓低点，就是指高程明显低于周围区域内所有点的高程的一些点。分离低点通常对数据点击进行初始化，寻找高程明显低于周围区域数据点集的高程的错误点。并以该点为中心在一定范围内与其他数据点的高程进行对比。若需要判断的点的高程明显低于目标中心点的高程，则将其归为低点。若遇到低点密度较大时使用单点分类会造成错误，这时需要采用数据分类的方法来分离低点，其具体操作过程及参数设置如下图所示。图3分离低点的参数设置界面(2)分离独立点。所谓独立点，是指在一定三维空间分布范围内分布比较稀疏的一些点，因此根据上述独立点的分布原则建立独立点的分离准则，其具体步骤为：首先寻找以空间独立点为目标点，以此点位中心建立一定范围的空间三维缓冲区，如果在该缓冲区内的点数少于规定的点数阀值，则该点就为独立点。图4分离独立点参数设置界面(3)分离空中点。所谓空中点，顾名思义就是指其高程明显高于数据点集周围所有点高程均值的一些点。分离空中点的基本原则同上述独立点的基本类似，首先设置需要判断的某点为目标点为中心点，设定一个三维搜索半径的缓冲区，包含在这个三维缓冲区内的点看做是目标点的邻近点。将目标点的高程与邻近点高程均值大小做比较，同时做差计算，若其差值达到高程标准偏差规定的倍数。就将其归类为空中点。图5分离空中点参数设置(4)分离植被点。在TerraScan中将植被信息根据高程值可将植被点归为三类：Lowvegetation、Mediumvegetation、Highvegetation,其中Lowvegetation高程值在0.5m～2m之间，Mediumvegetation在2m～5m之间，Highvegetation5m。分离植被点时，以从地面起的高度为准，先从地面以上的点分离出低的植被点，其基本的思想是首先在地面点建立临时的三角形模型，然后这个三角形模型的高程值与区域中的其他点高度比较，如果低于0.5m，就认为这个点属于低矮植被进行分离。然后根据相同的算法，从余下的点中分离出中等高度的植被点，最后分离出高植被点。图6TerraScan植被信息的提取操作界面(5)分离建筑物点。分离出的中、高植被中包含有建筑物点，因此从高植被点中分离出建筑物点。建筑物信息的提取应遵循的原则是建筑物上的点能形成一个平面。因此首先要完成地面点的分类，在此基础上再进行建筑物分类。图7建筑物分离操作界面(6)地面点集的获取建立地表初始不规则TIN网模型并通过不断的迭代计算来分类地面点。通过选择最大建筑物尺寸参数来确定初始点的选取。如一个建筑物的最大尺寸是50m，那么在50×50m的区域内最少有一个地面点存在，认为其中的最低点就是地面点。首先用已经选择好的最低点作为初始地面点建立不规则TIN网模型。然后通过不断的迭代计算、不停地加入新激光点来拟合初始的地面模型，通过不断地选择一些可以接受的新点到地面模型中的重复参数来更新地面模型使得模型更接近地表一些，最终得到一个近似地面模型。图8提取地面点集操作界面(7)提取Modelkeypoints。Modelkeypoints是提取DEM所用的点，为后面提取DEM方便，在分离Modelkeypoints的同时，要将点输出，选择macrostep-action–outputpoints,设置选择保存Modelkeypoints的文件路径，并设置好提取建立DEM关键点的参数，对其进行提取输出。图9Modelkeypoints的提取与文件输出提取Modelkeypoints是在已经获取了地面点集的基础上进行的，根据精度要求建立地面TIN网。地面TIN网模型精度受到参数AboveModel和BelowModel影响。前者是确定激光高程超过模型的高度阀值，后者是确定激光高程能地于模型的高度阀值。建立初始的地面TIN网后，通过不断的迭代寻找Modelkeypoints直至满足要求为止。最后根据分类数据获得的结果进行滤波处理，过滤掉数据点集中的非地面点，由地面点集中提取出来的ModelKeyPoints进行数据建模，生成数字高程模型，即DEM。图10TerraScan分类滤波后模型2结论地面三维激光扫描测量技术以其高效率、高采样率、高精度、高密度、非接触性、主动性、实时性等优点，在测绘科学与技术等领域有广泛的应用前景。随着技术的不断成熟和发展，测量获取的点云数据量会越来越庞大，因此点云数据的处理技术越来越受到人们的关注，本文以徕卡地面三维激光扫描系统获取的点云数据为对象，对点云数据的处理做了一些研究，得到了一系列结果，总结如下。对于点云数据的滤波处理，本文利用cyclone软件结合TerraScan两种软件，对获取的点云数据进行分类滤波处理。TerraScan中提供的分类方法都是基于面的分类理论，即假设分类过程目标点的计算值与周边一定范围内的点都有关系。这种假设比较适合实际点的分布情况，因为空间点不是独立存在的，与周边点相互关联，从而揭示了空间点间结构上的联系