Variational image reconstruction from arbitrarily spaced samples: a fast multiresolution spline solution

We propose a novel method for image reconstruction from nonuniform samples with no constraints on their locations. We adopt a variational approach where the reconstruction is formulated as the minimizer of a cost that is a weighted sum of two terms: 1) the sum of squared errors at the specified points and 2) a quadratic functional that penalizes the lack of smoothness. We search for a solution that is a uniform spline and show how it can be determined by solving a large, sparse system of linear equations. We interpret the solution of our approach as an approximation of the analytical solution that involves radial basis functions and demonstrate the computational advantages of our approach. Using the two-scale relation for B-splines, we derive an algebraic relation that links together the linear systems of equations specifying reconstructions at different levels of resolution. We use this relation to develop a fast multigrid algorithm. We demonstrate the effectiveness of our approach on some image reconstruction examples.

[1]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[2]  Richard Szeliski,et al.  Fast Surface Interpolation Using Hierarchical Basis Functions , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[3]  Jelena Kovacevic,et al.  Wavelets and Subband Coding , 2013, Prentice Hall Signal Processing Series.

[4]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[5]  Hans G. Feichtinger,et al.  Fast iterative and non-iterative reconstruction of band-limited functions from irregular sampling values , 1991, [Proceedings] ICASSP 91: 1991 International Conference on Acoustics, Speech, and Signal Processing.

[6]  Thomas Strohmer,et al.  Computationally attractive reconstruction of bandlimited images from irregular samples , 1997, IEEE Trans. Image Process..

[7]  Y. Meyer,et al.  Wavelets and Filter Banks , 1991 .

[8]  Richard K. Beatson,et al.  Fast Solution of the Radial Basis Function Interpolation Equations: Domain Decomposition Methods , 2000, SIAM J. Sci. Comput..

[9]  Richard K. Beatson,et al.  Surface interpolation with radial basis functions for medical imaging , 1997, IEEE Transactions on Medical Imaging.

[10]  Thierry Blu,et al.  Least-squares image resizing using finite differences , 2001, IEEE Trans. Image Process..

[11]  W. Bowen,et al.  Philadelphia , 1892 .

[12]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[13]  Fred L. Bookstein,et al.  Morphometric Tools for Landmark Data. , 1998 .

[14]  Jean Duchon,et al.  Splines minimizing rotation-invariant semi-norms in Sobolev spaces , 1976, Constructive Theory of Functions of Several Variables.

[15]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[16]  Michael Unser,et al.  Splines: a perfect fit for signal and image processing , 1999, IEEE Signal Process. Mag..

[17]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[18]  Yonina C. Eldar,et al.  Filterbank reconstruction of bandlimited signals from nonuniform and generalized samples , 2000, IEEE Trans. Signal Process..

[19]  Thierry Blu,et al.  MOMS: maximal-order interpolation of minimal support , 2001, IEEE Trans. Image Process..

[20]  Michael Unser,et al.  Multigrid adaptive image processing , 1995, Proceedings., International Conference on Image Processing.

[21]  Wen-Thong Chang,et al.  Wavelet transform in scattering data interpolation , 1993 .

[22]  A. Jacobson,et al.  Morphometric tools for landmark data , 1993 .

[23]  Akram Aldroubi,et al.  Nonuniform Sampling and Reconstruction in Shift-Invariant Spaces , 2001, SIAM Rev..

[24]  Michael Unser,et al.  Cardinal spline filters: Stability and convergence to the ideal sinc interpolator , 1992, Signal Process..

[25]  G. Wahba.  Spline models for observational data , 1990 .

[26]  Mark Berman.  Automated Smoothing of Image and Other Regularly Spaced Data , 1994, IEEE Trans. Pattern Anal. Mach. Intell..

[27]  Saad T. Bakir,et al.  Nonparametric Regression and Spline Smoothing , 2000, Technometrics.

[28]  T. Strohmer,et al.  Fast scattered data approximation with Neumann and other boundary conditions , 2004 .

[29]  C. Micchelli.  Interpolation of scattered data: Distance matrices and conditionally positive definite functions , 1986 .

[30]  Richard Szeliski,et al.  Fast surface interpolation using hierarchical basis functions , 1989, Proceedings CVPR '89: IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[31]  I J Schoenberg,et al.  SPLINE FUNCTIONS AND THE PROBLEM OF GRADUATION. , 1964, Proceedings of the National Academy of Sciences of the United States of America.

[32]  J. Meinguet.  Multivariate interpolation at arbitrary points made simple , 1979 .

[33]  Michael Unser.  Approximation power of biorthogonal wavelet expansions , 1996, IEEE Trans. Signal Process..