基于TM影像反演郑州市地表温度及城市绿化关系探究.docVIP

基于TM影像反演郑州市地表温度及城市绿化关系探究摘要:本文以Landsat 5 TM为主要遥感数据源,辅助以气象数据资料,通过NDVI计算植被覆盖度来反映郑州市不同下垫面的绿化情况,并利用单窗算法精确反演郑州市2009年6月25日的地表温度。研究表明,郑州市植被覆盖程度与地表温度呈明显的负相关,城市绿化对热岛效应具有明显的削弱作用。 关键词:Landsat 5 TM 植被覆盖度 单窗算法 热岛效应 城市绿化 城市地域的景观格局变化在不同尺度上对生态系统的结构与功能产生影响,其中植被覆盖情况是衡量区域地表温度（Land Surface Temperature,简称LST）的重要指标之一[1]。由于植被光合作用能够减弱太阳辐射热效应,因此植被覆盖能够有效降低地表温度。植被覆盖度（fegetation fraction,简称Fv）与LST间的关系一直是城市热岛（Urban Heat Island,以下简称UHI）研究的热点[1]。TM数据是目前环境研究中应用最多的卫星热红外遥感数据之一[2]。本文采用覃志豪[3]等提出的针对TM6地温反演的单窗算法,同时,利用归一化植被指数（NDVI）得到的Fv,对郑州市UHI的分布与城市绿化关系做初步探索研究,为郑州市环境监测评价、城市规划等提供借鉴。 1、研究区域介绍及数据源 1.1 研究区域概况 郑州市位于东经113°26’-113°51’和北纬34°35’-34°57’之间,属暖温带大陆性气候,四季分明,年平均气温14.4°C;7月最热,平均27.3°C;1月最冷,平均0.2°C;年平均降雨量640.9毫米。平均海拔110.4m,地势是西部高,东部低。 1.2 数据源 本文采用的Lansat5 TM影像从国际科学数据服务平台下载,卫星过境为2009年6月25日,轨道号/行号为124/36, 研究区无云,影像已经过系统几何校正,精度能够满足本文需要。气象数据是卫星过境当天的平均气温、相对湿度、水面温度等,来自中国气象科学数据共享服务网。遥感分析采用ENVI4.7软件。 2、研究要素的提取 2.1 基本原理 本研究采用Qin等[3]针对TM热红外波段反演实际地表温度Ts（K）的单窗算法,其公式为: Ts={a6(1-C6-D6)+[b6(1-C6-D6)+C6+D6]T6-D6Ta}/C6(1) C6=τε (2) D6=(1-τ)[1+τ(1-ε)] (3) T6为TM6的亮度温度(℃);Ta为大气平均作用温度（℃）;τ为大气透过率;ε为地表辐射率;a6、b6为常量,一般情况下（地表温度在0-70℃之间）,a6=-67.355351,b6=0.458606;C6和D6是中间变量,分别通过(2)、(3)式表示。 因此,只要知道参数Ta、τ和ε,即可用上述单窗算法反演任何像元的实际地表温度。其中Ta、τ可由覃志豪等[4]提出的单窗算法的大气参数估计方法求得。 2.2 亮度温度T6计算 对于TM,可用（4）式求算卫星高度的像元亮度温度[3] T6=1260.56/ln[1+60.776/(0.1238+0N)] (4) 式中,DN为TM6的像元值;T6单位是K。 2.3 地表辐射率ε的估计 首先由覃志豪等[5]提出的公式(5)求得植被覆盖度Fv。 Fv=(NDVI-NDVIs)/（NDVIv-NDVIs) (5) 其中NDVIs为裸土的NDVI值,NDVIv为植被的NDVI值,结合实验区具体情况,本研究取NDVIv为0.7,NDVIs为0。 根据覃志豪等[5]提出的估计方法,本研究将遥感影像分为水体、城镇和自然表面3种类型,水体像元比辐射率取0.995,自然表面ε1和城镇ε2像元的比辐射率根据(6)、(7)获得: ε1=0.9625+0.0614Fv-0.0461Fv (6) ε2=0.9589+0.086Fv-0.0671Fv (7) 进而可利用ENVI[6]软件由Fv、水体像元比辐射率、ε1、ε2,获得研究区每个像元的地表辐射率ε。 3、分析 3.1 城市地表温度分布特征分析与结果 地表温度反演结果（当日实测气温:最高41.9℃,最低23.6℃,平均33.3℃）如图1。 图1 郑州市地表温度反演结果图(单位:℃) 结果表明,研究区地表平均温度为44.5℃,最低温为22.3℃（黄河）,最高温度为59.4℃。市中心温度大多在45℃以上（除水体、绿地外）,郊区除西部的工业区和西南部的交通要道外,温度都在45℃以下,而且40℃以下的区域占有相当大的比例,市中心与郊区温度相差显著。 研究区内的高温区域主要有二七塔附近的商业街与车站,铁路沿线附近,西郊工