blocks|key|48925|text|根据文档|type|unstyled|depth|inlineStyleRanges|entityRanges|offset|length|data|48926|二进制操作可以在两个GeoSeries之间应用，在这种情况下，操作是按元素执行的。这两个系列将通过匹配的指数对齐。|blockquote|48927|你的例子是行不通的。因此，如果您想测试每个点是否位于一个多边形中，则必须这样做：|48928|poly+=+GeoSeries(Polygon([(0,0),+(0,2),+(2,2),+(2,0)]))
g1.intersects(poly.ix[0])+|code-block|syntax|javascript|48929|产出：|48930|++++0++++True
++++1++++True
++++2++++True
++++dtype:+bool|48931|或者，如果您想测试特定GeoSeries中的所有几何图形：|48932|points.intersects(poly.unary_union)|48933|土工品依靠形状来做几何工作。直接使用它有时很有用(而且更容易阅读)。下列代码也如所宣传的那样工作：|48934|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))|48935|您还可以查看形状四舍五入误差的处理中可能出现的边界点问题。|48936|entityMap|0|LINK|mutability|MUTABLE|url|http://geopandas.org/data_structures.html#geoseries|1|https://stackoverflow.com/questions/28028910/how-to-deal-with-rounding-errors-in-shapely^0|2|2|0|0|0|0|0|0|0|0|0|0|0|6|B|1|0^^$0|@$1|2|3|4|5|6|7|1B|8|@]|9|@$A|1C|B|1D|1|1E]]|C|$]]|$1|D|3|E|5|F|7|1F|8|@]|9|@]|C|$]]|$1|G|3|H|5|6|7|1G|8|@]|9|@]|C|$]]|$1|I|3|J|5|K|7|1H|8|@]|9|@]|C|$L|M]]|$1|N|3|O|5|6|7|1I|8|@]|9|@]|C|$]]|$1|P|3|Q|5|K|7|1J|8|@]|9|@]|C|$L|M]]|$1|R|3|S|5|6|7|1K|8|@]|9|@]|C|$]]|$1|T|3|U|5|K|7|1L|8|@]|9|@]|C|$L|M]]|$1|V|3|W|5|6|7|1M|8|@]|9|@]|C|$]]|$1|X|3|Y|5|K|7|1N|8|@]|9|@]|C|$L|M]]|$1|Z|3|10|5|6|7|1O|8|@]|9|@$A|1P|B|1Q|1|1R]]|C|$]]|$1|11|3|-4|5|6|7|1S|8|@]|9|@]|C|$]]]|12|$13|$5|14|15|16|C|$17|18]]|19|$5|14|15|16|C|$17|1A]]]]

According to the <a href="http://geopandas.org/data_structures.html#geoseries" rel="noreferrer">documentation</a>:

<blockquote>
 Binary operations can be applied between two GeoSeries, in which case
 the operation is carried out elementwise. The two series will be
 aligned by matching indices.
</blockquote>

Your examples are not supposed to work. So if you want to test for each point to be in a single polygon you will have to do:

<pre><code>poly = GeoSeries(Polygon([(0,0), (0,2), (2,2), (2,0)]))
g1.intersects(poly.ix[0]) 
</code></pre>

Outputs:

<pre><code> 0 True
 1 True
 2 True
 dtype: bool
</code></pre>

Or if you want to test for all geometries in a specific GeoSeries:

<pre><code>points.intersects(poly.unary_union)
</code></pre>

Geopandas relies on Shapely for the geometrical work. It is sometimes useful (and easier to read) to use it directly. The following code also works as advertised:

<pre><code>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))
</code></pre>

You might also have a look at
<a href="https://stackoverflow.com/questions/28028910/how-to-deal-with-rounding-errors-in-shapely">How to deal with rounding errors in Shapely</a> for issues that may arise with points on boundaries.

blocks|key|885135|text|解决这一问题的一种方法似乎是要么获得一个特定的条目(这对我的应用程序不起作用，但可能对其他人的：|type|unstyled|depth|inlineStyleRanges|entityRanges|data|885136|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])|code-block|syntax|javascript|885137|另一种方法(对我的应用程序更有用)是与第二层功能的一元合并相交：|885138|points.intersects(poly.unary_union)|885139|entityMap^0|0|0|0|0^^$0|@$1|2|3|4|5|6|7|M|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|N|8|@]|9|@]|A|$E|F]]|$1|G|3|H|5|6|7|O|8|@]|9|@]|A|$]]|$1|I|3|J|5|D|7|P|8|@]|9|@]|A|$E|F]]|$1|K|3|-4|5|6|7|Q|8|@]|9|@]|A|$]]]|L|$]]

One way around this seems to be to either get a particular entry (which doesn't work for my application, but may work for someone else's:

<pre><code>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])
</code></pre>

Another way (more useful for my application) is to intersect with a unary union of the features for the second layer:

<pre><code>points.intersects(poly.unary_union)
</code></pre>

blocks|key|1926835|text|使用下面这个简单的函数，您可以很容易地检查多边形内有哪些点：|type|unstyled|depth|inlineStyleRanges|entityRanges|data|1926836|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))|code-block|syntax|javascript|1926837|entityMap^0|0|0^^$0|@$1|2|3|4|5|6|7|I|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|J|8|@]|9|@]|A|$E|F]]|$1|G|3|-4|5|6|7|K|8|@]|9|@]|A|$]]]|H|$]]

You can easily check which points lies inside a polygon with this simple function below:

<pre><code>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))
</code></pre>

blocks|key|885203|text|由于地质公园最近经历了许多提高性能的变化，所以这里的答案已经过时了。Geopandas+0.8引入了许多更改，使得处理大型数据集的速度要快得多。|type|unstyled|depth|inlineStyleRanges|entityRanges|data|885204|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')|code-block|syntax|javascript|885205|entityMap^0|0|0^^$0|@$1|2|3|4|5|6|7|I|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|J|8|@]|9|@]|A|$E|F]]|$1|G|3|-4|5|6|7|K|8|@]|9|@]|A|$]]]|H|$]]

Since geopandas underwent many performance-enhancing changes recently, answers here are outdated. Geopandas 0.8 introduced many changes that makes handling large datasets a lot faster.
<pre><code>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')
</code></pre>

blocks|key|49096|text|我认为最快的方法是使用geopandas.sjoin。|type|unstyled|depth|inlineStyleRanges|offset|length|style|CODE|entityRanges|data|49097|import+geopandas+as+gpd

gpd.sjoin(pts,+poly,+how='left',+op='intersects')|code-block|syntax|javascript|49098|检查示例：链接|49099|entityMap|0|LINK|mutability|MUTABLE|url|https://gis.stackexchange.com/questions/287243/intersect-points-and-polygons-in-shapely#^0|B|F|0|0|5|2|0|0^^$0|@$1|2|3|4|5|6|7|U|8|@$9|V|A|W|B|C]]|D|@]|E|$]]|$1|F|3|G|5|H|7|X|8|@]|D|@]|E|$I|J]]|$1|K|3|L|5|6|7|Y|8|@]|D|@$9|Z|A|10|1|11]]|E|$]]|$1|M|3|-4|5|6|7|12|8|@]|D|@]|E|$]]]|N|$O|$5|P|Q|R|E|$S|T]]]]

I think the fastest way to do it is by using <code>geopandas.sjoin</code>.
<pre><code>import geopandas as gpd

gpd.sjoin(pts, poly, how='left', op='intersects')
</code></pre>
Check example: <a href="https://gis.stackexchange.com/questions/287243/intersect-points-and-polygons-in-shapely#">link</a>

I've been trying to use the "intersects" feature on a geodataframe, looking to see which points lie inside a polygon. However, only the first feature in the frame will return as true. What am I doing wrong?

<pre><code>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)
</code></pre>

How to find which points intersect with a polygon in geopandas?

翻译质量差，导致语言生硬或混乱。

没有提供实际的解决方法或示例。

解答不清晰，无法理解或解决问题。

页面排版不美观，阅读体验差。

文章

问答

视频

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

我一直在尝试使用地理数据图上的“交叉”特性，看看哪个点位于多边形内。但是，只有框架中的第一个功能才会返回为true。我做错了什么？from geopandas.geoseries import *p1 = Point(.5,.5)p2 = Point(.5,1)p3 = Point(1,1)g1 = GeoSeries...

问如何在地质公园中找到与多边形相交的点？
EN

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问如何在地质公园中找到与多边形相交的点？EN