Trouvez les points d'intersection de 2 polygones en Python

Question

J'essaie d'obtenir les pointcoordonnées croisées de 2 polygonsen Python.

Par exemple, les points d'intersection des deux polygones ( polys2). À des fins d'illustration, j'ai marqué les points d'intersection en vert.

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);

Quelle est la meilleure façon de trouver ces points d'intersection en Python?

Answer 1

Utilisez l'intersection de la limite des polygones. Dans votre cas, il renvoie les points d'intersection sous forme de multipoint.

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 ...>]