如何計算R中沿線的兩個點之間的地理距離？

Question

輸入

我有兩個導入到R中的shapefile，所以最終得到了。 包含公交路線的空間線數據框。 包含公交車站的空間點數據幀。

用停靠點繪制給定路線看起來像這樣。

樣本數據

此鏈接包含兩個shapefile，以兩個示例的路徑作為zip下載。

目標

我的目標是計算每對站點之間的地理距離，以米為單位：跨越基礎公交路線的長度的Stop 1至Stop 2，Stop 2至Stop 3等。

我發現的大多數方法都是計算歐幾里得距離的，也就是“烏鴉飛過”的。 在這里不起作用。

這篇文章提到了PBSmapping ，它具有calcLength函數，可以很好地計算路線的總距離，但是我找不到將其與停靠站情況相匹配的方法，也找不到真正的子集shapefile通過其屬性。

riverdist軟件包同樣有趣，但是針對河流進行了高度優化，我找不到應用它的方法。

Answer 1

從rgeos包嘗試gProject ：

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