获取与区域接壤的区域

Question

我想获得与特定区域接壤的区域的几何图形。

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

这给出了每行中d0的几何形状。但是我想要每行右表的几何形状。是否有可能获得左右表格几何？

Answer 1

你可以使用join您sjoin

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

GDF后得到右表中的几何图形将有一列/一系列名为index_right我们可以利用

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

不知道如何geopandas将处理两个几何。我猜所有的操作都会利用原来的操作d0