Растворение полигонов в Python приводит к отсутствию некоторых полигонов. ⇐ Python

[Anonymous]

У меня есть большой код, который по сути векторизует слой Tiff в шейп-файл многоугольника. последний шаг кода — растворение полигонов в зависимости от их значения (в диапазоне от 1 до 5). Я пробовал использовать функцию «dissolve» Geopandas, Shapely и «unary_union» Фионы, я пробовал функцию ogr2ogr GDAl, и во всех из них я получаю одну и ту же проблему.
Некоторые полигоны нижнего уровня подвергаются ' Проглочено высшими уровнями. почему-то на некоторых участках полигоны меняют свои границы и выходят за нижние уровни. это означает, что в конечном результате мне не хватает некоторых полигонов.
Тот же эффект происходит, когда я использую инструмент растворения в 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}")


Подробнее здесь: https://stackoverflow.com/questions/786 ... ns-missing