北京地区地表温度和植被指数及热岛效应相关性研究.docVIP

北京地区地表温度和植被指数及热岛效应相关性研究 　　[摘要]地表温度（LST）是研究区域地表能量平衡和资源环境变化重要的参数之一，也是区域生态环境过程的主要影响因素，同时也是地球科学中很重要的组成部分。本文选取2009年8月北京的Landsat TM影像对北京地区的地面温度进行反演。并通过城市热场变异指数（HI）对北京地区的热岛效应进行定量化分析，结果表明HI能够明显的反应出城市热岛现象。 　　[关键词]地表温度 植被指数 单窗算法 热岛效应 　　[中图分类号] P412.11 [文献码] B [文章编号] 1000-405X（2014）-6-243-3 　　0引言 　　地表温度[1]（LST，Land Surface Temperature）是研究区域能量平衡和资源环境变化的重要参数之一，也是区域生态环境过程的主要因素。LST作为地球科学非常重要的组成的部分，一方面决定了地表向大气的长波辐射能力，另一方面LST也取决于地表许多参数的状态，如地表湿度、植被覆盖和长势。因此掌握LST的空间分布状态及其时间变化态势，对于准备模拟大气和地表、地表地物之间的能量交换是十分必要的。 　　同时，陆地卫星（Landsat）的TM遥感影像的天顶视角下的像元地面分辨率为120m×120m[1]，这一地面分辨率远远比气象卫星NOAA-AVHRR[2， 3]遥感数据的地面分辨率要高。而且TM6热波段的波长区间为10.40～12.50μm [1]，可以很好的用来分析地球表面的热辐射和温度的区域性差异。因此，对于要求精确的区域温度分析来说，TM数据是很好的选择。 　　本次研究将利用Landsat-5 TM3和TM4波段进行计算NDVI，结合TM6波段数据对北京地区进行反演地表温度，研究北京地区的地表温度空间差异，并通过城市热场变异指数对北京地区的热岛效应进行分析，进一步认识北京地区的地表热力场空间差异[4]，为协调北京的资源生态环境提供参考。 　　1研究区概况 　　本次研究选择的研究区域选择北京市，位于39°50N～41°10N，115°51E～116°53E之间，面积约为15000平方千米。北京位于华北平原西北边缘，全市面积一万六千多平方公里，辖16区2县。北京为暖温带半湿润大陆性季风气候，夏季炎热多雨，冬季寒冷干燥，春、秋短促，年平均气温10～12摄氏度。并且，此区域位于平原地区，适合多种农作物生长，水陆差异明显，云量不大[5]，易于遥感信息的获取与判读。 　　2研究方法 　　2.1遥感数据 　　本次实验所用数据是Landsat-5 TM影像，拍摄地区为北京市，拍摄时间为2009年8月。椭球体采用KRASOVSKY1940，地图投影为高斯投影[4]。图像及其地理位置如图1，其中假彩色图像是由TM4、TM3和TM2波段分别赋予r、g、b波段合成。 　　2.2单窗算法 　　对于北京地区的地表温度反演是通过Qin等人提出的“单窗算法”[6]，所谓单窗，就是仅仅利用某个遥感热通道数据估算地表温度。基本思想是通过TM6热红外波段的DN值计算亮度温度，然后将亮度温度作为自变量构造和LST的关系函数，同时考虑地表热辐射受大气吸收作用后的衰减和大气的向下热辐射被地表反射的部分，其公式形式如下： 　　其中， 是要求计算的地表温度（单位：开尔文）；a和b是常量，a=-60.3263，b=0.43436； 是传感器在卫星高度测得的热辐射强度所对应的温度，即亮度温度[6，7]（单位：开尔文），对于Landsat TM，计算公式如下： 　　在（2）式中，DN6是TM6波段的灰度值，取值范围[0， 255]。 　　Ta是大气平均作用温度（单位：开尔文），由大气剖面气温分布和大气状态决定。而对于8月份的北京地区，采用覃志豪提出的经验公式对中纬度夏季平均大气进行估计，可以得到： 　　其中，T0是地面2米高度处的平均温度（单位：摄氏度）。因为缺少观测数据，所以T0采用中国气象局提供的北京地区1990年～2008年的8月份的平均气温t=24℃。 　　C和D是中间变量，分别按下式计算： 　　其中，ε是地表比辐射率[6]，τ是大气透射率[7]。在本次研究中， 采用覃志豪等提出的TM6波段范围内地表辐射率的估计方法来确定。该方法提出在像元尺度上，地表辐射率ε可以通过以下公式得到： 　　北京地区地势比较平坦，所以dε≈0。植被和裸地的比辐射率可按如下取EV=0.970，ES=0.986[8]。RV和RS分别是植被和裸地的温度比率，可以通过以下经验公式计算： 　　在公式（2-6）～（2-8）中出现的PV为植被覆盖度可用图下经验公式计算： 　　其中，NDVI是归一化植被指数，NDVIV和NDVIS分别是植被和裸地的NDVI数，在本次研究中，分别取0