def compute_centers(clusters):
    """
    Compute centers of each cluster

    Parameters
    ----------
    clusters : pandas.DataFrame
        RPM Clusters DataFrame

    Returns
    -------
    centers : ndarray
        n x 2 array of (lat, lon) centers
    """
    centers = clusters.groupby('cluster_id')[['latitude', 'longitude']].mean()

    return centers.values