BIOMEDICAL IMAGING GROUP (BIG)
Laboratoire d'imagerie biomédicale (LIB)

Fast Multi-Level Reconstruction of Biomedical Images Using Wavelet Sparsity Constraints

M. Unser, C. Vonesch, M. Guerquin-Kern, D. Van De Ville

Approximation and Optimization in Image Restoration and Reconstruction (AOIRR'09), Île de Porquerolles, French Republic, June 8-12, 2009.

Wavelet-domain ℓ1-regularization is a powerful approach for solving inverse problems. In their 2004 landmark paper, Daubechies et al. proved that one could solve such linear inverse problems by means of a "thresholded Landweber" (TL) algorithm [1]. While this iterative procedure is simple to implement, it is known to converge slowly. Here, we present a multilevel version of the algorithm that is inspired from the multigrid techniques used for solving PDEs, but with one important difference: instead of cycling through coarser versions of the problem (REDUCE part of multigrid), the multilevel algorithm cycles through the successive wavelet subspaces. The method works with arbitrary wavelet representations; it typically yields a 10-fold speed increase over the standard TL algorithm, while providing the same restoration quality. We illustrate the applicability of the method to three biomedical image reconstruction problems: the deconvolution of 3D fluorescence micrographs [2], the global reconstruction of dynamic PET from time measurements [3], and the reconstruction of magnetic resonance images from arbitrary (non-uniform) k-space trajectories. We present experimental results with real data sets in all three cases.

References

  1. I. Daubechies, M. Defrise, C. De Mol, "An Iterative Thresholding Algorithm for Linear Inverse Problems with a Sparsity Constraint," Communications on Pure and Applied Mathematics, vol. 57, no. 11, pp. 1413-1457, November 2004.

  2. C. Vonesch, M. Unser, "A Fast Multilevel Algorithm for Wavelet-Regularized Image Restoration," IEEE Transactions on Image Processing, vol. 18, no. 3, pp. 509-523, March 2009.

  3. J. Verhaeghe, D. Van De Ville, I. Khalidov, Y. D'Asseler, I. Lemahieu, M. Unser, "Dynamic PET Reconstruction Using Wavelet Regularization with Adapted Basis Functions," IEEE Transactions on Medical Imaging, vol. 27, no. 7, pp. 943-959, July 2008. 

@INPROCEEDINGS(http://bigwww.epfl.ch/publications/unser0909.html,
AUTHOR="Unser, M. and Vonesch, C. and Guerquin-Kern, M. and Van De
	Ville, D.",
TITLE="Fast Multi-Level Reconstruction of Biomedical Images Using
	Wavelet Sparsity Constraints",
BOOKTITLE="Approximation and Optimization in Image Restoration and
	Reconstruction ({AOIRR'09})",
YEAR="2009",
editor="",
volume="",
series="",
pages="",
address="{\^{I}}le de Porquerolles, French Republic",
month="June 8-12,",
organization="",
publisher="",
note="")
© 2009 ANR. 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 ANR. 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.