对于一个数据集中的每个点，计算到第二个数据集中最近点的距离

Question

试图为SpatialPointsDataFrame 中的每个点找到一秒钟SpatialPointsDataFrame 到最近点的距离（相当于ArcGIS 中两个SpatialPointDataFrames 的“最近”工具）。

我可以通过使用gDistance 计算所有成对距离并采用min (like answer 1 here) 来进行简单的实现，但我有一些庞大的数据集，并且正在寻找更有效的方法。

例如，这里是trick with knearneigh for points in same dataset。

在r-sig-geo上交叉发布

Answer 1

SearchTrees 包提供了一种解决方案。引用其文档，它“提供了 QuadTree 数据结构的实现 [它] 用于在二维中实现快速 k-最近邻 [...] 查找。”

您可以使用它来快速查找SpatialPoints 对象b 中的每个点，第二个SpatialPoints 对象B 中最近的两个点

library(sp)
library(SearchTrees)

## Example data
set.seed(1)
A <- SpatialPoints(cbind(x=rnorm(100), y=rnorm(100)))
B <- SpatialPoints(cbind(x=c(-1, 0, 1), y=c(1, 0, -1)))

## Find indices of the two nearest points in A to each of the points in B
tree <- createTree(coordinates(A))
inds <- knnLookup(tree, newdat=coordinates(B), k=2)

## Show that it worked
plot(A, pch=1, cex=1.2)
points(B, col=c("blue", "red", "green"), pch=17, cex=1.5)
## Plot two nearest neigbors
points(A[inds[1,],], pch=16, col=adjustcolor("blue", alpha=0.7))
points(A[inds[2,],], pch=16, col=adjustcolor("red", alpha=0.7))
points(A[inds[3,],], pch=16, col=adjustcolor("green", alpha=0.7))

Answer 2

R-Sig-Geo 的另一个建议是 nabor 库中的 knn 函数。