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

使用for循环估算r中一个参考列与几个列之间的距离？

1 个答案: