Geopandas对带Z坐标的LineString Geoseries缓冲丢失Z坐标如何解决

Geopandas的buffer方法底层依赖GEOS库的二维缓冲实现，GEOS库的原生二维缓冲操作默认会丢弃Z、M维度坐标，这不是Geopandas或Fiona的配置遗漏，属于底层依赖的能力限制，当前版本没有内置参数可以保留Z值，也不原生支持三维空间的缓冲计算。

根据你的实际需求可以选择两种处理方案：

方案1：XY平面缓冲后保留原始Z值 #

如果仅需要缓冲后的几何携带Z值，不需要Z方向的缓冲效果，可以手动将原始线要素的Z值赋值给缓冲结果的坐标点，参考实现代码：

import geopandas as gpd
import shapely

def buffer_preserve_z(geo_series, buffer_dist):
    # 先执行原生二维缓冲
    buffered = geo_series.buffer(buffer_dist)
    # 逐要素补全Z值
    for idx, raw_geom in enumerate(geo_series):
        # 可根据业务需求调整Z值逻辑，示例取原始线所有节点的平均Z
        z_list = [coord[2] for coord in raw_geom.coords]
        target_z = sum(z_list) / len(z_list)
        # 为缓冲后的所有坐标点补充Z值
        buffered_geom = shapely.transform(
            buffered.iloc[idx],
            lambda coords: [(x, y, target_z) for x, y in coords]
        )
        buffered.iloc[idx] = buffered_geom
    return buffered

# 调用示例
myGeoseries_with_z = buffer_preserve_z(myGeoseries, 2.5)

该方案输出的是带统一Z值的POLYGON Z类型几何，可直接在Geopandas中正常使用。

方案2：执行真正的三维空间缓冲 #

如果需要同时对X、Y、Z三个方向做缓冲生成三维实体，可使用GDAL/OGR库提供的Buffer3D方法实现，参考代码：

from osgeo import ogr
import shapely

def buffer_3d(geom, xy_buffer, z_buffer):
    # 转OGR几何对象
    ogr_geom = ogr.CreateGeometryFromWkt(geom.wkt)
    # 执行三维缓冲，第二个参数为XY平面缓冲距离，第三个参数为Z方向缓冲距离
    buffered_3d = ogr_geom.Buffer3D(xy_buffer, z_buffer)
    # 转回shapely几何对象
    return shapely.from_wkt(buffered_3d.ExportToWkt())

# 调用示例，XY方向缓冲2.5，Z方向缓冲2.5，可按需调整两个方向的缓冲距离
myGeoseries_3d = myGeoseries.apply(lambda x: buffer_3d(x, 2.5, 2.5))

该方案输出的是三维多面体几何，需注意当前版本Geopandas对三维多面体的可视化、空间计算支持有限，仅适合需要存储、导出三维缓冲结果的场景使用。

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