我有一个栅格列表。我想根据每个栅格的范围创建多边形,并将所有新创建的多边形组合成一个对象。
我尝试了以下方法:
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
不接受列表。
如何将多边形列表组合成多多边形对象?
最佳答案
你可以做到
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)
关于r - 将多边形列表组合成单个对象,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/61999131/