GeoSeries.
intersection
Returns a GeoSeries of the intersection of points in each aligned geometry with other.
GeoSeries
The operation works on a 1-to-1 row-wise manner:
The Geoseries (elementwise) or geometric object to find the intersection with.
If True, automatically aligns GeoSeries based on their indices. If False, the order of elements is preserved.
See also
GeoSeries.difference
GeoSeries.symmetric_difference
GeoSeries.union
Examples
>>> 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.00000 0.00000, 2.00000 2.00000, 0.... 1 POLYGON ((0.00000 0.00000, 2.00000 2.00000, 0.... 2 LINESTRING (0.00000 0.00000, 2.00000 2.00000) 3 LINESTRING (2.00000 0.00000, 0.00000 2.00000) 4 POINT (0.00000 1.00000) dtype: geometry
>>> s2 1 POLYGON ((0.00000 0.00000, 1.00000 1.00000, 0.... 2 LINESTRING (1.00000 0.00000, 1.00000 3.00000) 3 LINESTRING (2.00000 0.00000, 0.00000 2.00000) 4 POINT (1.00000 1.00000) 5 POINT (0.00000 1.00000) dtype: geometry
We can also do intersection of each geometry and a single shapely geometry:
>>> s.intersection(Polygon([(0, 0), (1, 1), (0, 1)])) 0 POLYGON ((1.00000 1.00000, 0.00000 0.00000, 0.... 1 POLYGON ((1.00000 1.00000, 0.00000 0.00000, 0.... 2 LINESTRING (0.00000 0.00000, 1.00000 1.00000) 3 POINT (1.00000 1.00000) 4 POINT (0.00000 1.00000) 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:
align=True
align=False
>>> s.intersection(s2, align=True) 0 None 1 POLYGON ((1.00000 1.00000, 0.00000 0.00000, 0.... 2 POINT (1.00000 1.00000) 3 LINESTRING (2.00000 0.00000, 0.00000 2.00000) 4 POINT EMPTY 5 None dtype: geometry
>>> s.intersection(s2, align=False) 0 POLYGON ((1.00000 1.00000, 0.00000 0.00000, 0.... 1 LINESTRING (1.00000 1.00000, 1.00000 2.00000) 2 POINT (1.00000 1.00000) 3 POINT (1.00000 1.00000) 4 POINT (0.00000 1.00000) dtype: geometry