1985~2025 年各省市区县不同土地覆盖类型的土地面积

今天给大家分享一份 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年各省份不同土地覆盖类型的土地面积_km2.dta
  • 1985~2025年各城市不同土地覆盖类型的土地面积_km2.dta
  • 1985~2025年各区县不同土地覆盖类型的土地面积_km2.dta

各文件预览如下:

图表展示

下图展示了 2025 年各城市土地覆盖类型分布:

下图展示了历年各省份耕地面积变化趋势:

下图展示了 2025 年各区县不透水区域面积占比:

处理代码

该数据基于 CLCD 土地覆盖栅格数据处理得到,主要使用 R 语言的 terra、sf 等包进行栅格数据裁剪与面积汇总。处理流程如下:

  1. 读取 CLCD 土地覆盖栅格数据和行政区划边界(省、市、县)
  2. 使用 terra 包进行栅格裁剪(crop)和掩膜处理(mask)
  3. 计算各土地类型的斑块面积(cellSize + zonal)
  4. 按行政区划代码合并汇总
  5. 转换为长格式并与行政区划属性表关联
  6. 输出为 .dta 面板数据

关键代码示例:

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

# 读取栅格数据和行政区划边界
rst <- rast("CLCD_v01_2025_albert.tif")
county <- read_sf("县.shp") %>% st_transform(crs(rst))

# 栅格裁剪与面积计算
rsttemp <- terra::crop(rst, countytemp) %>% terra::mask(countytemp)
cellSize(rsttemp, unit = "km") -> newrst
zonal(newrst, rsttemp, fun = sum, na.rm = T) %>%
set_names("class", "area_km2")

# 并行计算加速
library(parallel)
cl <- makeCluster(20)
parLapply(cl, 1:nrow(citydf), function(x){
# ... 逐个行政区划计算
})

数据引用格式

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

RStata 数据中心: 1985~2025年各省市区县不同土地覆盖类型的土地面积. 2025. https://tidyfriday.cn/rsdb2/

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

RStata Data Center: Panel Data on Land Area by Land Cover Type and Province/City/District/County, 1985–2025. 2025. https://tidyfriday.cn/rsdb2/

关联课程/数据推荐

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

栅格数据几何运算:https://rstata.duanshu.com/#/course/b579b775ee4f449283615a06f9cdcebb

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

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

2001~2021 年各省市区县不同土地覆被类型的区域面积:https://rstata.duanshu.com/#/course/9f6a61dfdc0948b1b30b1e272de7e766

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

点击这里跳转到 RStata 短书平台获取附件:1985~2025 年各省市区县不同土地覆盖类型的土地面积

评论