Encontre os pontos de intersecção de 2 polígonos em Python

Question

Estou tentando obter as pointcoordenadas de interseção de 2 polygonsem Python.

Por exemplo, pontos de intersecção dos dois polígonos ( polys2). Para fins de ilustração, marquei os pontos de intersecção em verde.

from shapely.geometry import Polygon,Point
import geopandas
polys1 = geopandas.GeoSeries([Polygon([(0,0), (2,0), (2,2), (0,2)]),
                               Polygon([(1,1), (3,1), (3,3), (1,3)])])
 
polys2 = geopandas.GeoSeries([Point([(1,2)]),Point([(2,1)]) ])
                                  
df1 = geopandas.GeoDataFrame({'geometry': polys1, 'df1':[1,2]})
df2 = geopandas.GeoDataFrame({'geometry': polys2, 'df2':[1,2]})

ax = df1.plot(color='red', alpha = 0.5)
df2.plot(ax=ax, color='green', alpha=0.5);

Qual é a melhor maneira de encontrar esses pontos de intersecção em Python?

Answer 1

Use a interseção do limite dos polígonos. No seu caso, ele retorna pontos de intersecção como multiponto.

geom1 = df1.geometry[0]
geom2 = df1.geometry[1]

points = geom1.boundary.intersection(geom2.boundary) # multipoint
points = list(points) # point list

print(points)
#[<shapely.geometry.point.Point at ...>,
# <shapely.geometry.point.Point at ...>]