def contrast_rms(data, *kwargs):
    """ Compute RMS contrast norm of an image
    """
    av = np.average(data, *kwargs)
    mal = 1 / (data.shape[0] * data.shape[1])
    return np.sqrt(mal * np.sum(np.square(data - av)))