在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