Question

我一直试图使用＆＃34;交叉＆＃34;地理数据框上的要素，用于查看多边形内的哪些点。但是，只有帧中的第一个特征才会返回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)

Answer 1

根据documentation：

二进制操作可以在两个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))

您可能也会看一下 How to deal with rounding errors in Shapely针对边界点上可能出现的问题。

Answer 2

解决这个问题的方法似乎是要么得到一个特定的条目（这对我的应用程序不起作用，但可能适用于其他人的：

from geopandas.geoseries import *

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)])])

points.intersects(poly.ix[0])

另一种方式（对我的应用程序更有用）是与第二层的一元特征的一元并集相交：

points.intersects(poly.unary_union)

Answer 3

由于大熊猫最近经历了许多提高性能的变化，因此这里的答案已经过时了。 Geopandas 0.8引入了许多更改，这些更改使处理大型数据集的速度大大加快。

import geopandas
from shapely 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')

Answer 4

您可以使用以下简单功能轻松检查哪些点位于多边形内：

import geopandas
from shapely.geometry import *

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

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

def point_inside_shape(point, shape):
    #point of type Point
    #shape of type Polygon
    pnt = geopandas.GeoDataFrame(geometry=[point], index=['A'])
    return(pnt.within(shape).iloc[0])

for p in [p1, p2, p3]:
    print(point_inside_shape(p, g))

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

4 个答案: