114130168宋国俊实验五.doc

本科学生实验报告 学号 114130168 姓名 宋国俊 学院 旅地学院 专业、班级 11地信 实验课程名称 GIS分析 实验名称: 提取地表温度 教师及职称 洪亮 开课学期 2014 至 2015 学年 上 学期 填报时间 2014 年 11 月 10 日 云南师范大学教务处编印 实验序号 2 实验名称 提取地表温度 实验时间 2014年11月15 实验室 GIS实验室 一．实验预习 实验材料 昆明2000年左右和2014年左右两个时期（尽量同一天）的遥感影像 实验目的和要求 2、 提取地表温度(LST)，制作地表温度专题图。 3、 提取绿地、城市不透水表面（包括：建筑物、道路等）、水体3类，3类专 题图。 4、 分析地表温度与城市绿地、城市不透水表面、水体分布的时空关系。 实验方法步骤及注意事项 一：提取地表温度(LST)，制作地表温度专题图 (1)对昆明影像数据进行辐射定标： Ⅰ、启动ENVI File→Open External File→Landsat→Geo TIFF with metadata→Enter Landsat MetaData Filenames（输入元数据） Ⅱ、Spectral→Preprocessing→Calibration utilities→ Landsatcalibration→Landsat calibration input file→输入第一步的元数据 （2）将辐射定标后的数据转化为BIL格式 在ENVI主菜单 Basic Tools→Convert Data(BSQ、BIL、BIP)→convert File input File→convert File parameter (3)大气校正：Spectral→FLAASH 得到的结果 (4)植被指数（NDVI）计算，打开Basic tool 打开band math ，比值植被指数RVI： RVI=(b1 le 0)*0+(b2 le 0)*0+((b1 gt 0 or b2 gt 0)*(float(b1)/b2)) B1为热红外波段band4，B2为可见光红波段band3 经辐射定标和大气校正之后的滇池流域 (5) 植被覆盖度计算 计算植被覆盖度Fv采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下： FV = (NDVI- NDVIS)／(NDVIV - NDVIS) （2） 其中，NDVI为归一化差异植被指数，取NDVIV = 0.70和NDVIS = 0.00，且有，当某个像元的NDVI大于0.70时，FV取值为1；当NDVI小于0.00，FV取值为0。 利用ENVI主菜单-Basic Tools-Band Math，在公式输入栏中输入： （b1 gt 0.7）*1+（b1 lt 0.）*0+（b1 ge 0 and b1 le 0.7）*（(b1-0.0)/(0.7-0.0)） b1:选择NDVI图像 得到植被盖度图像。 (6) 地表比辐射率计算 根据前人的研究，将遥感影像分为水体、城镇和自然表面3种类型。本专题采取以下方法计算研究区地表比辐射率：水体像元的比辐射率赋值为0.995，自然表面和城镇像元的比辐射率估算则分别根据下式（3）（4）进行计算： εsurface = 0.9625 + 0.0614FV - 0.0461FV2 （3） εbuilding = 0.9589 + 0.086FV - 0.0671FV2 （4 利用ENVI主菜单-Basic Tools-Band Math，在公式输入栏中输入： （b1 le 0）*0.995+（b1 gt 0 and b1 lt 0.7）*(0.9589 + 0.086*b2 - 0.0671*b2^2)+(b1 ge 0.7)*(0.9625 + 0.0614*b2 - 0.0461*b2^2) b1：NDVI值； b2：植被覆盖度值。 （7）计算相同温度下黑体的辐射亮度值 利用ENVI主菜单-Basic Tools-Band Math，在公式输入栏中输入： （b2-3.39-0.6*（1-b1）*5.12）/（0.6*b1） b1：60m分辨率的地表比辐射率值； b2：表示热红外波段的辐射定标值。 得到了温度为T的黑体在热红外波段的辐射亮度值 （8）反演地表温度 在获取温度为TS的黑体在热红外波段的辐