Question

我计算球体上大量位置（5000）之间的距离矩阵（使用Haversine距离函数）。

这是我的代码：

require(geosphere)
x=rnorm(5000)
y=rnorm(5000)
xy1=cbind(x,y)

计算距离矩阵所需的时间是

 system.time( outer(1:nrow(xy1), 1:nrow(xy1), function(i,j) distHaversine(xy1[i,1:2],xy1[j,1:2])))

执行此程序所需的时间很长。任何建议如何降低时间消耗来完成这项工作！谢谢。

Answer 1

尝试obj包中的内置功能？

geosphere

z <- distm( xy1 )的默认距离函数 - 计算一组点之间的距离矩阵 - 是Haversine（distm()）公式，但您可以使用"distHaversine"参数指定另一个

在我的2.6GHz酷睿i7 rMBP上，这需要5秒钟，持续5,000点。

Answer 2

我在下面添加了使用spatialrisk软件包的解决方案。此软件包中的关键功能是用C ++（Rcpp）编写的，因此非常快。

library(geosphere)
library(spatialrisk)
library(data.table)

x=rnorm(5000)
y=rnorm(5000)
xy1 = data.table(x,y)

# Cross join two data tables
coordinates_dt <- optiRum::CJ.dt(xy1, xy1)

system.time({
  z <- distm( xy1 )
})
# user  system elapsed 
# 14.163   3.700  19.072 

system.time({
  distances_m <- coordinates_dt[, dist_m := spatialrisk::haversine(y, x, i.y, i.x)]
})
# user  system elapsed 
# 2.027   0.848   2.913

R：计算大距离矩阵的更快方法

2 个答案: