循环中一个多边形的栅格裁剪中的新栅格 R

Question

从 19 个气候栅格中，我尝试仅裁剪多边形中的栅格区域（在法国）。 为此，我尝试生成一个循环，用多边形逐一裁剪每个栅格，如下所示：

# my list of raster
rlist <- list.files(path=mydir, pattern='tif$', full.names=FALSE)
# loop to generate the 19 raster
for(i in rlist) {assign(unlist(strsplit(i, "[.]"))[1], raster(i)) }

# import my polygon (France)
France <- readOGR("./QGIS/France/Shapefile/France.shp")

# crop and mask one raster by polygon France
c1 <- crop(raster1, extent(France))
c2 <- mask(c1, France)

它有效，但现在我试图找到一个循环来一次裁剪和遮罩 19 个栅格，并从我会命名的这个裁剪和遮罩中创建新栅格。 有什么建议么？

我不知道是否清楚，如果您需要更多信息，请告诉我。 谢谢。

Answer 1

我建议你切换到terra 。 它与raster库来自同一作者，但速度更快：

library(terra)

files <- list.files(path=mydir, pattern='tif$', full.names=TRUE)
v <- vect("./QGIS/France/Shapefile/France.shp")

l <- list()

for(i in seq_along(files)){
  rtemp <- rast(files[i])
  l[[i]] <- crop(rtemp, v, mask=TRUE)
}

s <- rast(l)

Answer 2

不需要循环

library(terra)
ff <- list.files(path=mydir, pattern='tif$', full.names=FALSE)
x <- rast(ff)
france <- vect("./QGIS/France/Shapefile/France.shp")
xf <- crop(x, france, mask=TRUE)

如果xf有 19 层并且您想将每个层写入一个单独的文件（如您在评论中指出的那样），您可以执行类似

fout <- paste0("bio_", 1:19)
writeRaster(xf, fout, filetype="GTiff")

此外，您永远不应该使用assign 。 在这种情况下，您可以制作一个列表，以便稍后可以迭代其项目（使用for或lapply 。您正在做的事情可以写成

x <- lapply(ff, rast)
y <- lapply(x, \(i) crop(i, France, mask=TRUE))
z <- rast(y)

但这仍然很笨拙