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

土地利用解译操作流程准备工作：在空间数据云平台（/）网站上下载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字段信息，方便变化前信息与变化后信息对比，查找变化信息。（注意：前期可以跳过1113步的获取土地转移矩阵部分，直接做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.打开感兴趣区，并建立5个新的