土地利用解译操作流程准备工作：在空间数据云平台（）网站上下载1992年-2012年遥感影像数据（注意：为确保郑州市在图幅范围内，检索时候输入与郑州市相近的经纬度坐标，控制云量大小，确保数据时间在5-10月份之间）考虑到下载的数据质量比较好，数据都自带坐标且云量基本都为0，省去配准和大气校正步骤，直接加载已知1995年的土地分类矢量图，并对遥感图像进行裁剪。1.由于下载的数据都是分波段的，因此先进行波段合成，具体操作如下：打开，打开File->Openexternalfile->Landsat->GeoTIFF，选择7个TIF格式文件打开，选择RGB显示，勾选742波段，LoadRGB，打开BasicTools->LayerStacking点击ImportFile，输入要合成的波段，选择742波段进行合成并输出一个文件（名为1992）2．加载1995年郑州市土地利用矢量图，具体操作如下：打开File->OpenVetorFilel浏览到要打开的SHP格式文件，点击打开，点OK。选择刚刚打开的SHP图层并加载到之前打开的RGB图层，如下图选择Display#1:得到如下结果：3.根据加载的矢量图对遥感影像进行不规则裁剪，具体操作如下：打开图像感兴趣区域（ROI）如下所图示：点击File->SubsetDataviaROIs，对图像进行裁剪输出文件名为1992裁剪，如图所示：加载裁剪图像，显示裁剪结果：4.按照同样步骤，将矢量图加载到裁剪后图像上（参考第二步）略5．将剪裁图像另存为.img格式文件，命名为1992裁剪.img。6.将1992裁剪.img和1992土地利用矢量文件加载到ArcMap中，如下图所示：7.新建Shpfile图层，定义坐标系（与遥感影像坐标系保持一致），打开编辑和捕捉工具，选择editor->snapping->snappingtoolbar如下图所示：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->AnalysisTools->Overlay->Union分析工具如上图所示，将B-1994年变化图层和图层zhengzhou1995添加到InputFeatures中，输出结果命名为1994，得到结果如下图所示：说明：对于合并后属性信息，删除矢量化过程中造成的明显拓扑错误信息，删除无用字段并添加1994_class新字段信息来表明1994年新增加的地类信息，保留1995-1994年变化信息和1995_class字段信息，方便变化前信息与变化后信息对比，查找变化信息。（注意：前期可以跳过11——13步的获取土地转移矩阵部分，直接做14步）11.分别对1994年、1995年土地利用变化图层做数据融合操作打开ArcToolbox，选择DataManagementTools->Generalization->Dissolve工具。InputFeature选择要融合的图层，OutputFeatureClass选择输出结果存储的位置及名称，DissolveField(s)选择土地利用类型字段（如Type1995），然后勾选Creatmultipartfeatures选项，点击OK完成，如下图所示：重复此过程，对1995年数据也进行融合，查看属性表和融合后图层变化如下图：12.对融合后数据进行叠置分析在ArcMap中打开两个时相融合后的数据，在ArcToolbox中选择AnalysisTools->Overlay->Intersect工具，InputFeatures选择两个时相的图层，OutputFeatureClass选择叠加结果存储的位置，命名为A，其余选项可以忽略，单击【OK】完成，如图所...
