def METHOD_NAME(dss, bbox=False, simplify=True, tolerance=0.001):
    from datacube.testutils.geom import epsg4326
    from datacube.utils.geometry import bbox_union
    geoms = [ds.extent.to_crs(epsg4326) for ds in dss]
    if simplify:
        geoms = [g.simplify(tolerance) for g in geoms]
    polygons = [g.__geo_interface__ for g in geoms]
    if bbox:
        return polygons, bbox_union(g.boundingbox for g in geoms)
    return polygons