[英]How to calculate distance between 2 points with leaflet, geosphere, and R shiny
[英]Calculate minimum distance between GPS points using geosphere gives inaccuate values
我正在尝试计算 R 中一系列 GPS 点的最近邻点之间的平均距离。我找到了两个代码来获取这些值。 但它似乎没有给出以米为单位的正确距离。 当我对照谷歌地图检查时,它已经离开了。
我找到了这个答案: R - 在给定半径内查找最近的相邻点和邻居数量,坐标经纬度
library(geosphere)
sp.groups <- groups
coordinates(sp.groups) <- ~Long+Lat
class(sp.groups)
d<- distm(sp.groups)
min.d<- apply(d, 1, function(x) order(x, decreasing=F)[2])
min.d
mean(min.d)
groupdist<- cbind(groups, groups[min.d,], apply(d, 1, function(x) order(x, decreasing=F)[2]))
colnames(groupdist)<- c(colnames(groups), 'neighbor', 'n.lat', 'n.long','dist')
这里使用 package rgeos,但它给出了相同的结果: Calculate the distance between points of two datasets (nearest neighbor)
library(rgeos)
sp.groups <- groups
coordinates(sp.groups) <- ~Long+Lat
proj4string(sp.groups) <-CRS("+proj=utm +datum=WGS84")
class(sp.groups)
d<- gDistance(sp.groups, byid=TRUE)
min.d<- apply(d, 1, function(x) order(x, decreasing=F)[2])
min.d
mean(min.d)
groupdist<- cbind(groups, groups[min.d,], apply(d, 1, function(x) order(x, decreasing=F)[2]))
colnames(groupdist)<- c(colnames(groups), 'neighbor', 'n.lat', 'n.long','dist')
当我拨打 go 并在 Google 地球上查看时,距离可能会很远。 它甚至为 160-200 米的最近邻居提供了不同的值。 还有一些最近的邻居没有相同的距离值,请参见 K11 和 K3,然后是 K3 和 K11。 这是我得到的结果,我添加了来自 Google 地图的预期值:
Group Lat Long neighbor n.lat n.long dist GMaps
K1 -26.96538 21.80965 K34 -26.96503 21.80940 27 44
K10 -26.96575 21.81132 K1 -26.96538 21.80965 1 172
K11 -26.96249 21.81120 K3 -26.96387 21.81053 22 166
K24 -26.96033 21.81090 K11 -26.96249 21.81120 3 240
K3 -26.96387 21.81053 K11 -26.96249 21.81120 3 166
K34 -26.96503 21.80940 K1 -26.96538 21.80965 1 44
怎么了?
我的资料
groups<-data.frame(Group = c('K1', 'K10', 'K11', 'K24', 'K3', 'K34'),
Lat = c(-26.96538, -26.96575, -26.96249, -26.96033, -26.96387, -25.96503),
Longitude = c(21.80965, 21.81132, 21.81120, 21.80190, 21.81053, 21.80940))
我不确定你在这里试图计算什么,但列distance
只是指距离最小的点的位置/编号。 我添加了一个实际最小距离的数字,看起来像您预期的结果。
library(geodist, include.only = NULL)
library(sp, include.only = NULL)
groups <- data.frame(Group = c('K1', 'K10', 'K11', 'K24', 'K3', 'K34'),
Lat = c(-26.96538, -26.96575, -26.96249, -26.96033, -26.96387, -25.96503),
Long = c(21.80965, 21.81132, 21.81120, 21.80190, 21.81053, 21.80940))
sp.groups <- groups
sp::coordinates(sp.groups) <- ~Long+Lat
# mindistance Matrix
d <- geodist::geodist(groups, measure = "cheap")
# position of minimum distance
diag(d) <- Inf
min.d <- max.col(-d)
min.d
#> [1] 2 1 5 5 3 4
groupdist <- cbind(groups, groups[min.d,], min.d)
colnames(groupdist) <- c(colnames(groups), 'neighbor', 'n.lat', 'n.long','closest_to')
# get minimum distance for each pair of coordinates
groupdist$distance <- d[cbind(seq_along(min.d), min.d)]
groupdist
#> Group Lat Long neighbor n.lat n.long closest_to distance
#> 2 K1 -26.96538 21.80965 K10 -26.96575 21.81132 2 171.1554
#> 1 K10 -26.96575 21.81132 K1 -26.96538 21.80965 1 171.1554
#> 5 K11 -26.96249 21.81120 K3 -26.96387 21.81053 5 167.2225
#> 5.1 K24 -26.96033 21.80190 K3 -26.96387 21.81053 5 944.4119
#> 3 K3 -26.96387 21.81053 K11 -26.96249 21.81120 3 167.2225
#> 4 K34 -25.96503 21.80940 K24 -26.96033 21.80190 4 110613.1362
由 reprex package (v2.0.1) 创建于 2021-08-22
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.