在一些土地利用变化的分析展示中,土地利用转移矩阵是最直观、常见的。之前的文章中我们分享了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、然后就到了如下的表格:

第一行和第一列表示的是各类土地利用对应的编号。

这就是土地利用转移的面积矩阵,整饬如下

我们还可以基于此进一步计算土地利用转移的百分比矩阵

有奖
同学们可以将转移面积和百分比的矩阵结果截图,在GIS点滴公众号回复给我们,我们会随机抽取一名送出GIS百宝盘哦