Interpolation with Splines in Tension: A Green's Function Approach
Source: Mathematical Geology, Volume 30, Number 1, January 1998 , pp. 77-93(17)
Interpolation and gridding of data are procedures in the physical sciences and are accomplished typically using an averaging or finite difference scheme on an equidistant grid. Cubic splines are popular because of their smooth appearances; however, these functions can have undesirable oscillations between data points. Adding tension to the spline overcomes this deficiency. Here, we derive a technique for interpolation and gridding in one, two, and three dimensions using Green's functions for splines in tension and examine some of the properties of these functions. For moderate amounts of data, the Green's function technique is superior to conventional finite-difference methods because (1) both data values and directional gradients can be used to constrain the model surface, (2) noise can be suppressed easily by seeking a least-squares fit rather than exact interpolation, and (3) the model can be evaluated at arbitrary locations rather than only on a rectangular grid. We also show that the inclusion of tension greatly improves the stability of the method relative to gridding without tension. Moreover, the one-dimensional situation can be extended easily to handle parametric curve fitting in the plane and in space. Finally, we demonstrate the new method on both synthetic and real data and discuss the merits and drawbacks of the Green's function technique.
Document Type: Research Article
Affiliations: 1: Department of Geology & Geophysics, School of Ocean and Earth Science and Technology, University of Hawaii at Manoa, 1680 East-West Road, Honolulu, Hawaii 96822. firstname.lastname@example.org 2: Department of Geology & Geophysics, School of Ocean and Earth Science and Technology, University of Hawaii at Manoa, 1680 East-West Road, Honolulu, Hawaii 96822
Publication date: January 1, 1998