r语言 - 计算多个记录的两点之间的距离以匹配行 - 循环两个矩阵的行

Question

我有两个带有坐标的矩阵，我正在尝试计算匹配行中点之间的距离，即第一个矩阵中的第 1 行和第二个矩阵中的第 1 行之间的距离。

我得到的是计算第 1 行和所有其他行之间的距离。这会产生内存问题，因为我有 800,000 行。有谁知道如何要求？

我正在使用

dist1 <- distm(FareStageMatrix[1:25000,], LSOACentroidMatrix[1:25000,], fun=distHaversine)

我正在尝试创建这样的东西，但似乎不起作用

for(i in 1:nrow(FareStageMatrix)) {
for(j in 1:nrow(LSOACentroidMatrix)) {
my_matrix[i] <- my_matrix[distm(FareStageMatrix[i], LSOACentroidMatrix[i], fun=distHaversine)]
}
}

改为

for (i in 1:nrow(FareStageMatrix)){
for (i in 1:nrow(LSOACentroidMatrix)){
r1<-FareStageMatrix[i,]
r2<-LSOACentroidMatrix[i,]
results[i]<-distm(r1, r2, fun=distHaversine)  
}
}

这是应该起作用的东西吗？

Answer 1

似乎我已经设法找到了解决方案：

results<-matrix(NA,nrow(FareStageMatrix))
for (i in 1:nrow(FareStageMatrix)){
for (i in 1:nrow(LSOACentroidMatrix)){
r1<-FareStageMatrix[i,]
r2<-LSOACentroidMatrix[i,]
results[i]<-distm(r1, r2, fun=distHaversine)  ## Example function
}
}

其中 FareStageMatrix 和 LSOACentroidMatrix 是带有坐标的矩阵

它似乎已经为给定的一对点计算了一个距离

Answer 2

为此，我已经调整了geosphere的distGeo函数(测地线距离)。

library(geosphere)
source("https://raw.githubusercontent.com/RomanAbashin/distGeo_v/master/distGeo_v.R")

数据

set.seed(1702)
m1 <- matrix(runif(20000, -10, 10), ncol = 2) 
m2 <- matrix(runif(20000, -10, 10), ncol = 2)

法典

result <- distGeo_v(m1[, 1], m1[, 2], 
m2[, 1], m2[, 2])

结果

> head(m1)
[,1]      [,2]
[1,]  8.087152  9.227607
[2,]  9.528334  9.103403
[3,]  5.637921 -2.213228
[4,] -2.473758 -9.812986
[5,] -2.844036 -5.245779
[6,] -4.824615 -4.330890
> head(m2)
[,1]       [,2]
[1,]  0.1673027  0.6483745
[2,] -2.5033184  0.1386050
[3,]  4.8589785  5.1996968
[4,]  8.3239454 -8.9810949
[5,]  0.8280422 -7.8272613
[6,] -6.2633738 -5.8725562
> head(result)
[1] 1292351.3 1661739.3  824260.0 1189476.4  496403.2  233480.2