def METHOD_NAME(fc):
    """Compute the lon/lat bounding box for all transects and regions"""
    bounds = shapely.geometry.GeometryCollection()
    for feature in fc.features:
        shape = shapely.geometry.shape(feature['geometry'])
        shape_bounds = shapely.geometry.box(*shape.bounds)
        bounds = shapely.geometry.box(*bounds.union(shape_bounds).bounds)
    return bounds.bounds