及时、有效地获取城市中的建成区信息对进行合理的城市规划和管理有十分重要的意义。遥感数据具有来源广、覆盖范围大、实时采集的特点,包含丰富的地面信息,是进行城市信息采集和监测的有效手段。针对现有的利用遥感影像进行建成区提取中存在的裸土和建成区混分问题,根据建成区的光谱特征,本文提出了一种基于Landsat影像的修正的归一化差分建成区指数(Modified Normalized Difference Built-up Index,MNDBI),并针对建成区的指数提取结果提出了相应的后处理框架,主要包括聚类去裸土和形态学偏重建两步。后处理中,第一步可以去除建成区提取结果中混杂的裸土,提高建成区提取的正确率,第二部可以减少提取结果中建成区的漏分,降低漏检率。实验结果显示,MNDBI相对于现有的建成区指数,可以更好的区分裸土和建成区,提高建成区的指数提取精度。指数提取后处理步骤可以有效的降低建成区指数提取的漏检率,减少虚警,进一步提高提取精度。 As remotely sensed data can provide consistent and substantially information of the earth, it has become a common tool to explore the land surface timely and effectively. Specifically, it is widely applied in monitoring the built-up surfaces in urban area, which is significant for urban planning and management. To alleviate the confusion between built-up area and bare soil which is one of the major difficulties in the extraction of built-up area, this article proposes a novel Modified Normalized Difference Built-up Index (MNDBI). Moreover, a two-step post processing scheme, including clustering and morphological partial reconstruction, is also presented follow the built-up extraction. During the two-step post processing, missed built-up pixels can be filled up and the rest bare soil pixels can also be excluded. The proposed MNDBI index demonstrates better results and accuracies comparing to the existing built-up indices in our experiments. In addition, the two-step post processing can further alleviate false alarm and missed alarm and improve the extraction accuracy of built-up area.
胡月瑶
武汉大学测绘遥感信息工程国家重点实验室,湖北 武汉
收稿日期:2017年6月1日;录用日期:2017年6月19日;发布日期:2017年6月21日
及时、有效地获取城市中的建成区信息对进行合理的城市规划和管理有十分重要的意义。遥感数据具有来源广、覆盖范围大、实时采集的特点,包含丰富的地面信息,是进行城市信息采集和监测的有效手段。针对现有的利用遥感影像进行建成区提取中存在的裸土和建成区混分问题,根据建成区的光谱特征,本文提出了一种基于Landsat影像的修正的归一化差分建成区指数(Modified Normalized Difference Built-up Index,MNDBI),并针对建成区的指数提取结果提出了相应的后处理框架,主要包括聚类去裸土和形态学偏重建两步。后处理中,第一步可以去除建成区提取结果中混杂的裸土,提高建成区提取的正确率,第二部可以减少提取结果中建成区的漏分,降低漏检率。实验结果显示,MNDBI相对于现有的建成区指数,可以更好的区分裸土和建成区,提高建成区的指数提取精度。指数提取后处理步骤可以有效的降低建成区指数提取的漏检率,减少虚警,进一步提高提取精度。
关键词 :Landsat,建成区提取,指数,后处理
Copyright © 2017 by author and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
从上世纪70年代改革开放至今将近40年的过程中,中国城市数量和面积都发生了翻天覆地的变化,城市化进程在全国甚至全球都在不断推进,近几年仍然有着愈演愈烈的趋势。及时、准确地获取城市区域空间信息可以对城市的发展进程进行宏观的控制,建成区作为城市区域的一个重要组成,是人类活动的主要载体,提取和监测城市中的建成区,对于进行合理的城市规划、提升城市生活品质有着重要意义。遥感技术是从20世纪60年代发展起来的一门新兴的技术,从1972年美国发射第一颗陆地卫星(Landsat系列卫星)开始快速融入人类社会发展的浪潮中,逐渐成为对地观测的主要方式之一,为人类提供宏观、大范围、持续的地表信息。其中Landsat系列影像具有较为丰富的空间、光谱信息,单景影像覆盖范围较大,易于获得,可以高效的监测城市尺度、区域尺度地表现状,具有实现大尺度城市区域实时监测和更新的潜力 [
利用Landsat影像进行城市中的建成区的提取,目前常见的方法包括监督分类 [
针对建成区和裸土的混分问题,本文使用Landsat作为数据源,主要创新成果包括:
1) 在现有的建成区指数NDBI基础上,提出了一种修正的归一化差分建成区指数(Modified Normalized Difference Built-up Index, MNDBI),可以有效的抑制裸土背景,提高提取精度。
2) 针对指数提取的不确定性,提出两步后处理框架,包括聚类去裸土和形态学偏重建两步,进一步剔除建成区提取结果中可能混淆的裸土像素并降低漏检率。
指数方法提取建成区具有便捷、快速的特点,可以进行快速城市建成区提取,为后续的应用分析提供基础。本文在归一化差分指数构造模式(Normalized Difference X Index, NDXI)下 [
由不同地物类型在Landsat不同波段上的光谱特征可知(图2),相对于其他几类地物类型,建成区在短波红外1波段(Swir1)反射率较强,近红外波段(Nir)反射率较低,在这两个波段上的差值相对其他几类地物类型较大,构造指数时可以通过这两个波段作差,增强建成区的响应。基于建成区在Swir1和Nir两个波段上的光谱特征,利用NDXI模型(公式(1)),查勇 [
图1. 建成区提取流程图
图2. 地物光谱特性曲线(Landsat 8)
公式(1)中,
观察图2中建成区和裸土在Nir、Swir1、Swir2三个波段上的光谱特征,可以看出,相对于森林、耕地和水体,建成区和裸土在短波红外2波段(Swir2)和短波红外1波段都具有较高的反射率,Swir2与Nir的差值和Swir1与Nir差值都可以增强建成区和裸土的响应值,将这两类地物与其他地物类型区分开。而其中,建成区和裸土在Swir2与Nir两个波段上的差值的差异性要明显大于在Swir1与Nir两个波段上的差值的差异性。基于建成区和裸土的光谱特征,本文提出了一种修正的归一化差分建成区指数(MNDBI,公式(3)),在NDBI的基础上,将短波红外1波段替换成了短波红外2波段,更有效的抑制裸土背景对于建成区提取的影响。
裸土与建成区一样具有比较亮的亮度,且裸土没有形态多变,形状不统一,分布规律性较弱,在分类时往往会和建成区混淆,引起错分。实验发现,对于指数提取建成区来说,对提取精度影响最大的因素就是裸土的错分。
Ward [
这一组矿物指数包括:
1) 黏土矿物质指数:
2) 氧化铁指数:
3) 铁类金属指数:
本文将黏土矿物质指数、氧化铁指数、铁类金属指数三种矿物指数影像进行叠加,作为输入数据,利用K-means聚类,识别并剔除混淆在提取结果中的裸土。
聚类去裸土可以降低建成区提取的虚警率,提高建成区提取的正确性。针对指数提取结果中可能存在的建成区漏分现象,如建成区中一些亮度较暗的建筑被漏分,我们对指数提取结果进行第二步的后处理,从保证建成区提取的完整性。
数学形态学已经被广泛应用于图像处理中 [
实验主要包括两个部分:第一部分是对本文提出的修正的归一化差分建成区指数的提取结果进行评价,并与现有的建成区指数方法进行对比;第二部分是比较经过后处理的建成区提取结果和指数提取结果,分析后处理各步骤对建成区提取效果产生的影响。本文对各个实验区域分别绘制了地面参考样本,在参考样本中随机抽取2000个像素点作为测试样本,其中包括:1000个建成区样本点、250个水体样本点、250个裸土样本点、250个草地样本点、250个森林样本点。利用混淆矩阵、总体精度OA、Kappa系数来衡量利用不同方法进行建成区提取的精度。
研究数据包括武汉市2014年Landsat-8 OLI数据、北京市1992年Landsat-5 TM数据、广州市2002年Landsat-7 ETM+数据,本文所使用的Landsat数据均为 L1T等级数据,可以从USGS上免费获取。L1T等级的Landsat数据已经经过了基本的地形数据参与的几何校正,一般情况下可以直接使用而不需要做几何校正。为了利用其丰富的波段光谱信息,获得地表反射率,我们对原始数据进行了辐射定标和FLAASH大气校正(ENVI 5.1)。本文对Landsat影像预处理步骤主要包括:
1) 统一坐标系;
2) 辐射定标;
3) 大气校正;
4) 拼接、裁剪。
统一所有数据的坐标系为WGS-84坐标系,使用UTM投影。辐射定标是将传感器直接获得的地表DN(Digital Number)值转换为辐射亮度值;大气校正则是将辐射亮度转换为地表真实反射率。本论文直接使用ENVI 5.1辐射定标和FLAASH大气校正相关工具对所有Landsat影像进行这两步处理。对于一个城市的多景影像,由于需要得到市级行政区域的Landsat影像,需要将多景影像进行拼接,再利用行政区矢量图(Database of Global Administrative Areas, GADM),对拼接得到的影像进行裁剪(这两步也可直接通过ENVI 5.1 Seamless Mosaic和Subset Data by ROIs两个工具实现),最终得到经过预处理的,武汉市、北京市、广州市的影像数据。
本节对比了归一化差分建成区指数NDBI、基于指数计算的建成区指数IBI,以及本文提出的修正的归一化差分建成区指数MNDBI在三个实验区域(武汉市、北京市、广州市)上的建成区提取结果及精度。
图3给出了各个指数提取建成区的结果图,图4给出广州市局部放大区域的提取细节。其中,IBI的提取结果在非主城区区域存在很多细碎点,经检查可知,这些细碎点大部分为土壤虚警,被错误的标记为建成区,NDBI的提取结果中同样存在一定的虚警。MNDBI提取结果中的虚警明显较少,同时建成区也被比较完整的提取出来。表1中的定量精度评定结果也展现出了类似的规律,总体上三个指数都可以较完整的提取实验区域的建成区,MNDBI在三个实验区域的建成区提取中获得了最高的精度,而IBI
图3. 建成区指数提取结果对比
研究区域 | NDBI | IBI | MNDBI | |||
---|---|---|---|---|---|---|
总体精度% | Kappa系数 | 总体精度% | Kappa系数 | 总体精度% | Kappa系数 | |
武汉2014 | 84.84 | 0.6970 | 85.05 | 0.7010 | 86.75 | 0.7350 |
北京1992 | 85.25 | 0.7050 | 81.95 | 0.6390 | 86.30 | 0.7260 |
广州2002 | 92.55 | 0.8510 | 86.10 | 0.7220 | 94.00 | 0.8800 |
表1. 建成区指数提取精度评定
图4. 广州市局部放大区域提取结果
和NDBI的精度相对略低。尤其是在广州市,MNDBI的提取精度(总体精度:94.00%;Kappa系数:0.8800)显著高于IBI(总体精度:86.10%;Kappa系数:0.7220),略高于NDBI(总体精度:92.55%;Kappa系数:0.8510)。以上的实验结果验证了MNDBI的有效性和优越性。
根据本文提出的建成区指数提取后处理框架,我们需要对MNDWI指数提取出的建成区进行两步后处理:
1) K-means聚类去除裸土;
2) 形态学偏重建减少漏分、去噪。
需要注意的是,为了提高效率,本文对于指数提取错分率小于10/1000的影像省略了后处理步骤1),直接进行后处理步骤2)。在本实验中,广州2002年的影像符合这一条件,只进行了第二步后处理。
分别对武汉、北京、广州三个城市的指数提取结果进行后处理,每一步的城市提取实验结果见图5,精度评定结果见表2。从各个步骤的处理结果看,聚类去裸土后非中心城市区域离散点显著减少,裸土区域被进一步去除,而经过形态学偏重建,建成区提取结果变得更加聚集,内部空洞被适当填补,中心城市区域建成区漏分被缓解。从每一步的精度评定混淆矩阵中也可以看出,经过聚类去裸土,被分到建成区的非建成区样本点显著减少,虚警率显著降低。经过形态学偏重建,被分到非建成区的建成区样本点数量明显减少,建成区提取的漏检率降低,建成区提取的更加正确和完整。随着每一步的后处理,总体精度和Kappa系数也在逐渐升高,三个实验区域的建成区提取精度在后处理之后均在90%以上,Kappa系数均在0.86以上,武汉市的总体精度达到了97.25% (Kappa: 0.9450),总体精度比指数提取结果提高了10.5% (Kappa系数提高了0.21)。在三个实验区域的建成区提取结果说明了指数提取后处理方法的有效性。
本文针对Landsat影像建成区提取问题,首先基于建成区的光谱特征提出了一种修正的归一化差分建成区指数MNDBI。MNDBI在NDBI基础上,将Swir1波段替换为Swir2,从而拉开建成区与裸土的光谱差异,抑制裸土背景,提高建成区提取的正确性。指数目标提取具有计算简便、易于实施的特点,但
实验区域 | 方法 | 类别 | 地表真实样本 | 生产者精度(%) | 使用者精度(%) | ||
---|---|---|---|---|---|---|---|
建成区 | 非建成区 | 总计 | |||||
武汉2014 | 指数 提取 | 建成区 | 980 | 245 | 1225 | 98.00 | 80.00 |
非建成区 | 20 | 755 | 775 | 75.50 | 97.42 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):86.75 | Kappa系数:0.7350 | ||||||
聚类 去裸土 | 建成区 | 821 | 44 | 865 | 82.10 | 94.91 | |
非建成区 | 179 | 956 | 1135 | 95.60 | 84.23 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):88.85 | Kappa系数:0.7770 | ||||||
形态学 偏重建 | 建成区 | 988 | 43 | 1031 | 98.80 | 95.83 | |
非建成区 | 12 | 957 | 969 | 95.70 | 98.76 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):97.25 | Kappa系数:0.9450 | ||||||
北京1992 | 指数 提取 | 建成区 | 782 | 56 | 838 | 78.20 | 93.32 |
非建成区 | 218 | 944 | 1162 | 94.40 | 81.34 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):86.30 | Kappa系数:0.7260 | ||||||
聚类 去裸土 | 建成区 | 739 | 38 | 777 | 73.90 | 95.11 | |
非建成区 | 261 | 962 | 1223 | 96.20 | 78.66 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):85.05 | Kappa系数:0.7010 | ||||||
形态学 偏重建 | 建成区 | 920 | 26 | 946 | 92.00 | 97.25 | |
非建成区 | 80 | 974 | 1054 | 97.40 | 92.41 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):94.70 | Kappa系数:0.8940 | ||||||
广州1992 | 指数 提取 | 建成区 | 824 | 6 | 830 | 82.40 | 99.28 |
非建成区 | 176 | 994 | 1170 | 99.40 | 84.96 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):90.90 | Kappa系数:0.8180 | ||||||
形态学 偏重建 | 建成区 | 870 | 2 | 872 | 87.00 | 99.77 | |
非建成区 | 130 | 998 | 1128 | 99.80 | 88.48 | ||
总计 | 1000 | 1000 | 2000 | ||||
总体精度(%):93.40 | Kappa系数:0.8680 |
表2. 建成区提取后处理精度评定
图5. 建成区提取后处理结果
对于建成区这样具有很高光谱多样性和复杂性的地物类型,指数提取结果往往依然存在很多的不确定性。针对建成区指数提取结果的不确定性,本文提出了一个建成区指数提取后处理框架,包括聚类去裸土和形态学偏重建两步。聚类去裸土可以有效剔除建成区提取结果中可能混淆的裸土像素,形态学偏重建可以填补提取结果中的空洞,降低漏检率,去除噪声点,使提取结果更加连续和完整。在三个实验区的建成区提取中,MNDBI可以获得最高的建成区提取精度,且建成区的提取精度随着后处理逐渐提高,提取结果中的虚警减少,漏检降低。实验结果验证了MNDBI进行建成区提取的优越性,同时验证了本文所提出的后处理框架对于提高建成区提取精度的有效性。
胡月瑶. 修正的归一化差分建成区指数及建成区提取后处理方法研究 Investigation of a Modified Normalized Built-Up Index and a Post Processing Scheme for BUILT-UP Extraction in Urban Area[J]. 测绘科学技术, 2017, 05(03): 83-92. http://dx.doi.org/10.12677/GST.2017.53011