摘要

Due to practical limitations, geophysical acquisitions for natural resource exploration are always performed at scattered points in space. To derive more rational and accurate interpretations about physical properties of earth forming materials, it is necessary to interpolate geophysical data in a dense and regular grid as a preliminary step. This article introduces a method for interpolation of scattered data using multilevel B-splines to estimate field values on a dense and regular grid based on original sparse and irregular data. The algorithm applies an effective B-splines approximation technique to a hierarchy of control lattices to generate a sequence of functions whose sum approaches the desired interpolation function. This method is implemented successfully in interpolation on both synthetic and real magnetic anomaly data sets. In comparison, it achieves better results than inverse distance weighing, lcriging and minimum curvature methods.