英国国家电网中两点间距离的计算

Question

我想以公里为单位计算数据集中每一行中两组英国国家网格坐标之间的距离(即我的示例中每一行的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]    

Answer 1

当我问我的问题时，我发现了答案。然而，我想我还是会贴出这个问题和我的答案，以防它在将来对其他人有帮助。

Points_sf <- st_as_sf(Points, coords = c("AreaA_X","AreaA_Y"),crs=27700)            
Points_sf$geometry2 <- st_as_sf(Points[,c("AreaB_X","AreaB_Y")], coords = c("AreaB_X","AreaB_Y"),crs=27700) 
Points_sf$Distance_M <- st_distance(Points_sf$geometry, Points_sf$geometry2, by_element = TRUE)