将多边形列表组合成一个对象

Question

我有一个栅格列表。我想根据每个栅格的范围创建多边形，并将所有新制作的多边形组合成一个对象。

我尝试了以下方法：

library(tidyverse)
library(raster)
library(sf)

lst(
  raster(ext = extent(20, 21, 10, 11)),
  raster(ext = extent(25, 26, 15, 16))
) %>% 
  map(
    ~ .x %>% 
      pluck("extent") %>% 
      as("SpatialPolygons")
  ) %>% 
  st_union()

但它会引发错误，因为 st_union 或 st_combine 不接受列表。

如何将多边形列表组合成一个多多边形对象？

Answer 1

你可以的

library(raster)

x <- list(
  raster(ext = extent(20, 21, 10, 11)),
  raster(ext = extent(25, 26, 15, 16))
) 

y <- lapply(x, function(i) as(extent(i), "SpatialPolygons"))
bind(y)