Get Bezirke am Rande des gegebenen Bezirks

Question

Ich möchte die Geometrien der Bezirke an einem bestimmten Bezirk.

districts 
d0 = districts[0] 
gpd.sjoin(d0, districts, op='intersects') 

Dies gibt die Geometrie von d0 in jeder Zeile. Aber ich möchte die Geometrie der richtigen Tabelle in jeder Zeile. Ist es möglich, sowohl linke als auch rechte Tabellengeometrien zu erhalten?

Answer 1

Sie könnten join verwenden, um die Geometrie von der rechten Tabelle nach dem sjoin

gdf = gpd.sjoin(d0, districts, op='intersects') 

GDF zu bekommen eine Spalte/Reihe namens index_right die wir

gdf.join(districts['geometry'], on='index_right', lsuffix='', rsuffix='_districts') 

Nicht sicher, wie geopandas nutzen können behandelt zwei Geometrien. Ich nehme an, dass alle Operationen das Original von d0

nutzen werden