ENVI下土地利用类型转换矩阵与ETM 数据反演地表.docVIP

ENVI下土地利用类型转换矩阵与ETM 数据反演地表 土地利用类型转换(变化)矩阵 查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少(或者面积)分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图(掩膜图像)，它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像(ENVI Classification)。 2、计算转换矩阵 打开两个土地利用分类结果。 (1)在主菜单中，选择Basic Tools→Change Detection→Change Detection Statistics。 (2)分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 (3)在Define Equivalent Classes对话框中(图1)，如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 (4)选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位(Report Type)：像素(Pixels)、百分比(Percent)和面积(Area)。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 (5)单击OK，执行土地利用类型转换矩阵计算过程。 图1 Define Equivalent Classes对话框 图2选择数据参数 3、查看结果 (1)如图3为得到的土地利用类型转换矩阵结果。横字段表示前一时间段(Initial State)的土地利用类别，纵字段为后一时间段(Final State)的土地利用类别。横字段和纵字段交叉处表示变化值，如有平方米林地用地变化为草地。 图3土地利用类型转换矩阵 (2)还可以为每一个地类生成一个变换掩膜图像，图4所示为其中一个地类的掩膜图像。掩膜图像的灰度值表示变化类型，如这里的2{草地}表示林地变化为草地的像元。 图4变化掩膜图像 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