计算数据帧中两个长纬度坐标之间的距离

Question

我想计算几个 GPS 点之间的距离。我试过

distm(c(lon1,lat1), c(lon2,lat2), fun = distHaversine)

这仅适用于一点，但不适用于我的数据框中的列。

所以我按照这里的建议尝试：

计算 2 lat longs 之间的距离

但是对于这两个计算，我确实得到了不同的结果：

df <- read.table(sep=",", col.names=c("lat1", "lon1", "lat2", "lon2"),text="
7.348687,53.36575,7.348940,53.36507 
7.348940, 53.36507,7.350939,53.36484")


# as recommended in the link above
distHaversine(df[,2:1], df[,4:3])

[1]  80.18433 223.97181

# with distm
distm(c(7.348687,53.36575), c(7.348940,53.36507), fun = distHaversine)

         [,1]
[1,] 77.54033

distm(c(7.348940, 53.36507), c(7.350939,53.36484), fun = distHaversine)

         [,1]
 [1,] 135.2317

那么如何计算数据帧列中两个 GPS 点之间的正确距离（即 distm(c(lon1,lat1), c(lon2,lat2), fun = distHaversine)）？我仔细检查了这么远的距离，我知道我通过这种方式获得了正确的距离。

提前致谢。

Answer 1

Given that the output you want to store in the new column is this:

77.54033 135.23165

Try this

df$distance<-distHaversine(df[,1:2], df[,3:4])

Which should return

> df
      lat1     lon1     lat2     lon2  distance
1 7.348687 53.36575 7.348940 53.36507  77.54033
2 7.348940 53.36507 7.350939 53.36484 135.23165

Answer 2

问题到底是什么？您还没有达到所需的所有距离吗distHaversine()？

您想将距离添加为数据框中的列吗？干得好：

f$dist <- distm(x = df[, c('lon1', 'lat1')], 
                y = df[, c('lon2', 'lat2')],
                fun = distHaversine
                )