ArcGIS：土地利用变化模式图绘制，土地利用转移矩阵

我们经常想知道土地利用在不同年份的变化过程，然而ArcGIS里面并没有现成的工具可以使用。因此，我制作了2010年到2020年的土地利用变化模式图。虽然只是两年的变化过程，但是我觉得挺有意思。所以记录下来，以便提供给你们一些解决思路。结果图如下：

每一个标签代表一种变化模式。比如11表示从耕地到耕地的变化，43表示水域到草地的变化。

---

方法一：
利用ArcGIS的 栅格计算器 ，输入公式 2010年的TIFF*10+2020的TIFF 。查看属性表就是土地利用转移矩阵。

---

方法二：
首先，让我们明确一下每一类土地利用类型对应的标签：

土地利用类型	标签
耕地	1
林地	2
草地	3
水域	4
建设用地	5
未利用地	6

---

然后，第一步我们对2010年的数据做重分类，重分类的依据是每一类标签之间的间隔恰好能隔土地利用种类的数值大小（比如本实验采用7类，以防有空值的转换加了一个间隔）。

土地利用类型	标签	重分类标签
耕地	1	10
林地	2	18
草地	3	26
水域	4	34
建设用地	5	42
未利用地	6	50

---

第三步，使用ArcGIS里面的栅格计算器，我们使用相加运算。
将重分类以后的2010年的数据和2020年的数据相加得到栅格数据B。
这样的结果就是他们每个像素的和是唯一的，每一个和对应一种变化模式，这样将求和以后的栅格数据B做一次重分类，对应回来正确的变化模式，就可以得到第一幅变化模式图的结果了。

土地利用类型	2010年重分类以后的标签	2020年标签	求和以后的标签	变化模式对应的标签
耕地	10	1…7	11…17	11…16
林地	18	1…7	19…25	21…26
草地	26	1…7	27…33	31…36
水域	34	1…7	35…41	41…46
建设用地	42	1…7	43…49	51…56
未利用地	50	1…7	51…57	61…66

---

注意：① 本实验使用的土地利用数据是从中科院（宫鹏）网站（http://www.globallandcover.com/defaults.html?src=/Scripts/map/defaults/download.html&head=download&type=data）上下载的。
② 本实验中的地类较少，为6类。所以如果地类比较多的话，可以尝试将前一年数据重分类之间的标签值之间的距离设置大一点。
③ 本实验只考虑了两年的变化模式，你们可以尝试三年，四年。。。的变化模式，当然这种变化模式比较复杂，后期的变化模式种类也比较多。但是这是一个好的思路，可以在时空两个维度来了解土地利用类型变化的过程，对于了解人地关系和土地利用类型的变化，解释它们之间的相互作用关系有很重要的意义。

---

第二个作用，生成土地利用转移矩阵。上面的步骤已经介绍了土地利用变化的过程了，还难做土地利用转移矩阵吗？
打开属性表（count计算的是像素数量，面积需要字段计算器新建一个字段，像素转换面积百度其他博客）。面积单位是平方公里。

土地利用类型转移矩阵（平方公里）	1	2	3	4	5	6
1	83814.5	3648	5168.75	2057.5	12924	299
2	2536.5	35657	5481.75	163.25	716.75	28.5
3	5429.25	5588.75	22380	369.5	1153.25	200.25
4	1594.25	295.75	428.75	2791.5	702.75	466
5	5668.25	235.75	361	1346.75	12316.5	54.5
6	695.5	59	326.5	142.75	143.75	624.25