Interpolation in the Presence of Noise 
Investigators: Michael Unser 

Summary: We derive splinefitting algorithms that are optimized for noisy data. 

Interpolation is a crucial operation in a variety of medical imaging tasks (e.g., image registration, tomographic reconstruction, and more). Presently, most techniques are optimized for a noisefree scenario. Our goal here is to optimize the prefilter that yields the Bspline expansion coefficients so that it can handle noisy data. 

We are investigating extensions of splinetype interpolation algorithms for the nonideal case where the signal samples are corrupted by noise.
The first aspect is that the spline space itself can be optimized depending on the smoothness properties or the statistics of the class of signal of interest. When the samples are corrupted by noise, it is appropriate to apply a smoothingspline estimator with a regularization parameter that is set inversely proportional to the signaltonoise ratio.
The second aspect is the design of the "optimal" prefilter that yields the expansion coefficient of the signal in a shiftinvariant (or spline) subspace spanned by the integer shifts of a given generating function. In our formulation, we treat in a unified way the interpolation problem from ideal, noisy samples, and the deconvolution problem in which the signal is filtered prior to sampling. We have proposed several alternative approaches to designing the correction filter, which differ in their assumptions on the signal and noise. In particular, we have adapted the classical deconvolution solutions (leastsquares, Tikhonov, and Wiener) to our particular situation and also proposed new methods that are optimal in a minimax sense. The solutions often have a similar structure and can be computed simply and efficiently by digital filtering. 

Collaborations: Prof. Yonina Eldar (Technion) 

Past Investigators: Sathish Ramani, Thierry Blu 


Funding: Grant 200020101821 from the Swiss Science Foundation, European HASSIP Network 


[2]  S. Ramani, D. Van De Ville, M. Unser, "NonIdeal Sampling and Adapted Reconstruction Using the Stochastic Matérn Model," Best student paper award, Proceedings of the IEEE ThirtyFirst International Conference on Acoustics, Speech, and Signal Processing (ICASSP'06), Toulouse, France, May 1419, 2006, pp. II73II76.


