如何在地质公园中找到与多边形相交的点？
如何在地质公园中找到与多边形相交的点？
提问于 2015-05-22 20:49:06
我一直在尝试使用地理数据图上的“交叉”特性，看看哪个点位于多边形内。但是，只有框架中的第一个功能才会返回为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)回答 5
Stack Overflow用户
回答已采纳
发布于 2015-05-25 14:15:53
根据文档
二进制操作可以在两个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)土工品依靠形状来做几何工作。直接使用它有时很有用(而且更容易阅读)。下列代码也如所宣传的那样工作：
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))您还可以查看形状四舍五入误差的处理中可能出现的边界点问题。
Stack Overflow用户
发布于 2020-11-23 13:45:41
由于地质公园最近经历了许多提高性能的变化，所以这里的答案已经过时了。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')Stack Overflow用户
发布于 2015-05-27 02:54:56
解决这一问题的一种方法似乎是要么获得一个特定的条目(这对我的应用程序不起作用，但可能对其他人的：
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)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/30405652
复制相关文章
相似问题
腾讯云开发者
