Reproducing kernel Hilbert space method for optimal interpolation of potential field data