无法在R中运行osmar导航演示。可能是因为演示需要igraph0，这是不推荐使用的

Question

我正在尝试在R中运行osmar导航演示。此演示将使用osmar和igraph绘制基于openstreetmap数据的慕尼黑市中心周围的交通路线。

我在Lubuntu上使用R版本3.1.1

这里将详细介绍演示和osmar库http://journal.r-project.org/archive/2013-1/eugster-schlesinger.pdf

要运行我输入的演示，

library("osmar")
library("igraph")     # The demo tries to call igraph0, but this is 
                      # no-longer available in my version of R, so I
                      # have changed it to "igraph"
demo("navigator")

演示完美运行，直到它到达igraph部分。

gr_muc<-as_igraph(hways_muc)      # make a graph of the highways from openstreetmap
summary(gr_muc)

这应该返回

Vertices: 2381
Edges: 2888
Directed: TRUE
No graph attributes.
Vertex attributes: name.
Edge attributes: weight, name.

但对我来说它会返回

IGRAPH DNW-2385 2894 --
attr: name (v/c), weight (e/n), name (e/n)

我知道gr_muc是图表，因为命令E(gr_muc)和V(gr_muc)会返回边和顶点列表。

然后演示

route <- get.shortest.paths(gr_muc,from = as.character(hway_start_node),to = as.character(hway_end_node))[[1]]

并返回错误

At structural_properties.c:4482 :Couldn't reach some vertices

这意味着它无法链接起始和结束顶点。然后脚本失败了。

我要更改什么才能运行演示脚本，为什么它不起作用？

Answer 1

有一些单行道阻止开始和结束节点的连接。这里有一些简洁的代码，用于绘制从hway_start节点可到达的节点：

plot(hways_muc)

hway_start # osmar object
plot_nodes(hway_start, add = TRUE, col = "red", pch = 19, cex = 2)

# get reachable nodes, return graph Vertex indexes
n1 = neighborhood(gr_muc,200,V(gr_muc)[as.character(hway_start_node)], mode="out")

# get graph subset
allpts = V(gr_muc)[n1[[1]]]

# use names to subset OSM nodes:
nds = subset(muc, node(allpts$name))
plot_nodes(nds, add=TRUE,col="green",pch=19,cex=1)

请注意，您无法到达右上角，这是演示中终端节点的位置。

如果您不介意通过使图表无向而在错误的单向街道上行驶，您可以解决这个问题：

gru_muc=as.undirected(gr_muc)
route <- get.shortest.paths(gru_muc,
                             from = as.character(hway_start_node),
                             to = as.character(hway_end_node))[[1]]

现在你有了一条路线：

> route
[[1]]
  [2] 1444  491 2055  334  331  481  478  479  [etc]

但是get_shortest_paths的返回是一个列表列表，因此您需要获取route的第一个组件才能继续演示代码：

route=route[[1]]
route_nodes <- as.numeric(V(gr_muc)[route]$name)

然后绘图：

所以我认为首先，开始和结束节点没有按照指示的方式连接，其次是演示代码中的一个错误，因此没有从{{1获得正确的返回元素}}。它与get_shortest_paths无关。