R语言landscapemetrics包批量迭代处理多栅格计算景观指标

所有栅格分辨率、CRS一致仅范围不同的场景下，不需要做额外栅格预处理，直接按以下步骤实现批量计算即可，不需要重复修改单文件代码。

前置准备 #

把208个待处理的土地覆盖栅格统一放到一个独立空文件夹中，不要混入其他无关文件，避免读取错误。

步骤1：加载依赖包 #

# 读写栅格用terra，处理速度更快，landscapemetrics原生支持terra对象
library(terra)
library(landscapemetrics)

步骤2：批量获取栅格文件路径 #

# 替换为你本地存放栅格的文件夹绝对路径
raster_dir <- "D:/your_workspace/landcover_rasters"

# 批量匹配路径，pattern参数要匹配你自己的栅格后缀，比如.tif文件就写"\\.tif$"，.img文件就写"\\.img$"
raster_paths <- list.files(
  path = raster_dir,
  pattern = "\\.tif$",
  full.names = TRUE,
  recursive = FALSE
)

# 运行length(raster_paths)检查读取到的文件数量，确认是208个再往下跑

实现方法1：for循环写法（适配已掌握的循环语法） #

提前创建空对象存储结果，避免循环中动态赋值降低运行效率，每个结果追加来源文件名标识，防止后续混淆指标对应的栅格。

# 初始化空结果表
all_metrics <- data.frame()

for (i in seq_along(raster_paths)) {
  # 读取当前顺序的栅格
  current_rast <- rast(raster_paths[i])
  # 提取当前栅格的文件名作为标识
  current_file <- basename(raster_paths[i])
  
  # 替换为你自己需要计算的景观指标即可，和你单文件运行的计算逻辑完全一致
  current_res <- calculate_lsm(
    landscape = current_rast,
    # what参数填你需要的指标名，不指定就默认计算所有指标，指定后运行速度更快
    what = c("lsm_l_pd", "lsm_l_ai", "lsm_c_pland"), 
    verbose = FALSE # 关闭逐行运行提示，避免批量计算时控制台刷屏
  )
  
  # 给当前结果加来源列
  current_res$source_raster <- current_file
  # 将当前结果追加到总结果表
  all_metrics <- rbind(all_metrics, current_res)
  
  # 打印运行进度
  cat("进度：", i, "/", length(raster_paths), " 已完成文件：", current_file, "\n")
}

# 计算完成后可直接导出为csv存本地
write.csv(all_metrics, "landscape_metrics_batch_result.csv", row.names = FALSE, fileEncoding = "UTF-8")

实现方法2：lapply简化写法 #

不需要手动写rbind追加逻辑，代码更简洁，适合对运行稳定性要求更高的场景：

metrics_list <- lapply(seq_along(raster_paths), function(i){
  current_rast <- rast(raster_paths[i])
  current_file <- basename(raster_paths[i])
  
  current_res <- calculate_lsm(
    landscape = current_rast,
    what = c("lsm_l_pd", "lsm_l_ai", "lsm_c_pland"), # 替换为你需要的指标
    verbose = FALSE
  )
  
  current_res$source_raster <- current_file
  cat("进度：", i, "/", length(raster_paths), " 已完成文件：", current_file, "\n")
  return(current_res)
})

# 一次性合并所有栅格的计算结果
all_metrics <- do.call(rbind, metrics_list)

实用优化技巧 #

• 容错处理：如果担心个别栅格损坏导致整个批量任务中断，可以把计算部分用tryCatch包裹，出错时自动跳过并提示，不会终止整体运行：

current_res <- tryCatch({
  calculate_lsm(current_rast, what = c("lsm_l_pd", "lsm_l_ai", "lsm_c_pland"), verbose = FALSE)
}, error = function(e){
  cat("文件", current_file, "计算出错：", e$message, "\n")
  return(NULL)
})

• 结果校验：全部计算完成后，运行table(all_metrics$source_raster)可以检查每个栅格对应的指标行数是否一致，快速排查漏算问题。
• 只要所有栅格的土地覆盖分类编码规则统一，不同空间范围的栅格计算结果可以直接对比，不需要做范围裁剪、拼接等预处理。

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