如何使用Geopandas提取LineString起止XY坐标并生成对应新列

首先确认你的运行环境已经安装geopandas，操作全程不需要额外依赖其他库，针对标准LineString类型的线要素可以直接按如下流程处理：

• 读入目标shapefile数据

  import geopandas as gpd
  
  # 括号内替换为你本地的折线shapefile路径
  line_gdf = gpd.read_file("your_polyline_layer.shp")
  

• 批量计算四个坐标字段
  线要素的geometry对象自带的coords属性会按绘制顺序存储所有折点的坐标元组，索引0对应线的起点，索引-1对应线的终点，直接通过字段赋值生成新列即可：

  # 写入起点坐标
  line_gdf["StartX"] = line_gdf.geometry.apply(lambda x: x.coords[0][0])
  line_gdf["StartY"] = line_gdf.geometry.apply(lambda x: x.coords[0][1])
  
  # 写入终点坐标
  line_gdf["EndX"] = line_gdf.geometry.apply(lambda x: x.coords[-1][0])
  line_gdf["EndY"] = line_gdf.geometry.apply(lambda x: x.coords[-1][1])
  

注意：如果你的数据集里包含MultiLineString（多段线）类型要素，直接运行上述代码会触发属性报错，需要先把多段线拆分为单段线再做提取，拆分代码如下：

# 拆分多线要素为独立单线要素
line_gdf = line_gdf.explode(index_parts=True).reset_index(drop=True)

• 结果导出
  字段添加完成后，可以直接将处理后的数据导出为新的矢量文件：

  line_gdf.to_file("polyline_with_start_end_coords.shp", encoding="utf-8")
  

处理完成后打开属性表就能看到新增的StartX、StartY、EndX、EndY四个字段，对应值就是每条线的起终点平面坐标，和数据本身的坐标系保持一致。

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