def intersect_bounds(
    bbox0: tuple[float, float, float, float],
    bbox1: tuple[float, float, float, float]
) -> tuple[float, float, float, float]:
    """Get the intersection in w s e n

    Parameters:
        bbox0: bounding coordinates
        bbox1: bounding coordinates

    Returns:
        intersection coordinate bounds (w, s, e, n)

    """
    w = max(bbox0[0], bbox1[0])
    s = max(bbox0[1], bbox1[1])
    e = min(bbox0[2], bbox1[2])
    n = min(bbox0[3], bbox1[3])
    bounds = (w, s, e, n)

    return bounds