Question

我试图在R中使用distHavrsine函数，在循环内计算几百行的某些纬度和经度坐标之间的距离。在我的循环中，我有这个代码：

<select>

之后，如果距离小于50米，我希望它记录这些行，并且它所引用的纬度和经度坐标如下所示：

Category_Id

和

if ((distHaversine(c(file[i,"long"], file[i,"lat"]),
                   c(file[j,"long"], file[j,"lat"]))) < 50 )

但我收到此错误

.pointsToMatrix（p1）中的错误：纬度＆gt; 90

Answer 1

我收到此错误“.pointsToMatrix（p1）：纬度＆gt; 90”时出错。能够有人解释为什么以及如何解决？

错误告诉您纬度值大于90，这超出了范围：

library(geosphere)
distHaversine(c(4,52), c(13,52))
# [1] 616422
distHaversine(c(4,52), c(1,91))
# Error in .pointsToMatrix(p2) : latitude > 90

您可以通过仅使用可接受范围内的坐标提供distHaversine来解决此问题。

我试图在循环中使用R中的distHavrsine函数计算一些纬度和经度坐标之间的距离几百行。（...）如果距离小于50 米我希望它记录那些行

查看distm函数，它可以轻松计算几百行的距离矩阵（即没有循环）。它默认使用distHaversine。例如，要获得接近650000米的数据帧行：

df <- read.table(sep=",", col.names=c("lon", "lat"), text="
4,52
13,52 
116,39")
(d <- distm(df))
#         [,1]    [,2]    [,3]
# [1,]       0  616422 7963562
# [2,]  616422       0 7475370
# [3,] 7963562 7475370       0

d[upper.tri(d, T)] <- NA
( idx <- which(d < 650000, arr.ind = T) )
#      row col
# [1,]   2   1
cbind(df[idx[, 1], ], df[idx[, 2], ])
#   lon lat lon lat
# 2  13  52   4  52

如何使用distHaversine功能？

1 个答案: