3D kriging of the ionosphere based on maximum likelihood and restricted maximum likelihood estimation of a non-stationary covariance model