Объединение полигонов в Python, GeoPandas или shapely (в единую геометрию)

Я пытаюсь найти объединение двух полигонов в Геопандах и вывести единую геометрию, которая включает точки из обоих полигонов в качестве своих вершин. Функция geopandas.overlay дает мне полигоны для каждого отдельного объединения, но я хотел бы иметь один полигон.

Для контекста я использую это, чтобы объединить две административные области вместе в одну область (т. е. включить городской округ в пределах страны).

Следующий пример взят с веб-сайта geopandas и иллюстрирует, что Я бы хотел:

from matplotlib import pyplot as plt

import geopandas as gpd

from shapely.geometry import Polygon



polys1 = gpd.GeoSeries([Polygon([(0,0), (2,0), (2,2), (0,2)]),

                         Polygon([(2,2), (4,2), (4,4), (2,4)])])



polys2 = gpd.GeoSeries([Polygon([(1,1), (3,1), (3,3), (1,3)]),

                         Polygon([(3,3), (5,3), (5,5), (3,5)])])



df1 = gpd.GeoDataFrame({'geometry': polys1, 'df1':[1,2]})

df2 = gpd.GeoDataFrame({'geometry': polys2, 'df2':[1,2]})



res_union = gpd.overlay(df1, df2, how='union')

res_union.plot()

Ни одна из выходных геометрий не является тем, что я ожидал, а именно:

poly_union = gpd.GeoSeries([Polygon([(0,0), (0,2), (1,2), (1,3), 

    (2,3), (2,4), (3, 4), (3, 5), (5, 5), (5, 3), (4, 3), (4, 2), 

    (3,2), (3,1), (2, 1), (2, 0), (0, 0)])])



poly_union.plot(color = 'red')

plt.show()

Во-первых, как вывести вышеупомянутый полигон (poly_union) из входных полигонов (df1, df2) с помощью Геопанд или фигуристой?

Во-вторых, какова правильная номенклатура, связанная с геометрией (poly_union), которую я пытаюсь найти? Я бы назвал это "Союзом", но каждый пример, который я нахожу, относится к "союзам". не выводит эту геометрию.

Примечание: этот пример , похоже, также не выводит один полигон:

poly1 = df1['geometry']; poly2 = df2['geometry']

mergedpoly = poly1.union(poly2)

mergedpoly.plot()