pub fn st_intersection(a: &Geometry, b: &Geometry) -> Geometry
ST_Intersection(a, b) — return the geometry representing the shared area.
Returns GeometryCollection with an empty geometries vec if there is no intersection.