Usando el paquete R sf
, ¿cómo se combinan los sfc
objetos? Por ejemplo, dado el siguiente código, ¿cómo se puede crear un solo sfc
objeto sfc12
que incluya las geometrías de ambos sfc1
y sfc2
? ( length(sfc12)
debe ser 2.)
library(sf)
pt1 = st_point(c(0,1))
pt2 = st_point(c(1,1))
sfc1 = st_sfc(pt1) # An sfc object
sfc2 = st_sfc(pt2) # Another sfc object
# sfc12 = ?
Algunos enfoques que no funcionan:
sf_sfc(sfc1, sfc2)
# Error in vapply(lst, class, rep("", 3)) : values must be length 3,
# but FUN(X[[1]]) result is length 2
sfc1 + sfc2 # Seems to add the points coordinate-wise.
# Geometry set for 1 feature
# geometry type: POINT
# dimension: XY
# bbox: xmin: 1 ymin: 2 xmax: 1 ymax: 2
# epsg (SRID): NA
# proj4string: NA
# POINT(1 2)
rbind(sfc1, sfc2)
# [,1]
# sfc1 Numeric,2
# sfc2 Numeric,2
do.call(c, list(sfc1, sfc2))
trabajó para mí.sfc
objetos usandoReduce(c, sfcs)
opurrr::reduce(sfcs , c)
Basándose en la respuesta de @ TimSalabim , si sus
sfc
objetos están en el mismo CRS, puede usardo.call(rbind, list(sfc1, sfc2))
.fuente