如何在R语言中将expanse函数结果合并到单个数据框并补充信息

可以借助purrr包的批量处理函数，结合自定义逻辑高效完成需求，以下是具体实现步骤：

1. 整理栅格为命名列表 #

先把所有需要计算的栅格整理成命名列表，列表名称可作为后续识别各栅格的标识：

library(raster)
library(purrr)
library(dplyr)

# 示例：将你的栅格对象放入命名列表
raster_list <- list(
  "耕地_2020" = r1,
  "林地_2010" = r2,
  "建设用地_2015" = r3
)

2. 自定义处理函数 #

写一个函数，同时完成面积计算和元数据提取（可根据需求调整要添加的信息列）：

calc_area_with_meta <- function(rast) {
  # 计算面积
  area_result <- expanse(rast, unit = "ha")
  # 添加栅格额外信息：示例添加投影、分辨率
  area_result %>%
    mutate(
      投影 = crs(rast, asText = TRUE),
      分辨率_X = res(rast)[1],
      分辨率_Y = res(rast)[2]
    )
}

3. 批量处理并合并结果 #

用map_dfr()批量处理列表中的每个栅格，自动合并为单数据框，同时保留列表名称作为栅格列：

final_result <- map_dfr(raster_list, calc_area_with_meta, .id = "栅格名称")

# 查看最终结果
print(final_result)

简化版（仅添加栅格名称） #

如果只需要添加栅格识别名称，不需要额外元数据，可直接简化：

final_simple <- map_dfr(raster_list, ~expanse(.x, unit = "ha"), .id = "栅格名称")

结果说明 #

最终的数据框会包含：

• 原expanse输出的layer和area列
• 自定义添加的栅格名称、投影、分辨率_X/Y等信息列
• 若某栅格为多层栅格，layer列会自动对应各层面积，同时关联所属栅格的元数据

内容的提问来源于stack exchange，提问作者Mike