Application of universal Kriging interpolation in geomagnetic map