Biomedical Imaging GroupSTI
English only   BIG > Publications > Least Squares vs. Interpolation

 Home Page
 News & Events
 Tutorials and Reviews
 Download Algorithms

 All BibTeX References

On the Approximation Power of Convolution-Based Least Squares versus Interpolation

M. Unser, I. Daubechies

IEEE Transactions on Signal Processing, vol. 45, no. 7, pp. 1697-1711, July 1997.

There are many signal processing tasks for which convolution-based continuous signal representations such as splines and wavelets provide an interesting and practical alternative to the more traditional sine-based methods. The coefficients of the corresponding signal approximations are typically obtained by direct sampling (interpolation or quasi-interpolation) or by using least squares techniques that apply a prefilter prior to sampling. We compare the performance of these approaches and provide quantitative error estimates that can be used for the appropriate selection of the sampling step h. Specifically, we review several results in approximation theory with a special emphasis on the Strang-Fix (1971) conditions, which relate the general O(hL) behavior of the error to the ability of the representation to reproduce polynomials of degree n = L - 1. We use this theory to derive pointwise error estimates for the various algorithms and to obtain the asymptotic limit of the L2-error as h tends to zero. We also propose a new improved L2-error bound for the least squares case. In the process, we provide all the relevant bound constants for polynomial splines. Some of our results suggest the existence of an intermediate range of sampling steps where the least squares method is roughly equivalent to an interpolator with twice the order. We present experimental examples that illustrate the theory and confirm the adequacy of our various bound and limit determinations.

AUTHOR="Unser, M. and Daubechies, I.",
TITLE="On the Approximation Power of Convolution-Based Least Squares
        {\textit{versus}} Interpolation",
JOURNAL="{IEEE} Transactions on Signal Processing",

© 1997 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from IEEE.
This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.