【GIS教程】土地利用转移矩阵

随着科技社会的不断进步，人类活动对地理环境的影响与塑造日益明显，土地不断的侵蚀与改变也导致一系列的环境问题日益突出。土地利用/覆盖（LUCC）作为全球环境变化研究的重点问题为越来越多的国际研究机构所重视，研究它的变化既是对已有的工业化、城市化过程的一个较为直观的反映，也是对于未来土地变化进行预测和模拟的一个不可或缺的基础。对于一个地区土地利用变化的研究，有利于对该地区进行更好的生态环境综合评价及合理的国土空间规划。

而在土地利用变化的研究中经常会用到土地利用转移矩阵，用以有效表达俩个不同时相地类之间面积的转换情况。那么什么是土地利用转移矩阵，土地利用转移矩阵又如何制作呢？让我们一起看看吧！

示例数据下载链接：数据下载链接

本文主要分为以下几个部分：

一.概念

二.含义

三.数据准备

四.两期转移矩阵制作

五.三期转移矩阵制作

六.参考文献

一、概念

土地利用转移矩阵是马尔科夫模型在土地利用变化方面的应用。马尔科夫模型不仅可以定量地表明不同土地利用类型之间的转化情况，还可以揭示不同土地利用类型间的转移速率。

土地利用转移矩阵来源于系统分析中对系统状态与状态转移的定量描述。通常的土地利用转移矩阵中，行表示T1时点土地利用类型，列表示T2时点土地利用类型。

二、含义

土地利用转移矩阵以矩阵的形式将俩个不同时期的土地覆盖类型之间互相转换的数量关系展现出来，可以全面的体现出一个区域土地覆盖类型的数值和转移方向。一般形式如下(下方解释说明可忽略)：

转移矩阵示意图1

其中，TL表示上一时相，TN表示下一时相。C1至Cn表示n种不同的土地覆盖类型。假设行Cn=Ci（上一时相土地类型），列Cn=Cj（下一时相土地类型）。Sij表示上一时相Ci转变为下一时相Cj的面积量，Si*表示上一时相Ci的土地覆盖类型面积的总和，S*j表示下一时相Cj的土地覆盖类型面积的总和。Si*-Sii为Ci土地类型的流出量，即上一时相Ci类型土地中转移为下一时相其它土地类型的面积总和。S*j-Sjj表示Cj土地类型的流入量，即下一时相Cj类型土地中由上一时相其它类型土地转变而来的面积总和。

上面的Sij也可以统一替换为Pij，用来展示转移量占整体的比例，能够体现整体与部分的关系，而避免因原某一土地类型总面积较小而得到转移量也较小，从而面积不能很好地表现此土地类型变化程度的情况。

转移矩阵示意图2

三、数据准备

土地利用数据

全国省市行政边界

数据来源引用：地理遥感生态网www.gisrs.cn

练习数据也在文末提供给了大家。

四、两期转移矩阵制作

我们先用2000年和2020年土地覆盖数据制作土地利用转移矩阵。

1.先将土地覆盖数据Tiff图幅按研究所需区域边界裁剪出来，此过程较为简单。此过程中，由于行政边界数据为GCS2000的地理坐标系，global30土地覆盖数据为WGS1984的UTM投影坐标系，我们将边界数据也转换为WGS1984的UTM投影。

2.将地类代码excel表格保存为97-2003版本。

3.栅格数据创建属性表后，在2000年土地覆盖数据属性表中添加文本类型字段“2000”，此处填写为具体年份而不是只写地类，以方便制作转移矩阵时区分时相。

4.由于土地覆盖栅格数据已进行了分类分析，将该数据创建属性表后，其中的VALUE字段代表的即是定义土地类型所用的像元值。我们将土地覆盖数据的图层属性表基于value字段和excel表格里的代码列字段连接。

5.打开字段计算器，将excel表中地类字段值赋予属性表的“2000”字段。

6.将excel连接移除，同理将2020年份的数据进行连接赋值。此外，此操作也可直接通过在excel里更改.dbf文件达到更改属性表的目的。

7.通过工具“栅格转面”将tiff格式数据转为shp矢量。字段选择“2000”，取消简化面勾选，同理处理2020年份数据。

8.通过“融合”工具，将每个年份相同地类属性要素融合为一个要素，方便计算。

此时，属性表如下。

9.将2000年和2020年的融合后数据进行相交处理。

10.在相交后的属性表中添加“面积”字段，右键字段名称选择计算几何计算出面积，单位选择平方公顷。

11.选择“表转excel”工具，转出属性表。

12.打开保存的.xls文件，选择插入-数据透视表，以2000为行，2020为列，面积求和项为值，得到下表，此即为初步的转移矩阵表。

13.将空值以0填充，代表没有发生此俩种地类的转移。修改表头，行列和总计项，最终得到转移矩阵如下，单位为平方公顷。

2000-2020年上海市土地利用转移矩阵

五、三期转移矩阵制作

现大部分涉及多时相的土地利用变化研究时，制作的转移矩阵一般按照时间节点分开绘制。所以涉及多个时相如三个时相时，我们需要绘制两个分开的转移矩阵和一个总的转移矩阵图。

如我们在上一部分所做的2000-2020年土地转移矩阵即为总矩阵图，还应增加2000-2010年和2010年-2020年的转移矩阵图，并将三图并列如下。

2000-2010年上海市土地利用转移矩阵

2010-2020年上海市土地利用转移矩阵

2000-2020年上海市土地利用转移矩阵

除了运用土地利用转移矩阵之外，结合计算各时段土地利用覆盖比例、土地利用变化幅度、土地利用动态度等信息，可以帮助我们分析研究区域时间段内土地利用变化状况，为国土空间规划双评估以及其它相关专题提供更多的土地利用变化决策参考。

值得一提的是，土地利用转移矩阵只是土地利用变化分析众多方法中的一种，是一种纯数字的表达，不够直观。而且只能展现出时间上转移量变化的数值关系，通常需要结合其它方法和图表用以展现各种不同维度的数量关系，或者来表现空间上的变化。

六、参考文献

[1]王徐凡. 土地覆被时空变化的可视化表达研究与实现[D].中国地质大学(北京),2019.

[2]鲁亚楠,姚顺波,邓元杰,丁振民,侯孟阳,郑雪,李雅男. 陕北地区土地利用及景观格局变化对生态服务价值的影响——基于退耕还林(草)背景[J]. 中国农业资源与区划,2019,40(11):180-192.

[3]潘登. 基于复杂网络的苏州市土地利用/覆被变化研究[D].浙江师范大学,2014.

[4]苗立志. 基于GIS的阿克苏地区土地利用趋势分析与研究[D].山东科技大学,2005.

[5]冯文靖,韦燕飞,童新华. 基于Landsat影像的贵港市港北区土地利用分类及变化分析[J]. 广东土地科学,2019,18(05):22-28.

未完待续。。。

数据来源引用：地理遥感生态网www.gisrs.cn

同时，地理遥感生态网www.gisrs.cn上分享了很多地理遥感领域的科学数据（土地利用数据、npp净初级生产力数据数据、NDVI数据、气象数据（降雨量、气温、蒸散量、辐射、湿度、日照时数、风速、水汽压数据）、径流量数据、夜间灯光数据、统计年鉴、道路网、POI兴趣点数据、GDP分布、人口密度分布、三级流域矢量边界、地质灾害分布数据、土壤类型、土壤质地、土壤有机质、土壤PH值、土壤质地、土壤侵蚀、植被类型、自然保护区分布、建筑轮廓分布等等地理数据，以及关于gis、遥感从方面的操作教程）。

【GIS教程】土地利用转移矩阵相关推荐

ArcGIS中实现土地利用转移矩阵
土地利用转移矩阵,就是根据同一地区不同时相的土地覆盖现状的变化关系,求得一个二维矩阵.通过对得到的转移矩阵进行分析,能够得到2个时相,不同的地类之间相互转化的情况,它描述了不同的土地利用的类型在不同年 ...
波段合成，去除黑边并制作土地利用转移矩阵（附练习数据下载）
这次呢,我们来做一做第七届全国大学生GIS技能大赛试题A(上午),因为这次的大赛其实有蛮多套试题,所以我们可能会多写一点类似的文章哦~~试题A告诉我们:土地利用每年都在发生变化,在统计一些大区域的土地 ...
土地利用转移矩阵图怎么做_土地利用转移矩阵生成的几种方法
土地利用转移矩阵生成的几种方法土地利用类型转换矩阵为相同区域在不同时间段内的土地利用类型的相互转换关系,一般用二维表来表达,从二维表中可以快速查看各个地类间相互转化的具体情况.比如某一类别的土地有 ...
转载《Python与开源GIS教程》随书源码网址
卜坤(著)<Python与开源GIS教程>的内容 https://www.osgeo.cn/pygis/ 说明:随书源代码托管在GitHub: https://github.com/buk ...
python与开源gis_转载《Python与开源GIS教程》随书源码网址
卜坤(著)<Python与开源GIS教程>的内容说明:随书源代码托管在GitHub: 正在撰写中,所以会有不完善之处,并会有大量的调整. aptitude install python3 ...
arcgis做土地利用转移矩阵（直接用栅格数据）
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档文章目录本次实验主要使用arcgis的combine工具如果你已经做好了两张土地利用类型的数据,不需要将栅格转为矢量数据→使用相交 ...
土地利用转移矩阵、土壤侵蚀模数变化等此类用地类型变化的Arcgis操作步骤
土地利用转移矩阵.土壤侵蚀模数变化等此类用地类型变化的Arcgis操作步骤总的步骤:两期栅格转面-融合-相交 1.两期土地利用栅格或矢量图 2.栅格转面有时会出错,是因为栅格数据不是整型或字符串 ...
土地利用转移矩阵_ArcGIS_土地利用转移矩阵（对比2010年-2020年土地利用变化）...
1.数据准备 30米全球地表覆盖数据GlobeLand30.http://www.globallandcover.com/ 2.具体操作 ArcGIS中--打开数据--查看字段--添加字段--根据分类 ...
[GIS教程] 7 空间数据查询与空间度量
在知识传播途中,向涉及到的相关著作权人谨致谢意! 文章目录 1 空间查询 1.1 空间属性查询 1.2 空间图形查询 1.3 空间关系查询 1.3.1 拓扑关系查询 1.3.2 缓冲区查询 2 空间度 ...
GIS教程之将栅格数据 raster data发布到 Web 的 3 个简单步骤
本教程向您展示如何创建您自己的栅格文件的 web 地图.为简单起见,本教程将重点关注栅格数据,但类似的方法也适用于矢量数据. 创建 web 地图变得相当容易.例如,使用传单和 Flask,您只需几行代 ...

【GIS教程】土地利用转移矩阵

【GIS教程】土地利用转移矩阵相关推荐

最新文章

热门文章