土地利用解译及精度评价操作流程最终版docWord文档格式.docx
- 文档编号:17802713
- 上传时间:2022-12-10
- 格式:DOCX
- 页数:11
- 大小:18.56KB
土地利用解译及精度评价操作流程最终版docWord文档格式.docx
《土地利用解译及精度评价操作流程最终版docWord文档格式.docx》由会员分享,可在线阅读,更多相关《土地利用解译及精度评价操作流程最终版docWord文档格式.docx(11页珍藏版)》请在冰豆网上搜索。
选择RGB显示,勾选742波段,LoadRGB,打开BasicTools->
LayerStacking
点击ImportFile,输入要合成的波段,选择742波段进行合成并输出一个文件(名为1992)
2.加载1995年郑州市土地利用矢量图,具体操作如下:
打开File->
OpenVetorFilel浏览到要打开的SHP格式文件,点击打开,点OK。
1
选择刚刚打开的SHP图层并加载到之前打开的RGB图层,如下图选择Display#1:
得到如下结果:
2
3.根据加载的矢量图对遥感影像进行不规则裁剪,具体操作如下:
3
打开图像感兴趣区域(ROI)如下所图示:
4
点击File->
SubsetDataviaROIs,对图像进行裁剪
5
输出文件名为1992裁剪,如图所示:
6
加载裁剪图像,显示裁剪结果:
4.按照同样步骤,将矢量图加载到裁剪后图像上(参考第二步)略
7
5.将剪裁图像另存为.img格式文件,命名为1992裁剪.img。
8
6.将1992裁剪.img和1992土地利用矢量文件加载到ArcMap中,如下图所示:
7.新建Shpfile图层,定义坐标系(与遥感影像坐标系保持一致),打开编辑和捕捉工具,
选择editor->
snapping->
snappingtoolbar如下图所示:
9
8.比对1994年与1995年的土地利用变化,将增加变化类型提取矢量化,并编辑属性表,打
开新建图层属性表,并添加字段信息,如下图所示1-3代表该增加部分是由类型1用地转化
成类型3用地的。
10
9.增加部分数字化提取结果如下图所示:
注意:
11
(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分析工具
12
如上图所示,将B-1994年变化图层和图层zhengzhou1995添加到InputFeatures中,输出
结果命名为1994,得到结果如下图所示:
13
对于合并后属性信息,删除矢量化过程中造成的明显拓扑错误信息,删除无用字段并添加
1994_class新字段信息来表明1994年新增加的地类信息,保留1995-1994年变化信息和
1995_class字段信息,方便变化前信息与变化后信息对比,查找变化信息。
(注意:
前期可以跳过11——13步的获取土地转移矩阵部分,直接做14步)
11.分别对1994年、1995年土地利用变化图层做数据融合操作
14
打开ArcToolbox,选择DataManagementTools->
Generalization->
Dissolve工具。
Input
Feature选择要融合的图层,OutputFeatureClass选择输出结果存储的位置及名称,
DissolveField(s)选择土地利用类型字段(如Type1995),然后勾选Creatmultipart
features选项,点击OK完成,如下图所示:
重复此过程,对1995年数据也进行融合,查看属性表和融合后图层变化如下图:
15
12.对融合后数据进行叠置分析
在ArcMap中打开两个时相融合后的数据,在ArcToolbox中选择AnalysisTools->
Overlay
->
Intersect工具,InputFeatures选择两个时相的图层,OutputFeatureClass选择叠
加结果存储的位置,命名为A,其余选项可以忽略,单击【OK】完成,如图所示:
16
查看A图层属性表信息,删除无用字段,添加面积字段Area,并在属性表中Area字段上单
击右键选择CalculateGeometry,在打开的CalculateGeometry对话框中,Property
选择Area,Units选择要使用的面积单位,单击【OK】完成图斑面积计算,属性信息如图所
示:
13.制作转移矩阵
保存ARCGIS中所得结果,用Excel打开文件,并另存为.xls格式文件。
17
在Excel中选中所有数据(不要点左上角,只选择有效数据),点击【插入】选项卡,选择【数据透视表】|【数据透视表】,点击【确定】。
18
在打开的数据透视表中按图示将字段拖入相应区域。
19
得到结果如下图所示:
20
14.遥感解译结果精度验证
(1)对1999_B数据进行融合处理
Feature选择要融合的图层1999_B,OutputFeatureClass选择输出结果存储的位置及名
称,DissolveField(s)选择土地利用类型字段(如B_class),然后勾选Creatmultipart
完成后查看新生成文件1999_B_Dis的属性表,发现全部融合为5类,如下图所示:
(2)在ArcGIS中生成随机点
如下图所示在工具箱里打开数据管理工具->
要素类->
创建随机点
21
如下图填写输出路径、输出名称、约束要素类选择上一步融合处理结果,长整型改成
50,
下边的线性单位可设为0m。
(3)添加XY坐标
如下图所示打开工具箱中的数据管理工具->
要素->
添加XY坐标
22
选择上一步生成的图层1999_B_R,并点击确定。
23
结果如下图所示:
打开1999_B_R属性表可以发现增加了X、Y坐标值如下图:
24
(4)结合谷歌地球进行校验并调整地类
将上一步得到的格式文件用Excel方式打开,可以得到每种地类20个随机点的X、Y坐标值。
打开谷歌地球,打开工具->
选项->
3D视图,勾选“通用横轴墨卡托投影”
25
(5)对比校验,调整地类
点击添加地标按钮分别将不同地类的X、Y坐标如下图所示输入Google地球中,并放大每一
个标记后点,查看标记后点所处地类是否与之前随机点地类相一致,如果有变化则在Excel
表格中对地类信息进行调整编辑(无变化不用调整,且只调整地类,如1调为2,坐标值无
需调整)
26
全部调整完毕,将新变化地类全部进行重新归类,如之前2地类变为1地类,则将其直接移
至1地类部分,每个地类单独按如下图格式粘贴到记事本中,存为.TXT格式:
(6)建立感兴趣区域
打开前期处理好的遥感影像,通过Overlay->
RegionofInsterest...打开感兴趣区,
并建立5个新的地类,如下图所示:
27
每个地类选择“ROI_Type->
inputfromASCII”,把保存成TXT格式的随机点一一加载进去,
如下图所示X选择2,Y选择3。
点击确定。
28
建立好后,保存ROI。
29
(7)精度检验
如下图所示打开Classification->
CreateClassImagefromROIS
30
31
点击OK得到最终精度评价结果:
32
33
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 土地利用 解译 精度 评价 操作 流程 最终版 doc