遥感实验7.docVIP

实验7 利用ETM+数据反演地表温度 1 实验目的与要求 实习目的： 利用成都市Landsat 7/8遥感影像中的热红外波段数据，计算两个时期内地表亮度温度，并分析地表温度的空间特征。 实习要求： ⑴ 掌握利用遥感影像进行地表温度反演的原理与方法。 ⑵ 几何校正。由于热红外波段影像丢失空间参考信息，故需要进行几何校正（配准）处理，实习中建议使用多项式法。 ⑶ 地表亮温空间特征分析。对计算结果进行分析，重点分析地表温度的空间分异特征，并分析其与地表覆被变化之间的关系。 ⑷ 采用ArcGIS、ENVI或Erdas软件进行地表温度。 ⑸ 详细记录操作过程，比较各种算法的优缺点并整理成实验报告上交。 2 实验数据 ⑴ 成都市Landsat 7/8遥感影像中的热红外波段数据。 ⑵ 矢量数据：成都市行政边界。 3 实验分组与时间 ⑴ 本次实验采取分组方式进行，每组4-5人，设组长一名，组长负责协调并记录每人的操作内容（工作量）。 ⑵ 每组交一份实验报告。 ⑶ 时间为两次实验课，如果可能完成不了，每组利用课后时间完成。 4 实验性质 本实验属于软件操作练习。 5 参考材料 5.1 计算公式 ⑴ 辐射定标运算。即由像元DN值计算辐射亮度L： ⑵ 亮度温度计算。即由辐射亮度L计算亮度温度： 其中，K1的取值为666.09W/(m2srμm)，K2的取值为1282.71K。 5.2 ENVI下利用ETM+数据反演地表温度 地表温度作为地球环境分析的重要指标，而遥感技术作为现代重要的对地观测手段，使得基于遥感图像的地表温度反演的研究越来越多。主要的地表温度反演方法有：大气校正法，单窗算法，单通道法等等。本文介绍用辐射传输方程法对地表温度进行反演。 技术流程： 例子数据为2002年9月2日的襄樊市Landsat ETM+数据。根据数据的特点以及地表温度反演研究的技术要求，采用的技术路线为：先对Landsat ETM+数据进行预处理：数据读取、辐射定标、大气校正、襄樊区域裁剪，利用大气校正，即：辐射传输方程法对其影像热红外波段数据进行操作反演，实现襄樊市地区的地表真实温度的反演研究。具体的处理流程如下： 具体的实现步骤如下： 第一步：准备数据 热红外数据使用的是Landsat的第六波段，已经做了传感器定标、几何校正、工程区裁剪，详细流程参考上面的流程图。文件为TM6-rad-subset-jz-xiangfan.img。 由TM影像（已经过大气校正）生成的NDVI数据，已经利用主菜单-Basic Tools-Resize Data（SFatial/SFectral）重采样为60米分辨率，与TMi6数据保持一致，文件名为：TM-NDVI-60m.img。 第二步：地表比辐射率计算 物体的比辐射率是物体向外辐射电磁波的能力表征。它不仅依赖于地表物体的组成，而且与物体的表面状态(表面粗糙度等)及物理性质(介电常数、含水量等)有关，并随着所测定的波长和观测角度等因素有关。在大尺度上对比辐射率精确测量的难度很大，目前只是基于某些假设获得比辐射率的相对值，本文主要根据可见光和近红外光谱信息来估计比辐射率。 (一) 植被覆盖度计算 计算植被覆盖度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图像 得到植被盖度图像。 (二) 地表比辐射率计算 根据前人的研究，将遥感影像分为水体、城镇和自然表面3种类型。本专题采取以下方法计算研究区地表比辐射率：水体像元的比辐射率赋值为0.995，自然表面和城镇像元的比辐射率估算则分别根据下式（3）（4）进行计算： εsurface?= 0.9625 + 0.0614FV?- 0.0461FV2? ?（3） εbuilding?= 0.9589 + 0.086FV?- 0.0671FV2??? （4） 式中，εsurface和εbuilding分别代表自然表面像元和城镇像元的比辐射率。 利用ENVI主菜单-Basic