Я хотел бы рассчитать расстояние между двумя наборами координат Briti sh National Grid в каждой строке в наборе данных (то есть расстояние между AreaA и AreaB для каждой строки в моем примере) в километрах.
library("sf")
BBox <- st_bbox(c(xmin = -10000, xmax = 15000, ymax = 20000, ymin = -5000), crs = st_crs(27700)) %>% st_as_sfc()
P1 <- st_sample(BBox, 50, exact = TRUE)
P2 <- st_sample(BBox, 50, exact = TRUE)
Points <- data.frame(AreaA=sample(110:190,50,replace=FALSE),AreaB=sample(10:90,50,replace=FALSE))
Points$AreaA_X <- data.frame(st_coordinates(P1))[,1]
Points$AreaA_Y <- data.frame(st_coordinates(P1))[,2]
Points$AreaB_X <- data.frame(st_coordinates(P2))[,1]
Points$AreaB_Y <- data.frame(st_coordinates(P2))[,2]