Question

我有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）$坐标中的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为单位的顶点数。所以我需要提取对应于el [，1]和el [，2]的V（g）$坐标。建议将不胜感激。

Answer 1

这里存在一个问题，split正在返回一个数据框，我们可以通过以下方式修复：

V(g)$coordinate <- lapply(split(postcode_df, 1:nrow(postcode_df)), unlist)

然后你基本上需要迭代两个列表，即每个顶点的坐标。

来自map2的{{1}}很容易。

purrr

计算iGraph中边缘的地理距离

1 个答案: