使用 for 循环估计 r 中的一个参考列和几列之间的距离?

Question

我有一个具有以下结构的数据框(df(：

ID   Long_Ref   Lat_Ref    Long_1   Lat_1    Long_2   Lat_2
A    -71.69     -33.39    -70.29   -33.39   -71.69   -34.19
B    -72.39     -34.34    -74.29   -31.19   -70.54   -33.38
C    -70.14     -32.38    -70.62   -32.37   -69.76   -32.22
D    -70.39     -33.54    -70.42   -34.99   -68.69   -32.33

我正在尝试使用 distHaversine 将新列附加到我现有的数据框中，其中包含 Ref 和 1、Ref 和 2 等之间的距离。

这是我到目前为止完成的代码：

for (i in 1:2){
for (j in 1:nrow(df)){
df[,paste("Dist_Mts_",i,sep="")][j]<-distHaversine(matrix(c(df$Long_Ref,df$Lat_Ref),ncol=2),
matrix(c(as.matrix(df[,paste("Long_",i,sep="")]),
as.matrix(df[,paste("Lat_",i,sep="")])),ncol=2))
}
}

但是，我不断收到以下错误：

Error in .pointsToMatrix(p2) * toRad : 
non-numeric argument to binary operator
In addition: Warning message:
In .pointsToMatrix(p2) : NAs introduced by coercion

我将不胜感激任何帮助。谢谢。

Answer 1

假设您有许多点可能大于 2 想要计算距离，您可以执行以下操作(不使用for循环(：

point.no <- 1:2 # depending on how many points you have
library(tidyverse)
point.no %>% 
map(~ paste0(c('Long_', 'Lat_'), .x)) %>% 
setNames(paste0('Dist_Mts_', point.no)) %>% 
map_df(~ distHaversine(df[c('Long_Ref', 'Lat_Ref')], df[.x]))

如果你只需要计算到这两个点的距离，那么你可以这样写：

list(Dist_Mts_1=c('Long_1', 'Lat_1'), Dist_Mts_2=c('Long_2', 'Lat_2')) %>% 
map_df(~ distHaversine(df[c('Long_Ref', 'Lat_Ref')], df[.x]))

输出：

# A tibble: 4 x 2
Dist_Mts_1 Dist_Mts_2
<dbl>      <dbl>
1     130123      89056
2     393159     201653
3      45141      39946
4     161437     208248