Abstract
The accuracy of many schemes for interpolating scattered data with radial basis functions depends on a shape parameter c of the radial basis function. In this paper we study the effect of c on the quality of fit of the multiquadric, inverse multiquadric and Gaussian interpolants. We show, numerically, that the value of the optimal c (the value of c that minimizes the interpolation error) depends on the number and distribution of data points, on the data vector, and on the precision of the computation. We present an algorithm for selecting a good value for c that implicitly takes all the above considerations into account. The algorithm selects c by minimizing a cost function that imitates the error between the radial interpolant and the (unknown) function from which the data vector was sampled. The cost function is defined by taking some norm of the error vector E = (E 1, ... , EN)T where E k = Ek = fk - Sk xk) and S k is the interpolant to a reduced data set obtained by removing the point x k and the corresponding data value f k from the original data set. The cost function can be defined for any radial basis function and any dimension. We present the results of many numerical experiments involving interpolation of two dimensional data sets by the multiquadric, inverse multiquadric and Gaussian interpolants and we show that our algorithm consistently produces good values for the parameter c.
Similar content being viewed by others
References
R.E. Carlson and T.A. Foley, The parameter R 2 in multiquadric interpolation, Comput. Math. Appl. 21 (1991) 29-42.
T.A. Foley, Interpolation and approximation of 3-D and 4-D scattered data, Comput. Math. Appl. 13 (1987) 711-740.
T.A. Foley, Near optimal parameter selection for multiquadric interpolation, J. Appl. Sci. Comput. 1 (1991) 54-69.
R. Franke, Scattered data interpolation: tests of some methods, Math. Comp. 38 (1982) 181-200.
R.L. Hardy, Multiquadric equations of topography and other irregular surfaces, J. Geophys. Res. 76 (1971) 1905-1915.
R.L. Hardy, Theory and applications of the multiquadric-biharmonic method, Comput. Math. Appl. 19 (1990) 163-208.
E.J. Kansa and R.E. Carlson, Improved accuracy of multiquadric interpolation using variable shape parameters, Preprint, Lawrence Livermore National Laboratory.
T. Lyche and K. Morken, Knot removal for parametric B-spline curves and surfaces, Comput. Aided Geom. Design 4 (1987) 217-230.
C.A. Micchelli, Interpolation of scattered data: distance matrices and conditionally positive definite functions, Constr. Approx. 2 (1986) 11-22.
T. Poggio and F. Girosi, Networks for approximation and learning, Proceedings of the IEEE 78 (1990) 1481-1497.
M.J.D. Powell, The theory of radial basis function approximation in 1990, in: Advances in Numerical Analysis, Vol. II: Wavelets, Subdivision Algorithms and Radial Functions, ed. W. Light (Oxford University Press, Oxford, UK, 1991) pp. 105-210.
W.H. Press, B.P. Flannery, S.A. Teukolsky and W.T. Vetterling, Numerical Recipes in C. The Art of Scientific Computing, 2nd ed. (Cambridge University Press, Cambridge, 1992).
S.I.M. Ritchie, Surface representation by finite elements, Ms. thesis, University of Calgary (1978).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Rippa, S. An algorithm for selecting a good value for the parameter c in radial basis function interpolation. Advances in Computational Mathematics 11, 193–210 (1999). https://doi.org/10.1023/A:1018975909870
Issue Date:
DOI: https://doi.org/10.1023/A:1018975909870