【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

相关文章

全网最全Tomcat完全卸载重装教程小结

《全网最全Tomcat完全卸载重装教程小结》windows系统卸载Tomcat重新通过ZIP方式安装Tomcat,优点是灵活可控,适合开发者自定义配置,手动配置环境变量后,可通过命令行快速启动和管理... 目录一、完全卸载Tomcat1. 停止Tomcat服务2. 通过控制面板卸载3. 手动删除残留文件4.

Python的pandas库基础知识超详细教程

《Python的pandas库基础知识超详细教程》Pandas是Python数据处理核心库,提供Series和DataFrame结构,支持CSV/Excel/SQL等数据源导入及清洗、合并、统计等功能... 目录一、配置环境二、序列和数据表2.1 初始化2.2  获取数值2.3 获取索引2.4 索引取内容2

python依赖管理工具UV的安装和使用教程

《python依赖管理工具UV的安装和使用教程》UV是一个用Rust编写的Python包安装和依赖管理工具,比传统工具(如pip)有着更快、更高效的体验,:本文主要介绍python依赖管理工具UV... 目录前言一、命令安装uv二、手动编译安装2.1在archlinux安装uv的依赖工具2.2从github

C#实现SHP文件读取与地图显示的完整教程

《C#实现SHP文件读取与地图显示的完整教程》在地理信息系统(GIS)开发中,SHP文件是一种常见的矢量数据格式,本文将详细介绍如何使用C#读取SHP文件并实现地图显示功能,包括坐标转换、图形渲染、平... 目录概述功能特点核心代码解析1. 文件读取与初始化2. 坐标转换3. 图形绘制4. 地图交互功能缩放

SpringBoot集成redisson实现延时队列教程

《SpringBoot集成redisson实现延时队列教程》文章介绍了使用Redisson实现延迟队列的完整步骤,包括依赖导入、Redis配置、工具类封装、业务枚举定义、执行器实现、Bean创建、消费... 目录1、先给项目导入Redisson依赖2、配置redis3、创建 RedissonConfig 配

基于C#实现PDF转图片的详细教程

《基于C#实现PDF转图片的详细教程》在数字化办公场景中,PDF文件的可视化处理需求日益增长,本文将围绕Spire.PDFfor.NET这一工具,详解如何通过C#将PDF转换为JPG、PNG等主流图片... 目录引言一、组件部署二、快速入门:PDF 转图片的核心 C# 代码三、分辨率设置 - 清晰度的决定因

Java Scanner类解析与实战教程

《JavaScanner类解析与实战教程》JavaScanner类(java.util包)是文本输入解析工具,支持基本类型和字符串读取,基于Readable接口与正则分隔符实现,适用于控制台、文件输... 目录一、核心设计与工作原理1.底层依赖2.解析机制A.核心逻辑基于分隔符(delimiter)和模式匹

spring AMQP代码生成rabbitmq的exchange and queue教程

《springAMQP代码生成rabbitmq的exchangeandqueue教程》使用SpringAMQP代码直接创建RabbitMQexchange和queue,并确保绑定关系自动成立,简... 目录spring AMQP代码生成rabbitmq的exchange and 编程queue执行结果总结s

python使用Akshare与Streamlit实现股票估值分析教程(图文代码)

《python使用Akshare与Streamlit实现股票估值分析教程(图文代码)》入职测试中的一道题,要求:从Akshare下载某一个股票近十年的财务报表包括,资产负债表,利润表,现金流量表,保存... 目录一、前言二、核心知识点梳理1、Akshare数据获取2、Pandas数据处理3、Matplotl

Python pandas库自学超详细教程

《Pythonpandas库自学超详细教程》文章介绍了Pandas库的基本功能、安装方法及核心操作,涵盖数据导入(CSV/Excel等)、数据结构(Series、DataFrame)、数据清洗、转换... 目录一、什么是Pandas库(1)、Pandas 应用(2)、Pandas 功能(3)、数据结构二、安