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

2024-06-07 07:44

本文主要是介绍【GIS教程】土地利用转移矩阵,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

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

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

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

一.概念

二.含义

三.数据准备

四.两期转移矩阵制作

五.三期转移矩阵制作

六.参考文献

一、概念

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

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

二、含义

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

图片

转移矩阵示意图1

其中,TL表示上一时相,TN表示下一时相。C1Cn表示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)人员,限定为关注小编的用户。

(2)各类项目(包括各类科研项目)申请本数据扔享受免费政策,但需向本号捐赠一定数量的硬盘才能获取。

(3)捐赠硬盘可免留言获取数据。

这篇关于【GIS教程】土地利用转移矩阵的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1038560

相关文章

Linux搭建Mysql主从同步的教程

《Linux搭建Mysql主从同步的教程》:本文主要介绍Linux搭建Mysql主从同步的教程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux搭建mysql主从同步1.启动mysql服务2.修改Mysql主库配置文件/etc/my.cnf3.重启主库my

SpringBoot操作MaxComputer方式(保姆级教程)

《SpringBoot操作MaxComputer方式(保姆级教程)》:本文主要介绍SpringBoot操作MaxComputer方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的... 目录引言uqNqjoe一、引入依赖二、配置文件 application.properties(信息用自己

Tomcat的下载安装与使用教程

《Tomcat的下载安装与使用教程》本文介绍了Tomcat的下载、安装和使用方法,包括在本机和云服务器上部署Tomcat的过程,以及解决启动失败问题的方法... 目录Tomcat的下载安装与使用Tomcat的下载与安装Tomcat在本机运行使用Tomcat在php云服务器上的使用总结Tomcat的下载安装与

SpringBoot基于沙箱环境实现支付宝支付教程

《SpringBoot基于沙箱环境实现支付宝支付教程》本文介绍了如何使用支付宝沙箱环境进行开发测试,包括沙箱环境的介绍、准备步骤、在SpringBoot项目中结合支付宝沙箱进行支付接口的实现与测试... 目录一、支付宝沙箱环境介绍二、沙箱环境准备2.1 注册入驻支付宝开放平台2.2 配置沙箱环境2.3 沙箱

IDEA接入Deepseek的图文教程

《IDEA接入Deepseek的图文教程》在本篇文章中,我们将详细介绍如何在JetBrainsIDEA中使用Continue插件接入DeepSeek,让你的AI编程助手更智能,提高开发效率,感兴趣的小... 目录一、前置准备二、安装 Continue 插件三、配置 Continue 连接 DeepSeek四

jdk21下载、安装详细教程(Windows、Linux、macOS)

《jdk21下载、安装详细教程(Windows、Linux、macOS)》本文介绍了OpenJDK21的下载地址和安装步骤,包括Windows、Linux和macOS平台,下载后解压并设置环境变量,最... 目录1、官网2、下载openjdk3、安装4、验证1、官网官网地址:OpenJDK下载地址:Ar

Java8需要知道的4个函数式接口简单教程

《Java8需要知道的4个函数式接口简单教程》:本文主要介绍Java8中引入的函数式接口,包括Consumer、Supplier、Predicate和Function,以及它们的用法和特点,文中... 目录什么是函数是接口?Consumer接口定义核心特点注意事项常见用法1.基本用法2.结合andThen链

C#集成DeepSeek模型实现AI私有化的流程步骤(本地部署与API调用教程)

《C#集成DeepSeek模型实现AI私有化的流程步骤(本地部署与API调用教程)》本文主要介绍了C#集成DeepSeek模型实现AI私有化的方法,包括搭建基础环境,如安装Ollama和下载DeepS... 目录前言搭建基础环境1、安装 Ollama2、下载 DeepSeek R1 模型客户端 ChatBo

JAVA集成本地部署的DeepSeek的图文教程

《JAVA集成本地部署的DeepSeek的图文教程》本文主要介绍了JAVA集成本地部署的DeepSeek的图文教程,包含配置环境变量及下载DeepSeek-R1模型并启动,具有一定的参考价值,感兴趣的... 目录一、下载部署DeepSeek1.下载ollama2.下载DeepSeek-R1模型并启动 二、J

Vue项目的甘特图组件之dhtmlx-gantt使用教程和实现效果展示(推荐)

《Vue项目的甘特图组件之dhtmlx-gantt使用教程和实现效果展示(推荐)》文章介绍了如何使用dhtmlx-gantt组件来实现公司的甘特图需求,并提供了一个简单的Vue组件示例,文章还分享了一... 目录一、首先 npm 安装插件二、创建一个vue组件三、业务页面内 引用自定义组件:四、dhtmlx