erdas 北京地表温度反演_实习报告

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反演温度的建模

该步骤使用简化的公式Ts?

1260.56

607.76ln(?1)B?2.6得到了地表温度的结果，然而显示的确实全白。

查

看

数

据

结

果

：

有大面积的非法数据无法显示，因此必须想办法滤除不能显示的数据。 2.8滤除非法数据的建模

和同学一切讨论，一起上网查资料，我选择了分类函数CONDITIONAL(),上面的公式的意思是如果n1_wendu1的栅格值大于0，则返回该栅格值，否则返回0. 3. 结果展示：

拉一条直线看其温度变化

本文来源：https://www.bwwdw.com/article/n1gg.html