实验报告被列为实验课程考核之中-课程中心-广州大学.doc

地表温度反演实验 开课学院及实验室：地理科学学院 理学实验楼625 2014 年6 月 10 日 学 院 地理科学学院　 年级、专业、班 地理信息系统121 姓名 　 　 实验课程名称 遥感数字影像处理实验　 　 　　 教 师 　 一、实验目的 根据一定的条件和要求，运用有关的原理和实验技能提高实验素质和科学研究能力，进行创造性能力的培养。 二、实验原理 地表温度对理解地气之间的能量交换十分重要。地表辐射温度通过TM和ETM+的热波段辐射值计算，地表温度通过陆地卫星热红外波段，用简化的普朗克函数获得（公式1）(Markham and Barker, 1986)： 式中，是地表温度（单位：K）；是波段6的光谱辐射值；是与NDVI相关的地表比辐射率（公式4）；是校准常数1，是校准常数2，对陆地卫星Landsat 5TM，、分别为60.776 mWcm-2 sr-1 μm-1、1260.56K（对陆地卫星Landsat 7ETM+，、分别为666.09 mWcm-2 sr-1 μm-1、1282.71K）。 对Landsat 5 TM，是DN值的线性函数： 式中，=0.0056322 mWcm-2 sr-1 μm-1，=0.1238 mWcm-2 sr-1 μm-1。 地表比辐射率（同温度下地表辐射能与黑体辐射能的比率）用于通过Landsat的热波段计算地表温度。地表比辐射率可用NDVI作参数，以经验公式进行估计(Bastiaanssen, 2000)： 当时（如水体）地表比辐射率取常数1。 三、使用仪器、材料 装有ERDAS 9.2的PC一台 四、实验步骤及过程记录 1、登入/或并下载相关研究区域遥感影像，我选的是云南大理，条带号为131-42。 2、由于使用的是landsat8的遥感影像，温度反演需要将波段4、5、10进行融合 并利用大理市shpfile对融合后的遥感影像进行裁剪。 Interpreter-Utilities-Layer Stack（波段融合） Interpreter-Utilities-Mask（研究区域裁剪） 3、打开ERDAS 9.2软件，选 择MODELER工具，单击Model?Maker打开工具面板，开始创建模型，生成ndvi.img，数据类型是float 4、从下载的遥感影像头文件中找出计算需要的相关数据值 RADIANCE_MULT_BAND_10 = 3.3420E-04…………………..m RADIANCE_ADD_BAND_10 = 0.10000………………………...d K1_CONSTANT_BAND_10 = 774.89 ……………………………K1 K2_CONSTANT_BAND_10 = 1321.08 …………………………..K2 5、继续运行Model?Maker创建模型，计算值，其中ndvi使用步骤3所得。 6、继续运行Model?Maker创建模型，计算R值，其中m、d值使用步骤4所得。 6、运行Model?Maker创建计算LST模型。所需已知值从上述步骤得知及调用。并输出为img文件 7、将LST计算后输出的img文件用arcmap打开，使用raster to polygon进行转化 8、去除30度（即300K）以上、-5（即268K）度以下和655K的记录 其中30度（即300K）以上、-5（即268K）度以下为误差较大值 655K为遥感图像边界 再符号化各温度段、利用fieldCalculator转化为摄氏度并进行输出 五、实验结果及分析 分析： 1、图像中央白色部分为洱海：洱海古代文献中曾称为叶榆泽、昆弥川、西洱河、西二河等，位于云南大理郊区，为云南省第二大淡水湖，洱海北起洱源，长约42.58公里，东西最大宽度9.0公里，湖面面积256.5平方公里，平均湖深10米，最大湖深达20米。洱海有两个出水口：1.在下关镇附近，经西洱河流出；2.“引洱入宾”。洱海是大理“风花雪月”四景之一“洱海月”之所在。据说因形状像一个耳朵而取名为“洱海”。洱海水质优良，水产资源丰富，同时也是一个有着迤逦风光的风景区。金梭岛在大理市海东乡西南，洱海之中，距洱海东岸250-600米，西与七里桥乡罗久邑隔海相望。 岛呈东南至西北走向，长约420米，宽约200-500米，高出海面260米。岛形两端较高阔，中部较狭，尤如漂浮海面的葫芦，更像一只织布用的梭子，故名金梭岛。 并且通过分析地表温度的分布能得出地形