计算R中多个坐标的距离

Question

正如标题中所写，我想使用 package osrm 计算从我家乡的所有加油站到这里附近的两个高速公路入口的距离。

station_ms 包含加油站的纬度和经度，highway_ms.df 包含高速公路入口的纬度和经度。

只计算我的数据集的一行的距离是没有问题的，但我无法创建一个循环/函数，它对每一行都执行此操作。

这是我的代码：

route4 <-  osrmRoute(src = c(stations_ms$longitude[1], stations_ms$latitude[1]), 
                    dst = highway_ms.df[1,],
                    overview = "FALSE")

for (i in 1:nrow(stations_ms)) {
 route[i] <- osrmRoute(src = c(stations_ms$longitude[i], stations_ms$latitude[i]),
                       dst = highway_ms.df[1,],
                       overwiew = "FALSE")
}
```

Maybe someone can help me :)

Answer 1

这是一个可行的示例，可能会有所帮助。

osrmRoute中的overview有以下选项：

“完整”、“简化”或 FALSE。 使用“full”返回详细几何，使用“simplified”返回简化几何，使用 FALSE 仅返回时间和距离。

如果您只想要时间和距离，使用FALSE应该可以正常工作。 我的评论是关于拼写（用“w”而不是“v”）。

我做了一些示例数据：

my_points <- data.frame(
  id = 1:3,
  longitude = c(13.4, 13.5, 13.3),
  latitude = c(52.4, 52.5, 52.3)
)

并且想找到到柏林药房的距离（使用osrm软件包附带的apotheke.df ）。 你可以这样做：

library(osrm)

route <- list()
for (i in 1:nrow(my_points)) {
  route[[i]] <- osrmRoute(src = c(my_points$longitude[i], my_points$latitude[i]),
                          dst = apotheke.df[1,],
                          overview = FALSE)
}

这从一个名为route的空列表开始。 然后，我们用时间和持续时间填充每个列表元素。 最终结果是以下列表：

R> route

[[1]]
duration distance 
   20.56    11.77 

[[2]]
duration distance 
   17.38     7.63 

[[3]]
duration distance 
   33.12    27.45

可以转换为矩阵或数据框（在这种情况下，我制作了一个矩阵）：

R> do.call(rbind, route)

     duration distance
[1,]    20.56    11.77
[2,]    17.38     7.63
[3,]    33.12    27.45