ArcGIS 计算土地利用转移矩阵

“All roads lead to Rome”

本文以中国2000年和2010年的欧空局300米分辨率LUCC数据为例，运用 ArcGIS + Excel 获取栅格影响的土地利用转移矩阵。

分三步(如果已分类完成，则可直接跳过第一步；如果仅想得到转移矩阵结果，可以跳过第二步，直接利用Tabulate Area工具，依次添加初始年份、终止年份信息，导出的结果中，左侧为转出结果，即初始年份情况，右侧为转入结果，为变化情况)。

1. 重分类土地利用类型。

2. 一个像元的属性值表征多期土地利用类型信息。

3. 结果导出。

---

1. 重分类土地利用类型

首先加载两期栅格数据：

 

 

 

根据需要，将若干类别，划分为几大类别：

在Excel中新建表格，将原地类编码和重分类后编码一一对应，这里分被设置属性名为【Value】和【RecCode】,保存为【.csv】或者【.txt】后加载进当前工作空间。

 

 

利用 【Reclass by Table】工具进行重分类

 

 

 

 得到处理后的影像。

---

 

2. 利用栅格计算器合并多时段属性信息。

 此处利用栅格计算器合并属性信息。

 思路是：单个像元从2000年的耕地【Value=1】转变为2010年的建设用地【Value=5】，

则可以用Value₂₀₀₀*10 +Value₂₀₁₀ 得到新值为15，（当然，地类数较多时乘100也是没问题的）

即公式为【"Reclass2000.tif"*10+"Reclass2010.tif"】

打开栅格计算器（可以批处理）

 

 

 得到结果的属性即包含了土地利用的变化情况。

打开属性表，所见如下：

 

 

 第二步即完成。

---

 

3.结果导出。

打开【Tabulate Area】面积制表工具，批处理，得到属性表。

 

 

 

 

 

 可以看出，有187.891651的耕地未发生变化，1.546543的1（耕地）转变为5（建设用地）

 

 

 然后用【Table to Excel】导出至EXCEL表中即可。

由于导出结果是单列数据，并且其中可能有部分地类转移未出现，需要首先检核，并手动补齐空行。（如上图仅出现55，需要手动补齐51、52、53、54、56、57）。

在此之后可以在Excel中利用【数据透视表】或者Offset函数等调整格式，转为转移矩阵，这里不做具体讲解。

 

 

---

 

希望这篇文章能对你有所帮助~