如何计算R中沿一条线的两点之间的地理距离？

Question

输入

我有两个形状文件，我输入到R，所以我结束了。一种包含公共汽车路线的空间数据。一个包含公交车站的空间点。

划出一个给定的路线和它的停站如下所示。

​

​

样本数据

此链接包含两个表单文件，作为一个zip下载，其中包含两个路径示例。

目标

我的目标是计算每一对车站之间的地理距离(以米为单位)：停1到停2，停2到停3，等等。

我发现的大多数计算欧几里德距离的方法，或者说是“乌鸦苍蝇”，在这里都是行不通的。

这个职位提到了PBSmapping，它有一个calcLength函数，它可以很好地计算路由的总距离，但是我无法找到一种方法来匹配它与停止对的情况，也无法找到一种方法来根据它的属性对shapefile进行子集。

riverdist包同样有趣，但对我无法应用的河流进行了高度优化。

Answer 1

尝试使用来自gProject包的rgeos：

library("rgdal")
library("rgeos")

# read shapefile and transfrom to UTM zone 36N (distances in m)
route <- spTransform(readOGR(".", "Trips"), CRS("+init=epsg:32636"))
stops <- spTransform(readOGR(".", "Stops"), CRS("+init=epsg:32636"))

l <- subset(route, route_id==1137)
p <- subset(stops, grepl("^1137_", UID))

plot(l, axes=TRUE, col="orange")
plot(p, add=TRUE, pch=19, cex=0.1)
text(p)

​

​

# distance along route
d <- sort(gProject(l, p))
d
# [1]     0  3051  3057  7221 10379 15657 20326 20326 22141 24262

# distance between stops
diff(d)
#[1] 3050.9166    5.9720 4164.2480 3157.7702 5278.5812 4668.1810    0.5878
#[8] 1814.9612 2120.8470