如何用geopandas计算LineString与Point之间的米/千米单位距离？

核心原理 #

你当前计算得到的距离单位为度，原因是输入几何使用的是WGS84地理坐标系（EPSG:4326，坐标为经纬度），只有转换为以米为单位的投影坐标系，或直接计算椭球面测地距离，才能得到米/千米单位的结果。

方案1：转换投影坐标系计算（推荐，性能更高） #

• 步骤1：确认当前CRS
  先验证几何的坐标系是否为WGS84：
  print(line.crs)
  若输出为EPSG:4326则符合预期。

• 步骤2：转换为适配研究区域的投影坐标系
  你提供的坐标属于欧洲区域，适配的UTM投影为EPSG:32632（UTM 32N，单位为米），转换代码如下：

  # 线几何转投影
  line_proj = line.to_crs(epsg=32632)
  # 点几何转投影
  point_proj = point.to_crs(epsg=32632)
  

  注意：若研究区域为其他地区，请替换为对应区域的官方/等距投影代码，避免计算误差

• 步骤3：计算距离并转换单位

  # 计算得到的距离单位为米
  distance_m = line_proj.distance(point_proj)
  # 转换为千米
  distance_km = distance_m / 1000
  

方案2：直接计算测地距离（无需选投影，适合小批量计算） #

如果不想手动匹配投影，可以直接计算WGS84椭球面上的真实距离：

from shapely.ops import nearest_points
from geopy.distance import geodesic

# 先找到线到点的最近点坐标
nearest_line_point = nearest_points(line, point)[0]
# 生成(纬度, 经度)格式的坐标对
p1 = (nearest_line_point.y, nearest_line_point.x)
p2 = (point.y, point.x)
# 直接得到千米单位的距离
distance_km = geodesic(p1, p2).km
# 转米的话取geodesic(p1,p2).m即可

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