遍历栅格波段计算NDVI遇报错，寻求非堆叠式实现方案

错误原因 #

• 错误地将list.files()的结果用list()再次包裹，导致red和nir变成了双层嵌套列表（仅含1个元素，该元素才是包含9个文件路径的向量），遍历1:9时取到的nir[[i]]和red[[i]]并非有效文件路径。
• 直接将文件路径字符串传入计算函数，未读取栅格数据。字符串属于非数值类型，无法执行加减乘除运算，触发报错。

修正方案 #

使用terra包（替代维护较少的旧版raster包，高效且功能完善）读取栅格，全程保留元数据，以下是完整修正代码：

# 加载依赖包
library(terra)

# 设置目标文件夹路径
subfolders <- '~/Desktop/hls project/fieldszn'

# 正确获取红波段、近红外波段的文件路径（无需嵌套list()）
red <- list.files(subfolders, pattern = "red", full.names = TRUE, recursive = TRUE)
nir <- list.files(subfolders, pattern = "nir", full.names = TRUE, recursive = TRUE)

# 定义NDVI计算函数：读取栅格+计算+保留元数据
calculate_ndvi <- function(nir_path, red_path){
  # 读取栅格文件，生成带元数据的SpatRaster对象
  nir_rast <- rast(nir_path)
  red_rast <- rast(red_path)
  
  # 执行NDVI计算，自动继承原栅格的投影、分辨率等元数据
  ndvi <- (nir_rast - red_rast) / (nir_rast + red_rast)
  
  return(ndvi)
}

# 遍历计算对应波段的NDVI
ndvi_list <- list()
for (i in seq_along(red)){ # 用seq_along替代固定1:9，适配波段数量变化
  ndvi_list[[i]] <- calculate_ndvi(nir[i], red[i])
}

# 可选：将每个NDVI结果保存为tif文件，保留元数据
for (i in seq_along(ndvi_list)){
  writeRaster(ndvi_list[[i]], 
              filename = paste0("ndvi_output_", i, ".tif"),
              overwrite = TRUE)
}

关键说明 #

• 移除list()对list.files()结果的包裹，让red和nir直接成为包含文件路径的向量，长度与波段数一致。
• 用rast()读取栅格后，生成的SpatRaster对象支持直接数值运算，且完整保留原栅格的投影、坐标系、分辨率等元数据，不会出现堆叠导致的元数据丢失问题。
• 使用seq_along(red)遍历，避免硬编码波段数量，代码灵活性更高。

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