有几种方法可以简化代码。由于对每个VOI执行相同的操作，因此可以将它们存储在一个向量中。注意，这里的voi包含所有5个VOI：

base_url <- "https://files.isric.org/soilgrids/latest/data/" # Path to the webDAV data.
voi <- c('sand', 'clay', 'silt', 'phh2o', 'soc')
depth <- "5-15cm"
quantile <- "mean"

我们可以使用sprintf来更简洁地创建fetch URL：

urls <- sprintf('%s%s/%s_%s_%s.vrt', base_url, voi, voi, depth, quantile)
names(urls) <- voi
                                                                       sand 
  "https://files.isric.org/soilgrids/latest/data/sand/sand_5-15cm_mean.vrt" 
                                                                       clay 
  "https://files.isric.org/soilgrids/latest/data/clay/clay_5-15cm_mean.vrt" 
                                                                       silt 
  "https://files.isric.org/soilgrids/latest/data/silt/silt_5-15cm_mean.vrt" 
                                                                      phh2o 
"https://files.isric.org/soilgrids/latest/data/phh2o/phh2o_5-15cm_mean.vrt" 
                                                                        soc 
    "https://files.isric.org/soilgrids/latest/data/soc/soc_5-15cm_mean.vrt"

最后，我们可以使用lapply对每个构造的URL执行raster()函数（lapply是一个循环函数）。

rasters <- lapply(urls, raster)

请注意，由于我们已经为上述urls的元素提供了名称，因此这些名称将保留在rasters结果中，因此您可以使用数字索引（rasters[1]）或命名索引（rasters['sand']）来检索各个栅格。

另一个变化

url = "https://files.isric.org/soilgrids/latest/data/"
voi = c("sand", "clay", "silt", "phh2o", "soc")
depth = "5-15cm"
quantile = "mean"

f <- paste0(url, voi, "/", paste(voi, depth, quantile, sep="_"), '.vrt')

“raster”已过时，我们改用“terra”并创建SpatRasterDataset

library(terra)
x <- sds(f)

或SpatRasters列表

y <- lapply(f, rast)

另请参阅geodata::soil_world和geodata::soil_world_vsi以获取这些数据。

一般来说，你可以创建一个函数来创建一个对象，在一个for循环中调用它，然后使用assign来为它们赋值变量名。我推荐使用“glue”包中的glue函数，因为它比paste更直观。我不能让raster函数在我的电脑上工作，因为我得到了this error。然而，这就是我的想法。

library(glue)

raster_object <- function(url, voi, depth, quantile){
  return(raster(glue("{url}/{voi}/{voi}_{depth}_{quantile}.vrt")))
}

for (voi in c("sand", "clay")){
  assign(glue("raster_{voi}"),  raster_object(url = "https://files.isric.org/soilgrids/latest/data", voi, depth = "5-15cm", quantile = "mean"))
}