R-链接两组不同的坐标

R-链接两组不同的坐标,r,coordinates,R,Coordinates,我有两个数据集——让我们称它们为“绘图”(734行)和“站点”(62行)。不久前,我发现这段代码应该允许我根据坐标将每个“绘图”链接到最近的“站点” 数据集有点像这样-(但没有Long和Lat的标题) 绘图-长纬度站-长纬度 13.2 60.5 14.6 55.4 15.4 62.6 15.5 62.9 15.6 62.4 16.4 58.9 16.5 58.7

我有两个数据集——让我们称它们为“绘图”(734行)和“站点”(62行)。不久前,我发现这段代码应该允许我根据坐标将每个“绘图”链接到最近的“站点”

数据集有点像这样-(但没有Long和Lat的标题)

绘图-长纬度站-长纬度
13.2  60.5                   14.6  55.4
15.4  62.6                   15.5  62.9
15.6  62.4                   16.4  58.9
16.5  58.7                   19.3  64.0
16.5  58.5
#将结果打印到“results.csv”
接收器(“results.csv”)
#识别每个数据集的long+lat坐标

p_坐标你的坐标参考系是什么?这些地点在斯堪的纳维亚半岛吗? 无论如何,您可以使用geosphere软件包并使用
distHaversine
distVincentyEllipsoid
(更精确)获取距离:

plot <- data.frame(Lon = c(13.2,15.4,15.6,16.5,16.5),
                   Lat = c(60.5,62.6,62.4,58.7,58.5))

stations <- data.frame(Lon = c(14.6,15.5,16.4,19.3),
                       Lat = c(55.4,62.9,58.9,64))

p_coord <- SpatialPoints(plot[,c(1,2)])  
s_coord <- SpatialPoints(stations[,c(1,2)])

library(geosphere)
apply(p_coord@coords, 1, function(x) {
  which.min(distHaversine(p1 = x, p2 = s_coord@coords))
})

这意味着地块1靠近桩号3,地块2链接到桩号2,依此类推。

您的坐标参考系是什么?这些地点在斯堪的纳维亚半岛吗? 无论如何,您可以使用geosphere软件包并使用
distHaversine
distVincentyEllipsoid
(更精确)获取距离:

plot <- data.frame(Lon = c(13.2,15.4,15.6,16.5,16.5),
                   Lat = c(60.5,62.6,62.4,58.7,58.5))

stations <- data.frame(Lon = c(14.6,15.5,16.4,19.3),
                       Lat = c(55.4,62.9,58.9,64))

p_coord <- SpatialPoints(plot[,c(1,2)])  
s_coord <- SpatialPoints(stations[,c(1,2)])

library(geosphere)
apply(p_coord@coords, 1, function(x) {
  which.min(distHaversine(p1 = x, p2 = s_coord@coords))
})
这意味着1号地块靠近3号站,2号地块与2号站相连,依此类推。

“4号和5号地块与4号站相连。”。。不应该是3号站吗?“4号和5号地块与4号站相连。”。。不应该是3号站吗?
[1] 3 2 2 3 3