如何在Geopandas数据框中按经度排序几何列构建LineString

步骤说明与代码实现 #

• 提取经度值到新列：从POINT Z类型的geometry列中提取经度（即点的x坐标），为排序做准备：

  gdf['lon'] = gdf.geometry.x
  

• 按经度排序数据框：使用sort_values方法对数据框按经度升序（从西到东）或降序（从东到西）排序，确保点按经度顺序排列，避免构建线串时出现乱序折线：

  # 按经度升序排序（从西到东）
  gdf_sorted = gdf.sort_values(by='lon', ascending=True)
  
  # 若需要从东到西排序，设置ascending=False
  # gdf_sorted = gdf.sort_values(by='lon', ascending=False)
  

• 基于排序后的点构建LINESTRING：

  • 仅保留XY坐标的线串：

    from shapely.geometry import LineString
    lineStringObj = LineString([[point.x, point.y] for point in gdf_sorted.geometry.values])
    

  • 保留Z坐标的带Z线串（若后续需用到高程信息）：

    from shapely.geometry import LineStringZ
    lineStringObjZ = LineStringZ([[point.x, point.y, point.z] for point in gdf_sorted.geometry.values])
    

• 缓冲生成POLYGON：对构建好的线串执行缓冲操作，替换your_buffer_distance为实际需要的缓冲距离（单位需与数据坐标系一致）：

  polygonObj = lineStringObj.buffer(your_buffer_distance)
  

额外提示 #

如果排序后仍出现异常折线，需检查数据中是否存在经度跳变的异常点（比如跨180°经线的情况），这类场景需要额外处理坐标转换或数据清洗。

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