如何在R中使用for循环将多个光栅裁剪到相同的程度

Question

我有4个光栅，我想裁剪到相同的程度。在这个脚本的未来迭代中，我将有超过4个，所以我试图编写一个循环，以相同的程度裁剪目录中的所有光栅。光栅是下载的Sentinel-2产品，其中至少包含4个波段，这些波段已使用sen2r((库转换为GeoTIFF。我已经尝试过在这里发布类似问题的答案，但在这个过程中不知何故丢失了波段，我稍后需要这些波段来做一些光栅计算。

到目前为止的代码：

raster_files <- list.files(here::here("data", "s2_rasters")) #dir with 4 rasters 
raster_paths <- paste0(here::here("data", "s2_rasters", raster_files))
wp_shp <- readOGR(here::here("data", "wp_boundary.shp"))
e <- extent(wp_shp)
n <- length(raster_paths)
for (i in 1:n) {
m <- raster_paths[i]
crop(x = m, y = e) 
}

编辑：我意识到我的循环没有意义。我是新手，知道我在做什么。直到脚本中的这一点，我一直在使用文件的路径来做一些事情(构建虚拟光栅，应用大气校正等(

这是我为一种作物做的一个例子，效果很好。

extent <- extent(802331.9, 802503.7, 9884986, 9885133)
ras_crop <- stack(here::here("data", "s2_rasters", "sample_raster.tif")) %>% 
crop(extent) %>%
writeRaster(filename=file.path(here::here("data", "s2_rasters"), "raster1_crop.tif"))

Answer 1

据我所知，您应该能够执行类似的操作

# input filenames
inf <- list.files("data/s2_rasters", pattern="tif$", full.names=TRUE)
# create output filenames and folder
outf <- gsub("data/s2_rasters", "output", inf)
dir.create("output", FALSE, FALSE)
library(raster)
wp_shp <- shapefile("data/wp_boundary.shp")
e <- extent(wp_shp)
for (i in 1:length(inf)) {
b <- brick(inf[i])
crop(b, e, filename=outf[i]) 
}