R Studio中如何批量计算数据框点位到固定位置的距离

报错原因 #

c(df$longitude, df$latitude)会将经纬度两列的所有值拼接为长度为2倍行数的一维长向量，不符合geosphere距离函数的输入格式要求：geosphere的距离计算函数要求点位输入为经度在前、纬度在后的两列结构（矩阵或数据框），每一行对应一组独立点位坐标，因此无法识别传入的长向量，无法完成逐行计算。

注意：不要误用R基础包自带的dist()函数，该函数仅用于计算平面欧式距离，不适用于球面经纬度坐标的地理距离计算，算出的结果没有实际地理意义。

正确实现代码 #

geosphere内置的distGeo()（高精度椭球面距离）、distHaversine()（快速Haversine公式大圆距离）原生支持批量点位计算，不需要写for循环或者apply逐行迭代，直接传入坐标矩阵即可一次性算出所有点位到中心点的距离，直接追加为数据框新列：

# 先加载geosphere包，未安装的话先运行install.packages("geosphere")
library(geosphere)

# 核心计算逻辑：把经纬度列转成矩阵传入即可
# 注意c("longitude", "latitude")要替换成你数据框里实际的经纬度列名
df$dist_to_center <- distGeo(
  p1 = citycenter,
  p2 = as.matrix(df[, c("longitude", "latitude")])
)

函数默认返回的距离单位为米，如果需要公里为单位，直接对计算结果做单位转换即可：

df$dist_to_center_km <- df$dist_to_center / 1000

常见踩坑提示 #

• 坐标顺序不能搞反：geosphere所有距离函数都要求参数顺序为「经度、纬度」，顺序写反会得到完全错误的计算结果
• 坐标系要匹配：如果你的经纬度不是WGS84坐标系（比如国内高德/谷歌地图用的GCJ02、百度地图用的BD09），需要先将坐标统一转换为WGS84再计算，否则距离结果会存在几十到几百米的偏差
• 不要做低效逐行计算：很多人遇到批量计算第一反应是用apply按行传参，实际上distGeo原生支持矩阵批量输入，计算效率比逐行迭代高几个数量级，十万级点位也能秒出结果

内容的提问来源于stack exchange，提问作者Arthur Lo