我有两个空间特征:
library(sf)
points1 <- data.frame(foo = seq(15, 75, 15),
long = c(-85, -80, -78, -75, -82),
lat = c(34, 36, 37, 38, 35)) %>%
st_as_sf(coords = c('long', 'lat'), crs = 4326)
points2 <- data.frame(bar = seq(15, 75, 15),
long = c(85, 80, 78, 75, 82),
lat = c(30, 32, 34, 36, 38)) %>%
st_as_sf(coords = c('long', 'lat'), crs = 4326)
cbind(points1, points2) -> df
这样做会产生以下结果:
foo bar geometry geometry.1
1 15 15 POINT (-85 34) POINT (85 30)
2 30 30 POINT (-80 36) POINT (80 32)
3 45 45 POINT (-78 37) POINT (78 34)
4 60 60 POINT (-75 38) POINT (75 36)
5 75 75 POINT (-82 35) POINT (82 38)
我想在df
中点对之间绘制一条线-从geometry
中的一个点到geometry.1
中的一个点。我尝试将这些点转换为一个LINESTRING
,方法如下:
df %>% summarise(do_union=F) %>% st_cast("LINESTRING") %>% plot()
但是这似乎不起作用。我得到的是一条连续的线,而我想要的是五个单独的行。