土地利用转移矩阵是马尔科夫模型在土地利用变化方面的应用。本期分享如何用ArcGIS 制作土地利用转移矩阵,GIS版本为10.1。

1

数据准备

本期所用数据为资源环境科学与数据中心下载的2000及2020年贵州省土地利用类型数据,格式为adf,其数据生产制作是以各期Landsat TM/ETM遥感影像为主要数据源,通过人工目视解译生成,土地利用类型包括耕地、林地、草地、水域、居民地和未利用土地6个一级类型以及25个二级类型,具体分类系统如下:

2

数据重分类

(1) 二级类型重分类

a. 打开ArcToolbox工具/Spatial Analyst工具/重分类,依次对两期数据进行重分类;

b. 输入栅格,重分类字段,旧值开头数字是几,新值就改为几,输出栅格,确定:

c. 重分类结果如下:

3

栅格转矢量

(1) 栅格转面:

a. 打开转换工具/由栅格转出/栅格转面,依次对两期数据进行转面;

b. 输入栅格,字段,选择输出路径,勾选掉简化面,确定:

c. 结果如下:

4

同一地类融合

(1) 融合

a. 打开数据管理工具/制图综合/融合,依次对两期数据进行融合;

b. 输入要素,选择输出路径,融合字段选择GRIDCODE字段,其余默认,确定:

c. 结果展示

打开融合后图层的属性表,可以看出,所有面融合为了6类:

(2) 添加字段:

a. 依次在两期融合后的数据里,打开属性表,添加字段,名称为2000/2020类型,文本型,长度默认,确定:

b. 点击编辑,在新字段里填写对应类型,停止编辑,保存编辑内容:

5

相交

(1) 相交工具

a. 打开分析工具/叠加分析/相交工具;

b. 输入要素,选择输出路径,其余默认,确定:

c. 结果如下:

拿标注这行为例,2000年是耕地的,2020年还是耕地:

(2) 添加字段:

为了看出2000至2020年土地利用类型的变化情况,在得到的相交图层属性里,添加字段,命名为变化,类型选择双精度,选择合适的投影坐标系,进行面积计算:

6

土地利用转移矩阵

(1) Excel打开及处理相交表

a. 把相交后的图层dbf文件复制一份,后缀改为xls,即可在Excel里打开:

b. 点击插入/数据透视表,操作如下:

(2) 土地利用转移矩阵:

调整格式及字体等,即可得到如下土地利用转移矩阵:

-----END-----

(请添加下方小助手微信)

来源:生态遥感笔记