GeoPandas按长度分割线要素且保留属性列的方法

你已有的单要素点分割逻辑可以直接复用，批量处理时核心要避免直接对geometry列做apply后直接返回，否则会丢失非空间属性。正确的实现思路是逐行处理要素，拆分得到子线段后直接复制原行的全部属性，最后重建成新的GeoDataFrame即可。

依赖导入 #

import geopandas as gpd
import numpy as np
from shapely.geometry import LineString
from shapely.ops import split, snap

核心函数实现 #

1. 单条线按长度阈值拆分 #

兼容你之前写的snap容差逻辑，自动沿线路径按固定间隔生成分割点，返回拆分后的子线段列表：

def split_line_by_length(line, max_length, tolerance: float=1.0e-12):
    # 长度未超阈值直接返回原要素
    if line.length <= max_length:
        return [line]
    # 沿线路径生成所有分割点
    split_points = []
    current_dist = max_length
    while current_dist < line.length:
        split_points.append(line.interpolate(current_dist))
        current_dist += max_length
    # 逐点迭代分割
    segments = [line]
    for pt in split_points:
        temp_segs = []
        for seg in segments:
            if seg.length > max_length:
                split_res = split(snap(seg, pt, tolerance), pt)
                temp_segs.extend(list(split_res.geoms))
            else:
                temp_segs.append(seg)
        segments = temp_segs
    # 过滤长度为0的无效碎片
    return [seg for seg in segments if seg.length > tolerance]

2. 批量处理GeoDataFrame、保留全属性 #

遍历原始表每一行，拆分后每个子线段都继承原行的所有字段值，最后重建为带坐标系的GeoDataFrame：

def split_gdf_lines(gdf, max_length, geom_col='geometry'):
    output_rows = []
    for _, row in gdf.iterrows():
        sub_lines = split_line_by_length(row[geom_col], max_length)
        for sub_line in sub_lines:
            new_row = row.copy()
            new_row[geom_col] = sub_line
            output_rows.append(new_row)
    return gpd.GeoDataFrame(output_rows, crs=gdf.crs, geometry=geom_col)

使用示例 #

针对你给出的测试数据，设置阈值为两点间距√2即可得到你要的拆分结果：

# 构造测试数据
test_gdf = gpd.GeoDataFrame([
    {"Feature":'Hi',"ID":1,"geometry":LineString([[1,1],[2,2],[3,3]])},
    {"Feature":'bye',"ID":2,"geometry":LineString([[10,10],[20,20],[30,30]])}
], crs="EPSG:3857")

# 执行拆分
result = split_gdf_lines(test_gdf, max_length=np.sqrt(2))

输出结果结构完全匹配需求：

• Feature、ID等所有原始非空间字段完整保留，子线段自动继承对应原线的属性值
• 拆分后每条子线段长度均不超过设定阈值
• 自带容差处理，不会因为浮点精度问题出现分割失败、多生成无效碎片的问题

注意：长度阈值的单位和你数据坐标系的单位保持一致，如果是经纬度地理坐标系，建议先投影到平面坐标系（单位为米）再做拆分，避免长度计算偏差。

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