This Article 
 Bibliographic References 
 Add to: 
Scattered Data Interpolation with Multilevel B-Splines
July-September 1997 (vol. 3 no. 3)
pp. 228-244

Abstract—This paper describes a fast algorithm for scattered data interpolation and approximation. Multilevel B-splines are introduced to compute a $C^2$-continuous surface through a set of irregularly spaced points. The algorithm makes use of a coarse-to-fine hierarchy of control lattices to generate a sequence of bicubic B-spline functions whose sum approaches the desired interpolation function. Large performance gains are realized by using B-spline refinement to reduce the sum of these functions into one equivalent B-spline function. Experimental results demonstrate that high-fidelity reconstruction is possible from a selected set of sparse and irregular samples.

[1] E. Arge, M. Dæhlen, and A. Tveito, "Approximation of Scattered Data Using Smooth Grid Functions," J. Computational and Applied Math., vol. 59, pp. 191-205, 1995.
[2] R. Barnhill, "Representation and Approximation of Surfaces," J.R. Rice, ed., Mathematical Software III, pp. 68-119.New York: Academic Press, 1977.
[3] W. Bohm and H. Prautzsch, "The Insertion Algorithm," Computer Aided Design, vol. 17, no. 2, pp. 58-59, 1985.
[4] A. Brandt, "Multi-Level Adaptive Solutions to Boundary-Value Problems," Mathematics of Computation, vol. 31, no. 138, pp. 333-390, 1977.
[5] W.L. Briggs, A Multigrid Tutorial.Lancaster, Penn.: SIAM, Lancaster Press, 1987.
[6] P.J. Burt, "Moment Images, Polynomial Fit Filters, and the Problem of Surface Interpolation," Proc. IEEE Conf. Computer Vision and Pattern Recognition, pp. 144-152, 1988.
[7] R. Clough and J. Tocher, "Finite Element Stiffness Matrices for Analysis of Plates in Bending," Proc. Conf. Matrix Methods in Structural Mechanics, pp. 515-545, 1965.
[8] E. Cohen, T. Lyche, and R. Riesenfeld, "Discrete B-Splines and Subdivision Techniques in Computer Aided Geometric Design and Computer Graphics," Computer Graphics and Image Processing, vol. 14, no. 2, pp. 87-111, 1980.
[9] I. Crain and K. Bhattacharyya, "Treatment of Nonequispaced Two-Dimensional Data with a Digital Computer," Geoexploration, vol. 5, pp. 173-194, 1967.
[10] G. Cressman, "An Operational Objective Analysis System," Mon. Wea. Rev., vol. 87, pp. 367-374, 1959.
[11] J. Duchon, "Splines Minimizing Rotation-Invariant Semi-Norms in Sobolev Spaces," C. Chui, L. Schumaker, and J. Ward, eds., Multivariate Approximation Theory, pp. 85-100.Basel, Switzerland: Birkhauser, 1975.
[12] N. Dyn, "Interpolation and Approximation by Radial and Related Function," C. Chui, L. Schumaker, and J. Ward, eds., Approximation Theory VI, pp. 211-234.San Diego, Calif.: Academic Press, 1989.
[13] N. Dyn, D. Levin, and S. Ruppa, "Numerical Procedures for Global Surface Fitting of Scattered Data by Radial Functions," SIAM J. Sci. Statis. Comput., vol. 7, pp. 639-659, 1986.
[14] J.D. Foley,A. van Dam,S.K. Feiner,, and J.F. Hughes,Computer Graphics: Principles and Practice,Menlo Park, Calif.: Addison-Wesley, 1990.
[15] T.A. Foley and G.M. Nielson, "Multivariate Interpolation to Scattered Data Using Delta Iteration," Approximation Theory III, E. Cheney, ed., pp. 419-424.New York: Academic Press, 1980.
[16] D.R. Forsey and R.H. Bartels, “Hierarchical B-Spline Refinement,” Computer Graphics (SIGGRAPH '88 Proc.), vol. 22, no. 4, pp. 205-212, Aug. 1988.
[17] D.R. Forsey and R.H. Bartels, "Tensor Products and Hierarchical Fitting," Curves and Surfaces in Computer Vision and Graphics II (Proc. SPIE 1610), pp. 88-96, 1991.
[18] D.R. Forsey and R.H. Bartels, "Surface Fitting with Hierarchical Splines," ACM Trans. Graphics, vol. 14, no. 2, pp. 134-161, 1995.
[19] R. Franke, "A Critical Comparison of Some Methods for Interpolation of Scattered Data," Technical Report NPS-53-79-003, Naval Postgraduate School, 1979.
[20] R. Franke, "Scattered Data Interpolation: Tests of Some Methods," Math. Comp., vol. 38, pp. 181-200, 1982.
[21] R. Franke and G.M. Nielson, "Scattered Data Interpolation of Large Sets of Scattered Data," Intl. J. Numerical Methods in Eng., vol. 15, pp. 1,691-1,704, 1980.
[22] R. Franke and G.M. Nielson, "Scattered Data Interpolation and Applications: A Tutorial and Survey," Geometric Modelling: Methods and Their Application, H. Hagen and D. Roller, eds., pp. 131-160.Berlin: Springer-Verlag, 1991.
[23] A. Glassner, Principles of Digital Image Synthesis.San Francisco: Morgan Kaufman, 1995.
[24] R.C. Gonzalez and R.E. Woods, Digital Image Processing, Addison-Wesley, New York, 1993.
[25] S.J. Gortler and M.F. Cohen, "Hierarchical and Variational Geometric Modeling with Wavelets," Symp. Interactive 3D Graphics, pp. 35-42, 1995.
[26] A. Goshtasby, “Piecewise Cubic Mapping Functions for Image Registration,” Pattern Recognition, Vol. 20, No. 5, 1987, pp. 525-533.
[27] R. Hardy, "Multiquadratic Equations of Topography and Other Irregular Surfaces," J. Geophysical Research, vol. 76, no. 8, pp. 1,905-1,915, 1971.
[28] J. Hoschek and D. Lasser, Computer Aided Geometric Design.Wellesley, Mass.: A.K. Peters, Ltd., 1993.
[29] W.M. Hsu,J.F. Hughes, and H. Kaufman,"Direct Manipulation of Free-Form Deformations," Computer Graphics, vol. 26, no. 2, pp. 177-184, 1992 (Siggraph 92).
[30] C. Lawson, "Software for$C^1$Surface Interpolation," Mathematical Software III, J.R. Rice, ed., pp. 161-194.New York: Academic Press, 1977.
[31] S. Lee, K.-Y. Chwa, J. Hahn, and S.Y. Shin, "Image Morphing Using Deformable Surfaces," Proc. Computer Animation '94, pp. 31-39, 1994.
[32] S. Lee, K.-Y. Chwa, J. Hahn, and S.Y. Shin, "Image Morphing Using Deformation Techniques. J. Visualization and Computer Animation, vol. 7, no. 1, pp. 3-23, 1996.
[33] S. Lee, K.-Y. Chwa, S.Y. Shin, and G. Wolberg, "Image Metamorphosis Using Snakes and Free-Form Deformations," Computer Graphics (Proc. SIGGRAPH '95), pp. 439-448, 1995.
[34] S. Lee, G. Wolberg, K.-Y. Chwa, and S.Y. Shin, "Image Metamorphosis with Scattered Feature Constraints," IEEE Trans. Visualization and Computer Graphics, vol. 2, no. 4, pp. 337-354, 1996.
[35] P. Litwinowicz and L. Williams, "Animating Images with Drawings," Computer Graphics (Proc. SIGGRAPH '94), pp. 409-412, 1994.
[36] T. Lyche, "Note on the Oslo Algorithm," Computer Aided Design, vol. 20, no. 6, pp. 353-355, 1988.
[37] T. Lyche, E. Cohen, and K. Morken, "Knot Line Refinement Algorithms for Tensor Product B-Spline Surfaces," Computer Aided Geometric Design, vol. 2, pp. 133-139, 1985.
[38] T. Lyche and K. Morken, "Making the Oslo Algorithm More Efficient," SIAM J. Numerical Analysis, vol. 23, no. 3, pp. 663-675, 1986.
[39] D. McLain, "Two Dimensional Interpolation from Random Data," Comp. J., vol. 19, pp. 178-181, 1976.
[40] D. Mitchell, "Generating Antialiased Images at Low Sampling Densities," Computer Graphics, vol. 21, no. 4, pp. 65-72, 1987.
[41] G.M. Nielson, “Scattered Data Modeling,” IEEE CG&A, Vol. 13, No. 1, Jan. 1993, pp. 60-70.
[42] W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipes in C. Cambridge Univ. Press, 1992.
[43] D. Ruprecht and H. Muller, "Image Warping with Scattered Data Interpolation," IEEE Computer Graphics and Applications, vol. 15, no. 2, Mar. 1995, pp. 37-43.
[44] L. Schumaker, "Fitting Surfaces to Scattered Data," Approximation Theory II, C. Chui, L. Schumaker, and G. Lorentz, eds., pp. 203-268.New York: Wiley, 1976.
[45] D. Shepard, “A Two-Dimensional Interpolation Function for Irregularly Spaced Data,” Proc. 23rd National ACM Conf., ACM Press, New York, 1968, pp. 517-524.
[46] R. Szeliski, "Fast Surface Interpolation Using Hierarchical Basis Functions," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 12, no. 6, pp. 513-528, June 1990.
[47] D. Terzopoulos, "Multilevel Computational Processes for Visual Surface Reconstruction," Computer Vision, Graphics, and Image Processing, vol. 24, pp. 52-96, 1983.
[48] D. Terzopoulos, "Image Analysis Using Multigrid Relaxation Methods," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 8, no. 2, pp. 129-139, Mar. 1986.
[49] D. Terzopoulos, "Regularization of Inverse Visual Problems Involving Discontinuities," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 8, no. 4, pp. 413-424, 1986.
[50] W. Welch and A. Witkin, “Variational Surface Modeling,” Proc. ACM SIGGRAPH, pp. 157-166, July 1992.
[51] G. Wolberg, Digital Image Warping, IEEE CS Press, 1990.

Index Terms:
Scattered data interpolation, multilevel B-splines, data approximation.
Seungyong Lee, George Wolberg, Sung Yong Shin, "Scattered Data Interpolation with Multilevel B-Splines," IEEE Transactions on Visualization and Computer Graphics, vol. 3, no. 3, pp. 228-244, July-Sept. 1997, doi:10.1109/2945.620490
Usage of this product signifies your acceptance of the Terms of Use.