GeoDataFrame.intersection

Returns a GeoSeries of the intersection of points in each aligned geometry with other. The operation works on a 1-to-1 row-wise manner. Note: Unlike most functions, intersection may return the unordered with respect to the index. If this is important to you, you may call sort_index() on the result.

Signature

intersection(other, align=None)

Parameters

Geoseries or geometric object

required

The Geoseries (elementwise) or geometric object to find the intersection with.

bool | None (default None)

default:"None"

If True, automatically aligns GeoSeries based on their indices. None defaults to True. If False, the order of elements is preserved.

Returns

GeoSeries

See also

GeoSeries.difference
GeoSeries.symmetric_difference
GeoSeries.union

Examples

>>> from sedona.spark.geopandas import GeoSeries
>>> from shapely.geometry import Polygon, LineString, Point
>>> s = 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 = 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
  POLYGON ((0 0, 2 2, 0 2, 0 0))
  POLYGON ((0 0, 2 2, 0 2, 0 0))
           LINESTRING (0 0, 2 2)
           LINESTRING (2 0, 0 2)
                     POINT (0 1)
dtype: geometry

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

We can also do intersection of each geometry and a single shapely geometry:

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

We can also check two GeoSeries against each other, row by row. The GeoSeries above have different indices. We can either align both GeoSeries based on index values and compare elements with the same index using align=True or ignore index and compare elements based on their matching order using align=False:

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

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

​Signature

​Parameters

​Returns

​Examples

Signature

Parameters

Returns

Examples