如何计算R数据集中坐标与参考点的距离并生成公里单位新列

常见失败原因说明 #

之前调用包失败大概率是经纬度参数顺序传反：

• geosphere 包所有距离计算函数要求输入的坐标顺序为 (经度, 纬度)
• geodist 包默认坐标顺序为 (纬度, 经度)，符合常规数据存储习惯

方案1：使用 geosphere 包实现 #

步骤代码 #

# 安装加载包
install.packages("geosphere")
library(geosphere)

# 构造示例数据集，直接运行即可生成你提供的测试数据df
df <- structure(list(Day = c("26", "05", "17", "18", "19", "19"), Month = c("07", 
"08", "08", "08", "08", "08"), Year = c("2021", "2021", "2021", 
"2021", "2021", "2021"), Location.Receiver = c("Den Oever Ijsselmeer", 
"Medemblik Ijsselmeer, gemaal", "Den Oever Ijsselmeer", "Den Oever Ijsselmeer", 
"Den Oever Ijsselmeer", "Den Oever Ijsselmeer"), Transmitter = c("A69-1602-59776", 
"A69-1602-59777", "A69-1602-59776", "A69-1602-59776", "A69-1602-59769", 
"A69-1602-59776"), Batch.location = c("Den Oever", "Den Oever", 
"Den Oever", "Den Oever", "Den Oever", "Den Oever"), BatchNr = c(8, 
9, 8, 8, 1, 8), Latitude = c(52.92343, 52.76098, 52.92343, 52.92343, 
52.92343, 52.92343), Longitude = c(5.04127, 5.12172, 5.04127, 
5.04127, 5.04127, 5.04127), Date = structure(c(18834, 18844, 
18856, 18857, 18858, 18858), class = "Date")), row.names = c(1095729L, 
1180267L, 1072657L, 1092667L, 716601L, 1077415L), class = "data.frame")

# 定义参考点 注意顺序是(经度, 纬度)
ref_point <- c(5.04127, 52.92343)

# 计算距离（Haversine算法适用于球面距离计算，返回值单位为米，除以1000转公里）
df$distance_km <- distHaversine(df[, c("Longitude", "Latitude")], ref_point) / 1000

结果验证 #

运行后第二行和参考点的距离约为18.88公里，其余和参考点坐标完全一致的行距离为0，符合预期。

方案2：使用 geodist 包实现（代码更简洁） #

步骤代码 #

# 安装加载包
install.packages("geodist")
library(geodist)

# 直接计算，默认返回公里单位
df$distance_km <- geodist(
  x = df[, c("Latitude", "Longitude")],
  y = c(52.92343, 5.04127),
  measure = "haversine"
)

补充说明 #

如果需要更高精度的距离计算，可以把两种方案里的haversine算法替换为vincenty椭球算法，误差更小，适合短距离高精度场景。

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