如何使用Geopandas/Shapely将多边形转换为直线？

核心思路很直接：多边形的边界本身就是线性几何要素，我们只需要提取多边形的外边界和内洞（如果存在），再将其转换为LineString或MultiLineString类型即可。

具体实现步骤 #

1. 读取GeoJSON数据：用GeoPandas加载多边形数据集
2. 编写转换函数：处理单个多边形/多多边形，提取边界并转换为线要素
3. 批量转换几何列：将函数应用到整个GeoDataFrame的几何列
4. 验证并保存结果（可选）

完整代码示例 #

import geopandas as gpd
from shapely.geometry import LineString, MultiLineString

# 1. 读取输入的GeoJSON多边形文件
gdf = gpd.read_file("your_input_polygons.geojson")

# 2. 定义多边形转直线的处理函数
def polygon_to_line_geometry(geom):
    # 处理单个Polygon
    if geom.geom_type == "Polygon":
        # 把外边界的LinearRing转成LineString
        exterior_line = LineString(geom.exterior.coords)
        # 处理内洞（如果有）
        interior_lines = [LineString(hole.coords) for hole in geom.interiors]
        # 有内洞则返回MultiLineString，否则返回单个LineString
        if interior_lines:
            return MultiLineString([exterior_line] + interior_lines)
        return exterior_line
    # 处理MultiPolygon
    elif geom.geom_type == "MultiPolygon":
        all_lines = []
        for poly in geom.geoms:
            # 提取每个子多边形的外边界
            all_lines.append(LineString(poly.exterior.coords))
            # 提取每个子多边形的内洞
            all_lines.extend([LineString(hole.coords) for hole in poly.interiors])
        return MultiLineString(all_lines)
    # 非多边形类型直接返回原几何
    else:
        return geom

# 3. 批量转换几何列
gdf["geometry"] = gdf["geometry"].apply(polygon_to_line_geometry)

# 可选：验证转换后的几何类型
print("转换后的几何类型：", gdf.geom_type.unique())

# 4. 保存为新的GeoJSON线要素文件
gdf.to_file("output_line_features.geojson", driver="GeoJSON")

关键细节说明 #

• LinearRing vs LineString：多边形的边界默认是LinearRing（闭合环），转成LineString后依然保持闭合（首尾坐标一致）。如果需要开放的线，可以去掉最后一个重复坐标：LineString(geom.exterior.coords[:-1])
• 内洞处理：如果多边形包含内洞（比如环岛的中心区域），函数会自动提取内洞的边界并合并为MultiLineString
• 多多边形兼容：函数支持直接处理MultiPolygon类型，遍历所有子多边形提取边界后合并

内容的提问来源于stack exchange，提问作者Paulo Henrique PH