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

Landsat TM6 地表温度反演实习报告

实习目的：

1、TM Level 1 数据的热红外波段辐射定标：学会阅读头文件，找出所需定标参数；利用定标参数将TM图像热红外波段DN值转换为辐射亮度；

2、运用单通道法，反演地表温度反演

实习步骤：

1.加载图像：import—>选择需要加载的图像

2.分部建模

2.1求算NDVI的建模如图所示

第三第四波段的辐射定标运算的增益和偏移均来自数据的头文件。

2.2第六波段辐射定标计算

说明：第六波段的辐射定标计算的增益和偏移不能再用头文件中的增益和偏移，否则误差会很大造成无法显示，因此必须在老师给的表格中查找。

然后再用老师给的公式进行计算.

从表格中找出L min 和L max 输入下面公式

255

G min

max L L rescale -=

min B L rescale = rescale cal rescale B Q G L +⋅=λ

即可求出增益和偏移，然后进行辐射定标运算即可求出所需结果。

2.3求解地表反射率（λε）的建模运算。

所使用的经验公式根据Van 的经验公式：

)ln(047.00094.1NDVI +=λε

2.4求解)T (B s λ的建模。)T (B s λ为温度为s T 的黑体在热红外波段的辐射亮度。

使用的公式是

L L o

o o s ↓

↑

--

-=λ

λ

λ

λλλλλεετε1L )T (B ）（，其中

L

o ↓

λ

表示大气向下辐射亮度，模拟结果为1.68 Wm -2um -1Sr -1，L o

↑

λ表示大气向上辐射亮度，

模拟结果为1.74 Wm -2um -1Sr -1，λτo 为大气在热红外波段的透过率，模拟结果为0.77。λε为上一步求解的结果。

2.5反演温度的建模

该步骤使用简化的公式)

176.607ln(1260.56

T +=

λ

B s

2.6得到了地表温度的结果，然而显示的确实全白。

查看数据结果：

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

2.8滤除非法数据的建模

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

3.结果展示：

拉一条直线看其温度变化

根据上图可以看出四周黑色部分已经被滤除值为0，中间部分温度300K左右，即摄氏27度左右。

直方图：

可以看出比率最大的部分温度为297.456K，即温度为摄氏24度左右的区域面积最大。

如图所示最小值是0，最大值329.79，平均值277.272，中间值298.87