獲取與區域接壤的區域

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