土地利用解译及精度评价操作流程最终版剖析.docVIP

土地利用解译操作流程 准备工作：在空间数据云平台（/）网站上下载1992年-2012年遥感影像数据（注意：为确保郑州市在图幅范围内，检索时候输入与郑州市相近的经纬度坐标，控制云量大小，确保数据时间在5-10月份之间） 考虑到下载的数据质量比较好，数据都自带坐标且云量基本都为0，省去配准和大气校正步骤，直接加载已知1995年的土地分类矢量图，并对遥感图像进行裁剪。 由于下载的数据都是分波段的，因此先进行波段合成，具体操作如下： 打开ENVI4.8，打开File-Open external file-Landsat-GeoTIFF，选择7个TIF格式文件打开，选择RGB显示，勾选742波段，Load RGB，打开Basic Tools-Layer Stacking 点击Import File，输入要合成的波段，选择742波段进行合成并输出一个文件（名为1992） 2．加载1995年郑州市土地利用矢量图，具体操作如下： 打开File-Open Vetor Filel浏览到要打开的SHP格式文件，点击打开，点OK。 选择刚刚打开的SHP图层并加载到之前打开的RGB图层，如下图选择Display #1: 得到如下结果： 3.根据加载的矢量图对遥感影像进行不规则裁剪，具体操作如下： 打开图像感兴趣区域（ROI）如下所图示： 点击File-Subset Data via ROIs，对图像进行裁剪 输出文件名为1992裁剪，如图所示： 加载裁剪图像，显示裁剪结果： 4.按照同样步骤，将矢量图加载到裁剪后图像上（参考第二步）略 5．将剪裁图像另存为.img格式文件，命名为1992裁剪.img。 6.将1992裁剪.img和1992土地利用矢量文件加载到ArcMap中，如下图所示： 7.新建Shpfile图层，定义坐标系（与遥感影像坐标系保持一致），打开编辑和捕捉工具，选择editor-snapping-snapping toolbar如下图所示： 8.比对1994年与1995年的土地利用变化，将增加变化类型提取矢量化，并编辑属性表，打开新建图层属性表，并添加字段信息，如下图所示1-3代表该增加部分是由类型1用地转化成类型3用地的。 9.增加部分数字化提取结果如下图所示： 注意： 在提取增加信息的过程中，由于变化部分图斑面积普遍较小，矢量化过程中记得打开捕捉工具，避免两个斑块离得太近造成拓扑错误； 数字化过程中可以利用查询工具快速查找原本地类编码，便于数字化过程中属性表变化信息编辑； 将1995年土地利用属性表中的地类编码进行合并，如将123、113等全部合并为第1类，将51、52、53等全部合并为第5类，整理属性表信息，删除冗余部分，主要属性信息应包括ID、地类编码、土地面积、变化信息、变化面积等基本信息。 说明： 对增加部分进行提取时，是以1995年为基准进行推算提取的，对于1995年以前的年份所提取出来的增加信息实际上来讲是返减信息，如1-3类型变化，按照时间顺序来看，应该是3-1的变化，而对于1995年后的数据则无此问题。 10.对图层zhengzhou1995和B-1994年变化信息图层做合并操作，得到1994年土地利用图层打开ArcToolbox-Analysis Tools-Overlay-Union分析工具 如上图所示，将B-1994年变化图层和图层zhengzhou1995添加到Input Features中，输出结果命名为1994，得到结果如下图所示： 说明： 对于合并后属性信息，删除矢量化过程中造成的明显拓扑错误信息，删除无用字段并添加1994_class新字段信息来表明1994年新增加的地类信息，保留1995-1994年变化信息和1995_class字段信息，方便变化前信息与变化后信息对比，查找变化信息。 （注意：前期可以跳过11——13步的获取土地转移矩阵部分，直接做14步） 分别对1994年、1995年土地利用变化图层做数据融合操作 打开ArcToolbox，选择Data Management Tools -Generalization -Dissolve工具。Input Feature选择要融合的图层，Output Feature Class选择输出结果存储的位置及名称，Dissolve Field(s)选择土地利用类型字段（如Type1995），然后勾选Creat multipart features选项，点击OK完成，如下图所示： 重复此过程，对1995年数据也进行融合，查看属性表和融合后图层变化如下图： 对融合后数据进行叠置分析 在ArcMap中打开两个时相融合后的数据，在ArcToolbox中选择Analysis Tools -Overlay -Intersect工具，Input Features选择两个时