1985～2025 年各乡镇各类型土地覆盖面积面板数据

809 字 约 3 分钟 数据资料 面板数据土地

今天给大家分享一份 1985～2025年各乡镇各类型土地覆盖面积面板数据。原始数据为栅格数据：「The 30m annual land cover datasets and its dynamics in China from 1985 to 2025」，下载链接为：https://zenodo.org/records/15853565。其中，土地覆被类型包括 9 种，具体为：Barren（裸地）、Cropland（耕地）、Forest（森林）、Grassland（草地）、Impervious（不透水面）、Shrub（灌丛）、Snow or Ice（冰雪）、Water（水体）、Wetland（湿地）。

由于这份数据处理量非常大，这里采用逐个乡镇裁剪汇总的方式，使用并行计算加速处理。

数据概览

经过处理即可得到下面的文件：

  • 1985～2025年各乡镇各类型土地覆盖面积面板数据.dta
  • 乡镇代码匹配省市区县.dta（用于匹配乡镇与省市区县信息）

数据包含了乡镇名称和行政区划代码、各种覆被类型面积（km²）。数据预览如下：

图表展示

下图展示了 2025 年各乡镇耕地面积（km²）：

处理代码

该数据使用 R 语言的 terra 和 sf 包进行处理，采用逐乡镇裁剪汇总的方式。处理流程如下：

  1. 读取 CLCD 土地覆盖栅格数据和乡镇行政区划边界
  2. 对每个乡镇进行栅格裁剪（crop）和掩膜（mask）
  3. 使用 cellSize + zonal 统计每种土地覆盖类型的面积
  4. 并行计算（16 个核心）加速，循环全部 41 年数据
  5. 与历史数据合并输出为 .dta 面板数据

关键代码示例：

library(tidyverse)
library(sf)
library(terra)

# 读取乡镇边界
town <- read_rds("town.rds") %>%
  rename(乡镇 = Name, 乡镇代码 = code) %>%
  st_transform(crs(rst))

# 并行计算各乡镇各类型面积
library(parallel)
cl <- makeCluster(16)
parLapply(cl, 1:nrow(towndf), function(x){
  rast(towndf$rstname[x]) -> rst
  town %>% filter(乡镇代码 == towndf$towncode[x]) %>% vect() -> towntemp

  rst %>% terra::crop(towntemp) %>% terra::mask(towntemp) -> rsttemp
  cellSize(rsttemp, unit = "km") -> newrst
  zonal(newrst, rsttemp, fun = sum, na.rm = T) %>%
    set_names("class", "area_km2") %>%
    mutate(year = str_extract(towndf$rstname[x], "\\d{4}"),
           code = towntemp$乡镇代码) %>%
    write_rds(paste0("res/town_", x, ".rds"))
})

附件中也提供了该数据的处理代码供参考：

数据引用格式

由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：

RStata 数据中心: 1985～2025年各乡镇各类型土地覆盖面积面板数据. 2026. https://tidyfriday.cn/rsdb2/

英文文献可以使用下面的格式引用：

RStata Data Center: Panel Data on Land Cover Area by Township and Land Cover Type, 1985–2025. 2026. https://tidyfriday.cn/rsdb2/

关联课程/数据推荐

栅格数据的裁剪与面积汇总——基于 R 语言的方法：https://rstata.duanshu.com/#/course/d2c2710928f44ad7a7518b8a178562c0

RStata 定制｜栅格数据转面板数据或时间序列数据：https://rstata.duanshu.com/#/course/839afebb01d54457b42b94e87e7dc085

1985～2024 年各乡镇各类型土地覆盖面积面板数据：https://rstata.duanshu.com/#/course/890e946ea9454a66ad5215516588cc94

1985～2025年各省市区县不同土地覆盖类型的土地面积：https://rstata.duanshu.com/#/course/dfd9cc9ba8af439d8b196a390813f50d

如果有相关需要可以联系李老师付费定制。

点击这里跳转到 RStata 短书平台获取附件：1985～2025 年各乡镇各类型土地覆盖面积面板数据

相关文章

评论