基于tm影像的地表温度反演算法研究

《基于tm影像的地表温度反演算法研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、基于TM影像的地表温度反演算法研究摘要：介绍了LandsatTM5遥感数据地表温度反演的流程，并对辐射定标、大气校正等算法进行了详细的阐述。选择四川省绵阳市和德阳市为实验区，对上述流程和算法进行实验分析，结果表明流程合理，相关算法具有普适性，算法中的大部分参数容易获取，反演效果较好，具有很好的实际应用价值。中国8/vie　　关键词：遥感影像；辐射校正；亮度温度；地表温度　　中图分类号：TP393文献标识码：A：1009-3044（2017）02-0223-03　　1概述　　当前生态环境日益恶化，与人们对

2、生活舒适度要求不断提高形成鲜明对比，人居环境质量状况备受关注。然而，随着我国工业化进程的进一步加深，工厂、汽车等人为热能的排放量也在迅速增加，加之全球气候变暖的大背景下，城市热岛效应这一生态环境问题也日益突出。自1818年城市热岛概念被提出至今已经过去近200年，各国的学者一直致力于该问题的研究。进而形成了三种不同的研究手段即：地面观测法、遥感监测法和数值模拟法[1]，目前，随着航空、航天技术以及计算机等技术的迅猛发展，遥感监测法已成主要手段。利用遥感影像反演地表温度这一指标在城市热岛效应分析中就显得尤

3、为重要。　　2数据与方法　　2.1数据　　研究以四川省绵阳市主城区和四川省德阳市主城区两个区域为实验样区，绵阳市属四川盆地盆中丘陵区北部，主城区以丘陵为主，德阳市位于成都平原东北部，主城区地形以平原为主，两个研究区具有一定的代表性。选择美国陆地卫星LandsatTM5为数据源，TM遥感影像热红外波段空间分辨率是120m，其余波段空间分辨率是30m。两个区域的影像各两个时相，时间均为1988年5月1日和2007年5月6日（轨道号129/38），两个区域影像质均较好清晰无云，适宜做相关的反演和分析。　　2.

4、2地表温度反演流程　　地表温度作为城市热岛效应研究的重要指标，其准确可靠性在后续定量研究中的重要地位是不言而喻的，有关地表温度反演算法的研究也越来越多，尤其是随着Landsat8遥感数据的逐步应用，辐射传输方程法[2]、单窗算法[3，4]和IB算法[5]等算法是目前TM影像最常用的地表温度反演算法。反演的流程主要包括：辐射定标、大气校正、区域裁剪、亮度温度反演、NDVI计算、地表比辐射率估算和地表温度反演几个步骤，详细流程见图1。　　2.3辐射定标　　由于遥感成像的复杂性，大气、太阳高度角等对成像质量都

5、会产生一定的影响，由此，获取的数据必须进行前期处理工作，为减弱其影响需完成辐射定标等工作。对于TM数据而言，根据波段光谱辐射亮度进行辐射定标的模型如下式所示[6，7]：　　[Lλ=Lminλ+Qλ×（Lmaxλ-Lminλ）/Qmax]（1）　　其中：λ代表波段值；Lλ代表光谱辐射值；Qλ则是标定后的像元值；Qmax是DN的最大值理论；Qmin是DN的最小值理论。Lminλ和Lmaxλ可在数据头文件中获得。　　2.4大气校正　　针对TM数据进行大气校正的方法比较多，各种算法也相对成熟，主要包括：辐射传输

6、方程法、直方图最小值去除法和COST模型法[8]。COST模型法具有操作简单、所需参数容易获取，校正精度相对较高等特点，应用比较广泛。其原理见如下：　　[LIhaze=LImin-LI1%]（2）　　上式：[LIhaze]为光谱辐射值；[LImin]代表传最小光谱辐射值；[LI1%]代表反射率为1%的黑体辐射值。　　经COST模型完成TM影像反射波段大气校正工作。图2展示了第3波段辐射校正前、后效果。辐射校正完成后DN值发生了明显变化，整体普遍增加，但是增加的幅度各不相同。　　2.5亮度温度反演　　由于

7、大多数传感器探测到的是城市下垫面的辐射温度（亮度温度，简称亮温）辐射定标和大气校正工作完成后，接下来进行亮度温度反演，模型如下：　　T=K2/ln（1+k1/Lλ）（3）　　式中，T代表亮温（K），K1=60.766m影像分辨率又相对较低，所以研究时将城市下垫面粗略认为是建成区、水体、自然表面的混合体，总结前人的成果，本研究认为水体像元的比辐射率选择为0.995比较合理。自然表面（公式4）和建成区（公式5）地表比辐射率的计算采用下式进行[10]。　　[εsurface=0.9625+0.0614FV-0

8、.0461FV2]（4）　　[εbuild-up=0.9589+0.086FV-0.0671FV2]（5）　　其中：εsurface为自然表面像元的比辐射率；εbuild-up为建成区像元的比辐射率；FV为植被覆盖度，它由下式计算获得。　　[FV=（NDVI-NDVIS）/（NDVIV-NDVIS）]（6）　　其中：NDVI是归一化植被指数。　　[NDVI=（ρ4-ρ3）/（ρ4+ρ3）]（7）　　其中：ρ3是3波段地表反射率，和ρ4是4波