DocumentCode :
1135316
Title :
A Fast Multilevel Algorithm for Wavelet-Regularized Image Restoration
Author :
Vonesch, Cédric ; Unser, Michael
Author_Institution :
EPFL, Biomed. Imaging Group, Lausanne
Volume :
18
Issue :
3
fYear :
2009
fDate :
3/1/2009 12:00:00 AM
Firstpage :
509
Lastpage :
523
Abstract :
We present a multilevel extension of the popular ldquothresholded Landweberrdquo algorithm for wavelet-regularized image restoration that yields an order of magnitude speed improvement over the standard fixed-scale implementation. The method is generic and targeted towards large-scale linear inverse problems, such as 3-D deconvolution microscopy. The algorithm is derived within the framework of bound optimization. The key idea is to successively update the coefficients in the various wavelet channels using fixed, subband-adapted iteration parameters (step sizes and threshold levels). The optimization problem is solved efficiently via a proper chaining of basic iteration modules. The higher level description of the algorithm is similar to that of a multigrid solver for PDEs, but there is one fundamental difference: the latter iterates though a sequence of multiresolution versions of the original problem, while, in our case, we cycle through the wavelet subspaces corresponding to the difference between successive approximations. This strategy is motivated by the special structure of the problem and the preconditioning properties of the wavelet representation. We establish that the solution of the restoration problem corresponds to a fixed point of our multilevel optimizer. We also provide experimental evidence that the improvement in convergence rate is essentially determined by the (unconstrained) linear part of the algorithm, irrespective of the type of wavelet. Finally, we illustrate the technique with some image deconvolution examples, including some real 3-D fluorescence microscopy data.
Keywords :
image resolution; image restoration; image sequences; iterative methods; wavelet transforms; 3D deconvolution microscopy; 3D fluorescence microscopy data; fast multilevel algorithm; image sequence; magnitude speed improvement; multiresolution versions; standard fixed-scale implementation; subband-adapted iteration parameters; thresholded Landweber algorithm; wavelet-regularized image restoration; 3-D; $ell_{1}$ -regularization; Bound optimization; confocal; convergence acceleration; deconvolution; fast; fluorescence; inverse problems; majorize-minimize; microscopy; multigrid; multilevel; multiresolution; multiscale; nonlinear; optimization transfer; preconditioning; reconstruction; restoration; sparsity; surrogate optimization; variational; wavelets; widefield; Algorithms; Artificial Intelligence; Image Enhancement; Image Interpretation, Computer-Assisted; Imaging, Three-Dimensional; Microscopy, Fluorescence; Pattern Recognition, Automated; Reproducibility of Results; Sensitivity and Specificity; Signal Processing, Computer-Assisted;
fLanguage :
English
Journal_Title :
Image Processing, IEEE Transactions on
Publisher :
ieee
ISSN :
1057-7149
Type :
jour
DOI :
10.1109/TIP.2008.2008073
Filename :
4770145
Link To Document :
بازگشت