shapely_geometry_intersection_area_and_spatial_comparison.py

1from shapely.geometry import Point
2
3# Create a patch of a circle (a Point buffered by a distance)
4patch = Point(0.0, 0.0).buffer(10.0)
5
6# Create a square Polygon
7from shapely.geometry import Polygon
8square = Polygon([(0, 0), (1, 0), (1, 1), (0, 1)])
9
10# Calculate the intersection of the two objects
11intersection = patch.intersection(square)
12
13# Check the area of the intersection
14print(f"Area of intersection: {intersection.area}")
15
16# Check if the intersection is equal to the square (it should be in this case)
17print(f"Is intersection equal to square? {intersection.equals(square)}")