在2个单独的矩阵中查找坐标的最近邻点以及2之间的距离

Question

我有2个数据框，它们只是2维（纬/长）矩阵。这两个数据帧看起来像下面的输入：

latitude longitude 
27.78833 -82.28197 
27.79667 -82.29294 

我们称之为“dfref”和“dfnew”。我想在dfref中的每个点找到dfnew中的最近点以及以米为单位的2个点之间的距离。

输出应该是这样的：

dr.latitude dr.longitude dn.latitude dn.longitude dist 
27.78833  -82.28197  27.54345  -82.33233  162.34 
27.79667  -82.29294  27.56543  -82.12323  232.23 

我已经使用了类包的KNN功能和Searchtrees包试过，但我的脚本只发现在dfref矩阵最近的点，我不知道如何添加测量。

knn1(train=cbind(dfref), test=cbind(dfnew), cl=seq_len(nrow(dfnew))) 

是否有确实有效且我怎样才能到一个脚本这一个功能呢？

Answer 1

我不是地球数学专家，但似乎你可以像这样开始：

dfref <- read.table(text = 
"latitude longitude 
27.78833 -82.28197 
27.79667 -82.29294", header = T) 
dtref <- data.table(dfref) 

dfnew <- read.table(text = 
"latitude longitude 
27.54345  -82.33233", header = T) 
dtnew <- data.table(dfnew) 

# Make cartesian product of to tables. 
dtref$fake <- 1 
dtnew$fake <- 1 
dtall <- merge(dtref, dtnew, by = "fake", allow.cartesian = T) 

# Calculate distance. 
library(geosphere) 
dtall[, distance := distVincentyEllipsoid(c(longitude.x, latitude.x), c(longitude.y, latitude.y)), by = 1:nrow(dtall)] 

# Print results. 
dtall[, .(latitude.x, longitude.x, latitude.y, longitude.y, distance)] 

#  latitude.x longitude.x latitude.y longitude.y distance 
# 1: 27.78833 -82.28197 27.54345 -82.33233 27587.29 
# 2: 27.79667 -82.29294 27.54345 -82.33233 28328.19