不透水层作为一种典型的地表覆盖成分,是衡量城市化水平的标准之一,同时也是衡量城市的生态环境的重要指标。本研究以南京市为研究区,以2013~2018年Landsat8影像为数据源,提取不透水层并反演地表温度,分析了2013~2018年南京市不透水层的时空变化以及不透水层与地表温度之间的相关性。研究结果表明:1) 南京市11个行政区县中,鼓楼区不透水层所占的面积最大,六合区所占的面积最小。2) 南京市的不透水面积随时间的推移而增大,越接近市中心的区县,不透水层覆盖面积越大。3) 不透水层与地表温度的相关性较高。 Impervious layer, as a typical surface covering component, is one of the standards to measure the level of urbanization and an important indicator to measure the ecological environment of a city. In this study, Nanjing was taken as the study area, and Landsat8 images from 2013 to 2018 were taken as the data source to extract impermeable layer and invert surface temperature. The spatial and temporal changes of impervious layer and the correlation between impervious layer and surface temperature in Nanjing from 2013 to 2018 were analyzed. The results show that: 1) Among the 11 administrative districts and counties in Nanjing, Gulou district occupies the largest area of impervious layer, while Liuhe district occupies the smallest area. 2) The impervious layer of Nanjing increases with the passage of time. The closer to the downtown area, the greater the impervious layer coverage area. 3) Impervious layer has a high correlation with surface temperature.
不透水层作为一种典型的地表覆盖成分,是衡量城市化水平的标准之一,同时也是衡量城市的生态环境的重要指标。本研究以南京市为研究区,以2013~2018年Landsat8影像为数据源,提取不透水层并反演地表温度,分析了2013~2018年南京市不透水层的时空变化以及不透水层与地表温度之间的相关性。研究结果表明:1) 南京市11个行政区县中,鼓楼区不透水层所占的面积最大,六合区所占的面积最小。2) 南京市的不透水面积随时间的推移而增大,越接近市中心的区县,不透水层覆盖面积越大。3) 不透水层与地表温度的相关性较高。
不透水层,NDISI指数,地表温度反演,时空变化分析
Xunyun Li1,2
1Faculty of Geography, Yunnan Normal University, Kunming Yunnan
2GIS Technology Research Centre of Resource and Environment in Western China of Ministry of Education, Yunnan Normal University, Kunming Yunnan
Received: Mar. 1st, 2024; accepted: Apr. 19th, 2024; published: Apr. 30th, 2024
Impervious layer, as a typical surface covering component, is one of the standards to measure the level of urbanization and an important indicator to measure the ecological environment of a city. In this study, Nanjing was taken as the study area, and Landsat8 images from 2013 to 2018 were taken as the data source to extract impermeable layer and invert surface temperature. The spatial and temporal changes of impervious layer and the correlation between impervious layer and surface temperature in Nanjing from 2013 to 2018 were analyzed. The results show that: 1) Among the 11 administrative districts and counties in Nanjing, Gulou district occupies the largest area of impervious layer, while Liuhe district occupies the smallest area. 2) The impervious layer of Nanjing increases with the passage of time. The closer to the downtown area, the greater the impervious layer coverage area. 3) Impervious layer has a high correlation with surface temperature.
Keywords:Impervious Layer, NDISI Index, Inversion of Surface Temperature, Spatio-Temporal Change Analysis
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着城市的不断发展,我国城市化的水平逐渐提高,城市的范围不断扩大。在城市化的过程中,耕地不断减少、不透水层迅速扩张 [
遥感技术的特性对于提取不透水层具有一定的优势。因此本研究将利用遥感技术的优点,研究不透水层覆盖面积变化,探讨城市不透水层和地表温度的相关性,从而为城市规划及未来发展提供依据和决策支持 [
遥感提取城市不透水层的方法主要有指数方法,影像和阴影信息和线性光谱的方法,还有一些学者考虑到城市不透水层和城市热环境、热岛效应之间的关系 [
指数方法是根据不同地物的光谱特征,提出可以真实反映该地物的特征以及明显区分于其他地物的指数,使用该指数来进行该类地物的提取研究 [
利用影像纹理和阴影信息提取城市不透水面的研究也有很多。黄曦涛等 [
此外,还有一些研究基于线性光谱模型估算城市的不透水层。杨朝斌等 [
还有新兴的基于Landsat8数据的对于城市不透水面和相应的热环境的研究。陈婉佳等 [
南京位于中国东部、长江下游中部地区,地理坐标为北纬31˚14″至32˚37″,东经118˚22″至119˚14″。南京市共包括11个行政区,总面积6587 km2,2017年建成区面积1398.69 km2。是长三角辐射带动中西部地区发展的国家重要门户城市,也是“一带一路”战略与长江经济带战略交汇的节点城市,南京都市圈核心城市,南京经济区主席方城市 [
南京属宁镇扬丘陵地区,以低山缓岗为主,低山占土地总面积的3.5%,丘陵占4.3%,岗地占53%,平原、洼地及河流湖泊占39.2%。南京的平面位置南北长、东西窄,成正南北向;南北直线距离150公里,中部东西宽50至70公里,南北两端东西宽约30公里。南面是低山、岗地、河谷平原、滨湖平原和沿江河地等地形单元构成的地貌综合体。
本文采用的数据是2013~2018年南京市的Landsat8 OLI_TIRS卫星遥感数据,来源于地理空间数据云(http://www.gscloud.cn)。影像的具体信息如表1所示,其中,由于2015年没有符合实验要求的低云量影像,故研究中没有采用。此外,用于影像剪裁的南京市行政边界数据来源于国家基础地理信息中心(欢迎光临国家基础地理信息中心网站(ngcc.cn))。
图像采集时间 | 数据标识 | 条代号 | 行编号 | 云量 |
---|---|---|---|---|
2013-04-07 | LC81200382013097LGN02 | 120 | 38 | 0.02 |
2014-04-24 | LC81200382014114LGN01 | 120 | 38 | 9.79 |
2016-04-29 | LC81200382016120LGN01 | 120 | 38 | 1.68 |
2017-03-15 | LC81200382017074LGN00 | 120 | 38 | 6.56 |
2018-04-19 | LC81200382018109LGN00 | 120 | 38 | 0.31 |
表1. 研究采用的影像信息
本研究使用ENVI5.3对软件对遥感影像进行预处理,具体包括辐射定标、大气校正以及图像剪裁等。在此过程中,选用WGS_84_UTM_Zone_50N为投影坐标系,选用双线性内插法进行重采样,并运用南京市行政边界图对下载的遥感影像进行裁剪。
本研究选取归一化差值不透水面指数(Normalized Difference Impervious Surface Index, NDISI) [
NDISI = R T I R − M N D W I + R M I R + R N I R 3 R T I R + M N D W I + R M I R + R N I R 3 (3.1)
其中,RNIR、RMIR和RTIR分别为影像的近红外、中红外1和热红外波段band10,MNDWI为改进型归一化水体指数RGreen,Green为绿光波段。
MNDWI = R G − R M I R R G + R M I R (3.2)
NDISI具有归一化指数的特征,即指数值介于−1和+1之间;被增强的信息大于0,受抑制的信息普遍小于或等于0。这有利于所增强信息的快速自动提取。
然而,NDISI指数只是增强不透水层信息,为了提取不透水层范围,还需要进行阈值化处理。通过设定阈值,将像元分为不透水层、非不透水层以及水体三类,并通过目视解译进行精度验证。
本研究使用大气校正法进行地表温度的反演。大气校正法反演地表温度的原理是从卫星传感器所观测到的热辐射总量中减去大气对地表热辐射的影响,从而得到地表热辐射强度 [
L λ = [ ε B ( T S ) + ( 1 − ε ) L ↓ ] τ + L ↑ (3.3)
式中,ε为地表比辐射率,B(TS)为黑体热辐射亮度,TS为地表真实温度(K),τ为大气在热红外波段的透过率。则温度为T的黑体在热红外波段的辐射亮度B(TS)为:
B ( T S ) = [ L λ − L ↑ − τ ( 1 − ε ) L ↓ ] / τ ε (3.4)
Ts可以用普朗克公式的函数获取。
T S = K 2 / ln ( K 1 / B ( T S ) + 1 ) (3.5)
对于TIRS Band10, K 1 = 774.89 W / ( m 2 × μ m × s r ) ,K2=1321.08 K。
由上可知此类算法需要2个参数:大气剖面参数和地表比辐射率。大气剖面的参数在NASA提供的网站(http://atmcorr.gsfc.nasa.gov/)中,输入成影时间以及中心经纬度可以获取大气剖面参数。
由2013~2018年Landsat8遥感数据计算得到的南京市NDISI指数分布如图1所示。在此基础上,经反复试验,选取阈值0.467提取各个年份南京市的不透水层,结果如图2所示。
以谷歌高清影像作为参考数据对NDISI阈值化结果进行精度验证。在研究区生成150个随机点,提取每个点的分类结果(不透水层或、非透水层或水体),同时将这些随机点坐标导入对应年份的谷歌高清影像进行目视解译和比对判读,得到对应位置上的真实类型,以此来验证2013~2018年提取的NDISI指数在经过阈值化处理之后数据的精度。验证精度如表2所示,总精度均达到80%以上。
图1. 2013~2018年南京市NDISI指数图
图2. 2013~2018年南京市不透水层分布图
使用者精度 | 生产者精度 | 总精度 | Kappa系数 | |
---|---|---|---|---|
2013 | 88.9% | 87.0% | 88.0% | 0.762 |
2014 | 83.7% | 86.7% | 85.2% | 0.689 |
2016 | 86.7% | 91.5% | 89.1% | 0.779 |
2017 | 90.6% | 95.6% | 93.1% | 0.833 |
2018 | 82.9% | 81.1% | 82.0% | 0.655 |
表2. 精度验证结果
为了更直观详细的观察不透水层在空间上的分布,以2018年为例,结合南京市行政区划矢量图,将2018年不透水层分布图分为南京市的11个行政区,将结果按不同行政区来显示,分别显示南京市区范围内各个区的不透水层面积所占的百分比,如表3所示。
通过表3可以发现,鼓楼区和秦淮区不透水层所占的城市土地面积明显要高于其他行政区。南京市的市中心一共由四个区组成,分别是鼓楼区、秦淮区、玄武区以及建邺区。从表中可以看出,四个组成市中心区县的不透水率均排在11个区县的前五。通过观察和横向比对所获得的2018年不透水层的比例数据,可以发现在南京市的11个行政区县中,鼓楼区不透水层比例为61.16%,所占的面积最大;六合区的不透水层比例为20.71%,所占的面积最小。二者之间相差40.45%,鼓楼区所占的不透水层面积是六合区的2.95倍。由此可以看出,例如鼓楼区这样位于集中进行社会活动发展的区域不透水层覆盖的面积最广泛,而六合区位于江苏北部,被称为是国家东部地区现代工业基地,华东地区先进制造业聚集区和科技创新基地的地区,其主要发展的是农业和工业,总面积中利用类型为耕地的土地占43.12%,工矿用地占17.58%,其中用于商业建设的土地利用面积远远小于位于市中心的鼓楼区。
非不透水层(m2) | 不透水层(m2) | 水体(m2) | 不透水层比例 | |
---|---|---|---|---|
鼓楼区 | 11,205,900 | 32,749,200 | 9,587,700 | 61.16% |
秦淮区 | 14,519,700 | 28,708,200 | 5,283,000 | 59.18% |
雨花台区 | 43,590,600 | 40,721,400 | 7,200,000 | 44.50% |
玄武区 | 38,029,500 | 32,582,700 | 6,455,700 | 42.28% |
建邺区 | 30,859,200 | 35,105,400 | 24,891,300 | 38.64% |
栖霞区 | 178,870,500 | 120,994,200 | 72,656,100 | 32.48% |
浦口区 | 545,869,800 | 258,844,500 | 99,761,400 | 28.62% |
江宁区 | 1,013,364,900 | 446,454,000 | 155,887,200 | 27.63% |
高淳区 | 281,908,800 | 207,333,900 | 310,994,100 | 25.91% |
溧水区 | 626,411,700 | 260,538,300 | 173,397,600 | 24.57% |
六合区 | 1,021,345,200 | 304,758,000 | 145,668,600 | 20.71% |
表3. 南京市各个区的不透水层面积所占的百分比
年份 | 不透水层面积(m2) | 不透水率 |
---|---|---|
2013 | 774,394,200 | 11.75% |
2014 | 866,862,000 | 14.03% |
2016 | 1,307,335,500 | 19.84% |
2017 | 1,464,299,100 | 22.21% |
2018 | 1,769,500,800 | 26.85% |
表4. 2013~2018年南京市不透水层面积和所占比
依据2013~2018的不透水层提取结果,计算各年份城市不透水层的面积和所占比例(表4,图3)。结果表明,随着时间的推移,不透水层的覆盖的面积也随之扩大。2013年南京市的不透水层面积仅占总面积11.75%,而2018年不透水层面积已经达到26.85%。2018年的不透水层面积是2013年的2.29倍。结合城市经济建设发展的状况可以发现,城市的飞速发展和人口迁移促使了城市的扩张,更多的土地被利用来建设成住房、金融中心或者商业建筑,城市不透水层的覆盖面积也不断扩大。
图3. 2013~2018年南京市不透水层比例变化
年份 | 鼓楼区 | 六合区 | ||
---|---|---|---|---|
面积(m2) | 比例 | 面积(m2) | 比例 | |
2013 | 12,231,900 | 22.85% | 140,388,300 | 9.54% |
2014 | 18,741,600 | 35.50% | 157,131,000 | 12.99% |
2016 | 24,363,000 | 49.63% | 191,214,000 | 16.68% |
2017 | 26,569,800 | 53.95% | 297,325,800 | 20.20% |
2018 | 28,885,500 | 61.16% | 304,758,000 | 20.71% |
表5. 不同年份鼓楼区和六合区的不透水层比率
为了直观地了解不透水层在南京市分布的状况以及2013~2018年各行政区内的变化状况,结合对南京市2018年不透水层的研究状况,确定了不透水层覆盖面积最大的行政区是位于市中心的鼓楼区,而覆盖面积最小的是位于南京市最北部的工业基地利六合区。由表5可知,在研究区范围内,鼓楼区的不透水层盖度最高,比例达61.16%,六合区的不透水层盖度最低,比例只有20.71%。
按照时间序列纵向分析,鼓楼区从2013年到2018年不透水层覆盖面积增长了38.31%,其中增长的最快的区间为2013~2014年,增长了12.65%。随着时间的推移,不透水层增长的速度逐渐减缓,由此证明南京市市中心由前期的飞速发展到后期的逐渐平稳建设,是城市商业建筑、公共轨道交通逐渐发展完善的结果;六合区从2013到2018不透水层覆盖面积增长了11.17%,其中增长的最快的区间为2014~2016年,增长了3.69%。六合区的不透水层覆盖面积增长速度相比于鼓楼区缓慢了许多,而且到2017~2018年不透水层覆盖面增长比例仅有0.51%。
由此可以证明,城市不透水层覆盖面积随时间推移而增长,城市的建设发展就是促进不透水层覆盖面积扩大的最主要原因之一。位于南京市市中心的鼓楼区不透水面扩张速度远远大于位于南京市最北部的六合区,说明了城市市中心是不透水层最容易扩张的区域。
为分析不透水层和地表温度的相关性,将NDISI指数的计算结果归一化作为不透水层盖度,并以1为一级进行重分类,总共分为10级。统计每一级所对应的平均地表温度 [
图4. 地表温度反演结果
图5. 不透水层盖度分级
图6. 地表平均温度分布图
图7. 不透水层所占面积与地表温度的相关性分析图
将不透水层盖度和平均地表温度分别作为自变量和因变量进行相关性分析,结果如图7所示。可以看出,R2= 0.9744,证明线性拟合的效果较好,不透水面与地表温度之间存在着相关性,对地表温度的上升起着促进的作用,且在高不透水层覆盖的地区表现得比不透水面盖度的地区更为突出。
不透水层面积分分级 | 温度/℃ |
---|---|
≤10% | 33.01108 |
11%~20% | 34.41853 |
21%~30% | 36.58246 |
31%~40% | 36.69606 |
41%~50% | 39.08632 |
51%~60% | 41.0614 |
61%~70% | 41.26859 |
71%~80% | 45.9986 |
81%~90% | 46.76027 |
>90% | 48.16331 |
表6. 不透水面对应面积及其对应的地表温度
本文通过归一化差值不透水面指数(Normalized Difference Impervious Surface Index, NDISI)提取城市不透水层方法,对南京市2013~2018年的城市不透水层覆盖面积进行时空变化分析及与地表温度的相关性分析,得到以下结论:
1) 对不透水层总体特征分析显示,2018年南京市11个行政区县中,鼓楼区的不透水层覆盖面积最大,所占比例为61.16%,六合区的不透水层覆盖面积最小,所占比例为20.71%,二者相差40.45%。
2) 对不透水层时间变化分析显示,2013~2018年南京市的不透水层随时间的推移而逐渐扩张,2018年的不透水层面积是2013年的2.29倍。由此可见近年来南京市城市化发展迅速,整体呈上升趋势。
3) 对不透水层空间变化分析显示,越靠近城市市中心的区县,不透水层所占面积比例越高;距离城市市中心越远的区县,不透水层所占的面积比例越低。在南京市11个行政区县中,组成市中心的秦淮区、鼓楼区、玄武区、建邺区均排在城市不透水层覆盖面积前五,而距离市中心最远的南京市最北部的工业中心六合区不透水层覆盖的面积最少。
4) 南京市在2013~2018的6年间,不透水层面积飞速增长,增加了近2.29倍。不透水面的增加主要通过大面积侵占植被和耕地来实现。根据相关性分析结果显示,不透水面与地表温度之间存在着相关性,对地表温度的上升起着促进的作用,且在高不透水层覆盖的地区表现得比低不透水面盖度的地区更为突出。
李巽云. 基于Landsat遥感数据城市不透水层信息提取与分析Extraction and Analysis of Urban Impervious Layer Information Based on Landsat Remote Sensing Data[J]. 地理科学研究, 2024, 13(02): 434-445. https://doi.org/10.12677/gser.2024.132042
https://doi.org/10.1080/02757250109532436
https://doi.org/10.1016/j.rse.2006.09.003
https://doi.org/10.1175/2008JAMC1983.1
https://doi.org/10.1016/j.rse.2009.10.008