我有iGraph对象中每个顶点的地理坐标。现在我想计算现有边之间的距离。
library(igraph)
library(ggmap)
library(geosphere)
g <- graph.ring(6)
V(grph)$postcode <- c("Johannesburg 2017",
"Rondebosch 8000",
"Durban 4001",
"Pietermaritzburg 3201",
"Jeffreys Bay 6330",
"Pretoria 0001" )
postcode_df <- geocode(V(g)$postcode, sensor = FALSE,
output = "latlon", source = "google")
V(g)$coordinate <- split(postcode_df, 1:nrow(postcode_df))
V(g)$coordinate[1]
[[1]]
lon lat
1 28.03837 -26.18825
我想使用以下通用方法计算距离：
el <- get.edgelist(g, names=FALSE)
E(g)$distance <- distHaversine(V(g)$coordinate[[el[,1]]],V(g)$coordinate[[el[,2]]])
问题是 V(g)$coordinate 中的 lon 和 lat 不能以这种方式引用。我在级别 3 处获得递归索引失败。显然，我无法将一个数据帧的索引嵌套在另一个数据框中。
str(V(g)$coordinate)
List of 6
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 28
..$ lat: num -26.2
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 28.3
..$ lat: num -25.8
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 31
..$ lat: num -29.8
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 30.4
..$ lat: num -29.7
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 24.9
..$ lat: num -34.1
$ :'data.frame': 1 obs. of 2 variables:
..$ lon: num 28.2
..$ lat: num -25.7
计算两点之间距离的一般方法是
distHaversine(p1, p2, r=6378137).
P1 由 EL[，1] 定义，P2 由 EL[，2] 定义。 EL[，1：2] 是指以 G 为单位的顶点数。所以我需要提取 V(g)$coordinate对应于 el[，1] 和 el[，2]。建议将不胜感激。
这里有一个问题，因为split返回一个数据帧，我们可以通过以下方式修复：
V(g)$coordinate <- lapply(split(postcode_df, 1:nrow(postcode_df)), unlist)
然后，您基本上需要迭代两个列表，即每个顶点的坐标。
这很容易从purrrmap2。
library(purrr)
el <- get.edgelist(g, names=FALSE)
E(g)$distance <- unlist(map2(V(g)$coordinate[el[,1]], V(g)$coordinate[el[,2]], distHaversine))