本站消息

  本站每日ip已达10000,出租广告位,位置价格可谈,需要合作请联系站长


+关注
已关注

分类  

暂无分类

标签  

暂无标签

日期归档  

暂无数据

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

发布于2019-09-11 19:42     阅读(713)     评论(0)     点赞(11)     收藏(5)


我一直在尝试在地理数据框架上使用“相交”功能,以查看哪些点位于多边形内。但是,只有帧中的第一个特征才会返回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)

解决方案


根据文件

可以在两个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中的舍入错误,以了解边界上的点可能出现的问题。



所属网站分类: 技术文章 > 问答

作者:黑洞官方问答小能手

链接: https://www.pythonheidong.com/blog/article/108051/

来源: python黑洞网

任何形式的转载都请注明出处,如有侵权 一经发现 必将追究其法律责任

11 0
收藏该文
已收藏

评论内容:(最多支持255个字符)