如何将含属性的Pandas点序列转为Shapely LineString并计算衍生属性

针对你的需求，推荐用Pandas分组+Shift向量化操作来处理，既避免了逐行遍历的低效，也解决了单纯groupby处理跨行属性计算的不便，代码简洁且性能优异。

核心思路 #

1. 确保每个Path分组内的点按locIdx有序（点序列的顺序是生成线的基础）
2. 用shift(-1)获取同一分组内下一个点的几何信息和Arr时间
3. 过滤掉每个分组的最后一行（无后续点，无法生成线）
4. 用当前点和下一个点的坐标构造LineString
5. 计算衍生属性timeDiff，并保留当前点的原有属性

代码实现（基于GeoPandas/Pandas） #

假设你的原始数据是GeoDataFrame（或普通DataFrame，需提前解析POINT几何）：

import pandas as pd
import geopandas as gpd
from shapely.geometry import LineString

# 1. 按Path分组并按locIdx排序（确保点序列顺序正确）
sorted_df = df.groupby('Path').apply(lambda grp: grp.sort_values('locIdx')).reset_index(drop=True)

# 2. 获取下一个点的几何和Arr时间
sorted_df['next_geom'] = sorted_df.groupby('Path')['geometry'].shift(-1)
sorted_df['next_Arr'] = sorted_df.groupby('Path')['Arr'].shift(-1)

# 3. 过滤掉无后续点的行（每个Path的最后一行）
line_data = sorted_df[sorted_df['next_geom'].notna()].copy()

# 4. 生成LineString几何列
line_data['geometry'] = line_data.apply(
    lambda row: LineString([row['geometry'].coords[0], row['next_geom'].coords[0]]),
    axis=1
)

# 5. 计算衍生属性并整理最终结果
line_data['timeDiff'] = line_data['next_Arr'] - line_data['Dep']
# 保留需要的属性列，可根据需求调整
final_lines = line_data[['Path', 'locIdx', 'PostLength', 'timeDiff', 'geometry']]

方案优势 #

• 性能高效：Pandas的向量化操作远快于逐行遍历，数据量越大优势越明显
• 逻辑清晰：用shift直接获取跨行属性，解决了groupby难以处理的跨行计算问题
• 可维护性强：代码符合Pandas惯用写法，后续调整属性或逻辑更便捷

如果你的原始数据不是GeoDataFrame，只需先把geometry列解析为Shapely POINT对象即可，比如：

from shapely.wkt import loads
df['geometry'] = df['geometry'].apply(loads)

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