发布于2019-09-11 19:42 阅读(809) 评论(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/0fef75f03a5265e1fd7c/
来源:python黑洞网
任何形式的转载都请注明出处,如有侵权 一经发现 必将追究其法律责任
昵称:
评论内容:(最多支持255个字符)
---无人问津也好,技不如人也罢,你都要试着安静下来,去做自己该做的事,而不是让内心的烦躁、焦虑,坏掉你本来就不多的热情和定力
Copyright © 2018-2021 python黑洞网 All Rights Reserved 版权所有,并保留所有权利。 京ICP备18063182号-1
投诉与举报,广告合作请联系vgs_info@163.com或QQ3083709327
免责声明:网站文章均由用户上传,仅供读者学习交流使用,禁止用做商业用途。若文章涉及色情,反动,侵权等违法信息,请向我们举报,一经核实我们会立即删除!