pub fn polygon_intersection(
poly1: &ArrayView2<'_, f64>,
poly2: &ArrayView2<'_, f64>,
) -> SpatialResult<Array2<f64>>pub fn polygon_intersection(
poly1: &ArrayView2<'_, f64>,
poly2: &ArrayView2<'_, f64>,
) -> SpatialResult<Array2<f64>>