我想知道r中是否有一个特定函数来计算步长，是两个连续的GPS GPS时间注册位置之间的距离（以米为单位）。
我的数据集如下：
> head(datagps)
Date & Time [Local] Latitude Longitude
1: 2018-06-18 03:01:00 -2.434901 34.85359
2: 2018-06-18 03:06:00 -2.434598 34.85387
3: 2018-06-18 03:08:00 -2.434726 34.85382
4: 2018-06-18 03:12:00 -2.434816 34.85371
5: 2018-06-18 03:16:00 -2.434613 34.85372
6: 2018-06-18 03:20:00 -2.434511 34.85376
并希望创建一个Step列，该列将进行上述操作。也许geosphere软件包具有这样的功能？如果不是，那么最紧凑的方法是什么？
任何输入都将受到赞赏！
您可以使用geosphere并计算两个坐标之间的distHaversine：
library(geosphere)
distm(c(lon1, lat1), c(lon2, lat2), fun = distHaversine)
使用mutate添加Step字段
i在https://github.com/michaelmalick/r-malick/blob/master/master/r/hhaversine.r
haversine <- function(lon1, lat1, lon2, lat2, r = 6378137) {
if(!is.numeric(c(lon1, lat1, lon2, lat2)))
stop("Inputs are not numeric")
# Convert degrees to radians
lon1 <- lon1 * pi / 180
lat1 <- lat1 * pi / 180
lon2 <- lon2 * pi / 180
lat2 <- lat2 * pi / 180
delta.lon <- (lon2 - lon1)
delta.lat <- (lat2 - lat1)
a <- sin(delta.lat/2)^2 + cos(lat1) * cos(lat2) *
sin(delta.lon/2)^2
c <- 2 * asin(min(1,sqrt(a)))
d <- r * c
return(d) # Distance
}
vectorized_haversine <- Vectorize(haversine, vectorize.args = c("lon1", "lat1", "lon2", "lat2"))
接下来，我使用了dplyr函数"滞后"one_answers"突变"，以及矢量化的haversine函数以获得连续点之间的距离（使用'tribble'函数来重新创建head(datagps)）。
library(dplyr)
tribble(
~`Date & Time [Local]`, ~Latitude, ~Longitude,
"2018-06-18 03:000", -2.434901, 34.85359,
"2018-06-18 03:06:00", -2.434598, 34.85387,
"2018-06-18 03:08:00", -2.434726, 34.85382,
"2018-06-18 03:12:00", -2.434816, 34.85371,
"2018-06-18 03:16:00", -2.434613, 34.85372,
"2018-06-18 03:20:00", -2.434511, 34.85376
) %>%
mutate(Step =
vectorized_haversine(Longitude, Latitude, lag(Longitude), lag(Latitude)))
Date & Time [Local] Latitude Longitude Step
1 2018-06-18 03:000 -2.434901 34.85359 NA
2 2018-06-18 03:06:00 -2.434598 34.85387 45.90731
3 2018-06-18 03:08:00 -2.434726 34.85382 15.29559
4 2018-06-18 03:12:00 -2.434816 34.85371 15.81292
5 2018-06-18 03:16:00 -2.434613 34.85372 22.62521
6 2018-06-18 03:20:00 -2.434511 34.85376 12.19500
查看： gmapdistance
注意：包含起点描述的字符串或矢量。应该在Quoutes内部（"）。如果使用了同一位置多个单词，则应通过一个加号分开，例如"波哥大 哥伦比亚"。，只要可以通过Google Maps识别，以纬度格式的坐标格式也是有效的输入。
https://cran.r-project.org/web/packages/gmapsdistance/gmapsdistance.pdf
函数gmapsdistance使用Google Maps距离矩阵API，以计算两个点之间的距离和时间。为了能够使用该函数，您将需要一个API键，并启用Google Developers Console Console
中的距离矩阵API警告：
请注意，Google改变了其计费惯例。（欢乐）
您必须有一个有效的API密钥和一个计费帐户才能访问我们的API。当您启用帐单时，每月将获得200美元的地图，路线或地点的免费使用。根据当今使用我们的API的数百万用户，其中大多数可以继续免费使用Google Maps平台，并使用此信用。拥有一个计费帐户可以帮助我们更好地了解开发人员的需求，并使您无缝扩展。
有关如何获取钥匙的更多信息，请访问https://developers.google.com/maps/documentation/distancematrix/get-ap-papi-key#key
有关Google Maps距离矩阵API的更多信息，请访问https://developers.google.com/maps/documentation/distance-matrix/intro？hl = en
示例：
在原点和目的地之间的旅行时间和距离的列表以及状态
results = gmapsdistance(origin = "38.1621328+24.0029257",
destination = "37.9908372+23.7383394",
mode = "walking")