A method for level correction and interpolation of anisotropically sampled aeromagnetic data is described. The level-correction is based on a piecewise comparison of adjacent flight-lines. The interpolation scheme is based on cubic splines and preserves thin linear structures in anisotropic potential field data. The design of such a method is seen as an important building block in the design of an image-based information system for geologic applications