Large-scale magnetic field maps using structured kernel interpolation for Gaussian process regression