Я пытаюсь найти объединение двух полигонов в GeoPandas и вывести единую геометрию, которая включает точки из обоих полигонов в качестве вершин. Функция 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
) с помощью GeoPandas или shapely?
Во-вторых, какова правильная номенклатура, связанная с геометрией (poly_union
), которую я пытаюсь найти? Я бы назвал это «объединением», но каждый нахожий мне пример, относящийся к «объединениям», не выводит эту геометрию.
Примечание. Этот пример также не выводит ни одного многоугольника:
poly1 = df1['geometry']; poly2 = df2['geometry']
mergedpoly = poly1.union(poly2)
mergedpoly.plot()
overlay
не изменились (как я наблюдал с v 0.6.1). - person swatchai   schedule 21.11.2020