在一些土地利用变化的分析展示中,土地利用转移矩阵是最直观、常见的。之前的文章中我们分享了IDRISI这款软件【分享】土地利用变化模拟预测软件IDRISI、教程,从这篇文章开始,我们将利用ArcGIS软件并桥接IDRISI软件进行一些GIS的分析应用。
本篇文章比较简单,介绍了如何在ArcGIS和Excel中计算两期土地利用的转移矩阵。主要流程如下
一、土地利用编码分类;
二、在ArcGIS中计算各类土地利用转移面积;
三、在Excel中对处理结果进行整饬
大家可以下载数据,按照以下流程进行操作学习。
【数据下载链接】:GIS点滴公众号界面回复“土地利用示例数据”
本数据只适于学习交流使用,切勿用于商业或者学术研究中。
一、土地利用编码分类
我们下载到了“XX地”2000年和2005年的土地利用数据,部分展示如下:
1、实际原始的数据分类比较多,为了方便演示,选取土地的一级分类,新建landcode字段,使用calculate field(字段计算器)工具,截取出CODE字段的第一个字符即可
2、然后按照数字编号对应的中文名称,新建landname字段,
1---林地
2---草地
3---水域
4---耕地
5---建筑用地
6---未利用地
3、同样使用calculate field(字段计算器)工具,得到每一条记录对应的中文名称
这里写了个条件赋值函数,文本如下(按理说第2步本需不做,但是该场景中涉及到的方法比较有用,适宜练手):
def landTransfer(landcode):
if landcode==1:
land=u"林地"
elif landcode==2:
land=u"草地"
elif landcode==3:
land=u"水域"
elif landcode==4:
land=u"耕地"
elif landcode==5:
land=u"建筑用地"
elif landcode==6:
land=u"未利用地"
return land
完成如下
同上,将另一期也按上述流程进行处理。
二、在ArcGIS中计算各类土地利用转移面积
1、首先将两期土地利用分别使用融合(dissolve)工具合并,参数如下设置:
意合并的是landcode字段
结果如下
2、将dissolve后的两个图层新增一个字段,2000年这一期的增加landcode2000,2005年这一期的增加landcode2005,并将landcode的值赋过来
3、使用intersect工具,添加融合后的图层数据
4、得到结果图层的属性表如下,
5、导出该属性表为txt
三、在Excel中对处理结果进行整饬
1、使用excel或者wps直接打开该dbf文件(或者是打开excel,按Ctrl+O快捷键打开)
注意选择文件类型为txt
2、分隔符勾选逗号
3、然后下一步下一步就打开了
4、点击excel中数据选项卡下的透视表功能
5、将如下三个字段拖入相应的区域
6、然后就到了如下的表格:
第一行和第一列表示的是各类土地利用对应的编号。
这就是土地利用转移的面积矩阵,整饬如下
我们还可以基于此进一步计算土地利用转移的百分比矩阵
近期评论