docs: 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)))