使用R的sf库计算UTM坐标对间公里距离时遇报错求解

错误原因分析 #

1. 坐标系统设置错误：你的坐标是UTM 32N（单位为米），但代码里错误指定了crs="4326"（WGS84经纬度，单位为度），这会导致坐标解析错误，距离计算完全不准确。UTM 32N对应的EPSG代码是32632。
2. st_distance参数使用错误：by_element=TRUE需要传入两个sf对象（分别对应起点和终点），用来计算两个对象中对应位置元素的距离。你只传入了一个包含起点的sf对象，所以触发报错。

修正后的代码 #

library(sf)

# 原始数据框
df <- data.frame(
  id = c(1,2,3),
  x_origin = c(642683.2, 373775, 383881),
  y_origin = c(5082920, 4997274, 4994504),
  x_dest = c(642683.3, 1126050, 942763.9),
  y_dest = c(5082920, 4374481, 4534235)
)

# 创建起点sf对象，设置正确的UTM 32N坐标系（EPSG:32632）
points_origin <- st_as_sf(df, coords = c("x_origin", "y_origin"), crs = 32632)

# 创建终点sf对象，使用相同坐标系
points_dest <- st_as_sf(df, coords = c("x_dest", "y_dest"), crs = 32632)

# 计算每组起点到终点的距离，by_element=TRUE对应逐元素计算
distances_m <- st_distance(points_origin, points_dest, by_element = TRUE)

# 转换为公里
distances_km <- as.numeric(distances_m) / 1000

# 把结果合并到原始数据框
df$distance_km <- distances_km

print(df)

代码说明 #

• 分别创建起点和终点的sf对象，确保两者使用相同的UTM坐标系，这样距离计算的单位是米。
• st_distance传入两个sf对象并设置by_element=TRUE，会返回长度为3的距离向量，对应每组起点到终点的距离。
• 最后将米转换为公里，方便阅读。

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