如何从R中的栅格列表中提取相同范围的栅格

Question

我有一个由栅格组成的列表对象，就其范围而言，栅格有两种类型。如何提取具有相同范围的文件并将其另存为R中的新列表？不可能给出可复制的示例，但是下面的代码可以帮助您理解问题。

一个假设的例子：

library(raster)
list_raster # Suppose this is a list having 48 rasters.
names(list_raster) <- paste0("raster", seq(1:48)) # Assigning names makes it possible to use the dollar sign.
lapply(list_raster, extent) # Gives 48 results but only two unique raster extents.
#  I would like to know which of the rasters has == extent(list_rasters$raster1)

Answer 1

您可以使用sapply：

same_as_r1 <- sapply(list_raster, function(x) extent(x) == extent(list_raster$raster1))

并使用它来细分您的列表：

group1 <- list_raster[same_as_r1]
group2 <- list_raster[!same_as_r1]