Calibrated prediction regions for Gaussian random fields