A Tikhonov-based projection iteration for nonlinear Ill-posed problems with sparsity constraints

In this paper, we consider nonlinear inverse problems where the solution is assumed to have a sparse expansion with respect to a preassigned basis or frame. We develop a scheme which allows to minimize a Tikhonov functional where the usual quadratic regularization term is replaced by a one-homogeneous (typically weighted ℓp) penalty on the coefficients (or isometrically transformed coefficients) of such expansions. For (p < 2), the regularized solution will have a sparser expansion with respect to the basis or frame under consideration. The computation of the regularized solution amounts in our setting to a Landweber-fixed-point iteration with a projection applied in each fixed-point iteration step. The performance of the resulting numerical scheme is demonstrated by solving the nonlinear inverse single photon emission computerized tomography (SPECT) problem.

[1]  I. Daubechies.  Ten Lectures on Wavelets , 1992 .

[2]  Ingrid Daubechies,et al.  Wavelet-based image decomposition by variational functionals , 2004, SPIE Optics East.

[3]  Ronny Ramlau,et al.  On the use of fixed point iterations for the regularization of nonlinear ill-posed problems , 2005 .

[4]  Ingrid Daubechies,et al.  Variational image restoration by means of wavelets: simultaneous decomposition , 2005 .

[5]  Stanley Osher,et al.  Modeling Textures with Total Variation Minimization and Oscillating Patterns in Image Processing , 2003, J. Sci. Comput..

[6]  Ingrid Daubechies,et al.  Variational image restoration by means of wavelets , 2005 .

[7]  R. Clackdoyle,et al.  Accurate attenuation correction in SPECT imaging using optimization of bilinear functions and assuming an unknown spatially-varying attenuation distribution , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[8]  R. DeVore,et al.  Nonlinear Approximation and the Space BV(R2) , 1999 .

[9]  P. G. Ciarlet,et al.  Introduction to Numerical Linear Algebra and Optimisation , 1989 .

[10]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[11]  Emmanuel J. Candès,et al.  New multiscale transforms, minimum total variation synthesis: applications to edge-preserving image reconstruction , 2002, Signal Process..

[12]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[13]  Volker Dicken,et al.  Simultaneous Activity and Attenuation Reconstruction in Single Photon Emission Computed Tomography, a Nonlinear Ill-Posed Problem , 1997 .

[14]  S. H. Manglos,et al.  Constrained IntraSPECT reconstruction from SPECT projections , 1993, 1993 IEEE Conference Record Nuclear Science Symposium and Medical Imaging Conference.

[15]  Volker Dicken,et al.  A new approach towards simultaneous activity and attenuation reconstruction in emission tomography , 1999 .

[16]  Yves Meyer,et al.  Oscillating Patterns in Image Processing and Nonlinear Evolution Equations: The Fifteenth Dean Jacqueline B. Lewis Memorial Lectures , 2001 .

[17]  I. Daubechies,et al.  Harmonic analysis of the space BV. , 2003 .

[18]  Y. Censor,et al.  A New Approach to the Emission Computerized Tomography Problem: Simultaneous Calculation of Attenuation and Activity Coefficients , 1979, IEEE Transactions on Nuclear Science.

[19]  Grant T. Gullberg,et al.  Toward accurate attenuation correction in SPECT without transmission measurements , 1997, IEEE Transactions on Medical Imaging.

[20]  Y. Meyer.  Oscillating Patterns in Some Nonlinear Evolution Equations , 2006 .

[21]  R. Ramlau.  MOROZOV'S DISCREPANCY PRINCIPLE FOR TIKHONOV-REGULARIZATION OF NONLINEAR OPERATORS , 2002 .

[22]  G. Teschke,et al.  Tikhonov replacement functionals for iteratively solving nonlinear operator equations , 2005 .

[23]  Ronny Ramlau,et al.  A steepest descent algorithm for the global minimization of the Tikhonov functional , 2002 .

[24]  Ronny Ramlau,et al.  TIGRA—an iterative algorithm for regularizing nonlinear ill-posed problems , 2003 .