土地利用解译及精度评价操作流程(最终版)

土地利用解译操作流程土地利用解译操作流程 准备工作 在空间数据云平台 1992 年 2012 年遥感影 像数据 注意 为确保郑州市在图幅范围内 检索时候输入与郑州市相近的经纬度坐标 控制云量大小 确保数据时间在 5 10 月份之间 考虑到下载的数据质量比较好 数据都自带坐标且云量基本都为 0 省去配准和大气校正 步骤 直接加载已知 1995 年的土地分类矢量图 并对遥感图像进行裁剪 1 由于下载的数据都是分波段的 因此先进行波段合成 具体操作如下 打开 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 增加部分数字化提取结果如下图所示 注意 1 在提取增加信息的过程中 由于变化部分图斑面积普遍较小 矢量化过程中记得打 开捕捉工具 避免两个斑块离得太近造成拓扑错误 2 数字化过程中可以利用查询工具快速查找原本地类编码 便于数字化过程中属性表 变化信息编辑 3 将 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 步 步 11 分别对 1994 年 1995 年土地利用变化图层做数据融合操作 打开 ArcToolbox 选择 Data Management Tools Generalization Dissolve 工具 Input Feature 选择要融合的图层 Output Feature Class 选择输出结果存储的位置及名称 Dissolve Field s 选择土地利用类型字段 如 Type1995 然后勾选 Creat multipart features 选项 点 击 OK 完成 如下图所示 重复此过程 对 1995 年数据也进行融合 查看属性表和融合后图层变化如下图 12 对融合后数据进行叠置分析 在 ArcMap 中打开两个时相融合后的数据 在 ArcToolbox 中选择 Analysis Tools Overlay Intersect 工具 Input Features 选择两个时相的图层 Output Feature Class 选择叠加结果存 储的位置 命名为 A 其余选项可以忽略 单击 OK 完成 如图所示 查看 A 图层属性表信息 删除无用字段 添加面积字段 Area 并在属性表中 Area 字段上 单击右键选择 Calculate Geometry 在打开的 Calculate Geometry 对话框中 Property 选 择 Area Units 选择要使用的面积单位 单击 OK 完成图斑面积计算 属性信息如图所 示 13 制作转移矩阵 保存 ARCGIS 中所得结果 用 Excel 打开 A dbf 文件 并另存为 xls 格式文件 在 Excel 中选中所有数据 不要点左上角 只选择有效数据 点击 插入 选项卡 选 择 数据透视表 数据透视表 点击 确定 在打开的数据透视表中按图示将字段拖入相应区域 得到结果如下图所示 14 遥感解译结果精度验证 1 对 1999 B 数据进行融合处理 打开 ArcToolbox 选择 Data Management Tools Generalization Dissolve 工具 Input Feature 选择要融合的图层 1999 B Output Feature Class 选择输出结果存储的位 置及名称 Dissolve Field s 选择土地利用类型字段 如 B class 然后勾选 Creat multipart features 选项 点击 OK 完成 如下图所示 完成后查看新生成文件 1999 B Dis 的属性表 发现全部融合为 5 类 如下图所示 2 在 ArcGIS 中生成随机点 如下图所示在工具箱里打开数据管理工具 要素类 创建随机点 如下图填写输出路径 输出名称 约束要素类选择上一步融合处理结果 长整型改成 50 下边的线性单位可设为 0m 3 添加 XY 坐标 如下图所示打开工具箱中的数据管理工具 要素 添加 XY 坐标 选择上一步生成的图层 1999 B R 并点击确定 结果如下图所示 打开 1999 B R 属性表可以发现增加了 X Y 坐标值如下图 4 结合谷歌地球进行校验并调整地类 将上一步得到的 1999 B R dbf 格式文件用 Excel 方式打开 可以得到每种地类 20 个随机 点的 X Y 坐标值 打开谷歌地球 打开工具 选项 3D 视图 勾选 通用横轴墨卡托投影 5 对比校验 调整地类 点击添加地标按钮分别将不同地类的 X Y 坐标如下图所示输入 Google 地球中 并放大每 一个标记后点 查看标记后点所处地类是否与之前随机点地类相一致 如果有变化则在 Excel 表格中对地类信息进行调整编辑 无变化不用调整 且只调整地类 如 1 调为 2 坐 标值无需调整 全部调整完毕 将新变化地类全部进行重新归类 如之前 2 地类变为 1 地类 则将其直接 移至 1 地类部分 每个地类单独按如下图格式粘贴到记事本中 存为 TXT 格式 6 建立感兴趣区域 打开前期处理好的遥感影像 通过 Overlay Region of Insterest 打开感兴趣区 并建立