导出含LineString几何的GeoDataFrame至Shapefile时出现错误

问题分析 #

你遇到的'Series' object has no attribute 'has_z'或'geom_type'错误，本质是GeoDataFrame的geometry列不再是合法的GeoSeries类型，导致to_file方法无法识别几何属性。虽然Kepler能可视化（它支持解析WKT字符串），但GeoPandas导出Shapefile要求几何列必须是包含shapely几何对象的GeoSeries。

解决方案步骤 #

1. 排查几何列的类型

执行以下代码确认问题：

import geopandas as gpd

# 检查geometry列的类型
print(type(test_gdf['geometry']))
# 检查列内元素的类型
print(test_gdf['geometry'].apply(type).unique())

如果输出不是geopandas.geoseries.GeoSeries，或者元素类型不是shapely.geometry.linestring.LineString，说明几何列已被损坏。

2. 修复WKT格式的几何列

如果你的geometry列是类似LINESTRING (...)的WKT字符串，需要转换为shapely几何对象并重新初始化GeoDataFrame：

from shapely.wkt import loads

# 将WKT字符串转换为shapely几何对象
test_gdf['geometry'] = test_gdf['geometry'].apply(loads)
# 重新构造GeoDataFrame并指定几何列
test_gdf = gpd.GeoDataFrame(test_gdf, geometry='geometry')
# 设置正确的CRS（示例为WGS84，根据你的数据修改）
test_gdf.crs = "EPSG:4326"

3. 修复合并后丢失的GeoSeries属性

如果是合并多个GeoDataFrame后几何列丢失了GeoSeries特性，直接重新初始化即可：

test_gdf = gpd.GeoDataFrame(test_gdf, geometry='geometry')
# 确认CRS设置，若无则补充
if test_gdf.crs is None:
    test_gdf.crs = "EPSG:4326"  # 替换为你的数据坐标系

4. 验证修复结果

执行以下代码确认修复成功：

print(test_gdf.geom_type.unique())  # 应输出统一几何类型（如['LineString']）
print(test_gdf.crs)  # 应显示正确的坐标系

确认无误后再执行导出代码：

output_shapefile_path = r'C:\Users\User.Name\Desktop\Folder\combined_data.shp'
test_gdf.to_file(output_shapefile_path)

额外注意事项 #

• 合并多个Shapefile时，确保所有源数据的几何类型一致（Shapefile不支持混合几何类型）。
• 处理过程中优先使用GeoPandas的合并方法，避免用普通pd.concat直接合并，若用了则合并后重新指定几何列。

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