如何将给定经纬度转换为KM单位出行距离，distGeo函数是否适用？

结论 #

distGeo是非常合适的选择，完全匹配你的需求。

原因说明 #

• 它是R语言geosphere包中专门用于计算地表两点间大地距离的函数，基于WGS84椭球体模型实现，精度比传统球面模型（比如Haversine公式）更高，适合所有常规的经纬度距离计算场景。
• 函数默认返回值单位为米，只需将结果除以1000即可得到你需要的公里单位数据，转换成本极低。
• 支持批量向量化运算，无需循环即可直接处理整表的上万条出行路线数据，运行效率高。

示例代码 #

用你给出的示例坐标演示运算逻辑：

# 加载依赖包
library(geosphere)

# 单条路线计算示例
start_point <- c(-87.68000, 42.08000) # 注意顺序：先经度，后纬度
end_point <- c(-87.69000, 42.07000)
distance_km <- distGeo(start_point, end_point) / 1000
# 本次示例计算结果约为1.35公里

# 整表批量计算示例
# 假设你的数据存储在数据框df中
df$distance_km <- distGeo(
  df[, c("start_lng", "start_lat")],
  df[, c("end_lng", "end_lat")]
) / 1000

注意事项 #

• 一定要注意参数的传入顺序：distGeo要求每个坐标点的第一个值为经度，第二个为纬度，顺序传反会得到完全错误的结果。
• 如果你的经纬度数据不是WGS84坐标系，需要先完成坐标转换再调用该函数计算，否则会有明显误差。

内容的提问来源于stack exchange，提问作者Maria Phoonsin