st_network_paths генерирует только node_path с одиночными узлами

Question

Я пытаюсь создать список узлов на кратчайшем пути между двумя узлами, используя st_network_paths(). Однако я получаю только одно значение для индекса узла в node_path.

Он работает с данными игрушек, но не с данными реального мира. Что должно произойти, чтобы реальный мир, потоковая сеть, играли в мяч?

Данные линии доступны здесь

      library(sfnetworks)
library(sf)
library(tidygraph)
library(tidyverse)

ln <- st_read("river.gpkg")
    
net = as_sfnetwork(ln)

paths <- st_network_paths(net, 
                          from = 2,
                          to = 50)

# plot one path to check
node_path <- paths %>%
  slice(1) %>%
  pull(node_paths) %>%
  unlist()

node_path

plot(net, col = "grey")
plot(slice(activate(net, "nodes"), 2), 
     col = "blue", add = TRUE)
plot(slice(activate(net, "nodes"), 50), 
     col = "red", add = TRUE)
plot(slice(activate(net, "nodes"), node_path), 
     col = "green", add = TRUE) # just recreates the node_path

Answer 1

Проблема в том, что вы пытаетесь вычислить путь между двумя разными ветвями в направленной сети:

      # packages
library(sf) 
library(tidygraph)
library(sfnetworks)

# data
ln <- st_read("C:/Users/andre/Desktop/river.gpkg", quiet = TRUE)
net = as_sfnetwork(ln)

# plot
par(mar = rep(0, 4))
plot(net)
plot(net %N>% slice(2, 50), add = TRUE, cex = 2, col = c("red", "blue"))

^{Создано 15 декабря 2021 г. пакетом reprex (v2.0.1)}

Если желаемые пути также можно рассчитать, путешествуя «назад», вы можете установить аргумент directedв as_sfnetworkк FALSE(т.е. net = as_sfnetwork(ln, directed = FALSE)и запустите тот же код, что и раньше. Вы можете распознать этот тип проблем, поскольку st_network_paths()возвращает предупреждающее сообщение, например

      > paths <- st_network_paths(
+   net, 
+   from = 2, 
+   to = 50
+ )
Warning message:
In shortest_paths(x, from, to, weights = weights, output = "both",  :
  At structural_properties.c:4745 :Couldn't reach some vertices