一种改进顶帽变换与LBP高程纹理的城区建筑物提取算法

何曼芸，程英蕾，邱浪波，赵中阳

空军工程大学信息与导航学院，陕西 西安 710077

何曼芸，程英蕾，邱浪波,等.一种改进顶帽变换与LBP高程纹理的城区建筑物提取算法[J].测绘学报，2017,46(9)：1116-1122.

10.11947/j.AGCS.2017.20170158.

HE Manyun,CHENG Yinglei,QIU Langbo,et al.An Algorithm of Building Extraction in Urban Area Based on Improved Top-hat Transformations and LBP Elevation Texture[J]. Acta Geodaetica et Cartographica Sinica,2017,46(9):1116-1122. DOI:10.11947/j.AGCS.2017.20170158.

一种改进顶帽变换与LBP高程纹理的城区建筑物提取算法

何曼芸，程英蕾，邱浪波，赵中阳

空军工程大学信息与导航学院，陕西 西安 710077

利用LiDAR数据的建筑物提取存在植被点与建筑物点难以区分的问题，利用航空影像进行城区建筑物提取则无法有效剔除阴影区域植被。本文融合LiDAR和航空影像两种数据源，提出了改进顶帽变换及局部二进制模式(LBP)高程纹理分析的建筑物提取算法。首先将LiDAR数据进行规则格网化，通过改进顶帽变换提取地面数据点，然后根据航空影像计算归一化差值植被指数(NDVI)值进行植被粗提取，计算LBP高程纹理，精细区分植被点与建筑物点，最后利用形态学操作填充建筑物孔洞，以检测出的建筑物点为种子点进行区域生长，得到完整的建筑物点集合。试验基于ISPRS提供的Vaihingen数据集中复杂多植被城区场景，试验结果表明，本文算法能够有效区分植被与建筑物，实现建筑物准确提取。

LiDAR；数字航空影像；LBP；形态学；建筑物提取

近年来，基于机载传感器数据的城区目标自动检测已经成为研究的热点，机载LiDAR(light detection and ranging)系统不易受天气、时段的影响，在地形测绘、环境监测及武器精确制导的军民事领域均得到了广泛的应用。建筑物是城市区域的重要目标，近几年，基于LiDAR系统的建筑物自动检测更是研究的重点。

LiDAR点云数据含有丰富的高程信息，通过高程信息即可有效地区别地物点与地面点。但是LiDAR点云数据具有以下3个特点：一是离散性、不均匀性使得点云数据的邻域点查找困难，目标提取处理速度慢，常用规则格网[1-2]，三角网[3-4]，K-D树[5]等数据结构构建点云数据的邻接关系，如文献[6]使用原始LiDAR点云数据建立Deluanay三角网，结合区域生长算法进行建筑物提取；二是在目标边缘处信息空缺、边界不明显，文献[7]通过Hough变换规则化LiDAR数据建筑物边界，实现规则的建筑物提取；三是缺乏目标的颜色纹理信息，仅用高程信息难以对目标物点进行精确有效地区分，例如常使用目标物面积大小来区分建筑物与植物等其他目标物，但面积阈值的设定需要依赖先验知识，且不能剔除与建筑物紧挨的植被点，易造成错检，文献[7]和文献[8]分别提出基于规则点云数据的建筑物提取方法，定义了不同的植被点剔除准则，实现了建筑物的精确提取，但是分类准则中包含的多个阈值，需要通过目视判读确定。光学航空影像具有目标明确的形状、色彩、光谱等信息，但在建筑物提取中，建筑物与地面的色彩，甚至纹理较为相近，仅从光学影像中难以区分城区建筑物与地面区域，未充分利用建筑物的高程信息和几何信息[9]，而且光学影像数据易受光照、遮挡等外界环境因素的影响，因此，融合光学航空影像与LiDAR点云数据的目标提取与分类能够实现两类数据的优势互补，成为研究的热点与趋势。

融合航空影像与机载LiDAR点云数据进行建筑物提取主要有两类方法，一类根据点云及光谱特征制定建筑物判定规则，如文献[10]融合航空影像和LiDAR点云实现建筑物点云的精确分类以及建筑轮廓边的提取，分类正确率高；文献[11]利用机载LiDAR数据和高分辨率图像提取复杂城区建筑物，特别针对区域内阴影处植被与冷色屋顶，有效地减少了由阴影和不同屋顶特征所造成的错误识别，显著提高了建筑物提取精度。但是该类方法中判定规则需要设定多个分类阈值，算法自动化程度较低。另一类方法对两种数据分别提取相关特征，采用机器学习方法训练分类器，再进行分类，如文献[12]构建了17个特征构成的点云扩展特征向量，采用支持向量机模型对LiDAR点云进行树木识别，文献[13]则利用随机森林[14]提出了一种自适应特征选择方法，取得了较高的精度，但缺点是需要大量训练样本及较长的训练时间。

本文融合光学航空影像与LiDAR点云数据，充分利用两者优势，提出改进顶帽变换及LBP高程纹理分析的建筑物提取算法，减少分类提取过程中的阈值设定，自适应能力较强，实现建筑物的快速精确提取。

1 算法原理

本文采用规则格网的数据结构对LiDAR点云数据进行组织，得到LiDAR数据的距离图像，提出改进的顶帽变换进行地物点提取，利用光学航空影像归一化差值植被指数(normalized difference vegetation index,NDVI)[5]提取植被区域，提出LBP高程纹理分析的建筑物提取算法，并结合形态学操作与基于NDVI指数的区域生长算法，实现建筑物的精确提取。具体算法流程如图1所示。

图1 建筑物提取算法流程Fig.1 Workflow of building extraction algorithm

1.1 地物提取

本文采用规则格网索引来有效的组织LiDAR点云数据。规则格网根据点云数据在水平X，竖直Y方向上的最大值和最小值确定格网的位置，设置网格的大小，通过式(1)，得到点云数据与格网的对应关系

(1)

式中，Grow、Gcol为规则格网的总行、列数；xmax、xmin为区域内点云横坐标的最大值与最小值；ymax、ymin为区域内点云纵坐标的最大值与最小值；i、j为其中一点p的行、列数；点p的横纵坐标为(xp,yp)；Δx、Δy为格网的大小，可根据精度需求与LiDAR点云的平均点密度来确定，本文设定规则网格的网格大小为1×1 m，将LiDAR点云数据转化为格网矩阵，每个格网存储内插得到的高程值。采用最邻近插值[15]算法，减缓了由插值操作导致数据准确度下降的程度，最大限度保留了突起物在边缘处高程的突变特性。规则格网是将LiDAR数据处理转化为图像处理的重要桥梁，打破了LiDAR数据的离散性，使得对LiDAR数据点邻域点的查找变得快速便捷，可以运用成熟的数字图像处理方法实现更加有效的目标提取与分类。

在具有起伏地形的场景内，通过设定一个高程阈值不能有效区分LiDAR数据地面点与非地面点。本文采用形态学滤波[15-16]中多尺度滤波[15]的思想，提出改进的顶帽变换地物点提取方法。顶帽变换[17](top-hat transformations)是形态学的重要应用之一，在图像处理中解决图像的非均匀光照问题，将顶帽变换应用到LiDAR数据处理滤波与目标提取中，LiDAR数据点的高程即对应图像的灰度，起伏的地形即为高程不均匀的背景，因此能够解决地形起伏区域目标提取问题。顶帽变换定义为高程H与高程H的开运算之差，公式表示为

J=H-(H∘s)

(2)

式中，s为开运算的结构元素，结构元素的结构选取常用的圆形结构，结构元素的尺寸过小无法将面积较大的建筑物提取，结构元素过大则无法有效区分地面与地物点。本文通过统计不同结构元素大小顶帽变换的结果，确定地面点与地物点的分类阈值。具体步骤如下：

(2) 进行N=maxSize-minSize+1次顶帽变换；

(3) 统计每个LiDAR点云数据在N次顶帽变换中被标记为地物点的次数；

(5) 计算标记为i次地物点的点的频率pi=qi/pointnum，pointnum为LiDAR点云数据的总个数；

该算法可以降低在顶帽变换中对于高程阈值和结构元素大小的依赖，提升了算法的自适应性。

1.2 植被提取

归一化差值植被指数(NDVI)是遥感影像中常用的植被区域分析指标，其原理是植被生物物理机理：太阳光谱在红波段(Red,R)(0.69 nm)被植被叶绿素强烈吸收，在近红外(infrared,IR)处太阳光谱被叶细胞结构强烈反射[5]，其计算如式(3)所示

(3)

式中，NDVI的数值在[-1,1]之间[18]；植被区域IR分量数值大于R分量数值，其NDVI值大于0；建筑物及地面区域IR分量数值小于R分量数值，NDVI值小于0。设定判定植被的NDVI阈值为0，能够有效区分植被与非植被区域。但是NDVI值容易受到植冠背景的影响[5]，当植物处在阴影区域时，如图3(b)中高大的建筑物阴影遮住了部分植被，利用其NDVI值无法进行有效提取，造成分类错误。而且随着建筑物形式的不断丰富及屋顶绿化的推广，蓝绿色的建筑物屋顶也不再少见，其NDVI值与植被极为相似，仅使用NDVI值对植被与建筑物进行区分，会导致错分。因此针对阴影区域的植被以及蓝绿色屋顶建筑物需要进行进一步分析与分类。

1.3 LBP高程纹理

LiDAR数据地物点最大的特征是在局部具有高程变化形成高程起伏，这些高程变化可以理解为高程纹理。高程纹理可定性、定量地定义为局部区域的高程变化以及由此产生的对比度、均匀性等物理特性[5]，一般有高程差、地形坡度、高程变化、高程差粗糙度[19]等。

建筑物点与植被点的区别在于其与周围邻近点高程的变化规律不同。植物的树冠结构稀疏，反射激光点的位置、强度及次数难以确定，得到的LiDAR数据点与邻域点之间高程变化剧烈，可视作高程纹理较“粗糙”；建筑物邻域点之间高程变化较规则，比如平顶建筑物LiDAR点云数据的高程在理论上是相同的，即使由于探测时环境及设备造成的误差，屋顶LiDAR点云数据高程值在一小范围内波动，斜坡屋顶LiDAR点云数据高程变化率相近，构成倾斜平面，这些情况均可视作高程纹理较“光滑”。基于这一原理，文献[17]在进行城区LiDAR点云数据的树木提取时，采用了梯度分割方法，通过试验确定梯度阈值，实现树木与建筑物的区分。

局部二进制模式(local binary patterns,LBP)[12]通过逐行扫描图像，对图像的每一个像素点进行操作，以该点的灰度值作为判断阈值，对该点的3×3的8邻域进行二值化，按照规则将二值化结果形成一个8位的二进制数，再将二进制数转化为十进制数(0～255)，获得该点的LBP值。对于数字图像来说，LBP值包含图像局部的亮度结构，是用来描述图像纹理的有效方法。将LBP应用到LiDAR数据处理中，对图像灰度值的处理转化为对LiDAR点云数据高程的处理，是对局部高程变化的很好描述。LBP高程纹理分析具体步骤是：

(1) 设中心点为m，中心点8邻域内某一待判断点为n，hm为中心点高程，hn为待判断点高程；

(2) 对中心点m的3×3的8邻域内点的高程逐一进行判断，待判断点高程与中心点高程差值在阈值Th范围内，即标记为1，否则标记为0，如式(4)所示

(4)

(4) 将8位二进制码L转换为十进制数，成为中心点的最终LBP值，如公式(5)所示

(5)

在建筑物点处，8邻域点与中心点间的高程差值均小于阈值，被标记为1，其LBP值较大，而在建筑物边缘点或者植被点处，8邻域点与中心点间的高程差值仅有部分小于阈值，导致此时LBP值较小。为保证将植被点完全地剔除，特别是利用NDVI值检测不到的阴影区植被点，将8邻域内全部被标记为1的点判定为建筑物点，此时LBP值为255，即LBP高程纹理值等于255的点为建筑物点，LBP高程纹理值小于255的点为非建筑物点。由于建筑物边缘点与其邻域存在较大高程差值，此时会被检测为非建筑物点，后期会对提取出的建筑物点进行区域生长操作，进一步提取建筑物边缘点。本文提出的LBP高程纹理建筑物提取方法从理论上分析得到分类规则，无需反复试验确定分类阈值。

1.4 建筑物精提取

建筑物精确提取包括两个方面，一是蓝绿色屋顶建筑物提取，二是建筑物屋顶孔洞与边缘点提取。

2 试验仿真

本文试验是基于MATLAB2014a平台，采用ISPRS提供的德国Vaihingen城市测试数据集中的area1、area2和area3区域，文献[21]提供了该数据集的相关说明。area1区域是城市中心区域，包含许多屋顶结构复杂的建筑物以及植被，area2是高程较大的民用建筑物区域，建筑物周围有许多植被，area3是居民建筑区域，包含许多离散的小型建筑物。3个场景包含丰富的地理环境、城区环境及建筑物类型，能够充分地验证本文算法。Vaihingen数据集中提供的area1、area2和area3区域数字航空影像为彩色红外图像，如图2所示，其光谱组成为IR-R-G，可直接用来计算图像的NDVI值，其对应的LiDAR点云数据平均点距为4 points/m2。

图2 Vaihingen试验区3个场景彩色红外图像Fig.2 The pan-sharpened color infrared images of three test sites in Vaihingen

3个场景LiDAR点云数据通过规则格网化得到数字高程模型，与光学航空影像数据配准融合，进行建筑物提取实验。以area2区域的试验为例，对LiDAR点云数据进行规则格网化，并根据高程对点进行着色，结果如图3(a)所示，规则格网化的插值操作对建筑物边缘的高程变化影响不大，能够保留边缘点的突变特性；计算彩色红外图像的NDVI值，区分植被与非植被区域，结果如图3(b)所示，植被区域在图中表现为灰度较亮的区域，其余地面及建筑物为黑色；对规则格网化的LiDAR数据进行改进的顶帽变换，结果如图3(c)，图中黑色背景部分为起伏的地形，被完全去除，仅剩区域内的建筑物、植被等目标物；在提取出目标物的基础上，运用LBP高程纹理分析，对目标物进行进一步提取与分类，得到去除阴影区植被的不完整建筑物点，如图3(e)所示；最后进行形态学填充孔洞以及区域生长算法，提取完整的建筑物，结果如图3(f)所示。

3 试验结果分析

试验结果定量分析采用文献[22]提出的评价指标，非建筑物点检测为建筑物点定义为一类误差，即误检率，对应算法提取的正确率(Correctness)，建筑物点检测为非建筑物点定义为二类误差，即漏检率，对应算法提取的完整率(Completeness)，文献[23]中定义了提取算法的提取质量(Quality)的计算方法，如式(6)所示

(6)

图3 area2区域建筑物提取结果Fig.3 The extraction results of area2

对area1、area2和area3场景的建筑物提取结果分别进行目标级与像素级的准确度分析，其准确度计算结果见表1；与Grigillo[7]，Awrangjeb[8]和Niemeyer[23]的建筑物提取结果进行像素级比较，见表2，本文算法的提取质量有一定提升，3个区域像素级平均提取质量为90.07%，与ISPRS公布的融合LiDAR数据与光学影像方法的提取结果[24]中像素级最好提取质量89.7%相比，提升了0.37%。

表1 算法提取准确度分析

表2建筑物提取质量对比

Tab.2Thecomparisonofextractionqualitiesindifferentpapers(%)

建筑物提取结果如图4所示，图中红色的点为漏检的点，蓝色的点为错检的点。由试验结果可以看出，本文算法有效区分了紧挨建筑物的植被点，如图4(a)右上所示；成功剔除了受阴影影响的植被点，得到边界明显的建筑物点，如图4(b)右上和右中所示；通过规则格网化处理，在原始LiDAR数据稀疏的建筑物屋顶处也取得了较好的结果，如图4(c)右上所示。但是过于低矮建筑物(图4(b)右下、图4(c)右中)、面积较小且存在高程差的连续屋顶(如图4(a)右中)以及小平台式屋顶(图4(a)右下、图4(c)右下)处仍出现了漏检现象。原因是建筑物面积过小，在LBP高程纹理分析时均被检为边缘点，或者因高程较低在地物点与地面点分类时产生了错分类。本文结果与Grigillo、Awrangjeb和Niemeyer文中的提取结果相比，具有更小的漏检率。

图4 不同区域建筑物提取结果Fig.4 The extraction results in different areas

4 结 论

本文算法通过对点云规则格网化处理，充分利用航空影像的多光谱信息以及LiDAR点云数据的高程信息，提出改进的顶帽变换与LBP高程纹理分析建筑物提取方法，充分考虑了LiDAR数据点的与其邻域点高程的关联性。通过与文献[24]中融合LiDAR数据与光学影像建筑物提取方法对比分析，本文算法实现了复杂城区场景建筑物的有效提取。减少了分类规则中阈值设定的数量，提升了建筑物的提取精度与自适应性，具有一定的普适性。但是此算法在建筑物边界以及特异型建筑物处，仍存在少量错检和漏检现象，下一步可考虑结合建筑物规则边界[25]特征以及特征自适应选择进行研究解决。

[1] HOVAD J，KOMARKOVA J，SEDLAK P.Slope Based Grid Creation Using Interpolation of LiDAR Data Sets[C]∥Proceedings of the 8th International Joint Conference on Software Technologies.Reykjavík，Iceland: SCITE，2013：227-232.

[2] 屈亚运，程英蕾，邱浪波.改进的八邻域搜索提取建筑物立体特征方法[J].空军工程大学学报(自然科学版)，2015，16(4)：66-69.

QU Yayun，CHENG Yinglei，QIU Langbo.An Improved Method of Extracting 3D Features of Buildings Based on Eight Neighborhood Search Method[J].Journal of Air Force engineering University (Natural Science Edition)，2015，16(4)：66-69.

[3] 庞世燕，刘亚文，左志奇，等.结合区域增长法和TIN边缘分割的建筑物立面几何特征提取[J].武汉大学学报(信息科学版)，2015，40(1)：102-106.

PANG Shiyan，LIU Yawen，ZOU Zhiqi，et al.Combination of Region Growing and TIN Edge Segmentation for Extraction of Geometric Features on Building Facades[J].Geomatics and Information Science of Wuhan University，2015，40(1)：102-106.

[4] 吴军，李伟，彭智勇，等.融合形态学灰度重建与三角网分层加密的LiDAR点云滤波[J].武汉大学学报(信息科学版)，2014，39(11)：1298-1303.

WU Jun，LI Wei，PENG Zhiyong，et al.Integrating Morphological Grayscale Reconstruction and TIN Models for High-Quality Filtering of Airborne LiDAR Points[J].Geomatics and Information Science of Wuhan University，2014，39(11)：1298-1303.

[5] 管海燕.LiDAR与影像结合的地物分类及房屋重建研究[D].武汉：武汉大学，2009.

GUAN Haiyan.Research on Objects Classification and Simple Buildings Reconstruction Based on Digital Images and LiDAR[D].Wuhan：Wuhan University，2009.

[6] 何曼芸，程英蕾，聂玉泽，等.综合不规则三角网和区域生长的建筑物提取方法[J].空军工程大学学报(自然科学版)，2017，18(3)：54-59.

HE Manyun，CHENG Yinglei，NIE Yuze，et al.A Building Extraction Method on Comprehensive Delaunay TIN Models and Region Growing[J].Journal of Air Force engineering University (Natural Science Edition)，2017，18(3)：54-59.

[7] GRIGILLO D，KANJIR U.Urban Object Extraction from Digital Surface Model and Digital Aerial Images[C]∥ISPRS Annals of the Photogrammetry，Remote Sensing and Spatial Information Sciences.Melbourne，Australia：ISPRS，2012(I-3)：215-220.

[8] AWRANGJEB M，FRASER C S.Rule-based Segmentation of LIDAR Point Cloud for Automatic Extraction of Building Roof Planes[C]∥ISPRS Annals of the Photogrammetry，Remote Sensing and Spatial Information Sciences.Antalya，Turkey：ISPRS，2013(II-3/W3)：1-6.

[9] 巩翼龙，闫利.结合机载LiDAR数据与航空可见光影像多层次规则分类建筑物变化检测[J].光谱学与光谱分析，2015，35(5)：1325-1330.

GONG Yilong，YAN Li.Building Change Detection Based on Multi-level Rules Classification with Airborne LiDAR Data and Aerial Images[J].Spectroscopy and Spectral Analysis，2015，35(5)：1325-1330.

[10] 程效军，程小龙，胡敏捷，等.融合航空影像和LiDAR点云的建筑物探测及轮廓提取[J].中国激光，2016，43(5)：514002.

CHENG Xiaojun，CHENG Xiaolong，HU Minjie，et al.Buildings Detection and Contour Extraction by the Fusion of Aerial Images and LiDAR Point Cloud[J].Chinese Journal of Lasers，2016，43(5)：514002.

[11] 王雪，李培军，姜莎莎，等.利用机载LiDAR数据和高分辨率图像提取复杂城区建筑物[J].国土资源遥感，2016，28(2)：106-111.

WANG Xue，LI Peijun，JIANG Shasha，et al.Building Extraction Using Airborne LiDAR Data and Very High Resolution Imagery over a Complex Urban Area[J].Remote Sensing for Land & Resources，2016，28(2)：106-111.

[12] 李海亭，王厚之，李艳红，等.支持向量机的车载雷达点云目标识别[J].测绘科学，2016，41(5)：45-49.

LI Haiting，WANG Houzhi，LI Yanhong，et al.Object Recognition for Vehicle-borne LiDAR Point Clouds Based on SVM[J].Science of Surveying and Mapping，2016，41(5)：45-49.

[13] 张爱武，肖涛，段乙好.一种机载LiDAR点云分类的自适应特征选择方法[J].激光与光电子学进展，2016，53(8)：267-277.

ZHANG Aiwu，XIAO Tao，DUAN Yihao.A Method of Adaptive Feature Selection for Airborne LiDAR Point Cloud Classification[J].Laser & Optoelectronics Progress，2016，53(8)：267-277.

[14] BREIMAN L.Random Forests[J].Machine Learning，2001，45(1)：5-32.

[15] 黄先锋，李卉，王潇，等.机载LiDAR数据滤波方法评述[J].测绘学报，2009，38(5)：466-469.DOI：10.3321/j.issn:1001-1595.2009.05.014.

HUANG Xianfeng，LI Hui，WANG Xiao，et al.Filter Algorithms of Airborne LiDAR Data：Review and Prospects[J].Acta Geodaetica et Cartographica Sinica，2009，38(5)：466-469.DOI：10.3321/j.issn:1001-1595.2009.05.014.

[16] 隋立春，张熠斌，柳艳，等.基于改进的数学形态学算法的LiDAR点云数据滤波[J].测绘学报，2010，39(4)：390-396.

SUI Lichun，ZHANG Yibin，LIU Yan，et al.Filtering of Airborne LiDAR Point Cloud Data Based on the Adaptive Mathematical Morphology[J].Acta Geodaetica et Cartographica Sinica，2010，39(4)：390-396.

[17] GONZALEZ R C，WOODS R E.Digital Image Processing[M].3rd ed.Beijing：Publishing House of Electronics Industry，2010.

[18] LIN Y，LIN C，TSAI M D，et al.Object-based Analysis of LIDAR Geometric Features for Vegetation Detection in Shaded Areas[C]∥The International Archives of the Photogrammetry，Remote Sensing and Spatial Information Sciences.Prague，Czech Republic：ISPRS，2016(XLI-B1)：43-46.

[19] 马东洋，尚大帅，保永强.一种基于高程差粗糙度的建筑物激光脚点提取方法[J].测绘科学技术学报，2013，30(2)：158-161.

MA Dongyang，SHANG Dashuai，BAO Yongqiang.A Building Extraction Method Based on Roughness of Height Difference for Airborne LIDAR Point Cloud Data[J].Journal of Geomatics Science and Technology，2013，30(2)：158-161.

[20] 张齐勇，岑敏仪，周国清，等.城区LiDAR点云数据的树木提取[J].测绘学报，2009，38(4)：330-335.DOI：10.3321/j.issn:1001-1595.2009.04.008.

ZHANG Qiyong，CEN Minyi，ZHOU Guoqing，et al.Extracting Trees from LiDAR Data in Urban Region[J].Acta Geodaetica et Cartographica Sinica，2009，38(4)：330-335.DOI：10.3321/j.issn:1001-1595.2009.04.008.

[21] ROTTENSTEINER F，SOHN G，GERKE M，et al.ISPRS Test Project on Urban Classification，3D Building Reconstruction and Semantic Labeling[EB/OL].(2013-12-12)[2017-5-15].http:∥www2.isprs.org/commissions/comm3/wg4/tests.html.

[22] SITHOLE G，VOSSELMAN G.ISPRS Test on Extracting DEMs from Point Clouds：A Comparison of Existing Automatic Filters[R].Delft：ISPRS，2003.

[23] NIEMEYER J，ROTTENSTEINER F，SOERGEL U.Conditional Random Fields for LIDAR Point Cloud Classification in Complex Urban Areas[C]∥ISPRS Annals of the Photogrammetry，Remote Sensing and Spatial Information Sciences.Melbourne，Australia：ISPRS，2012(I-3)：263-268.

[24] ROTTENSTEINER F，SOHN G，GERKE M，et al.Results of the ISPRS Benchmark on Urban Object Detection and 3D Building Reconstruction[J].ISPRS Journal of Photogrammetry and Remote Sensing，2014(93)：256-271.

[25] ZHAO Zongze，DUAN Yansong，ZHANG Yongjun，et al.Extracting Buildings from and Regularizing Boundaries in Airborne Lidar Data Using Connected Operators[J].International Journal of Remote Sensing，2016，37(4)：889-912.

An Algorithm of Building Extraction in Urban Area Based on Improved Top-hat Transformations and LBP Elevation Texture

HE Manyun,CHENG Yinglei,QIU Langbo,ZHAO Zhongyang

Information and Navigation College，Air Force Engineering University，Xi’an 710077，China

Classification of building and vegetation is difficult solely by LiDAR data and vegetation in shadows can’t be eliminated only by aerial images. The improved top-hat transformations and local binary patterns (LBP) elevation texture analysis for building extraction are proposed based on the fusion of aerial images and LiDAR data. Firstly, LiDAR data is reorganized into grid cell, the algorithm removes ground points through top-hat transform. Then, the vegetation points are extracted by normalized difference vegetation index (NDVI). Thirdly, according to the elevation information of LiDAR points, LBP elevation texture is calculated and achieving precise elimination of vegetation in shadows or surrounding to the buildings. At last, morphological operations are used to fill the holes of building roofs, and region growing for complete building edges. The simulation is based on the complex urban area in Vaihingen benchmark provided by ISPRS, the results show that the algorithm affording higher classification accuracy.

LiDAR；digital aerial images；LBP；morphology；building extraction

The National Natural Science Foundation of Shaanxi Province of China (No. 2015JM6346)

HE Manyun(1993—),female，postgraduate，majors in LiDAR data processing.

CHENG Yinglei

P237

A

1001-1595(2017)09-1116-07

陕西省自然科学基金(2015JM6346)

(责任编辑：张艳玲)

2017-04-06

修回日期： 2017-07-06

何曼芸(1993—)，女，硕士生，研究方向为机载LiDAR点云数据处理。

E-mail： 18515373732@163.com

程英蕾

E-mail： ylcheng718@163.com