def METHOD_NAME(compound_region):
    """Compute center for a CompoundRegion
    The center of the compound region is defined here as the geometric median
    of the individual centers of the regions. The geometric median is defined
    as the point the minimises the distance to all other points.
    Parameters
    ----------
    compound_region : `CompoundRegion`
        Compound region
    Returns
    -------
    center : `~astropy.coordinates.SkyCoord`
        Geometric median of the positions of the individual regions
    """
    regions = compound_region_to_regions(compound_region)
    if len(regions) == 1:
        return regions[0].center
    positions = SkyCoord([region.center.icrs for region in regions])
    def f(x, coords):
        """Function to minimize"""
        lon, lat = x
        center = SkyCoord(lon * u.deg, lat * u.deg)
        return np.sum(center.separation(coords).deg)
    ra, dec = positions.ra.wrap_at("180d").deg, positions.dec.deg
    ub = np.array([np.max(ra), np.max(dec)])
    lb = np.array([np.min(ra), np.min(dec)])
    if np.all(ub == lb):
        bounds = None
    else:
        bounds = Bounds(ub=ub, lb=lb)
    result = minimize(
        f,
        x0=[np.mean(ra), np.mean(dec)],
        args=(positions,),
        bounds=bounds,
        method="L-BFGS-B",
    )
    return SkyCoord(result.x[0], result.x[1], frame="icrs", unit="deg")