基于Landsat TM影像的地表温度反演算法差异及区域发展分析.docVIP

基于Landsat TM影像的地表温度反演算法差异及区域发展分析 [摘 要]以徐州市作为研究区，利用基于影像的反演算法（IB），单通道法（SC）和单窗算法（MW）分别对2005年的Landsat TM影像进行地表温度反演，并将结果与研究区亮度温度做比较，与此同时研究了区内热环境及热岛效应，对区域发展进行了分析，结果表明：（1）三种算法总体差异不大，都比亮度温度稍微高出一点，IB算法高出1.414℃，SC算法高出2.186℃ ，MW高出2.3096℃；（2）徐州市表现出较强的城市热岛效应，城区内地表温度普遍高于郊区，今后必须环境与经济和谐发展，在加快区域经济发展的同时，缓和城市热岛效应。 [关键词]Landsat TM； 地表温度反演算法差异； 热岛效应 ；区域发展； 中图分类号：TP75 文献标识码：A 文章编号：1009-914X（2015）12-0110-02 近年来，在温室效应及城市快速发展背景下，徐州城市热岛效应越来越强。地表温度作为环境监测和评估热岛效应的重要参数，已被广泛应用。目前，国内外基于LandsatTM/ETM+的地表温度反演方法主要有以下几种：热辐射传导法，基于影像的反演算法，单窗算法，单通道法，劈窗算法，多通道算法。本文应用基于影像的反演算法，单窗算法，单通道法这三种方法对研究区进行地表温度反演，得到徐州市地表温度时空分布图像，对徐州市热岛效应作出初步探究，为同步提高生态环境与经济发展提供理论基础。同时，将三种算法结果与卫星高度上传感器所探测到的亮度温度对比，作差异分析。 1.研究区概况及数据处理 1.1 简介 徐州市位于东经116°22′～118°40′、北纬33°43′～34°58′之间。东西长约210公里，南北宽约140公里，总面积11258平方公里，占江苏省总面积的11%。属于华北平原的东南部，域内除中部和东部存在少数丘岗外，大部皆为平原。丘陵海拔一般在100-200米左右，丘陵山地面积约占全市9.4%。徐州属暖温带半湿润季风气候，年气温14℃，年日照时数为2284至2495小时，日照率52%至57%，年均无霜期200至220天，年均降水量800至930毫米，雨季降水量占全年的56%。徐州是全国重要的煤炭产地、华东地区的电力基地，拥有多种矿产，由于过度开采，近年来，矿产资源日渐枯竭，生态环境日渐脆弱。 1.2 数据处理 本研究使用的数据是美国陆地资源卫星LandsatTM第六波段数据，数据获取时间是2009年5月3日，影像稍有云层，可能会造成微量误差。利用ENVI4.8对其进行图像裁剪，辐射校正、大气校正、预处理之后，进行地表温度反演。 2.研究方法 2.1 基于影像的反演算法 基于影像的反演算法（Image-basedMethod，以下简称IB算法），将TM热红外波段的灰度值（DN值）转换成辐射值；生成亮度温度，进而反演实际温度。具计算过程： TB 1260.56/ln[1+607.76/（1.2378+0.055158×DNTM6）]. 其中TB是地面亮度温度，DNTM6为TM6的像元灰度值（DN值），0≤DNTM6≤255. 利用地物比辐射率进行校正，得出地面真实温度： 普朗克常数h 6.626×10-34Js；波耳兹曼常数σ 1.38×10-23J/K； ε为地表辐射率。 2.2 单通道法 单通道算法（Single-channel Method，以下简称SC算法），TM6数据来说，计算公式为： 式中，Ts是地表温度（K）；Lsensor是卫星高度上遥感传感器测得的辐射强度（W?m-2?sr-1?um-1），其计算公式为： 其中，DNTM6为TM6的像元DN值，0≤DNTM6≤255。此外，在上述公式中，ε是地表比辐射率；γ，ψ1，ψ2，ψ3都是中间变量。 2.3 单窗算法 单窗算法（Mono-windowAlgorithm，以下简称MW算法），是有覃志豪等人提出的，计算公式如下： 热带平均大气（北纬15°，年平均）： Ta 17.9769+0.91715T0； 中纬度夏季平均大气（北纬45°，7月）： Ta 16.0110+0.92621T0. 对于LandsatTM数据，大气透射率可从NASA官网，通过输入经纬度可查出区域内平均透射率为0.9。 2.4 地表辐射率的计算 三种算法都会用到地表辐射率，这里直接从地表辐射率库里查找每种地表类型的近似或者经验比辐射率值，我国典型的地物标准波普数据库http：// 。 3.结果分析 3.1 三种算法差异分析 （1）对比b，c，d影像图可知，颜色逐渐加深，温度范围为15-35℃左右，表明地表温度结果逐渐增大，但总体差异不大，从而可知3种算法反演的结果总体趋势比较接近，IB算法反演的结果低一些为26.1℃，SC算法的结果居中为26.875℃，MW算