基于辐射传输方程的Landsat数据地表温度反演教程
一、数据准备
Landsa 8遥感影像数据一景，本教程以重庆市2015年7月26日的=行列号为（128，049）影像（LC8LGN00）为例。 同时需提前查询影像的基本信息（详见下表） 标识 日期 采集时间 中心经度 中心纬度 】 LC8LGN00 3:26:56 2016/7/26 …… …… …… …… …… 注：基本信息在影像头文件中均可查询到，采集时间为格林尼治时间。 ;
二、地表温度反演的总体流程
三、具体步骤 1、辐射定标
地表温度反演主要包括两部分，一是对热红外数据，二是多光谱数据进行辐射定标。
（1）热红外数据辐射定标
选择Radiometric Correction/Radiometric Calibration。在File Selection对话框中，选择数据LC8LGN02_MTL_Thermal，单击Spectral Subset选择Thermal Infrared1（），打开Radiometric Calibration面板。
Scale factor 不能改变，否则后续计算会报错。保持默认1即可。
?
（2）多光谱数据辐射定标
选择要校正的多光谱数据“LC8LGN02_MTL_MultiSpectral” 进行辐射定标。
因为后续需要对多光谱数据进行大气校正，可直接单击Apply Flaash Settings，如下图。
注意与热红外数据辐射定标是的差别，设置后Scale factor值为。
2、大气校正
本教程选择Flaash 校正法。FLAASH Atmospheric Correction，双击此工具，打开辐射定标的数据，进行相关的参数设置进行大气校正。
注意：如果在多光谱数据辐射定标时Scale factor值忘记设置，可在本步骤中打开辐射定标数时设置single scale faceor 值为，若已设置，则默认值为1即可。
…
1) Input Radiance Image：打开辐射定标结果数据； 2) 设置输出反射率的路径，由于定标时候；
3) 设置输出FLAASH校正文件的路径，最优状态：路径所在磁盘空间足够大； 4) 中心点经纬度Scene Center Location：自动获取；
5) 选择传感器类型：Landsat-8 OLI；其对应的传感器高度以及影像数据的分辨率自动读取；
6) 设置研究区域的地面高程数据；
7) 影像生成时的飞行过境时间：在layer manager中的Lc8数据图层右键选择View Metadata，浏览time字段获取成像时间；
注：也可以从元文件“”中找到，具体名称：DATE_ACQUIRED = 2013-05-12；SCENE_CENTER_TIME = 02:55:；
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\\
把它重新指向：..\\Program Files\\Exelis\\ENVI51\\resource\\filterfuncs\\
注：这是因为版本的一个小bug，即Classic中的L8的波谱响应函数不正确，另外一个一劳永逸的方法是：将“\\Exelis\\ENVI51\\resource\\filterfuncs”中的 和 两个文件拷贝覆盖：“...\\ENVI51\\classic\\filt_func” 中的两个文件。 否则SWIR1波段大气校正后的结果全为0。见下图
%
12) 高级参数设置：根据内存大小设置Tile Size（Mb）：100(8g物理内存)，其他参数默认即可，详细见下图
经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响，植被的波谱曲线趋于正常。
3、地表比辐射率计算 （1）植被覆盖度计算
计算植被覆盖度Fv采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下： FV = (NDVI- NDVIS)／(NDVIV - NDVIS)
landsat 遥感影像地表温度反演教程(大气校正法)
