Landsat TM6 地表温度反演实习报告
实习目的:
1、TM Level 1 数据的热红外波段辐射定标:学会阅读头文件,找出所需定标参数;利用定标参数将TM图像热红外波段DN值转换为辐射亮度; 2、运用单通道法,反演地表温度反演
实习步骤:
1. 加载图像:import—>选择需要加载的图像
2. 分部建模
2.1求算NDVI的建模如图所示
第三第四波段的辐射定标运算的增益和偏移均来自数据的头文件。 2.2第六波段辐射定标计算
说明:第六波段的辐射定标计算的增益和偏移不能再用头文件中的增益和偏移,否则误差会很大造成无法显示,因此必须在老师给的表格中查找。
然后再用老师给的公式进行计算.
从表格中找出Lmin和Lmax 输入下面公式
Grescale?Lmax?Lmin
255Brescale?Lmin L??Grescale?Qcal?Brescale
即可求出增益和偏移,然后进行辐射定标运算即可求出所需结果。
2.3求解地表反射率(??)的建模运算。
所使用的经验公式根据Van 的经验公式:
???1.0094?0.047ln(NDVI)
2.4求解B?(Ts)的建模。B?(Ts)为温度为Ts的黑体在热红外波段的辐射亮度。
使用的公式是
?(L??Lo?)1???B?(Ts)??????o???L??o,其中
L表示大气向下辐射亮度,模拟结果为1.68 Wm-2um-1Sr-1,o?L?表示大气向上辐射亮度,
o?模拟结果为1.74 Wm-2um-1Sr-1,?o?为大气在热红外波段的透过率,模拟结果为0.77。??为上一步求解的结果。
2.5反演温度的建模