R中的sf和网络：从不同的线串和点sf对象创建一个边缘列表

Question

我有一个节点（站点）的 sf 对象和另一个线串（路线）。

# toy example 

## nodes
p1 = st_point(c(7, 51))
p2 = st_point(c(7, 52))
p3 = st_point(c(7, 53))
p4 = st_point(c(8, 52))
nodes = st_as_sf(st_sfc(p1, p2, p3,p4, crs = 4326))



## routes

e1 = st_cast(st_union(p1,p3), "LINESTRING")
e2 = st_cast(st_union(p1,p4), "LINESTRING")
e3 = st_cast(st_union(p3,p4), "LINESTRING")
lines = st_as_sf(st_sfc(e1, e2, e3, crs = 4326))

如何获得直接连接两个节点的路径段的边缘列表？

#Desired output

from | to
p1   | p2
p2   | p3
p1   | p4
p3   | p4

这里 p1 和 p3 之间有一条线，但 p2 在中间，所以 p1-p2 和 p2-p3 之间有 2 条边

我知道sfnetwork包可以用线串或点来构建网络，但是如何用线和点的交点来构建空间网络呢？

Answer 1

不确定您在寻找什么，但我会再试一次，作为您评论的后续！因此，请在下面找到使用sf、sfnetworks 和dplyr 库的完全编辑的reprex。

Reprex

• 第 1 步：使用您的点/节点 p1、p3 和 p4 创建一个网络

library(sf)
library(sfnetworks)
library(dplyr)

# toy example 
p1 = st_point(c(7, 51))
p3 = st_point(c(7, 53))
p4 = st_point(c(8, 52))
nodes = st_as_sf(st_sfc(p1, p3, p4, crs = 4326))
nodes$names <- c("p1", "p3", "p4") # add the name of nodes


## routes
e1 = st_cast(st_union(p1,p3), "LINESTRING")
e2 = st_cast(st_union(p1,p4), "LINESTRING")
e3 = st_cast(st_union(p3,p4), "LINESTRING")
lines = st_as_sf(st_sfc(e1, e2, e3, crs = 4326))
lines$from  <-  c("p1", "p1", "p3")
lines$to  <-  c("p3", "p4", "p4")


# Create the network
network <- sfnetwork(nodes, lines, node_key = "names")
#> Checking if spatial network structure is valid...
#> Spatial network structure is valid


# What the network looks like:
network
#> # A sfnetwork with 3 nodes and 3 edges
#> #
#> # CRS:  EPSG:4326 
#> #
#> # A directed acyclic simple graph with 1 component with spatially explicit edges
#> #
#> # Node Data:     3 x 2 (active)
#> # Geometry type: POINT
#> # Dimension:     XY
#> # Bounding box:  xmin: 7 ymin: 51 xmax: 8 ymax: 53
#>             x names
#>   <POINT [°]> <chr>
#> 1      (7 51) p1   
#> 2      (7 53) p3   
#> 3      (8 52) p4   
#> #
#> # Edge Data:     3 x 3
#> # Geometry type: LINESTRING
#> # Dimension:     XY
#> # Bounding box:  xmin: 7 ymin: 51 xmax: 8 ymax: 53
#>    from    to                x
#>   <int> <int> <LINESTRING [°]>
#> 1     1     2     (7 51, 7 53)
#> 2     1     3     (7 51, 8 52)
#> 3     2     3     (7 53, 8 52)

• 第 2 步：将补充节点（即 p2）添加到“网络”对象

# Create a 'sf' object for point 2 'p2'
p2 = st_point(c(7, 52))
p2_sf <- st_as_sf(st_sfc(p2, crs = 4326))
p2_sf$names <- "p2"

# Add 'p2' to the network and get the 'new_network' object
new_network <- st_network_blend(network, p2_sf)


# Cleaning the 'new_network' object
new_network %>% 
  activate("nodes") %>% 
  mutate(names = coalesce(names.x, names.y)) %>% 
  arrange(., names) %>%
  select(names,x) -> new_network

• 第 3 步：创建具有所需结果的 sf 对象

new_network %>% 
  activate("edges") %>% 
  st_as_sf() %>% 
  mutate(from = paste0("p", from), to = paste0("p", to))

#> Simple feature collection with 4 features and 2 fields
#> Geometry type: LINESTRING
#> Dimension:     XY
#> Bounding box:  xmin: 7 ymin: 51 xmax: 8 ymax: 53
#> Geodetic CRS:  WGS 84
#> # A tibble: 4 x 3
#>   from  to                   x
#> * <chr> <chr> <LINESTRING [°]>
#> 1 p1    p2        (7 51, 7 52)
#> 2 p2    p3        (7 52, 7 53)
#> 3 p1    p4        (7 51, 8 52)
#> 4 p3    p4        (7 53, 8 52)

^{由reprex package (v2.0.1) 于 2021-12-05 创建}