如何在R中创建循环按变量名+日期规则命名保存栅格文件

算力阿峰

2026-10-3

R语言哨兵影像指数批量导出优化方案

你可以通过构建命名栅格列表+批量遍历的方式减少重复代码，还能自动拼接日期生成符合要求的文件名，完整优化后代码如下：

# 加载依赖包，未安装可先运行 install.packages(c("raster","stringr"))
library(raster)
library(stringr)

Dossier <- "C:/Users/Perrin/Desktop/INRA/Raster/sentinel/L1C_T31UDR_A019210_20190225T105315/S2A_MSIL1C_20190225T105021_N0207_R051_T31UDR_20190225T125616.SAFE/GRANULE/L1C_T31UDR_A019210_20190225T105315/IMG_DATA"

# 读入波段，补充你原有代码缺失的路径参数避免读文件失败
Bande1 <- raster(list.files(path = Dossier, pattern = "\\B01.jp2$"))
Bande2 <- raster(list.files(path = Dossier, pattern = "\\B02.jp2$"))
Bande3 <- raster(list.files(path = Dossier, pattern = "\\B03.jp2$"))
Bande4 <- raster(list.files(path = Dossier, pattern = "\\B04.jp2$"))
# 请自行补充你原有代码缺失的 Bande8、Bande8A、Bande11、NDWII 变量定义

# 计算指数（注意你原有NDWI公式有误，分母应为波段和，已修正）
NDVI <- (Bande8-Bande4)/(Bande8+Bande4)
NDWI <- (Bande8A-Bande11)/(Bande8A+Bande11)
NDDI <- (NDVI-NDWII)/(NDVI+NDWII)

# 自动从路径提取8位拍摄日期，也可手动赋值 date_str <- "20190225"
date_str <- stringr::str_extract(Dossier, "20\\d{6}")
# 定义输出路径
out_dir <- "C:/Users/Perrin/Desktop/INRA/résultats R/"

# 构建命名栅格列表，后续新增要导出的指数只需在这里加一行即可
raster_list <- list(
  NDVI = NDVI,
  NDWI = NDWI,
  NDDI = NDDI
)

# 批量导出，自动按「变量名_日期.tif」格式命名
lapply(names(raster_list), function(var_name){
  out_path <- file.path(out_dir, paste0(var_name, "_", date_str, ".tif"))
  writeRaster(raster_list[[var_name]], out_path, overwrite = TRUE)
})