如何在geopandas中获取LINESTRING Z几何的Z坐标并计算最大值

性能最优方案（优先使用） #

如果你的环境满足 shapely >= 2.0 + geopandas >= 0.10，可以直接使用内置接口，不需要遍历嵌套列表，性能提升明显：

• 提取全表所有点的全局最大Z值：

import shapely

# 批量提取所有坐标的三维值，取第三列Z值直接求最大值
all_z = shapely.get_coordinates(line.geometry, include_z=True)[:, 2]
max_z = all_z.max()

• 为每行几何单独计算该行轨迹的最大Z值，新增列存储：

line['max_agl'] = line.geometry.apply(lambda g: max(p[2] for p in g.coords))

单条几何快速提取方案 #

如果仅需要处理单个LINESTRING Z对象，有两种更简洁的写法：

1. 直接遍历坐标取最大值（可读性最高）

max_z = max(p[2] for p in single_line.coords)

不需要额外导入itemgetter依赖，代码更简洁。
2. 调用三维边界接口（性能最高）

# always_2d=False会返回三维边界：(xmin, ymin, zmin, xmax, ymax, zmax)
three_d_bound = shapely.bounds(single_line, always_2d=False)
max_z = three_d_bound[5]

该方法不需要遍历坐标点，直接读取几何内置的三维边界计算结果，效率最高。

原方案优化说明 #

你原来的实现逻辑是正确的，只是可以通过生成器表达式替代列表推导式减少中间内存占用，同时去掉不必要的依赖导入，代码会更轻量。

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