使用Python减少多边形节点数量并保持拓扑一致性

嗨，这个需求用Python的geopandas和shapely库就能轻松搞定，我给你一步步拆解，代码里加了详细注释，哪怕编程基础弱也能跟着做～

1. 先安装必要的库 #

首先得装处理空间数据的两个核心库，打开终端/命令行跑这个命令：
pip install geopandas shapely

2. 核心思路拆解 #

你说的按角度删节点，本质是遍历多边形的每个中间节点，计算该节点前后两条线段的夹角，保留角度不在175°-180°之间的节点。同时要注意：

• 多边形是闭合的，首尾坐标重复，处理时要跳过重复的首尾
• 多边形可能有内洞（比如湖泊），所以外边界和内边界都要处理
• 处理后要保证拓扑一致性，避免出现自交等无效几何

3. 完整代码实现 #

import geopandas as gpd
from shapely.geometry import Polygon, LinearRing
import math

def calculate_angle(a, b, c):
    """计算三点a-b-c中，b点的夹角（单位：度）"""
    # 计算向量ba和bc
    vec_ba = (a[0] - b[0], a[1] - b[1])
    vec_bc = (c[0] - b[0], c[1] - b[1])
    
    # 用atan2计算夹角（弧度），再转成角度
    angle_rad = math.atan2(vec_bc[1], vec_bc[0]) - math.atan2(vec_ba[1], vec_ba[0])
    angle_deg = math.degrees(angle_rad)
    
    # 把角度调整到0-360范围，再取最小的那个角（比如350°等价于10°）
    angle_deg = abs(angle_deg) % 360
    if angle_deg > 180:
        angle_deg = 360 - angle_deg
    return angle_deg

def simplify_polygon_by_angle(polygon, min_angle=175, max_angle=180):
    """按角度简化多边形，删除夹角在[min_angle, max_angle]之间的节点"""
    # 处理外边界
    exterior_coords = list(polygon.exterior.coords)
    # 闭合多边形首尾坐标重复，先去掉最后一个重复的
    if exterior_coords[0] == exterior_coords[-1]:
        exterior_coords = exterior_coords[:-1]
    
    simplified_exterior = [exterior_coords[0]]  # 保留第一个节点
    for i in range(1, len(exterior_coords)-1):
        a = exterior_coords[i-1]
        b = exterior_coords[i]
        c = exterior_coords[i+1]
        angle = calculate_angle(a, b, c)
        # 如果角度不在要删除的范围，就保留该节点
        if not (min_angle <= angle <= max_angle):
            simplified_exterior.append(b)
    simplified_exterior.append(exterior_coords[-1])  # 保留最后一个节点，保证闭合
    
    # 处理内边界（洞）
    simplified_interiors = []
    for interior in polygon.interiors:
        interior_coords = list(interior.coords)
        if interior_coords[0] == interior_coords[-1]:
            interior_coords = interior_coords[:-1]
        
        simplified_interior = [interior_coords[0]]
        for i in range(1, len(interior_coords)-1):
            a = interior_coords[i-1]
            b = interior_coords[i]
            c = interior_coords[i+1]
            angle = calculate_angle(a, b, c)
            if not (min_angle <= angle <= max_angle):
                simplified_interior.append(b)
        simplified_interior.append(interior_coords[-1])
        simplified_interiors.append(LinearRing(simplified_interior))
    
    # 创建新的多边形
    simplified_poly = Polygon(simplified_exterior, simplified_interiors)
    
    # 修复可能的拓扑问题（比如自交）
    if not simplified_poly.is_valid:
        simplified_poly = simplified_poly.buffer(0)
    
    return simplified_poly

# 主流程：读取Shapefile -> 简化 -> 保存
if __name__ == "__main__":
    # 替换成你的Shapefile路径
    input_shp_path = "你的输入文件.shp"
    output_shp_path = "简化后的文件.shp"
    
    # 读取数据
    gdf = gpd.read_file(input_shp_path)
    
    # 对每个几何进行简化
    gdf["geometry"] = gdf["geometry"].apply(simplify_polygon_by_angle)
    
    # 保存结果
    gdf.to_file(output_shp_path)
    print(f"简化完成，结果已保存到 {output_shp_path}")

4. 关键说明 #

• 角度计算逻辑：函数calculate_angle会把夹角转换为0-180°的范围，这样355°的角会被视为5°，避免误删接近反向的节点
• 拓扑修复：用buffer(0)是shapely里常用的修复无效几何的小技巧，能自动处理自交、重叠等问题
• 备份原始数据：处理前一定要备份你的原始Shapefile，避免操作失误导致数据丢失
• 参数调整：如果觉得简化效果不够或过度，可以修改min_angle和max_angle的值，比如改成170-180°会删更多节点

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