k最近的邻居在球体上

Question

我正在尝试使用 R 找到球体上的 k 最近邻。

由于我要处理数百万个点，因此无法使用蛮力。

有谁知道我该怎么做？谢谢

Answer 1

使用 Ben Bolker 提供的 link 我设法解决了我的问题。

lonlat2xyz=function (lon, lat, r) 
{
lon = lon * pi/180
lat = lat * pi/180
if (missing(r)) 
    r <- 6378.1
x <- r * cos(lat) * cos(lon)
y <- r * cos(lat) * sin(lon)
z <- r * sin(lat)
return(cbind(x, y, z))
}

lon1=runif(100,-180,180);lon2=runif(100,-180,180);lat1=runif(100,-90,90);lat2=runif(100,-90,90)

xyz1=lonlat2xyz(lon1,lat1)
xyz2=lonlat2xyz(lon2,lat2)

library(nabor)

out=knn(data=xyz1,query = xyz2,k=20)

library(maps)

map()
points(lon1,lat1,pch=16,col="black")
points(lon2[1],lat2[1],pch=16,col="red")
points(lon1[out$nn.idx[1,]],lat1[out$nn.idx[1,]],pch=16,col="blue")

注意 knn 函数给出的距离不是地理距离！

k最近的邻居在球体上

k nearest neighbor on a sphere

r

nearest-neighbor

k最近的邻居**在球体上**

k nearest neighbor **on a sphere**

r

nearest-neighbor

k最近的邻居在球体上

k nearest neighbor on a sphere