Suchen Sie die Schnittpunkte von 2 Polygonen in Python

Question

Ich versuche, die Schnittkoordinaten pointvon 2 polygonsin Python zu erhalten.

Zum Beispiel Schnittpunkte der beiden Polygone ( polys2). Zur Veranschaulichung habe ich die Schnittpunkte grün markiert.

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

Was ist der beste Weg, um diese Schnittpunkte in Python zu finden?

Answer 1

Verwenden Sie den Schnittpunkt der Polygonengrenze. In Ihrem Fall werden Schnittpunkte als Mehrpunkt zurückgegeben.

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