在R中如何将栅格列表拆分为单独栅格并导出为GeoTiff

是否需要拆分栅格列表？ #

不需要单独拆分该栅格列表，直接遍历列表内的每个RasterLayer元素即可完成批量导出。

具体操作步骤 #

1. 前置准备 #

首先确保已安装并加载raster包，同时将unlist后的结果存储为列表对象方便调用：

# 未安装依赖包时先执行安装
install.packages("raster")
library(raster)

# 存储unlist后的栅格列表
ras_list <- unlist(MCD12Q1.wgs84.ras.01L)

2. 批量导出GeoTiff #

可以选择两种常用实现方式：

方式1：for循环实现（逻辑清晰易调整）

# 循环遍历列表内每个栅格导出
for (i in seq_along(ras_list)) {
  # 自定义导出文件名，按序号区分不同图层
  out_path <- paste0("MCD12Q1_layer_", i, ".tif")
  # 导出为GeoTiff格式
  writeRaster(
    x = ras_list[[i]],
    filename = out_path,
    format = "GTiff",
    overwrite = TRUE, # 已存在同名文件时覆盖，不需要可删除
    options = "COMPRESS=LZW" # 开启LZW压缩减少文件体积，可选参数
  )
}

方式2：lapply实现（代码更简洁）

lapply(seq_along(ras_list), function(i) {
  writeRaster(
    x = ras_list[[i]],
    filename = paste0("MCD12Q1_layer_", i, ".tif"),
    format = "GTiff",
    overwrite = TRUE,
    options = "COMPRESS=LZW"
  )
})

注意事项 #

• 若需要导出到指定文件夹，可在文件名前加上文件夹路径，例如paste0("./output/MCD12Q1_layer_", i, ".tif")，需提前创建对应文件夹
• 如果图层对应特定标识（如年份、产品序号），可替换文件名中的序号为对应标识，方便后续识别

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