在R中找到4个坐标的中点坐标

Question

首先感谢大家对我之前的问题提供的所有帮助。

我试图在下面寻找一个可以找到我的数据集中心坐标的函数：

df <- read.table(sep=",", col.names=c("Start_Latitude","Start_Longitude","End_Latitude", "End_Longitude"),text="43.9567343,  -78.8571382, 43.9399364, -78.8497342")

Start_Latitude    Start_Longitude   End_Latitude    End_Longitude      
  43.9567343       -78.8571382       43.9399364      -78.8497342

我试过这段代码，但它只计算从起始坐标到结束坐标的距离，我想找到中点坐标而不仅仅是距离。我的意思是我需要找到中点的经纬度。

#dd the distance as column in the dataframe
df1$dist <- distm(x = df[, c('Start_Longitude', 'Start_Longitude')], 
                y = df[, c('End_Longitude', 'End_Latitude')],
                fun = distHaversine
)

我有这个数据集：

 Start_Latitude    Start_Longitude   End_Latitude    End_Longitude       dist
      43.9567343       -78.8571382       43.9399364      -78.8497342     13669708

有没有办法找到中点的纬度和经度？另外，我如何计算以公里为单位的距离。

提前谢谢你！

Answer 1

这是geosphere::midPoint 的一种方法：

library(geosphere)
data.frame(df,
           dist = distHaversine(df[, c('Start_Longitude', 'Start_Latitude')],
                                df[, c('End_Longitude', 'End_Latitude')]) / 1000, 
           midPoint(p1 = df[, c('Start_Longitude', 'Start_Latitude')],
                    p2 = df[, c('End_Longitude', 'End_Latitude')]))
#  Start_Latitude Start_Longitude End_Latitude End_Longitude    dist       lon      lat
#1       43.95673       -78.85714     43.93994     -78.84973 1.96183 -78.85344 43.94834

Answer 2

可以通过平均纬度/经度组合来计算中点。使用原始帖子中的数据，我们使用dplyr::mutate() 来计算中点。

textFile <- "Start_Latitude    Start_Longitude   End_Latitude    End_Longitude      
  43.9567343       -78.8571382       43.9399364      -78.8497342
41.947460   -87.654730   41.890620   -87.624480"

data <- read.table(text=textFile,header=TRUE)

library(dplyr)
data %>% rowwise(.) %>% 
    mutate(midpoint_lat = sum(Start_Latitude,End_Latitude)/2,
           midpoint_long = sum(Start_Longitude,End_Longitude)/2) %>%
    as.data.frame(.) # print as df, not tibble

...和输出：

  Start_Latitude Start_Longitude End_Latitude End_Longitude id midpoint_lat
1       43.95673       -78.85714     43.93994     -78.84973  1     43.94834
2       41.94746       -87.65473     41.89062     -87.62448  2     41.91904
  midpoint_long
1     -78.85344
2     -87.63961

以公里为单位计算距离

正如geosphere 包中的distHaversine() 函数的帮助中所述，结果中的测量单位是根据为地球半径给出的参数返回的。默认值为r = 6378137 米，因此以米为单位返回测量单位。

要换算成公里，有两种选择：

1. 将结果除以 1,000（如本问题的其他答案所示），或
2. 将参数r更改为6378.137