利用机载激光点云数据生产DEM的关键技术分析

蒋桂美，聂 倩，陈小松

(宁波市测绘设计研究院，浙江 宁波 315042)



利用机载激光点云数据生产DEM的关键技术分析

蒋桂美，聂 倩，陈小松

(宁波市测绘设计研究院，浙江 宁波 315042)

着重分析了利用机载激光点云数据提取DEM的关键技术：预处理、点云滤波与分类、高精度DEM制作与质量评价技术；建立了基于机载激光点云的DEM生产技术流程，并将其应用于宁波市地理国情普查市情专项——高精度地表模型制作项目中；构建了宁波市建成区规则格网1∶2000比例尺的高精度DEM模型，为宁波市地理国情普查提供了数据基础。

机载LiDAR；激光点云；DEM

机载激光扫描技术(LiDAR)是20世纪80年代中期逐渐发展起来的一项高新技术，该技术通过位置、距离、角度等观测数据直接获取对象表面的三维点云信息，具有数据获取速度快、空间与时间分辨率高、动态探测范围大、主动性强等特点，广泛应用于城市大比例尺测图、林业调查、电力巡线等各个领域[1-2]。与传统航空摄影测量技术相比，机载LiDAR技术不受日照、天气影响，可全天候主动快速获取精确的高分辨率数字地面模型及地物的三维坐标，具备采集速度快、高程精度高、成图周期短等特点。因此，利用机载LiDAR数据快速完成数字高程模型(DEM)的大规模生产，已成为最近几年测绘生产及其他领域的一个研究和应用热点。

本文结合机载LiDAR技术的工作原理和点云数据特点，着重分析利用机载激光点云数据提取DEM的关键技术，即激光点云滤波与分类、点云数据内插、DEM构建与评价，建立基于机载LiDAR数据的DEM生产技术流程，并将其应用于宁波市地理国情普查市情专项——高精度地表模型制作项目中，构建宁波市建成区规则格网1∶2000比例尺的高精度DEM模型。项目试验结果表明，该技术流程与方法正确可行，成果质量优良，满足测绘产品的要求。

1 机载LiDAR技术工作原理

机载LiDAR技术是全三维测量模式，它以飞机为载体，集成激光扫描系统、GPS和惯性测量单元(inertial measurement unit，IMU)3个基本的数据采集工具及控制单元系统等其他设备于一体，从而获取地面的三维空间信息。其中，激光扫描系统以脉冲激光测距作为主要技术手段，以激光束扫描的工作方式测量从传感器到地物对象的激光照射点间的距离，即通过测量地面采样点激光回波脉冲相对于发射激光主波之间的时间延迟得到传感器到地面采样点之间的距离[3-6]，同时还可以获得反射率、激光脉冲回波次数等信息；GPS接收机用于确定激光发射点的空间位置；IMU测量激光发射瞬间激光的空间姿态参数。机载LiDAR系统通过激光扫描系统获得激光点的距离信息，结合由惯性导航系统输出的激光发射点空间位置和姿态信息，可以解算出激光脚点在当地水平坐标系下的三维坐标，即一系列离散的、空间分布不规则的三维点云数据，具体流程如图1所示。

图1 机载激光点云数据处理技术流程

2 关键技术研究

2.1 数据预处理

机载激光点云数据预处理主要是将航飞过程中获得的差分GPS数据和IMU数据进行联合解算获得定位定姿数据，然后利用联合定位定向数据、原始点云信息及系统所提供的各类参数进行解算，获得每一个激光点的在WGS-84坐标系下的三维坐标。根据项目成果需要，将LiDAR点云数据由WGS-84坐标系转换到地方坐标系或CGCS2000坐标系，所得结果为随机分布的带有高程、位置和强度信息的激光点云。此外，由于机载LiDAR系统获得的海量点云数据中包含了大量的粗差和系统误差，数据预处理中应进行粗差检测和剔除处理，即剔除极高点和极低点。

2.2 激光点云滤波与分类

机载LiDAR系统获取的原始点云数据是离散的孤立点，其主要的数据值为回波信号点的三维空间坐标及一些附带的属性信息(如强度、反射波次数等)，点与点之间不存在任何拓扑关系[7-9]。机载激光点云滤波技术是指从离散的点云数据中区分出地面点和非地面点的过程，其基本原理是基于邻近激光脚点间的高程突变(一般不是由地形的突然起伏变化所造成的)，进行滤波计算时需要设置一定的阈值，判断激光脚点是否位于地形表面[10-11]。利用机载激光点云数据制作DEM时滤波与分类具体流程可概况为：按照回波次数分类—地面点分类—水系分类。

(1) 按照激光回波次数分类。由于机载激光具有穿透性，导致在扫描过程中不同的地物有着不同的回波次数和强度信息，当激光脉冲照射到建筑物顶部或裸露的地表时只产生一次回波，而激光脉冲照射到植被时，由于激光信号可以穿透植被从而形成多次回波。因此，可对多次回波中的最后一次回波数据和仅一次回波的激光数据进行提取，从而滤除大量非地面点。

(2) 地面点自动分类。通过反复建立地表三角网模型的方式分离出地表点，并反复加入新的激光点扩展模型，确保每个加入的激光点使得模型更加贴近地表，从而最终得到一个近似的地表面[12-13]。具体分类方法为：首先设定测区内的最大房屋尺寸，选择少许局部激光低点，确保该局部激光低点间的距离不小于最大房屋尺寸；其次基于上述低点构建初始三角网，假定该三角网的最高顶点贴合地面，并通过不断迭代将该初始三角网上方的点加入以构建新模型，从而得到一个近似的地表面。

(3) 人机交互分类。自动分类后的数据可以获得概略地形地貌，但存在一些错分或误分激光点云数据，同时一些小的地形不连续部分也会被平滑或去掉，因此需要通过人机交互方式分类。人机交互分类方法是指参照高分辨率影像，通过人工判断的方式对自动分类的结果进行修正或对特定地物类别进行分类。可通过建立不规则三角网模型(TIN格网)，对区域内的整体点云进行分类，剔除坑点或突出点，细部可再根据横截面修正。由于建筑物规模和形态变化较为复杂，需以人机交互分类为主。水系一般无激光回波，但也存在浑浊水体稀疏回波现象，需要结合影像和周围地形进行判断。

2.3 DEM构建与评价

DEM主要基于TIN格网构建，即对于滤波分类的地面点云数据，利用ArcGIS软件创建TIN，构建约束三角网模型，并进一步生成规则格网Grid格式的DEM数据[14-15]。

3 试验与分析

3.1 试验数据

本文试验数据采用宁波市第一次地理国情普查市情专项——高精度地表模型制作项目所实施的机载激光雷达航飞数据，具体试验区为宁波慈城镇，如图2所示。该地区三面环山，地形北高南低，采用传统测量方法获取高精度DEM数据需要大量的人力物力投入。该区域LiDAR数据获取时间为2015年底，采用瑞士AHAB公司的蝙蝠机载激光雷达水文/陆地测量系统，以运-12固定翼飞机作为飞行平台，点云密度为每平方米不少于6个点，影像分辨率优于20 cm，要求制作1∶2000比例尺精度的DEM成果。图3为获取的机载激光点云数据。

图2 慈城镇区域

3.2 试验结果与分析

利用TerraSolid软件对试验区激光点云数据进行预处理、滤波与分类、点云数据内插，最后得到DEM。图4为所有激光点云形成的DSM，图5为点云滤波分类后地面点所形成的DEM产品。

对于DEM成果质量，内业采用人机交互检查的方法检查发现，点云数据存在的少数漏分或错分现象；外业采用实测并与1∶500地形图比对的方式检查DEM的精度情况，其精度统计见表1。检测结果表明：采用本文技术路线制作的DEM中误差小于10 cm，满足项目设计的要求。

图3 机载激光点云数据

图5 DEM结果

4 结 语

本文主要介绍了利用LiDAR激光点云数据制作高精度DEM的技术流程和方法。试验表明，本文的技术路线能够满足现阶段城市大比例尺DEM生产的需要。

表1 DEM成果精度统计

[1] 张小红.机载激光雷达测量技术理论与方法[M].武汉：武汉大学出版社，2007：9-14.

[2] HAMMOUDI K，DORNAIKA F，SOHEILIAN B，et al. Extracting Wire-frame Models of Street Facades from 3D Point Clouds and the Corresponding Cadastral Map[J].Remote Sensing and Spatial Information Sciences， 2010,38(3A)：91-96.

[3] 周淑芳，李增元，范文义，等.基于机载激光雷达数据的DEM获取及应用[J].遥感技术与应用，2007，22(3)：356-360.

[4] 林昀，吴敦，李丹农.基于机载激光雷达的高精度电力巡线测量[J].城市勘测，2011(5)：71-74.

[5] 杨晓云.一种基于形态学的激光雷达数据滤波算法的改进[J].测绘科学，2011，36(4):166-168.

[6] 王继承.机载激光雷达技术和摄影测量匹配技术在DEM生产中的应用[J].测绘与空间地理信息，2015(5)：128-129.

[7] 周晓明.机载激光雷达点云数据滤波算法的研究与应用[D].郑州：信息工程大学，2011.

[8] 郭波，黄先锋，张帆，等.顾及空间上下文关系的JointBoost点云分类及特征降维[J].测绘学报，2013，42(5)：715-721.

[9] 黄先锋，李卉，王潇，等.机载LiDAR数据滤波方法综述[J].测绘学报，2009，38(5)：466-469.

[10] 熊爱武，杨蒙蒙.机载LiDAR点云数据误差分析[J].测绘通报，2014(3):75-78.

[11] 胡永杰，程朋根，陈晓勇，等.机载激光雷达点云滤波算法分析与比较[J].测绘科学技术学报，2015(1):72-77.

[12] 曹广强，瞿娜.浅谈利用激光点云数据制作DEM[J].测绘标准化，2011，27(2)：44-46.

[13] 王庆，陈磊，韩文泉.基于机载LiDAR技术的高精度DEM生产研究[C]∥第九届长三角科技论坛(测绘分论坛)论文集.宁波：[s.n.]，2012：114-116.

[14] 段佳，孙敬宇，刘海飞，等.利用机载激光三维点云构建张家界地貌DEM[J].测绘通报，2015(10)：107-109.

[15] 李鹏程，王慧，刘志青，等.一种从机载LiDAR点云数据获取DEM的方法[J].测绘通报，2012(5):59-62.

Key Technologies Analysis of DEM Production Based on LiDAR Data

JIANG Guimei,NIE Qian,CHEN Xiaosong

(Ningbo Institution of Surveying and Mapping, Ningbo 315042, China)

The paper focuses on the analysis of the key steps of extracting DEM from airborne-LiDAR data， including point cloud data preprocessing, the data filtering and classification, high precision DEM manufacturing and quality evaluation. In this paper, the DEM production technology process based on airborne-LiDAR data is established, and will be used onacitysituationspecialsurveyofgeographicalconditions：highprecisionsurfacemodelingprojectinNingbo. A high precision DEM model of 1∶2000 scale rule grid of Ningbo built-up area is constructed which provides a data base for the survey of geographical conditions in Ningbo.

airborne-LiDAR; laser cloud points; DEM

蒋桂美，聂倩，陈小松.利用机载激光点云数据生产DEM的关键技术分析[J].测绘通报，2017(6)：90-93.

10.13474/j.cnki.11-2246.2017.0197.

2016-10-26

宁波市科协资助项目(X201601-3)

蒋桂美(1982—)，女，工程师，研究方向为摄影测量与遥感、地理信息系统。E-mail：xiaomeijgm@126.com

P23

A

0494-0911(2017)06-0090-04