Question

我有两个坐标列表， mapped_coords，unmapped_coords ，它们都是坐标列表。

我想取unmapped_coords并为每个元素返回mapped_coord中最小距离的点的索引。

> head(mapped_coords)
[[1]]
[1] -79.2939  43.8234

[[2]]
[1] -79.7598  43.4381

[[3]]
[1] -79.4569  43.6693

[[4]]
[1] -81.2472  42.9688

[[5]]
[1] -79.1649  43.8073

[[6]]
[1] -79.7388  43.6753

 str(mapped_coords)
List of 62815
 $ : num [1:2] -79.3 43.8
 $ : num [1:2] -79.8 43.4
 $ : num [1:2] -79.5 43.7

使用geosphere包我可以使用distHaversine计算一对的距离，但我不知道如何在整个列表中进行。

> distHaversine(unlist(unmapped_coords[1]), unlist(mapped_coords[1]))
[1] 100594.6

Answer 1

您可以使用geosphere::distm创建一个距离矩阵，您可以使用which.min找到最小列（除了对角线，这是无用的）：

l <- list(c(-79.2939, 43.8234), 
          c(-79.7598, 43.4381), 
          c(-79.4569, 43.6693), 
          c(-81.2472, 42.9688), 
          c(-79.1649, 43.8073), 
          c(-79.7388, 43.6753))

m <- geosphere::distm(do.call(rbind, l))
diag(m) <- NA

apply(m, 1, which.min)
#> [1] 5 6 1 2 1 3

如果您有第二个距离列表，请将其作为第二个参数传递给distm，使对角线有用。由于没有NA s，因此可以使用max.col(-m)计算最小列。

Answer 2

您可以将distHaversine作为输入提供一对坐标和一个坐标矩阵（带有2列），这将返回一个距离的向量，该距离的长度与矩阵中的行数相同。您可以使用lapply

遍历未映射坐标列表

数据：

mapped_coord = list(c(-79.29,43.82),c(-79.76,43.44))
[[1]]
[1] -79.29  43.82

[[2]]
[1] -79.76  43.44

unmapped_coord = list(c(-79.16,43.12),c(-80.52,42.95))
[[1]]
[1] -79.16  43.12

[[2]]
[1] -80.52  42.95

方法：

library(geosphere)
## Transform the list of mapped coordinates into a matrix
mat = do.call(rbind,mapped_coord)
      [,1]  [,2]
[1,] -79.29 43.82
[2,] -79.76 43.44
## Find the coordinates with the min distances
lapply(unmapped_coord,function(x) which.min(distHaversine(x,mat)))
[[1]]
[1] 2

[[2]]
[1] 2

计算列表上的距离

2 个答案: