doi: 10.12052/gdutxb.200090
陈炳杰，陈佩敏，黄泽鹏，吴希文，王华
（广东工业大学 土木与交通工程学院，广东 广州 510006）
反监听
摘要: 利用1995年~2019年间的15景的Landsat-8 OLI/TIRS及Landsat-5 TM影像的热红外波段, 采用大气校正法和单窗算法反演广州市的地表温度。通过对研究结果进行时空分析可知, 自1995年~2019年, 广州市的热岛面积先处于稳定阶段, 然后在2001年~2014年间整体呈快速增长趋势, 于2014年后以缓慢速度减少。结果表明, 热岛效应得到缓解的地区为旧城区, 其变化原因可能与工业区外迁、绿化改善等原因有关。热岛空间分布格局由中心集聚型, 逐渐向南北方向和东方向蔓延, 演变成带状主导、多中心围绕的空间分布特征, 其变化的趋势与城市发展模式基本一致。 关键词: 热岛；温度反演；广州；Landsat
中图分类号: TP79 文献标志码: A 文章编号: 1007–7162(2021)02–0053–07
A Spatiotemporal Analysis of Heat Island Effect in
Guangzhou Based on Landsat
Chen Bing-jie, Chen Pei-min, Huang Ze-peng, Ng Alex Hay-Man, Wang Hua (School of Civil and Transportation Engineering, Guangdong University of Technology, Guangzhou 510006, China) Abstract: 15 Landsat-8 OLI/TIRS and Landsat-5 TM images from 1995 to 2019 were used to retrieve the land surface temperature in Guangzhou by atmospheric correction and mono-window algorithm. According to the spatiotemporal analysis of the results, from 1995 to 2019, the area of Guangzhou's heat island began to be in a stable stage, then it showed a rapid growth trend from 2001 to 2014, and then decreased at a slow rate after 2014.
The results show that the areas where the heat island effect has been alleviated are old urban areas, and the reasons for the change may be related to the relocation of industrial areas and the improvement of greening. The spatial distribution pattern of heat islands was agglomerated from the center, gradually spreading to the north-south direction and the east direction, and then evolved into a strip-shaped dominant, multi-center spatial distribution characteristic. Its changing trend was basically consistent with the urban development mode.
Key words: heat island; temperature retrieval; Guangzhou; Landsat
城市热岛效应是指城市中的温度明显高于周边环境温度的现象[1]，已成为现代城市气候的主要特征之一。影响热岛效应的因素一般有人为热量的排放、下垫面性质以及城市结构等[2]。随着城市化进程的加快，大面积的自然地表被热容量更小、吸热更快的人造表面所替代，导致城市热岛现象的强度和影响范围越来越大，对城市的生态环境和居民生活质量造成了不良的影响。因此，研究热岛问题对城市的科学规划、生态城市的建设有着重要的参考意义。常用的研究方法有地面气象站观测[3]、空间模拟[4]和热红外遥感技术。相对于前两种依靠气象数据的技术，卫星热红外遥感有着能获取大范围的地表温度信息的优势，因此成为分析城市热岛的空间分布的主要手段[5]。
改革开放以来，广州经历了快速的城市扩张[6]，城市热岛问题日益凸显。已有学者利用Landsat卫星数据对广州区域进行热岛研究，但多数研究重点关注于广州热岛效应与土地利用类型之间的关联性[7-8]、单一时相的热岛空间格局[9]，而多时相、长时间序列
第 38 卷 第 2 期广东工业大学学报Vol. 38 No. 2 2021 年 3 月Journal of Guangdong University of Technology March 2021
收稿日期：2020-07-15
基金项目：广东省自然科学基金资助项目(2018A030310538)；广州市科技计划项目(201904010254)；广东工业大学大学生创新训练资助项目(xj202011845274)
作者简介：陈炳杰(1995–)，男，硕士研究生，主要研究方向为遥感影像处理与地表形变监测
通信作者：吴希文(1983–)，男，教授，博士生导师，主要研究方向为InSAR，E-mail：*********************
的遥感分析还比较少见。因此，本文基于1995年~ 2019年的Landsat卫星数据，利用温度反演算法对广州地区24年来的热岛空间分布变化及其原因进行分析，为减缓热岛效应，合理规划城市发展和生态环境建设提供参考。
1 研究区及数据源概况
1.1 研究区概况
广州市是广东省省会，由荔湾、越秀、海珠、天河、白云、黄埔、番禺、花都、南沙、从化、增城11个行政区组成，总面积为7 434.4 km2。该城市具有典型的亚热带季风气候，年平均气温21.8 ℃~22.8 ℃，最高月平均气温发生在7、8月份，为28 ℃~30.3 ℃[10]。改革开放以来，由于经济的快速增长，人口及城市规模迅速扩大。至2019年末，广州市常住人口已达1 530.59
万人，城镇化率达86.46%。本文将研究区划分成6个
图 1 研究区域及行政区位置示意图(广州)
Fig.1 Location map of study area and administrative region (Guangzhou)
1.2 数据来源
本次研究使用的卫星影像为美国地质调查局(United States Geological Survey, USGS)提供的6期Landsat-8 OLI/TIRS 及9期Landsat-5 TM影像，每期均包含行列号为122/044和122/043的两景影像，其热红外波段空间分辨率均重采样至30 m[11]。选择的影像为L1TP级，经过了地形校正，均为云量较少的影像，具体日期信息如表1所示。下载的影像再经过辐射校正[12]、大气校正[13]
巴布尔等预处理操作，消除大气造成的辐射误差。最后通过裁剪和镶嵌操作生成研究区的图像。
表 1 Landsat数据日期
Table 1 Date of Landsat images
卫星获取日期
Landsat-5
1995-12-30 1997-11-01 1999-12-25
2001-03-01 2003-12-04 2004-11-20
朝鲜改革开放
2005-11-23 2007-01-29 2008-12-17 Landsat-8
2009-11-02 2013-11-29 2014-10-15
2015-10-18 2017-10-23 2019-11-14
另外，本次研究使用的气象数据来自美国国家环境预报中心(National Centers for Environmental Prediction, NCEP)提供的全球地面观测数据以及欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)的全球陆地表面观测反馈记录报告，广州国家基本气象站编号为59287，观测的时间间隔为3 h。
2 基本原理与方法
2.1 土地利用分类
本文采用支持向量机(Support vector machines, SVM)的监督分类算法[14]来获取土地利用分类图，土地利用类型包括：耕地、草地、林地、建筑用地、水体、未利用地。采用分层随机抽样的方法，来抽取验证样本进行精度评定。最终所有分类图的总体分类精度均在85%以上。
2.2 地表温度反演
根据不同的数据源，目前地表温度的反演算法主要包括大气校正法、单窗算法[15]、劈窗算法[16]、单通道法[17]等。本研究中Landsat-8地表温度反演采用大气校正法[18]，由于2000年前缺少大气参数，Landsat-5采用广泛应用的单窗算法。
2.2.1 大气校正法
大气校正法又称辐射传导方程法，该方法主要使用大气模型来模拟大气对地表热辐射的影响，然后把这部分大气的影响从卫星传感器观测到的热辐射总量中减去来求解地表温度。该方法表达式为[19]
I sensor
式中，为卫星高度上传感器测得的热辐射强度
54广 东 工 业 大 学 学 报第 38 卷
W ·m −2·sr −1·μm −1εB (T s )T s τI ↑
I ↓
W ·m −2
·
sr −1
·μm −1
τI ↑
I ↓
()；为地表比辐射率；是由Plank 定律推导得到的黑体热辐射强度，其中为地表温度 (K)；为大气透射率，、分别为大气的上行
和下行热辐射强度()。、、可以通
过NASA 的大气校正参数计算器(atmcorr.v )获取。
B (T s )在求得地表比辐射率(见式(7)~式(10))后，便能得到，再通过式(2)求解地表温度
[20]
K 1=774.89(W ·m −2·sr −1·μm −1K 2=1321.08式中，对于Landsat-8波段10，)，
K 。
2.2.2 单窗算法
根据覃志豪提出的单窗算法，地表真实温度反演公式为[21]
福建江夏学院学报a 6=−60.3263
b 6=0.北京汉铭
43436T s T 6T a τ6ε6T a 式中，，，为地表温度(K)，
为卫星高度传感器探测的像元亮度温度(K)；为
大气平均作用温度，为大气透射率，为地表比辐射率。的估算，根据热带平均大气经验公式计算
T 0g /cm 2
τ6式(6)中，为气象站的实测气温数据(K)。当大气水分含量在0.4~3.0 变动区间时，的估计方程如表2所示[21]。
表 2 大气透射率估计方程
Table 2 Atmospheric transmittance estimation equation
大气剖面水分含量w /(g·cm –2)
大气透射率估计方程
高气温
0.4~1.6τ6=0.974290−0.08007w 1.6~3.0τ6=1.031412−0.11536w 低气温
0.4~1.6τ6=0
.982007−0.09611w 1.6~3.0
τ6=1.053710−0.14142w
根据已有的研究表明，大气水分含量同地面露
点温度存在相关性。地面露点温度可由气象站实测数据获得，再根据经验公式估计出大气水分含量[22]
。根据土地分类图可分成自然表面、水域和城镇表面，用以下公式计算地表比辐射率
时彩族
[23]
εwater εnature εbuilt −up P v
式中，、、分别代表水域、自然表面、城镇表面的比辐射率。通过公式(10)求植被覆盖度NDVI v NDVI s NDVI v =0.7NDVI s =0.05P v P v 式中，、分别为植被和裸土的归一化植被指数值，取经验值和。当某个像元的NDVI 大于0.7时，取值为1；某个像元的NDVI 小于0.05时，取值为0。
2.3 热岛等级划分
不同时相影像的天气条件都存在差异，反演的温度值不能直接进行比较，但温度的空间分布是相对稳定的状态。采用划分强度等级的方法，既能将不同时相的结果统一到一个标准，也能更直观描述研究区地表温度差异情况。本文采用均值−标准差的方法将地表温度划分成5个等级，如表3所示。其中，低
温区和较低温区统称冷岛区，较高温区和高温区统称热岛区。
表 3 地表温度划分标准1)
Table 3 Classification standard of surface temperature
温度等级取值范围
低温区T s <µ−2σ较低温区µ−2σ<T s <µ−σ无热岛区
µ−σ<T s <µ+σ较高温区µ+σ<T s <µ+2σ
高温区
T s >µ+2σ
µT s σ1) 为研究区地表温度的平均值，为标准差
3 结果
3.1 温度反演结果验证
本文所用Landsat 影像获取时间均在GMT(Green-wich Mean Time)1:50 ~3:00之间，而气象数据观测时间为GMT3：00。由于时间相近，可近似将气象站所在位置的地表温度反演值与气象站观测值比较，用以检验温度反演结果，如图2所示。结果表明，温度反演值均比气象观测值大，这是由于广州在卫星过境时间处于地表升温阶段，且气温值是在距地1.5 m 高处观测所得，地表温度比气温升温更快，温度更高。检验差值最大有3.8 ℃，大部分差值小于2 ℃，均方根误差为1.88，地气温差属于正常范围，温度反演结果可靠。
3.2 热岛动态变化分析
图3显示了1995年~2019年间共15期温度反演图
第 2 期陈炳杰，等：基于Landsat 的广州市热岛效应时空变化分析55
的热岛和冷岛面积的变化，结果表明，广州热岛区域的面积在1995年~2001年间变化不大，呈相对稳定的状态；2001年~2014年进入快速增长阶段，整体呈增长趋势，由2001年的722 km 2扩大至2014年的1
176 km 2，年均增长率为3.82%；自2014年起，热岛面积呈缓慢减少的趋势。而冷岛区域面积在2008年前处于波动状态，2008年后整体呈增长趋势，由2008年的733 km 2增长到2019年的1
139 km 2。
图4
显示了其中6年的地表温度等级分布图，可
看出多年以来温度区的空间分布变化。1995年，城市主要热岛区域集中在珠江沿岸的港口区域以及海珠区西部、新塘镇；北部的从化、花都的热岛出现在裸地和因季节性变化而裸露的耕地(见图4(a))。1997年，沿珠江形成了荔湾−海珠−天河−黄埔的强热岛带。除此以外，各个市辖区的中心镇也开始形成较高温区为主的热岛中心：花都的新华街道、从化的街口街道、增城的荔城街道、番禺的市桥街道(见图4(b))。1997年到2004年间，主城区中心区域热岛效应加剧，
图 2 温度反演结果检验
Fig.2 Verification of temperature inversion results
图 3 广州市热岛、冷岛面积变化图
Fig.3 Area variation of heat island and cold island in Guangzhou
(a) 1995 年
040 km 040 km 040 km 040 km 040 km 040 km (b) 1997 年
(c) 2004 年
N
(d) 2009 年
(e) 2014 年
(f) 2019 年
图 4 广州市不同年份的温度等级分布图
Fig.4 Temperature class distribution of Guangzhou in different years
56
广 东 工 业 大 学 学 报
第 38 卷
南边番禺、南沙热岛面积增加，北边的花都新华镇和白云机场形成两大热岛中心(见图4(c))。2004年到2014年间，热岛区以沿珠江强热岛带为基础，往北沿主要交通道路呈辐射状蔓延至花都，往东向新塘镇、永和经济技术开发区扩张，往南向番禺和南沙扩散(见图4(d)和图(e))。2014年~2019年间，热岛区域无明显的大范围扩张(见图4(f))。
km 2km 2图5为广州市各区热岛占比统计图。可以看出，番禺、花都、主城区等区热岛面积都先呈上升的趋势，并于2014年后有所下降；从化和增城由于其地形和高植被覆盖率的生态环境，热岛占比维
持在较低的水平；南沙则是持续增加的趋势。各区之间比较，主城区热岛面积最大，热岛占比最高可达37.9%，是热岛效应最严重的区域；番禺区仅次于主城区，在2014年可达31.2%；1997年~2014年间，主城区热岛面积增长速度最快，为12.7 /a ，其次是花都区的10.3 /a 。
图 5 广州市分区热岛占比统计图
Fig.5 Statistical chart of heat island proportion in sub-district of
Guangzhou
3.3 土地利用类型与温度区分布的关系通过计算不同温度水平下不同土地利用类型的面积比，可分析出不同土地利用类型对广州热岛效应的影响。表4计算了所有年份的土地利用类型的3种温度区的平均值，可以看出，热岛区占比最高的土地类型是建筑用地，有79.2%，其次是未利用地的8.29%。随着温度水平的降低，建筑用地所占的比例有所下降，林地的比例则开始增加。冷岛区占比最高是林地，高达9
4.18%，在比例上有着绝对的优势，其次是水域的4.69%。
结合谷歌地球和OpenStreetMap 地图进行叠加分析，可以发现热岛区域在建筑用地中的分布特征：容易依附于主要的交通网络分布，如许广高速和广园快速路；较高温区广泛分布于人密集的住宅区
和商业圈，如员村和石牌等地；高温区基本以斑块状分布在工业区、产业聚集地和经济开发区，代表性区域为广州经济开发区、狮岭镇等。而冷岛区域则大面积聚集在森林茂密的广州北部，特别是拥有大量水库、湖泊和森林的从化。
4 讨论
从20世纪90年代初开始，广州进入了快速的城市化建设时期，围绕着中心城区进行向心聚集型的城市发展。2000年后，《广州市城市总体规划(2001~2010)》中提出城市空间发展的基本策略为：南拓、北优、东进、西联，确定南部、东部为城市发展的主要方向。城市发展模式转变为多中心式、组团式，并将旧城区的传统产业向黄埔至新塘一线集中迁移，形成密集的产业发展带。从图6显示的广州城市的扩张可知，除了中心城区的建筑密度增大外，南北方向和
东方向的建筑密度均有明显的增加；外围离散的建筑均开始往中心城镇聚集，并沿交通网络向外延伸。这样的城市扩张的模式与政策建设的方向是一致的。与之对应的，图4显示的热岛空间分布格局呈现着相似的扩散模式。1995~2001年间，广州虽然处于快速的城市化建设进程中，但除了珠江沿岸的中心城区和番禺，其他地区的人口、建筑密度以及工业发展规模依然相对较小，因而还没有大范围的热岛面积增长，在图3中表现为变化不大的热岛曲线。
2001年以后，结合图5可以看出，花都、番禺热岛扩散速度最快，主城区也随着白云区的发展，热岛
面积也有很大的增长。沿珠江强热岛带工业区形成的高温区有所减少，永和镇和新塘镇的热岛效应愈发强烈，各区中心城镇也因人口和产业的聚集形成局部热岛并通过交通网络向外扩散。
为分析热岛效应减弱原因，将1995年的城市建筑用地(旧城区)与2010年后新建城区的热岛面积做比较，如表5所示。可以看出2014年后，在建筑用地面积持续增长的情况下，新建城区的热岛面积仅有微
表 4 广州市不同土地利用类型的温度区比例
Table 4 The proportion of temperature zone in different land
use types in Guangzhou
%
土地类型热岛区无热岛区冷岛区耕地 5.6914.380.16草地0.870.800.05林地 4.3152.3594.18建筑用地79.2016.150.58水域 1.6412.77 4.69未利用地
8.29
3.55
0.34
第 2 期陈炳杰，等：基于Landsat 的广州市热岛效应时空变化分析
57
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议