地表温度反演实验.docVIP

地表温度反演实验 实验室 理学实验楼 实验日期: 2014 年 6 月 5 日 姓名 鲁俊 学院 地理科学学院 学号 1101200031 年级、专业名称 地信111 实验平台 ERDAS IMAGINE 9.2 一、原理与方法 地表温度对理解地气之间的能量交换十分重要。地表辐射温度通过TM和ETM+的热波段辐射值计算，地表温度通过陆地卫星热红外波段，用简化的普朗克函数获得（公式1）(Markham and Barker, 1986)： 式中，是地表温度（单位：K）；是波段6的光谱辐射值；是与NDVI相关的地表比辐射率（公式4）；是校准常数1，是校准常数2，对陆地卫星Landsat 8 OLI_TIRS，、存放在影像的头文件中。 对Landsat 8 OLI_TIRS，是DN值的线性函数： 式中，=0.0056322 mWcm-2 sr-1 μm-1，=0.1238 mWcm-2 sr-1 μm-1。 地表比辐射率（同温度下地表辐射能与黑体辐射能的比率）用于通过Landsat的热波段计算地表温度。地表比辐射率可用NDVI作参数，以经验公式进行估计(Bastiaanssen, 2000)： 当时（如水体）地表比辐射率取常数1。 二、实验目的 根据一定的条件和要求，运用有关的原理和实验技能提高实验素质和科学研究能力，进行创造性能力的培养。 三、实验步骤 1、理论根据和实验原理； （1）理论依据：地表辐射温度通过OLI_TIRS的热红外波段辐射值计算，地表温度通过陆地卫星热红外波段，使用简化的普朗克函数获得公式： (Markham and Barker, 1986)： 式中，是地表温度（单位：K）；是波段6的光谱辐射值；是与NDVI相关的地表比辐射率（公式4）；是校准常数1，是校准常数2，对陆地卫星Landsat 8 OLI_TIRS，、存放在影像的头文件中。 对Landsat 8 OLI_TIRS，是DN值的线性函数： 地表比辐射率（同温度下地表辐射能与黑体辐射能的比率）用于通过Landsat的热波段计算地表温度。地表比辐射率可用NDVI作参数，以经验公式进行估计(Bastiaanssen, 2000)： 当时（如水体）地表比辐射率取常数1。 2、选择实验所必需的仪器和材料； 实验仪器：华硕笔记本 实验材料：1、佛山三水边界（sanshui.img） 2、Landsat 8 OLI_TIRS 卫星数字产品 Landsat 8 OLI_TIRS影像基本信息： 数据标识 LC81220442013333LGN00 卫星名称 LANDSAT8 数据类型 OLI_TIRS 传感器 OLI_TIRS 接收站 LGN 白天/晚上 DAY 条带号 122 行编号 44 太阳高度角 41.4377 太阳方位角 154.8008 日期信息 获取时间 2013-11-29 ? ? 开始时间 2013-11-29 02:53:28.0 结束时间 2013-11-29 02:53:59.0 3、实验操作过程： （1）将下载到的数据的波段合并，合并的波段为B4、B5、B10。然后对影像进行裁切，裁切的区域为本次研究的区域佛山山水。 使用Erdas的建模工具求出影像的NDVI值。 NDVI的计算公式说明：NDVI=(Band5-Bnad4)/(Band5+Bnad4)(针对Landsat8) 图5 使用Erdas的建模工具求出LST（地表温度）。 其中： RADIANCE_ADD_BAND_10 = 0.10000 （对应上式中的“d”） RADIANCE_MULT_BAND_10 = 3.3420E-04（对应上式中的“m”） K1=774.89 K2=1321.08 最终结果分析： 结果分析： 经过地表温度反演出来佛山三水的温度与真实的温度相似，但是还是存在一定的误差，这或许是与所用的参考值的准确性有关联，所以以后要收集更多的资料来确保数据的精确性。 地表温度反演非常具有可用性，可以用于很多方面，比如说城市热岛监测、计算土壤湿度指数等等。 地表温度反演有很多方法，我们可以学习掌握，比如说：大气校正法、单窗算法、辐射传输方程法等等。