我有一个包含芝加哥地址的 DataFrame，我已将其地理编码为纬度和经度值，然后编码为 Point 对象（使 DataFrame 成为 GeoDataFrame）。芝加哥以外的一小部分使用 LatLong 值进行了错误的地理编码。我有一个芝加哥边界的形状文件（GeoDataFrame），我想选择点位于芝加哥边界多边形之外的所有行。
选择多边形内的所有点（通过 geopandas sjoin 函数）很容易，但我还没有找到选择不在多边形内的点的好方法。有吗？
如果将芝加哥边界 GeoDataFrame 转换为单个多边形，例如：
chicago = df_chicago.geometry.unary_union
那么你可以使用布尔过滤within运营商选择芝加哥境内和境外的点：
within_chicago = df[df.geometry.within(chicago)]
outside_chicago = df[~df.geometry.within(chicago)]
using ~反转布尔条件。
或者，您可以使用disjoint空间谓词：
outside_chicago = df[df.geometry.disjoint(chicago)]
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
