landsat遥感影像地温度反演教程大气校正法

基于辐射传输方程的Landsat数据地表温度反演教程

一、数据准备

Landsa 8遥感影像数据一景，本教程以重庆市2015年7月26日的=行列号为（128，049）影像（LC81280402016208LGN00）为例。 同时需提前查询影像的基本信息（详见下表） 标识 日期 LC81280402016208LGN00 2016/7/26 …… …… 采集时间 中心经度 3:26:56 …… 中心纬度 106.11288 30.30647 …… …… 注：基本信息在影像头文件中均可查询到，采集时间为格林尼治时间。 二、地表温度反演的总体流程

三、具体步骤 1、辐射定标

地表温度反演主要包括两部分，一是对热红外数据，二是多光谱数据进行辐射定标。

（1）热红外数据辐射定标

选择Radiometric Correction/Radiometric Calibration。在File Selection对话框中，选择数据LC81230322013132LGN02_MTL_Thermal，单击Spectral Subset选择Thermal Infrared1（10.9），打开Radiometric Calibration面板。

Scale factor 不能改变，否则后续计算会报错。保持默认1即可。

（2）多光谱数据辐射定标

选择要校正的多光谱数据“LC81230322013132LGN02_MTL_MultiSpectral” 进行辐射定标。

因为后续需要对多光谱数据进行大气校正，可直接单击Apply Flaash Settings，如下图。

注意与热红外数据辐射定标是的差别，设置后Scale factor值为0.1。 2、大气校正 本教程选择Flaash 校正法。FLAASH Atmospheric Correction，双击此工具，打开辐射定标的数据，进行相关的参数设置进行大气校正。 注意：如果在多光谱数据辐射定标时Scale factor值忘记设置，可在本步骤中打开辐射定标数时设置single scale faceor 值为0.1，若已设置，则默认值为1即可。

1) Input Radiance Image：打开辐射定标结果数据； 2) 设置输出反射率的路径，由于定标时候；

3) 设置输出FLAASH校正文件的路径，最优状态：路径所在磁盘空间足够大；

4) 中心点经纬度Scene Center Location：自动获取；

5) 选择传感器类型：Landsat-8 OLI；其对应的传感器高度以及影像数据的分辨率自动读取；

6) 设置研究区域的地面高程数据；

7) 影像生成时的飞行过境时间：在layer manager中的Lc8数据图层右键选择View Metadata，浏览time字段获取成像时间；

注：也可以从元文件“LC81230322013132LGN02_MTL.txt”中找到，具体名称：DATE_ACQUIRED = 2013-05-12；SCENE_CENTER_TIME = 02:55:26.6336980Z； 8) 大气模型参数选择：Sub-Arctic Summer（根据成像时间和纬度信息选择）； 9) 气溶胶模型Aerosol Model：Urban，气溶胶反演方法Aerosol Retrieval：2-band（K-T）； 10) 其他参数按照默认设置即可。

11) 多光谱参数设置中，

K-T反演选择默认模式：Defaults->Over-Land Retrieval standard（600:2100） 波谱响应函数：默认指向..

\\Program Files\\Exelis\\ENVI51\\classic\\filt_func\\landsat8_oli.sli

把它重新指向：..\\Program Files\\Exelis\\ENVI51\\resource\\filterfuncs\\landsat8_oli.sli