Używając pakietu R sf
, jak łączyć sfc
obiekty? Na przykład, biorąc pod uwagę poniższy kod, jak stworzyć pojedynczy sfc
obiekt, sfc12
który zawiera geometrie z obu sfc1
i sfc2
? ( length(sfc12)
powinno być 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 = ?
Niektóre podejścia, które nie działają:
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))
pracował dla mnie.sfc
obiektów za pomocąReduce(c, sfcs)
lubpurrr::reduce(sfcs , c)
Czerpiąc z @ TimSalabim za odpowiedź , jeśli
sfc
obiekty są w tych samych CRS można użyćdo.call(rbind, list(sfc1, sfc2))
.źródło