Fast covariance parameter estimation of spatial Gaussian process models using neural networks