时间戳

首页 > 解决方案 > 在R中，如何求和行之间的距离（每一行都是一个GPS坐标）？

r - 在R中，如何求和行之间的距离（每一行都是一个GPS坐标）？

问题描述

我有很多 GPS 点，我想要的是给定日期内两个后续点（行）之间的距离总和，这样我就可以获得每日跟踪距离。

每天大约有 200 个 GPS 点。后续点是指两行，其中第一行早于第二行。因为我需要这些点之间的总距离，所以必须考虑给定日期内“时间”列（“日期”列）的行顺序。

谢谢！

我的表格看起来像这样：

date        time       lat        lon    
18-Jan-18 12:48:39 -24.061464 -47.99523
18-Jan-18 12:48:48 -24.06163  -47.995354
18-Jan-18 12:53:17 -24.06175  -47.995277

标签： rgps

解决方案

我们可以使用pointDistancefrom 包raster来计算距离。lagfromdplyr将有助于计算后续点。replace_nafromtidyr非常方便，但您可以使用自己喜欢的方式处理NA.

library(raster)
library(dplyr)
library(tidyr)
data %>% 
  mutate(Distance = pointDistance(cbind(lon,lat),cbind(lag(lon),lag(lat)),lonlat = TRUE)) %>%
  mutate(TotalDistance = cumsum(replace_na(Distance,0)))
#       date     time       lat       lon Distance TotalDistance
#1 18-Jan-18 12:48:39 -24.06146 -47.99523       NA       0.00000
#2 18-Jan-18 12:48:48 -24.06163 -47.99535 22.29547      22.29547
#3 18-Jan-18 12:53:17 -24.06175 -47.99528 15.42660      37.72207

数据

data <- structure(list(date = structure(c(1L, 1L, 1L), .Label = "18-Jan-18", class = "factor"), 
    time = structure(1:3, .Label = c("12:48:39", "12:48:48", 
    "12:53:17"), class = "factor"), lat = c(-24.061464, -24.06163, 
    -24.06175), lon = c(-47.99523, -47.995354, -47.995277)), class = "data.frame", row.names = c(NA, 
-3L))

推荐阅读