融合航空影像的震后机载LiDAR建筑物点云提取[PDF全文]
邓飞, 窦爱霞, 王晓青
摘要: 基于机载激光雷达LiDAR(Light Detection and Ranging)数据识别震后建筑物震害,其前提是快速准确地提取建筑物点云。通过分析地震灾区机载激光雷达点云中提取建筑物点云的诸多难点,已有的方法难以达到预期效果,因此提出融合同机航空影像数据的方法,实现了震后灾区建筑物点云的获取。该方法首先在数据预处理的基础上,利用布料模拟滤波CSF(Cloth Simulation Filtering)算法进行点云滤波,得到地面点云和非地面点云(主要是建筑物、植被和车辆行人等),并将航空影像红波段光谱信息赋予非地面点云;然后基于灰度直方图阈值分割的方法剔除植被点;最后对剩余激光脚点利用具有噪声的基于密度的空间聚类DBSCAN(Density-Based Spatial Clustering of Applications with Noise)算法进行聚类提取最终的建筑物点,并与参考建筑物点比对,进行精度验证,得到建筑物点云提取的漏检概率、虚警概率分别为15.61%、7.52%,总体精度可达84.39%。结果表明,在一定精度要求范围内,该方法能有效实现地震灾区建筑物点云的提取,可为震后机载LiDAR建筑物点云提取提供技术参考和方法借鉴,为建筑物震害识别做好基础工作。
关键词: 机载LiDAR     航空影像     地震灾区     建筑物点云     DBSCAN聚类算法    
DOI: 10.11834/jrs.20187191    
收稿日期: 2017-05-31
中图分类号: TP79    文献标识码: A    
作者简介: 邓飞,1991年生,男,硕士研究生,研究方向为震后建筑物提取与震害识别。E-mail:feid163@163.com
通信作者: 窦爱霞,1977年生,女,副研究员,研究方向为遥感在地震灾害中的应用。E-mail:axdothy@163.com.
基金项目: 国家自然科学基金(编号: 41404046);中国地震局地震预测研究所基本科研业务专项(编号:2015IES0403)
Fusion of aerial imagery and airborne LiDAR data for post-earthquake building point extraction
DENG Fei, DOU Aixia, WANG Xiaoqing
Institute of Earthquake Science China Earthquake Administration, Beijing 100036, China
Abstract: Devastating earthquakes are one of the most serious natural disasters that affect mankind. Although humans cannot predict or prevent earthquakes, we can change the way we respond to them. Earthquake-induced building damage is one of the main causes of casualties and property loss. Extracting the information of post-earthquake buildings quickly and accurately is crucial for identifying earthquake-induced building damage and initiating effective rescue of survivors. The geographical environment is complex in the seismic field, and remote sensing technology is a potential alternative to capturing post-event scene data. Numerous studies have focused on researching intact buildings. The identification of building seismic damage based on airborne Light Detection and Ranging (LiDAR) is studied only by using object-oriented analysis, and few papers use the characteristics of building points. This paper proposes a method that extracts building points in the earthquake zone from airborne LiDAR point cloud data. This paper treats the fast and reliable extraction of the information of a post-earthquake building as a challenging task in which many difficulties arise because of airborne LiDAR point cloud data, and pre-existing methods cannot easily achieve the desired effect. Therefore, this paper presents a method that combines aerial images and LiDAR data to extract post-event building points in the affected quake areas. The approach consists of several steps. We first extract terrain points from the LiDAR point cloud by using the cloth simulation filtering algorithm, and off-terrain points are removed. Then, we manually collect training samples of building and vegetation points and identify the red band as the best spectral characteristic information to maintain two isolated points. This procedure is based on the minimum probability of total error. Afterwards, we assign the gray value of the aerial image red band to the terrain points (which mainly include buildings, vegetation, pedestrians, and vehicle). The vegetation points are eliminated from off-terrain points by using gray-level histogram threshold segmentation. We used the density-based spatial clustering of applications with the noise (DBSCAN) algorithm to extract the final building points from the rest of the laser footprints, which contain discrete noise points. We can obtain intact and destructive building points from the LiDAR point cloud by using techniques such as point cloud filtering and threshold segmentation based on the gray intensity of the red band. However, some low-density noise points on the road, such as vehicles, pedestrians, and vegetation, can be obtained. The high-density building points remain nearly constant, and discrete low-density noise points are removed effectively by using the DBSCAN algorithm. The extracted point cloud of a building is compared with reference points to quantitatively verify the extraction accuracy. The false-dismissal probability and false-alarm probability are 15.61% and 7.52%, respectively. The two main causes of undetected error are as follows: (1) The damaged building points, which have a similar height as the ground points, are filtered when the point cloud filter is executed. (2) The building points in the shadow of tall buildings or trees are divided into vegetation points during threshold segmentation. False-alarm error is due to crowded vehicle points in local areas, and these points cannot be removed by using the DBSCAN clustering algorithm. The two errors are small, and the total extraction accuracy is as high as 84.39%, thereby indicating that the extraction of building points in quake-hit regions satisfies requirements in some bounds of precision with the application of the aforementioned method. Results show that the suggested method can effectively extract building points in the earthquake zone within a certain accuracy requirement and provide technical reference for the extraction of post-earthquake intact or destroyed buildings by using airborne LiDAR point cloud data. The proposed approach could establish the foundation for identifying seismic building disasters and completing post-earthquake building damage assessment. With the development of laser radar technology, an airborne LiDAR remote sensing system can receive dense point cloud data while acquiring multispectral or hyperspectral data. Hence, we will improve the ability of the proposed approach to distinguish between vegetation and buildings by using threshold segmentation based on normalized difference vegetation index or other band combinations to further improve the accuracy of building point extraction in the seismic field.
Key Words: airborne LiDAR    aerial imagery    the earthquake zone    building point cloud    DBSCAN (Density_Based Spatial Clustering of Applications with Noise) clustering algorithm    
1、引 言

机载激光雷达是一种新型的测量技术,是获取高精度数字地面模型DEM(Digital Elevation Model)的重要手段(Liu,2008马洪超 等,2008周淑芳 等,2007),也是建筑物提取和城市3维建模的有效技术途径(Forlani 等,2006)。建筑物点云提取,是进行建筑物3维识别与重建的前提工作,利用LiDAR点云提取建筑物,难点在于判别、分离与建筑物高度、面积相近的植被(Awrangjeb 等,2012)。

目前,从数据源角度,建筑物点云提取分为仅基于机载LiDAR点云的方法、机载LiDAR点云与遥感影像相结合的方法。基于LiDAR点云,已有文献(Cheng 等,2013李亮 等,2016张爱武 等,2016)一般是在点云滤波的基础上,利用高程特性剔除车辆、行人点,根据建筑物点与植被点在高程纹理(高程差、高程变化、坡度等)、法向量、回波次数和强度等特征表现的差异性,利用这些特征、特征组合或衍生特征,基于反向迭代数学形态学、阈值分割、随机森林等方法进行点云的分类或建筑物点云提取。还有文献(陈光洲 等,2016赖旭东和万幼川,2005),通过点云内插,生成DSM深度影像,利用边缘检测的方法提取建筑物区域,从而可提取建筑物点云。随着机载雷达技术的发展,在获取3维激光脚点的同时,可获取航空影像,众多学者结合航空影像或卫星遥感影像提取LiDAR建筑物点云。融合遥感影像,主要是利用光谱特征差异解决建筑物与植被难以分离的问题。Lodha等人(2007)利用点云的高程(Height)及其衍射特征、回波强度(Return Intensity)和影像灰度(Image Intensity)等特征参数,采用最大期望EM(Expectation-Maximization)算法进行点云分类,提取建筑物点云。文献(Awrangjeb 等,2012Hartfield 等,2011Lee 等,2008)也均利用了光谱特征信息,首先根据高度阈值滤除车辆和低矮植被点,然后融合某一波段或是植被指数等光谱信息,将建筑物点云和高大植被点云进行分离。

以上研究对建筑物点云提取都具有一定的意义和效果,但主要是针对规则、完好建筑物的情况。基于机载LiDAR数据的震后建筑物点云提取或震害识别研究相对较少(Elberink 等,2011)。由于往往缺少震前机载LiDAR数据,难以用变化检测的方法(Vetrivel 等,2015),多数方法仅限于利用面向对象的方法(Ural 等,2011于海洋 等,2011)。对于场景复杂的震后灾区,窦爱霞等人(2013)分析了基于机载LiDAR识别建筑物震害的难点。受到地震影响,建筑物单体往往可能呈现倒塌、局部倒塌、未倒塌等多种状态(王晓青 等,2015),建筑物高度、屋顶发生变化。倒塌或局部倒塌建筑物与植被在纹理特征、坡度变化、法向量以及回波次数等方面表现相似的特性,难以继续利用这些点云特征参数进行二者的分离。此外,房屋建筑种类繁多、形状各式各样、材质复杂(水泥、彩钢棚、青瓦等),这给利用激光回波强度信息带来了困难。

基于以上困难或不足,本文尝试提出了一种融合航空影像光谱特征信息提取震后机载LiDAR建筑物点云的方法。该方法主要特点是借助航空影像光谱特征信息有效剔除难于与建筑物点区分的植被点,探讨性地利用DBSCAN聚类算法对空间低密度、离散的噪声点(车辆、行人点以及未滤除的少量植被点)进行识别,从而最终实现建筑物点云的提取。

2、建筑物点云提取方法

融合航空影像的震后机载LiDAR建筑物点云提取方法,首先对点云滤波处理,分离地面点云和非地面点云(或称“地物点云”,包括建筑物、植被、车辆和行人等),将航空影像光谱特征信息赋予非地面点云;然后,基于阈值分割的方法,剔除植被点;最后,根据车辆、行人点等噪声点分散、密度小的特点,采用基于密度的DBSCAN算法进行聚类,去除密度过小的噪声点,实现建筑物点云的提取。建筑物点云提取主要步骤有:点云滤波、基于阈值的分割和基于DBSCAN聚类的噪声点去除,整个过程如图1所示。

图 1 建筑物点云提取流程图 Figure 1 Flow chart of building point cloud extraction
    (2.1) CSF算法滤波

点云滤波实现地面点云和非地面点云的分离,目前常用的滤波算法主要有基于不规则三角网(TIN)滤波法(Axelsson,2000)、基于坡度/斜率滤波法(Susaki,2012)、数学形态学滤波(Zhang 等,2003)、移动曲面拟合法(苏伟 等,2009张小红,2007)和一些其他的滤波算法等。本文采用Zhang等人(2016)提出的布料模拟滤波(CSF)方法。该滤波方法具有应用于各种地形效果较好、设置参数个数少等优点。算法基本原理来源于一个简单的物理模拟过程:假设一块足够柔软的布料放置在地形之上,由于重力作用,布料下沉,并逐渐吸附在地行表面,最后布料的形状即为数字表面模型(DSM);若将地形倒置,具有一定硬度的布料放置其上,布料的形状即为数字地形模型(DTM),如图2所示。同理,将原始点云倒转,布料覆盖在该倒置点云表面,通过分析布料质点与激光脚点的相互作用,最终确定布料的形状。该物理过程采用“质点—弹簧(Mass-Spring)模型”的“布料模拟(Cloth Simulation)”技术实现。

图 2 CSF基本原理(Zhang等人(2016)) Figure 2 Basic principle of CSF (Zhang et al, 2016)
    (2.2) 阈值分割

阈值分割用于剔除非地面点云中的植被点。根据绿色植被吸收红、蓝光,对绿光具有反射作用,在近红外波段范围有强反射作用的光谱特性(林文鹏 等,2010),对波段间进行数值运算,构建植被指数,用于植被识别或提取。常用且特征明显的归一化植被指数(NDVI)需近红外波段参与运算,而本文航空影像仅有红(R)、绿(G)、蓝(B)共3个波段,无法构建NDVI。因此本文考虑基于单波段灰度值的阈值分割。采集建筑物、植被点云训练样本,提取红、绿、蓝3个波段灰度值,确定训练样本的灰度分布范围,基于核密度估计绘制两者训练样本在各单波段中的密度分布图(图3)。

图3可看出,建筑物、植被点云训练样本的密度曲线在红、绿、蓝3个波段均有相交,都会出现混淆的情况。统计相交部分混淆点云,按照式(1)计算总错分概率:

图 3 建筑物、植被训练样本灰度值分布图 Figure 3 Grey value distribution of training samples of buildings and vegetation
$P\left( t \right) = \left( {1 - \theta } \right){P_{\rm{B}}}\left( t \right) + \theta {P_{\rm{V}}}\left( t \right)$ (1)

式中,P(t)为总错分概率,t为分割阈值,θ为建筑物点云占总训练样本的比例,PB(t)表示分割阈值为t时植被点错分为建筑物点云的概率,PV(t)表示分割阈值为t时建筑物点错分为植被点云的概率。

最佳波段选择时,取定各波段的t值,以建筑物点和植被点的错分概率之和最小为原则。当选定最佳波段后,分割阈值的确定,以该波段中建筑物点和植被点的错分概率之和最小为原则,即当P(t)最小时对应的t值,其解算表达式为

$P\left( t \right) \geqslant 2\sqrt {\left( {1 - \theta } \right){P_{\rm{B}}}\left( t \right)\theta {P_{\rm{V}}}\left( t \right)} $ (2)

当且仅当

$\left( {1 - \theta } \right){P_{\rm{B}}}\left( t \right) = \theta {P_{\rm{V}}}\left( t \right)$ (3)

时,P(t)取得最小值。

    (2.3) DBSCAN聚类算法

DBSCAN聚类算法用于去除车辆、行人和极少的植被点,提取最终的建筑物点。DBSCAN算法(Ester 等,1996) 进行建筑物点云空间密度聚类,其基本思想:点云空间2维平面内组成簇(或称为“类”)的每个激光脚点, 在给定其一定邻域ε范围内,包含的激光脚点个数应大于给定的阈值Tmin。该算法具有无需指定簇的数目、能有效处理噪声点和发现任意形状空间聚类的簇等优点。下面介绍DBSCAN聚类,首先给出算法中所用几个参数的定义:

定义1  ε邻域:给定一个激光脚点,其半径为ε的区域称为该脚点的ε邻域。

定义2 核心脚点:若一个脚点的ε邻域内,包含的脚点个数大于或等于Tmin,则该脚点为核心脚点。

定义3 直接密度可达:在点云集合D中,若核心脚点pε邻域内有一脚点q,则称脚点q是从脚点p出发直接密度可达的(图4(a))。

定义4 密度可达:在点云集合D中,若存在一个脚点链p1p2 $ \cdots $ pnpiD(1≤in),给定ε邻域和Tmin,脚点pi+1是脚点pi的直接密度可达,则脚点pi是从脚点p1的密度可达(图4(b))。

定义5 密度相连:给定邻域ε和半径阈值Tmin,若存在脚点oD,脚点o到脚点p和脚点q均是密度可达的,则脚点p和脚点q是关于邻域ε和半径阈值Tmin密度相连的(图4(c))。

图 4 直接密度可达、密度可达、密度相连 Figure 4 Directly density-reachable, density-reachable, density-connected

DBSCAN算法就是找到从核心脚点出发的最大密度相连脚点,这些最大集合脚点对象形成“簇”,不在簇中的脚点则为噪声点。算法的具体聚类过程如下:

(1) 在点云集合D中,随机选择未被经处理的激光脚点p,脚点pε邻域包含的脚点个数大于或等于Tmin,即选择的脚点对象为核心脚点。将脚点pε邻域内所有脚点形成一个簇,标记为“C”。

(2) 遍历C中所有未处理的激光脚点q,若脚点q也是给定邻域半径ε和阈值Tmin条件下的核心脚点,则将未被规到任何簇的脚点q加入C。若脚点q不是核心对象,则将脚点q暂时标记为“噪声”。

(3) 重复步骤(2),遍历C中未被处理的脚点,直到没有新的激光脚点加入C

(4) 重复步骤(1)到(3),直到点云集合D中所有激光脚点均被“处理”、没有脚点加入任何簇时,整个过程结束。最终未被加入簇的脚点,即为噪声脚点。

DBSCAN聚类算法,设置的参数较少,仅有邻域半径ε和最小脚点个数阈值Tmin。合适εTmin的选取十分重要,ε取值过大,导致很多非建筑物点归到同一个簇中,过小则导致建筑物点被错误地标记为噪声点;Tmin取值越大,簇的密度越高,建筑物边缘脚点将视作噪声点被去除,过小则导致大量核心脚点,从而将噪声归入簇。

3、数据处理与结果分析     (3.1) 实验数据及其预处理

本文研究区位于海地共和国首都太子港,2010年1月12日发生里氏7.0级地震,震中位于首都太子港以西约25 km处,震源深度约10 km(由USGS提供),地震造成大量房屋建筑损毁和人员伤亡。实验数据受世界银行全球减灾与灾后恢复机构资助GFDDR (Global Facility for Disaster Deduction and Recovery),罗切斯特理工学院和ImageCat公司采集完成,获取时间为2010年1月21日,主要包括机载LiDAR点云,平均点云密度约为2 points/m2;高分辨航空影像(仅3个波段),空间分辨率为0.13 m。本文利用总统府及附近数据进行实验测试,面积为260×307 m2,点云个数为163215个,该区域场景复杂,包含建筑物、高大植被、道路、草坪、车辆以及行人等地物(图5)。

图 5 实验数据 Figure 5 The data of study site

数据预处理主要包括噪声脚点去除和影像匹配两部分。点云去噪是点云滤波、分类的前提工作,是指剔除由空中飞鸟、激光多路径效应、飞行姿态以及扫描系统误差等引起的高程异常点云,包括极高点和极低点,可通过查看点云高程分布直方图来剔除。点云通过最邻近算法内插生成数字表面模型(DSM),利用DSM进行航空影像的正射校正,从而完成航空影像向点云数据的配准。点云内插,根据一个格网内至少包含一个点云的原则,采样间隔取 $1/\sqrt n $ n为单位面积内激光脚点数目。

    (3.2) 实验结果

基于上述方法和流程对建筑物点云进行提取。在数据预处理的基础上,利用CSF算法进行点云滤波,很好地实现了地面点云与非地面点云的分离(图6(a)图6(b))。

图 6 地震灾区建筑物点云提取 Figure 6 Extraction of building point cloud in earthquake disaster area

在确定最佳波段时,按照式(1)最小为原则,不妨把t的值设为各波段中两密度曲线相交点对应的灰度值,即

${f_{\rm{B}}}\left( t \right) = {f_{\rm{V}}}\left( t \right)$ (4)

式中,t为灰度值,其中fB(t)表示灰度值为t的建筑物点频率值,fV(t)表示灰度值为t的植被点频率值。计算各波段建筑物点和植被点的总错分概率分别为:0.92%、2.98%、1.96%,可见基于红波段灰度阈值分割的总错分概率最小、不易混淆,因此选择红波段灰度值作为植被点分离的光谱特征信息参数。确定红波段为最佳波段后,通过统计训练样本的红波段灰度值,根据式(3)计算,当分割阈值为209时,总错分概率达到最小。利用该阈值,基于红波段灰度阈值分割的方法,得到图6(c)所示结果,从图6(c)中可见非地面点云中的植被点被有效地剔除,剩下建筑物点、离散分布在道路上的车辆行人点以及极少的植被点(图6(d))。最后,采用基于密度的DBSCAN算法进行聚类,其主要是对半径ε邻域参数和最小点云个数阈值Tmin的设置。根据实验区建筑物面积和点云密度情况,对参数ε在3—8 m、Tmin在20—50进行试验。当εTmin分别取值为4和35时,漏检概率与虚警概率之和最小(图7)。经过聚类后,高密度的建筑物点几乎不变,离散的低密度噪声点(如道路上的车辆、行人以及未分割去除的植被点)被有效去除(图6(e))。

图 7 参数取值 Figure 7 Parameter setting
    (3.3) 精度分析

为定量评价提取效果,验证该方法对震后建筑物点云提取的正确率,通过与手工提取的建筑物点云结果进行比对。建筑物点云提取误差来源包括正确建筑物点未被准确提取和其他类别点云划分为建筑物点云,二者分别用漏检概率和虚警概率衡量。建筑物点云提取误差统计见表1

表 1 建筑物点云提取误差 Table 1 Extraction error of building points

表1表明,漏检概率和虚警概率分别为15.61%、7.52%。漏检主要原因有二:一是建筑物毁坏严重,在进行点云滤波时,把倒塌建筑物点当做地面点滤除;二是阴影区的建筑物点,在基于红波段阈值分割时错分为植被点。虚警点产生主要是局部地区车辆等地物点过于密集,基于DBSCAN算法聚类难以去除。漏检概率和虚警概率均较小,建筑物点云提取总精度达到84.39%,满足一定精度范围的要求。因此,该方法提取震后灾区建筑物点云是比较准确有效的,具有一定的实用性,为震后建筑物点云提取建立了一种可靠的方法。

4、结 论

地震灾区建筑物震害识别是建筑物损失评估的基础,利用机载LiDAR点云数据识别建筑物震害的前提是提取震后建筑物点云。建筑物点云提取难点在于其与植被点云分离,而地震灾区场景更加复杂,建筑物点云,特别是遭受毁坏的建筑物点云与植被点云在法向量、回波次数等特征表现更为相似。许多关于建筑物点云提取的研究多针对规则、完好的建筑物开展;此外,基于机载LiDAR点云识别建筑物震害的研究多仅用面向对象的方法,缺乏基于建筑物点云自身特征识别震害的研究。

因此,本文提出了融合航空遥感数据光谱信息提取震后机载LiDAR建筑物点云的方法。该方法实现简单、流程清晰,为地震灾区建筑物点云提取提供了一种思路。其充分利用了遥感数据资源,通过融合航空影像红波段信息去除难于与建筑物点分离的植被点,基于密度空间聚类的思想剔除低密度噪声点,最终实现灾区建筑物点云的提取,总体精度可达到84.39%。这表明在一定的精度要求范围内,该方法提取地震灾区建筑物点云是可行的。

本文航空影像仅有3个波段,提取建筑物点云的误差来源之一是基于红波段阈值分割时,建筑物点与植被点易混淆。随着激光雷达技术的发展,机载激光雷达遥感可获取多光谱甚至高光谱数据,从而探讨利用NDVI或其他波段组合进行阈值分割,提高分辨植被点云与建筑物点云的能力。此外,DBSCAN算法聚类参数设置有赖于人工试验确定,工作量较大,有待于改进DBSCAN算法使参数自适应确定。建筑物点云提取是建筑物震害识别的前提工作,下一步的研究是根据提取的建筑物点云自身特征进行建筑物震害的识别。

参考文献
[1] Awrangjeb M, Zhang C S and Fraser C S. Building detection in complex scenes thorough effective separation of buildings from trees[J]. Photogrammetric Engineering and Remote Sensing, 2012, 78 (7) : 729 –745. DOI: 10.14358/PERS.78.7.729
[2] Axelsson P. DEM generation from laser scanner data using adaptive tin models[J]. International Archives of Photogrammetry and Remote Sensing, 2000, 33 : 110 –117.
[3] 陈光洲, 岳东杰, 许颖, 郭珍珍. 基于DSM深度影像的机载LiDAR建筑物边界提取[J]. 勘察科学技术, 2016 (1) : 44 –46, 56. Chen G Z, Yue D J, Xu Y and Guo Z Z. Building boundary extraction by airborne LiDAR based on DSM deep image[J]. Site Investigation Science and Technology, 2016 (1) : 44 –46, 56. DOI: 10.3969/j.issn.1001-3946.2016.01.012
[4] Cheng L, Zhao W, Han P, Zhang W, Shan J, Liu Y X and Li M C. Building region derivation from LiDAR data using a reversed iterative mathematic morphological algorithm[J]. Optics Communications, 2013, 286 : 244 –250. DOI: 10.1016/j.optcom.2012.08.028
[5] 窦爱霞, 马宗晋, 黄文丽, 王晓青, 袁小祥. 基于机载LiDAR和多光谱图像的建筑物震害自动识别方法[J]. 遥感信息, 2013, 28 (4) : 103 –109. Dou A X, Ma Z J, Huang W L, Wang X Q and Yuan X X. Automatic identification approach of building damages caused by earthquake based on airborne LiDAR and multispectral imagery[J]. Remote Sensing Information, 2013, 28 (4) : 103 –109. DOI: 10.3969/j.issn.1000-3177.2013.04.018
[6] Elberink S O, Shoko M A, Fathi S A and Rutzinger M. 2011. Detection of collapsed buildings by classifying segmented airborne laser scanner data. ISPRS International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, XXXVIII-5/W12: 307–312 [DOI: 10.5194/isprsarchives-XXXVIII-5-W12-307-2011]
[7] Ester M, Kriegel H P, Sander J and Xu X W. 1996. A density-based algorithm for discovering clusters in large spatial databases with noise//Proceedings of the 2nd International Conference on Knowledge Discovery and Data Mining. Portland, Oregon: AAAI, 226–231
[8] Forlani G, Nardinocchi C, Scaioni M and Zingaretti P. Complete classification of raw LIDAR data and 3D reconstruction of buildings[J]. Pattern Analysis and Applications, 2006, 8 (4) : 357 –374. DOI: 10.1007/s10044-005-0018-2
[9] Hartfield K A, Landau K I and Leeuwen W J D V. Fusion of high resolution aerial multispectral and LiDAR data: land cover in the context of urban mosquito habitat[J]. Remote Sensing, 2011, 3 (12) : 2364 –2383. DOI: 10.3390/rs3112364
[10] 赖旭东, 万幼川. 机载激光雷达距离图像的边缘检测研究[J]. 激光与红外, 2005, 35 (6) : 444 –446. Lai X D and Wan Y C. The study of the edge detection to the depth-image of light detection and ranging[J]. Laser and Infrared, 2005, 35 (6) : 444 –446. DOI: 10.3969/j.issn.1001-5078.2005.06.019
[11] Lee D H, Lee K M and Lee S U. Fusion of lidar and imagery for reliable building extraction[J]. Photogrammetric Engineering and Remote Sensing, 2008, 74 (2) : 215 –225. DOI: 10.14358/PERS.74.2.215
[12] 李亮, 王成, 李世华, 习晓环. 基于机载LiDAR数据的建筑屋顶点云提取方法[J]. 中国科学院大学学报, 2016, 33 (4) : 537 –541. Li L, Wang C, Li S H and Xi X H. Building roof point extraction based on airborne LiDAR data[J]. Journal of University of Chinese Academy of Sciences, 2016, 33 (4) : 537 –541. DOI: 10.7523/j.issn.2095-6134.2016.04.015
[13] 林文鹏, 李厚增, 黄敬峰, 刘冬燕, 宗玮, 胡小猛. 上海城市植被光谱反射特征分析[J]. 光谱学与光谱分析, 2010, 30 (11) : 3111 –3114. Lin W P, Li H Z, Huang J F, Liu D Y, Zong W and Hu X M. Analysis on urban vegetations reflectance characteristics in Shanghai[J]. Spectroscopy and Spectral Analysis, 2010, 30 (11) : 3111 –3114. DOI: 10.3964/J.ISSN.1000-0593(2010)11-3111-04
[14] Liu X Y. Airborne LiDAR for DEM generation: some critical issues[J]. Progress in Physical Geography, 2008, 32 (1) : 31 –49. DOI: 10.1177/0309133308089496
[15] Lodha S K, Fitzpatrick D M and Helmbold D P. 2007. Aerial Lidar data classification using expectation-maximization//Proceedings of SPIE 6499, Vision Geometry XV. San Jose, CA, United States: SPIE, 6499: 64990L [DOI: 10.1117/12.714713]
[16] 马洪超, 姚春静, 张生德. 机载激光雷达在汶川地震应急响应中的若干关键问题探讨[J]. 遥感学报, 2008, 12 (6) : 925 –932. Ma H C, Yao C J and Zhang S D. Some technical issues of airborne lidar system applied to Wenchuan earthquake relief works[J]. Journal of Remote Sensing, 2008, 12 (6) : 925 –932. DOI: 10.11834/jrs.200806124
[17] 苏伟, 孙中平, 赵冬玲, 孙崇利, 张超, 杨建宇. 多级移动曲面拟合LIDAR数据滤波算法[J]. 遥感学报, 2009, 13 (5) : 827 –839. Su W, Sun Z P, Zhao D L, Sun C L, Zhang C and Yang J Y. Hierarchical moving curved fitting filtering method based on LIDAR data[J]. Journal of Remote Sensing, 2009, 13 (5) : 827 –839. DOI: 10.11834/jrs.20090506
[18] Susaki J. Adaptive slope filtering of airborne LiDAR data in urban areas for digital terrain model (DTM) generation[J]. Remote Sensing, 2012, 4 (12) : 1804 –1819. DOI: 10.3390/rs4061804
[19] Ural S, Hussain E, Kim K, Fu C S and Shan J. Building extraction and rubble mapping for city Port-au-Prince post-2010 earthquake with GeoEye-1 imagery and Lidar data[J]. Photogrammetric Engineering and Remote Sensing, 2011, 77 (10) : 1011 –1023. DOI: 10.14358/PERS.77.10.1011
[20] Vetrivel A, Gerke M, Kerle N and Vosselman G. Identification of damage in buildings based on gaps in 3D point clouds from very high resolution oblique airborne images[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2015, 105 : 61 –78. DOI: 10.1016/j.isprsjprs.2015.03.016
[21] 王晓青, 窦爱霞, 王龙, 袁小祥, 丁香, 张伟. 2013年四川芦山7.0级地震烈度遥感评估[J]. 地球物理学报, 2015, 58 (1) : 163 –171. Wang X Q, Dou A X, Wang L, Yuan X X, Ding X and Zhang W. RS-based assessment of seismic intensity of the 2013 Lushan, Sichuan, China MS7.0 earthquake [J]. Chinese Journal of Geophysics, 2015, 58 (1) : 163 –171. DOI: 10.6038/cjg20150114
[22] 于海洋, 程钢, 张育民, 卢小平. 基于LiDAR和航空影像的地震灾害倒塌建筑物信息提取[J]. 国土资源遥感, 2011, 23 (3) : 77 –81. Yu H Y, Cheng G, Zhang Y M and Lu X P. The detection of earthquake-caused collapsed building information from LiDAR data and aerophotograph[J]. Remote Sensing for Land and Resources, 2011, 23 (3) : 77 –81. DOI: 10.6046/gtzyyg.2011.03.14
[23] 张爱武, 李文宁, 段乙好, 孟宪刚, 王书民, 李含伦. 结合点特征直方图的点云分类方法[J]. 计算机辅助设计与图形学学报, 2016, 28 (5) : 795 –801. Zhang A W, Li W N, Duan Y H, Meng X G, Wang S M and Li H L. Point cloud classification based on point feature histogram[J]. Journal of Computer-Aided Design and Computer Graphics, 2016, 28 (5) : 795 –801. DOI: 10.3969/j.issn.1003-9775.2016.05.012
[24] Zhang K Q, Chen S C, Whitman D, Shyu M L, Yan J H and Zhang C C. A progressive morphological filter for removing nonground measurements from airborne LIDAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41 (4) : 872 –882. DOI: 10.1109/TGRS.2003.810682
[25] Zhang W M, Qi J B, Wan P, Wang H T, Xie D H, Wang X Y and Yan G J. An easy-to-use airborne LiDAR data filtering method based on cloth simulation[J]. Remote Sensing, 2016, 8 (6) : 501 . DOI: 10.3390/rs8060501
[26] 张小红. 2007. 机载激光雷达测量技术理论与方法. 武汉: 武汉大学出版社: 105–113 Zhang X H. 2007. Theories and Methods of Airborne Laser Radar Measurement Technology. Wuhan: Wuhan University PRESS: 105–113
[27] 周淑芳, 李增元, 范文义, 庞勇, 陈尔学, 李晓松. 基于机载激光雷达数据的DEM获取及应用[J]. 遥感技术与应用, 2007, 22 (3) : 356 –360. Zhou S F, Li Z Y, Fan W Y, Pang Y, Chen E X and Li X S. DEM extraction and its application based on airborne lidar data[J]. Remote Sensing Technology and Application, 2007, 22 (3) : 356 –360. DOI: 10.11873/j.issn.1004-0323.2007.3.356