使用 Geopandas，如何 select 所有不在多边形内的点？

Question

我有一个包含芝加哥地址的 DataFrame，我已将其地理编码为纬度和经度值，然后编码为 Point 对象（使 DataFrame 成为 GeoDataFrame）。一小部分已使用芝加哥以外的 LatLong 值进行了错误的地理编码。我有一个芝加哥边界的 shapefile (GeoDataFrame)，我想 select 点在芝加哥边界多边形之外的所有行。

select 多边形内的所有点（通过 geopandas sjoin 函数）很容易，但我还没有找到 select 不在多边形内的点的好方法。有吗？

Answer 1

如果将芝加哥边界 GeoDataFrame 转换为单个多边形，例如：

chicago = df_chicago.geometry.unary_union

然后您可以使用带有 within 运算符的布尔过滤来 select 芝加哥内外的点：

within_chicago = df[df.geometry.within(chicago)]
outside_chicago = df[~df.geometry.within(chicago)]

使用 ~ 反转布尔条件。

或者，您可以使用 disjoint 空间谓词：

outside_chicago = df[df.geometry.disjoint(chicago)]

使用 Geopandas，如何 select 所有不在多边形内的点？

Using Geopandas, how do I select all points not within a polygon?

python

pandas

shapely

geopandas