如何在R中循环遍历长度不同的栅格列表与阈值列表？

问题说明 #

我有两个列表对象：Threshold_List是长度为2的双精度数值列表，Raster_List是长度为10的栅格列表。需要实现：用第一个阈值处理前5个栅格，第二个阈值处理后5个栅格，现有代码仅在两列表长度相等时有效，需修改实现该重复逻辑。

修改后的完整代码 #

library(raster)
# 创建随机栅格列表
r1 <- raster(nrows=10,ncols=10,res = 10, xmn = -100, xmx = 100, ymn = -100, ymx = 100)
Raster_List <- lapply(1:10, function(i) setValues(r1,runif(ncell(r1))))
Raster_names<-c("a","b","c","d","e","f","g","h","i","j")
names(Raster_List)<-Raster_names
rm(r1)

# 创建阈值列表，并扩展为与栅格列长度匹配的重复序列
Threshold_List <- list(0.2, 0.9)
# 每个阈值重复5次，生成长度10的阈值序列
expanded_thresholds <- rep(Threshold_List, each = 5)

# 遍历栅格和对应的阈值
for (i in seq_along(Raster_List)) {
  tif <- Raster_List[[i]]
  thresh <- expanded_thresholds[[i]]
  name <- Raster_names[[i]]
  
  # 赋予坐标系
  crs(tif) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
  
  # 根据阈值掩膜
  tif[tif < thresh] <- NA
  
  # 保存输出
  tif_file <- paste0("Binary_", name)
  writeRaster(tif, tif_file, format="GTiff", overwrite=TRUE)
}

关键修改点 #

• 扩展阈值序列：使用rep(Threshold_List, each = 5)将长度为2的阈值列表扩展为长度10的序列，第一个阈值重复5次，第二个阈值重复5次，完美匹配栅格列的长度。
• 调整循环逻辑：去掉原有的嵌套循环，改为单循环遍历索引，同时取出对应位置的栅格和阈值，实现一一对应的处理逻辑，避免了原代码中每个栅格被所有阈值处理的错误。
• 简化阈值列表创建：直接用list(0.2, 0.9)创建阈值列表，比原代码的data.frame转换更简洁直观。

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