土地利用解译及精度评价与衡量操作流程最终版.docx
- 文档编号:8835737
- 上传时间:2023-05-15
- 格式:DOCX
- 页数:33
- 大小:3.93MB
土地利用解译及精度评价与衡量操作流程最终版.docx
《土地利用解译及精度评价与衡量操作流程最终版.docx》由会员分享,可在线阅读,更多相关《土地利用解译及精度评价与衡量操作流程最终版.docx(33页珍藏版)》请在冰点文库上搜索。
土地利用解译及精度评价与衡量操作流程最终版
土地利用解译操作流程
准备工作:
在空间数据云平台〔〕上下载1992年-2012年遥感影像数据〔注意:
为确保市在图幅围,检索时候输入与市相近的经纬度坐标,控制云量大小,确保数据时间在5-10月份之间〕
考虑到下载的数据质量比拟好,数据都自带坐标且云量根本都为0,省去配准和大气校正步骤,直接加载1995年的土地分类矢量图,并对遥感图像进展裁剪。
1.由于下载的数据都是分波段的,因此先进展波段合成,具体操作如下:
打开ENVI4.8,打开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】完成,如下列图:
查看A图层属性表信息,删除无用字段,添加面积字段Area,并在属性表中Area字段上单击右键选择CalculateGeometry…,在打开的CalculateGeometry对话框中,Property选择Area,Units选择要使用的面积单位,单击【OK】完成图斑面积计算,属性信息如下列图:
13.制作转移矩阵
保存ARCGIS中所得结果,用Excel打开A.dbf文件,并另存为.xls格式文件。
在Excel中选中所有数据〔不要点左上角,只选择有效数据〕,点击【插入】选项卡,选择【数据透视表】|【数据透视表】,点击【确定】。
在打开的数据透视表中按图示将字段拖入相应区域。
得到结果如如下图所示:
14.遥感解译结果精度验证
〔1〕对1999_B数据进展融合处理
打开ArcToolbox,选择DataManagementTools->Generalization->Dissolve工具。
InputFeature选择要融合的图层1999_B,OutputFeatureClass选择输出结果存储的位置与名称,DissolveField(s)选择土地利用类型字段〔如B_class〕,然后勾选Creatmultipartfeatures选项,点击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->RegionofInsterest...打开感兴趣区,
并建立5个新的地类,如如下图所示:
每个地类选择“ROI_Type->inputfromASCII〞,把保存成TXT格式的随机点一一加载进去,如如下图所示X选择2,Y选择3。
点击确定。
建立好后,保存ROI。
(7)精度检验
如如下图所示打开Classification->CreateClassImagefromROIS
点击OK得到最终精度评价结果:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 土地利用 解译 精度 评价 衡量 操作 流程 最终版