查看原文
其他

栅格数据土地覆被变化制图与土地利用转移矩阵

走天涯徐小洋 走天涯徐小洋地理数据科学 2022-05-17

土地覆被变化制图与土地利用转移矩阵(栅格版)

前面讲过如何利用矢量数据进行土地覆被变化制图和土地利用转移矩阵的制作,接下来讲讲栅格数据如何进行变化制图和转移矩阵制作。

土地覆被变化制图

先以京津冀土地覆被变化为例,看看栅格的土地覆被变化如何制图。

1990-2015京津冀土地覆被变化

以耕地的变化为例:

  1. 先从土地覆被类型数据中通过重分类等方法提取出单独的耕地地类。
提取耕地,耕地为1,其它地类0,界外Nodata
  1. 使用栅格计算器(Raster Calcutor)计算,2015-1990:
  • -1:耕地减少
  • 0:不变+其它
  • 1:耕地增加
计算耕地增加减少2015-1990

由于前面0值包含其它地类信息,为单独提取出耕地不变的区域,需要额外进行栅格的处理:

  • 保留耕地,其它NODATA
  • 对两期的耕地分类执行等于(Equal to)
  • 1:耕地不变
模型构建器,计算耕地不变的范围

把两次处理的栅格进行符号化,制图即可获得下面的耕地变化图。

1990-2015耕地变化图

土地利用转移矩阵

使用栅格数据制作土地覆被转移矩阵,原理和矢量一样,就是换了一些栅格处理所用的工具。

  1. 先使用合并(Combine)工具给两个年份的土地覆被栅格合并。
Combine
Combine得到的栅格属性表
  1. 由于栅格属性表中的Count是像元数,要想转变为地类面积,需要根据分辨率做个转换。本实验数据使用的是30米的土地覆被数据,那么一个像元实际面积为30*30=900平方米

  2. 那么实际的面积应为Count900,则新建一个Area字段,Area=Count900,面积单位:平方米。注意,Area要使用Double字段类型,否则可能字段长度不够,导致面积计算失败。

新建一个Area字段,用字段计算器,Count*900
  1. 将计算完的属性表导出,在EXCEL中使用数据透视表做土地利用转移矩阵。行1990土地覆被,列2015土地覆被,求和项Area
数据透视表
  1. 整理一下,就可以做成一个土地覆被转移矩阵表格啦
京津冀1990-2015土地覆被转移矩阵,面积:公顷

相关阅读

点击阅读原文查看视频课程讲解


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存