geopandas.GeoSeries.intersection

GeoSeries.intersection(other, align=None)

返回一个 GeoSeries，它是每个对齐几何中与 other 的点的交集。

该操作以一对一的行方式进行：

Parameters:

otherGeoseries or geometric object

要查找与之相交的 Geoseries（逐元素）或几何对象。

alignbool | None (default None)

如果为真，则根据其索引自动对齐GeoSeries。 如果为假，则保留元素的顺序。 None默认为真。

Returns:

GeoSeries

另请参阅

GeoSeries.difference

GeoSeries.symmetric_difference

GeoSeries.union

示例

>>> from shapely.geometry import Polygon, LineString, Point
>>> s = geopandas.GeoSeries(
...     [
...         Polygon([(0, 0), (2, 2), (0, 2)]),
...         Polygon([(0, 0), (2, 2), (0, 2)]),
...         LineString([(0, 0), (2, 2)]),
...         LineString([(2, 0), (0, 2)]),
...         Point(0, 1),
...     ],
... )
>>> s2 = geopandas.GeoSeries(
...     [
...         Polygon([(0, 0), (1, 1), (0, 1)]),
...         LineString([(1, 0), (1, 3)]),
...         LineString([(2, 0), (0, 2)]),
...         Point(1, 1),
...         Point(0, 1),
...     ],
...     index=range(1, 6),
... )

>>> s
0    POLYGON ((0 0, 2 2, 0 2, 0 0))
1    POLYGON ((0 0, 2 2, 0 2, 0 0))
2             LINESTRING (0 0, 2 2)
3             LINESTRING (2 0, 0 2)
4                       POINT (0 1)
dtype: geometry

>>> s2
1    POLYGON ((0 0, 1 1, 0 1, 0 0))
2             LINESTRING (1 0, 1 3)
3             LINESTRING (2 0, 0 2)
4                       POINT (1 1)
5                       POINT (0 1)
dtype: geometry

我们还可以计算每个几何图形与单个 shapely 几何图形的交集：

>>> s.intersection(Polygon([(0, 0), (1, 1), (0, 1)]))
0    POLYGON ((0 0, 0 1, 1 1, 0 0))
1    POLYGON ((0 0, 0 1, 1 1, 0 0))
2             LINESTRING (0 0, 1 1)
3                       POINT (1 1)
4                       POINT (0 1)
dtype: geometry

我们还可以逐行检查两个GeoSeries。 上面的GeoSeries具有不同的索引。我们可以基于索引值对两个GeoSeries进行对齐，并使用相同索引比较元素，使用 align=True，或者忽略索引并根据它们的匹配顺序比较元素，使用 align=False:

>>> s.intersection(s2, align=True)
0                              None
1    POLYGON ((0 0, 0 1, 1 1, 0 0))
2                       POINT (1 1)
3             LINESTRING (2 0, 0 2)
4                       POINT EMPTY
5                              None
dtype: geometry

>>> s.intersection(s2, align=False)
0    POLYGON ((0 0, 0 1, 1 1, 0 0))
1             LINESTRING (1 1, 1 2)
2                       POINT (1 1)
3                       POINT (1 1)
4                       POINT (0 1)
dtype: geometry