使用deparse/substitute转换栅格变量名时循环报错，求解决方案

错误原因 #

你遇到的报错是因为在函数内部使用deparse(substitute())时，substitute()并没有返回变量名，而是返回了变量对应的实际创建表达式（比如stack(list(r, r, r))的完整解析内容）。当你用parse()去解析这个长表达式时，因为内容包含换行和复杂对象结构，导致解析失败。

方法1：直接构建对象列表（最简洁高效） #

这是R中处理批量对象的标准方式，不需要把变量名转字符串，直接将创建好的栅格对象放进列表：

test <- function(resolution){
  
  ## Create RasterStack objects
  r <- raster(xmn= 0, ymn= 0, xmx = 100000, ymx = 100000, resolution = resolution, crs = '+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ')
  r <- setValues(r, sample(x = 0:1, size = ncell(r), replace = T))
  
  # 直接创建带命名的栅格对象列表
  raster_list <- list(
    environmental_layers_scenario_1 = stack(list(r, r, r)),
    environmental_layers_scenario_2 = stack(list(r, r, r)),
    environmental_layers_scenario_3 = stack(list(r, r, r))
  )
  
  return(raster_list)
}

调用后可直接遍历列表处理对象：

layer_list <- test(1000)
# 遍历示例
for (raster_obj in layer_list) {
  print(nlayers(raster_obj)) # 输出每个栅格栈的图层数量
}

方法2：从变量名字符串获取对象（适合已有变量的场景） #

如果已经创建了单独的变量，想要通过变量名字符串获取对象，用mget()替代parse()+eval()，避免解析错误：

test <- function(resolution){
  
  ## Create RasterStack objects
  r <- raster(xmn= 0, ymn= 0, xmx = 100000, ymx = 100000, resolution = resolution, crs = '+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ')
  r <- setValues(r, sample(x = 0:1, size = ncell(r), replace = T))
  
  environmental_layers_scenario_1 <- stack(list(r, r, r))
  environmental_layers_scenario_2 <- stack(list(r, r, r))
  environmental_layers_scenario_3 <- stack(list(r, r, r))
  
  # 直接定义变量名字符串
  layer_names <- c("environmental_layers_scenario_1",
                   "environmental_layers_scenario_2",
                   "environmental_layers_scenario_3")
  
  # 从当前函数环境中批量获取变量对应的对象
  raster_list <- mget(layer_names, envir = environment())
  
  return(raster_list)
}

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