如何找到哪些点与 geopandas 中的多边形相交？

Question

我一直在尝试在地理数据框上使用“相交”功能，以查看哪些点位于多边形内。但是，只有框架中的第一个特征会返回 true。我究竟做错了什么？

 from geopandas.geoseries import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

g1 = GeoSeries([p1,p2,p3])
g2 = GeoSeries([p2,p3])

g = GeoSeries([Polygon([(0,0), (0,2), (2,2), (2,0)])])

g1.intersects(g) # Flags the first point as inside, even though all are.
g2.intersects(g) # The second point gets picked up as inside (but not 3rd)

原文由 Thomas Pingel 发布，翻译遵循 CC BY-SA 4.0 许可协议

Answer 1

根据 文档：

可以在两个 GeoSeries 之间应用二进制操作，在这种情况下，操作是按元素执行的。这两个系列将通过匹配索引对齐。

你的例子不应该工作。所以如果你想测试每个点都在一个多边形中，你将不得不这样做：

 poly = GeoSeries(Polygon([(0,0), (0,2), (2,2), (2,0)]))
g1.intersects(poly.ix[0])

输出：

     0    True
    1    True
    2    True
    dtype: bool

或者，如果您想测试特定 GeoSeries 中的所有几何图形：

 points.intersects(poly.unary_union)

Geopandas 依赖 Shapely 进行几何工作。直接使用它有时很有用（也更容易阅读）。以下代码也像宣传的那样工作：

 from shapely.geometry import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

poly = Polygon([(0,0), (0,2), (2,2), (2,0)])

for p in [p1, p2, p3]:
    print(poly.intersects(p))

您还可以查看 如何处理 Shapely 中的舍入错误， 以了解边界上的点可能出现的问题。

原文由 Fabzi 发布，翻译遵循 CC BY-SA 4.0 许可协议

Answer 2

由于 geopandas 最近经历了许多性能增强的变化，这里的答案已经过时了。 Geopandas 0.8 引入了许多变化，可以更快地处理大型数据集。

 import geopandas
from shapely.geometry import Polygon

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

points = GeoSeries([p1,p2,p3])

poly = GeoSeries([Polygon([(0,0), (0,2), (2,2), (2,0)])])

geopandas.overlay(points, poly, how='intersection')

原文由 Magnus 发布，翻译遵循 CC BY-SA 4.0 许可协议