1985～2025 年各城市核心绿色生境景观破碎化指标面板数据

941 字 约 3 分钟 数据资料 面板数据生态环境

今天给大家分享一份 1985～2025年各城市核心绿色生境景观破碎化指标面板数据。该数据基于 CLCD 土地覆盖栅格数据处理得到，将森林（Forest）、灌丛（Shrub）、草地（Grassland）、水体（Water）、湿地（Wetland）合并定义为绿色生境，并使用 landscapemetrics 包计算 9 个景观破碎化指标。除绿色生境外，该数据还同步计算了耕地和不透水面的相关指标。

数据概览

数据包含 1985 至 2025 年各城市的景观破碎化指标，具体指标如下：

  • AI（聚集度指数）：取值 0～100，衡量生境斑块在空间上的聚集与连通程度
  • DIVISION（分离度指数）：取值 0～1，反映斑块被分割、相互隔离的程度
  • ED（边缘密度）：单位面积内的斑块边缘总长度，值越大边缘效应越强
  • FRAC_MN（平均分维数）：围绕 1 变化，描述斑块几何形状的复杂程度
  • LPI（最大斑块指数）：最大斑块占目标生境总面积的比例
  • NP（斑块数量）：研究区内目标生境的总斑块个数，直接反映破碎化程度
  • PD（斑块密度）：单位面积内的斑块数量
  • PLAND（景观比例）：目标生境面积占整个研究区总面积的百分比
  • SHAPE_MN（平均形状指数）：衡量斑块形状偏离正方形/圆形的程度

数据预览如下：

图表展示

下图展示了 2025 年各城市绿色生境景观分离度：

处理代码

该数据使用 R 语言的 landscapemetrics、terra 和 sf 包处理得到。处理流程如下：

  1. 读取 CLCD 土地覆盖栅格数据和城市行政区划边界
  2. 对每个城市进行栅格裁剪（crop）和掩膜（mask）
  3. 将目标生境类型（Forest/Shrub/Grassland/Water/Wetland）合并重分类
  4. 使用 landscapemetrics 包计算 9 个景观指数
  5. 循环所有年份（1985～2025）并行计算，输出为面板数据

关键代码示例：

library(landscapemetrics)
library(terra)
library(sf)

# 读取栅格数据和城市边界
rst <- rast("CLCD_v01_1985_albert.tif")
city <- read_sf("2021行政区划/市.shp") %>% st_transform(crs(rst))

# 并行计算各城市景观指数
library(parallel)
cl <- makeCluster(5)
parLapply(cl, 1:nrow(city), function(c){
  city %>% slice(c) -> testcity
  rst %>% terra::crop(vect(testcity)) %>%
    terra::mask(vect(testcity)) -> testrst

  # 合并生境类型
  testrst[testrst == 2 | testrst == 3 | testrst == 4 |
          testrst == 5 | testrst == 9] <- 10

  # 计算 9 个景观指数
  bind_rows(
    lsm_c_division(testrst, directions = 8),
    lsm_c_pd(testrst, directions = 8),
    lsm_c_np(testrst, directions = 8),
    lsm_c_ed(testrst, directions = 8),
    lsm_c_pland(testrst, directions = 8),
    lsm_c_lpi(testrst, directions = 8),
    lsm_c_ai(testrst),
    lsm_c_shape_mn(testrst, directions = 8),
    lsm_c_frac_mn(testrst, directions = 8)
  ) %>% mutate(year = str_extract(f, "\\d{4}"))
})

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

数据引用格式

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

RStata 数据中心: 1985～2025年各城市核心绿色生境景观破碎化指标面板数据. 2026. https://tidyfriday.cn/rsdb2/

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

RStata Data Center: Panel Data on Landscape Fragmentation Indices of Core Green Habitat by City, 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/dfd9cc9ba8af439d8b196a390813f50d

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

点击这里跳转到 RStata 短书平台获取附件：1985～2025 年各城市核心绿色生境景观破碎化指标面板数据

相关文章

评论