def get_bbox(geo_pd):
    """Gets the bounding box of a geopandas dataframe

    Args:
        geo_pd (geopandas): shape from which we can get the bounding box as a basis for the extend of the meshgrid

    Return:
        (tuple of floats): bbox (x_min, y_min, x_max, y_max)
    """
    return (geo_pd.total_bounds[0], geo_pd.total_bounds[1], geo_pd.total_bounds[2], geo_pd.total_bounds[3])