史上最简单的土地利用转移矩阵分析（附练习数据下载）
土地利用/覆被变化反映了土地利用的空间格局和各地类要素之间的空间动态变化，对于了解区域生态系统服务功能的整体状况，针对不同区域实施生态系统服务功能保护具有重要意义。
土地利用转移矩阵即表达不同土地利用在时间和空间上变化情况，作为研究土地利用变化的常用方法之一，亦作为CA-MarKov模型预测分析的前期基础。
本文就以陕西省延安市2000年和2020年监督分类得到的土地利用为基础进行土地利用转移矩阵的方法介绍。
数据下载
我们自己准备了一份练习数据，有需要的在文末有下载链接。
也可从下面这些网站上申请下载。
01.资源环境科学与数据中心
https://www.resdc.cn/Default.aspx
选择中国土地利用遥感监测数据即可
软件具体操作步骤如下：
加载数据
首先打开ArcMap10.5，添加数据延安市2000年和2020年土地利用基础数据。
右键打开该数据的属性表可以知道，每一类土地利用的代码，即栅格的Value值，1代表建设用地，以此类推。AREA表示每一类土地利用的面积，可以用Count*30*30/1000000计算得到，单位为平方千米。lulc表示每一类土地利用类型的名称。
注*：Count表示栅格数，因为该栅格的空间分辨率为30米，即一个栅格的面积为30*30平方米，有n个栅格，由此计算得到，公式后除以一百万是为了将平方米转换为平方千米。
土地利用转移矩阵
土地利用转移矩阵实质上是利用ArcGIS空间叠加分析的思想，即栅格计算器。打开空间分析工具——地图代数——栅格计算器。
在地图代数表达式中输入“LULC2000*10+LULC2020”，它表示将2000年的土地利用栅格属性Value值*10+2020年的土地利用类型栅格值，计算得到的新栅格数据表示2000年至2020年延安市土地利用的空间转移变化。
计算得到的新图层rastercalc12，其Value值变化从11~110。比如，2000年建设用地Value值为1，将其*10，加上2020年的每一类土地利用类型的Value值。即原则上存在10类变化情况，即建设用地转变为建设用地；建设用地转变为林地；建设用地转变为耕地等，以此类推。其面积计算方法如上所述，最终将该属性表导入至Excel进行统计即可得到延安市近20年土地利用类型变化的转移矩阵。
注*：利用栅格计算器进行叠加运算时，需注意进行栅格捕捉。
栅格捕捉的目的是防止栅格之间由于坐标系等多种因素发生的栅格错位，从而导致计算结果误差过大的情况。
设置栅格捕捉的具体方法如下：
进入本公众号后台
发送数字1535获取下载链接
- END -