Estimating the Rate-Distortion Function by Wasserstein Gradient Descent