Python多边形溶解后部分多边形丢失，ArcMap同现问题QGIS正常求解决

我将TIFF图层矢量化为多边形Shapefile后，需要基于SM_LV字段（取值1-5）执行溶解操作，但尝试Geopandas的dissolve方法、Shapely+Fiona的unary_union、GDAL的ogr2ogr工具均出现同一问题：部分低层级多边形被高层级多边形“吞没”，边界偏移覆盖导致输出丢失部分要素。该问题在ArcMap溶解工具中同样存在，但QGIS的溶解工具可正常输出。

当前出现问题的代码：

def dissolve_shapefiles(self, algo_class):
    """Dissolves merged shapefiles by SM_LV column."""
    try:
        dissolved = self.final_gdf.dissolve(by='SM_LV', as_index=False)
        algo_class.dissolved_output_path = os.path.abspath(os.path.join(self.output_directory, "..",
                                                                        self.project_name + '_Vectorized_Dissolved.shp'))
        dissolved.to_file(algo_class.dissolved_output_path)
        full_path = os.path.abspath(algo_class.dissolved_output_path)
        print(f"Dissolved shapefile saved as {algo_class.dissolved_output_path}")
        print(f"Full path: {full_path}")
        self.split_merged_shapefile(algo_class.dissolved_output_path)
    except Exception as e:
        print(f"Error dissolving shapefile: {e}") 

1. 按层级从低到高依次溶解并裁剪 #

核心逻辑是先保留低层级要素，再用低层级要素裁剪高层级要素，避免高层级直接覆盖低层级。修改后的代码如下：

def dissolve_shapefiles(self, algo_class):
    """按SM_LV层级从低到高溶解，避免高层级吞没低层级"""
    try:
        # 按层级升序排序，确保从最低层级开始处理
        sorted_gdf = self.final_gdf.sort_values(by='SM_LV')
        result_gdf = None
        
        # 遍历每个层级（1到5）
        for lv in sorted_gdf['SM_LV'].unique().sort():
            # 提取当前层级的所有多边形并溶解
            lv_gdf = sorted_gdf[sorted_gdf['SM_LV'] == lv]
            dissolved_lv = lv_gdf.dissolve(by='SM_LV', as_index=False)
            
            if result_gdf is None:
                # 第一个层级直接作为基础结果
                result_gdf = dissolved_lv
            else:
                # 用已有的低层级要素裁剪当前层级，去除重叠部分
                dissolved_lv['geometry'] = dissolved_lv['geometry'].difference(result_gdf.unary_union)
                # 合并到结果中
                result_gdf = pd.concat([result_gdf, dissolved_lv], ignore_index=True)
        
        # 保存最终结果
        algo_class.dissolved_output_path = os.path.abspath(os.path.join(self.output_directory, "..",
                                                                        self.project_name + '_Vectorized_Dissolved.shp'))
        result_gdf.to_file(algo_class.dissolved_output_path)
        full_path = os.path.abspath(algo_class.dissolved_output_path)
        print(f"Dissolved shapefile saved as {algo_class.dissolved_output_path}")
        print(f"Full path: {full_path}")
        self.split_merged_shapefile(algo_class.dissolved_output_path)
    except Exception as e:
        print(f"Error dissolving shapefile: {e}")

2. 先修复多边形拓扑错误再溶解 #

矢量化后的多边形可能存在微小重叠、缝隙或拓扑无效问题，这会导致溶解时出现异常。可以在溶解前添加拓扑修复步骤：

# 在dissolve前添加：修复多边形拓扑错误
self.final_gdf['geometry'] = self.final_gdf['geometry'].buffer(0)

buffer(0)是Geopandas中常用的拓扑修复方法，能自动修正多边形的自相交、无效环等问题。

3. 模拟QGIS溶解的拓扑处理逻辑 #

QGIS的溶解工具默认会处理要素的拓扑关系，避免重叠覆盖。可以查看QGIS溶解工具的参数（如是否勾选“保持边界”“处理重叠”），在代码中模拟：

• 对同层级多边形先执行unary_union生成单部件要素，再explode拆分回多部件
• 结合overlay工具处理不同层级的重叠关系，保留低层级要素的优先级

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