地理：使用环将每个多边形与GeodataFrame之前的每个多边形区分开来？

Question

我有一个带有许多多边形的GeodataFrame。我想通过循环自动化差异，以每个先前的每个和polygon进行此操作，然后将它们添加到新的gdf：

new_gdf = polygon [0：1] .overlay（polygon [1：2]，如何='差异'）

我一次可以做到这一点，但是我想找到一种更有效的方法，尤其是当我（某天）有很多多边形时。

有人会想知道我怎么做吗？我正在尝试使用iterrows（），但是在确切的语法上失败了。

先感谢您。

Answer 1

而不是覆盖（）您可能正在寻找 geoseries.differenies.difference（）

返回每个对齐几何形状中的点的地层
不在另一个。

可能是这样的：

from shapely.geometry import Polygon
import geopandas as gpd

# sample polygons
polys = gpd.GeoDataFrame(geometry=[
    Polygon([(0,0), (2,0), (2,2), (0,2)]),
    Polygon([(1,1), (3,1), (3,3), (1,3)]),
    Polygon([(2,2), (4,2), (4,4), (2,4)]),
    Polygon([(3,3), (5,3), (5,5), (3,5)])])

# polys.iloc[:-1] : rows 0 .. 2 (exclude last row as df sizes must match)
# polys.iloc[1:]  : rows 1 .. 3
# align=False : do not aligns GeoSeries based on their indices; 
# order of elements is preserved.
diffs = polys.iloc[:-1].difference(polys.iloc[1:], align=False)

示例polys geodataframe：

结果diffs geodataseries：