Landsat TM6 地表温度反演实习报告 实习目的: 1、TM Lev el 1 数据的热红外波段辐射定标:学会阅读头文件,找出所需定标参数;利用定标参数将TM 图像热红外波段DN 值转换为辐射亮度; 2、运用单通道法,反演地表温度反演 实习步骤: 1. 加载图像:import—>选择需要加载的图像 2. 分部建模 2.1 求算 NDVI 的建模如图所示 第三第四波段的辐射定标运算的增益和偏移均来自数据的头文件。 2 .2 第六波段辐射定标计算 说明:第六波段的辐射定标计算的增益和偏移不能再用头文件中的增益和偏移,否则误差会很大造成无法显示,因此必须在老师给的表格中查找。 然后再用老师给的公式进行计算. 从表格中找出 Lmin和Lmax 输入下面公式 255GminmaxLLrescale minBLrescale rescalecalrescaleBQGL 即可求出增益和偏移,然后进行辐射定标运算即可求出所需结果。 2.3 求解地表反射率( )的建模运算。 所使用的经验公式根据 Van 的经验公式: )ln (047.00094.1NDVI 2.4 求解)T(Bs的建模。)T(Bs为温度为sT 的黑体在热红外波段的辐射亮度。 使用的公式是LLooos1L)T(B)(,其中Lo 表示大气向下辐射亮度,模拟结果为1.68 Wm-2u m-1Sr-1,Lo 表示大气向上辐射亮度,模拟结果为1.74 Wm-2u m-1Sr-1,o 为大气在热红外波段的透过率,模拟结果为0.77。为上一步求解的结果。 2.5 反演温度的建模 该步骤使用简化的公式)176.607ln (1260.56TBs 2.6 得到了地表温度的结果,然而显示的确实全白。 查看数据结果: 有大面积的非法数据无法显示,因此必须想办法滤除不能显示的数据。 2 .8 滤除非法数据的建模 和同学一切讨论,一起上网查资料,我选择了分类函数 CONDITIONAL(),上面的公式的意思是如果 n1_w endu1 的栅格值大于 0,则返回该栅格值,否则返回 0. 3. 结果展示: 拉一条直线看其温度变化 根据上图可以看出四周黑色部分已经被滤除值为0,中间部分温度300K 左右,即摄氏27度左右。 直方图: 可以看出比率最大的部分温度为297.456K ,即温度为摄氏24 度左右的区域面积最大。 如图所示最小值是0 ,最大值3 2 9 .7 9 ,平均值2 7 7 .2 7 2 ,中间值2 9 8 .8 7