土地利用解译及精度评价操作流程

土地利用解译操作流程准备工作：在空间数据云平台（）网站上下载1992年-2012年遥感影像数据（注意：为确保郑州市在图幅范围内，检索时候输入与郑州市邻近的经纬度坐标，控制云量大小，保证数据时间在5-10月份之间）考虑到下载的数据质量比较好，数据都自带坐标且云量基本都为0，省去配准和大气校订步骤，直接加载已知1995年的土地分类矢量图，并对遥感图像进行裁剪。因为下载的数据都是分波段的，所以先进行波段合成，详细操作以下：翻开，翻开File->Openexternalfile->Landsat->GeoTIFF，选择7个TIF格式文件翻开，选择RGB显示，勾选742波段，LoadRGB，翻开BasicTools->LayerStacking点击ImportFile，输入要合成的波段，选择742波段进行合成并输出一个文件2．加载1995年郑州市土地利用矢量图，详细操作以下：翻开File->OpenVetorFilel阅读到要翻开的SHP格式文件，点击翻开，点

（名为OK。

1992）选择刚才翻开的SHP图层并加载到从前翻开的RGB图层，以下列图选择Display#1:获得以下结果：依据加载的矢量图对遥感影像进行不规则裁剪，详细操作以下：翻开图像感兴趣地区（ROI）以下所图示：点击File->SubsetDataviaROIs，对图像进行裁剪输出文件名为1992裁剪，以下图：加载裁剪图像，显示裁剪结果：依据相同步骤，将矢量图加载到裁剪后图像上（参照第二步）略5．将剪裁图像另存为.img格式文件，命名为1992裁剪.img。6.将1992裁剪.img和1992土地利用矢量文件加载到ArcMap中，以下列图所示：新建Shpfile图层，定义坐标系（与遥感影像坐标系保持一致），翻开编写和捕获工具，选择editor->snapping->snappingtoolbar以下列图所示：8.比对1994年与1995年的土地利用变化，将增添变化种类提取矢量化，并编写属性表，打开新建图层属性表，并增添字段信息，以下列图所示1-3代表该增添部分是由种类1用地转变成种类3用地的。增添部分数字化提取结果以下列图所示：注意：1）在提取增添信息的过程中，因为变化部分图斑面积广泛较小，矢量化过程中记得翻开捕获工具，防止两个斑块离得太近造成拓扑错误；2）数字化过程中能够利用查问工具迅速查找原当地类编码，便于数字化过程中属性表变化信息编写；（3）将1995年土地利用属性表中的地类编码进行归并，如将123、113等所有归并为第1类，将51、52、53等所有归并为第5类，整理属性表信息，删除冗余部分，主要属性信息应包含ID、地类编码、土地面积、变化信息、变化面积等基本信息。说明：对增添部分进行提取时，是以1995年为基准进行计算提取的，关于1995年从前的年份所提拿出来的增添信息实质上来讲是返减信息，如1-3种类变化，依据时间次序来看，应当是3-1的变化，而关于1995年后的数据则无此问题。10.对图层zhengzhou1995和B-1994年变化信息图层做归并操作，获得层翻开ArcToolbox->AnalysisTools->Overlay->Union剖析工具

1994年土地利用图如上图所示，将B-1994年变化图层和图层结果命名为1994，获得结果以下列图所示：

zhengzhou1995

增添到InputFeatures

中，输出说明：关于归并后属性信息，删除矢量化过程中造成的显然拓扑错误信息，删除无用字段并增添1994_class新字段信息来表示1994年新增添的地类信息，保存1995-1994年变化信息和1995_class字段信息，方便变化前信息与变化后信息对照，查找变化信息。（注意：先期能够跳过11——13步的获得土地转移矩阵部分，直接做14步）分别对1994年、1995年土地利用变化图层做数据交融操作翻开ArcToolbox，选择DataManagementTools->Generalization

->Dissolve

工具。InputFeature选摘要交融的图层，OutputFeatureClass选择输出结果储存的地点及名称，DissolveField(s)选择土地利用种类字段（如Type1995），而后勾选Creatmultipartfeatures选项，点击OK达成，以下列图所示：重复此过程，对1995年数据也进行交融，查察属性表和交融后图层变化以下列图：对交融后数据进行叠置剖析在ArcMap中翻开两个时相交融后的数据，在ArcToolbox中选择AnalysisTools->Intersect工具，InputFeatures选择两个时相的图层，OutputFeatureClass

->Overlay选择叠加结果储存的地点，命名为

A，其他选项能够忽视，单击【

OK】达成，以下图：查察A图层属性表信息，删除无用字段，增添面积字段Area，并在属性表中Area字段上单击右键选择CalculateGeometry，在翻开的CalculateGeometry对话框中，Property选择Area，Units选摘要使用的面积单位，单击【OK】达成图斑面积计算，属性信息如图所示：制作转移矩阵保存ARCGIS中所得结果，用Excel翻开文件，并另存为.xls格式文件。在Excel中选中所有数据（不重点左上角，只选择有效数据），点击【插入】选项卡，选择【数据透视表】|【数据透视表】，点击【确立】。在翻开的数据透视表中按图示将字段拖入相应地区。获得结果以下列图所示：遥感解译结果精度考证（1）对1999_B数据进行交融办理翻开ArcToolbox，选择DataManagementTools->Generalization

->Dissolve

工具。InputFeature选摘要交融的图层1999_B，OutputFeatureClass称，DissolveField(s)选择土地利用种类字段（如B_classfeatures选项，点击OK达成，以下列图所示：

选择输出结果储存的地点及名），而后勾选Creatmultipart达成后查察重生成文件1999_B_Dis的属性表，发现所有交融为5类，以下列图所示：2）在ArcGIS中生成随机点以下列图所示在工具箱里翻开数据管理工具->因素类->创立随机点以下列图填写输出路径、输有名称、拘束因素类选择上一步交融办理结果，长整型改成50，下面的线性单位可设为0m。（3）增添XY坐标以下列图所示翻动工具箱中的数据管理工具->因素->增添XY坐标选择上一步生成的图层1999_B_R，并点击确立。结果以下列图所示：翻开1999_B_R属性表能够发现增添了X、Y坐标值以下列图：（4）联合谷歌地球进行校验并调整地类将上一步获得的格式文件用Excel方式翻开，能够获得每种地类20个随机点的X、Y坐标值。翻开谷歌地球，翻动工具->选项->3D视图，勾选“通用横轴墨卡托投影”（5）对照校验，调整地类点击增添地标按钮分别将不一样地类的X、Y坐标以下列图所示输入Google地球中，并放大每一个标志后点，查察标志后点所处地类能否与从前随机点地类相一致，假如有变化则在Excel表格中对地类信息进行调整编写（无变化不用调整，且只调整地类，如1调为2，坐标值无需调整）所有调整完成，将新变化地类所有进行从头归类，如从前2地类变成1地类，则将其直接移至1地类部分，每个地类独自按以下列图格式粘贴到记事本中，存为.TXT格式：（6）成立感兴趣地区翻开先期办理好的遥感影像，经过Overlay->RegionofInsterest...翻开感兴趣区，并成立5个