Tikhonov replacement functionals for iteratively solving nonlinear operator equations

We shall be concerned with the construction of Tikhonov-based iteration schemes for solving nonlinear operator equations. In particular, we are interested in algorithms for the computation of a minimizer of the Tikhonov functional. To this end, we introduce a replacement functional, that has much better properties than the classical Tikhonov functional with nonlinear operator. Namely, the replacement functional is globally convex and can effectively be minimized by a fixed point iteration. On the basis of the minimizers of the replacement functional, we introduce an iterative algorithm that converges towards a critical point of the Tikhonov functional, and under additional assumptions for the nonlinear operator F, to a global minimizer. Combining our iterative strategy with an appropriate parameter selection rule, we obtain convergence and convergence rates. The performance of the resulting numerical scheme is demonstrated by solving the nonlinear inverse SPECT (single photon emission computerized tomography) problem.

[1]  J. Timmer,et al.  Asymptotic scaling laws for precision of parameter estimates in dynamical systems , 2003 .

[2]  Til Aach,et al.  Estimation of multiple orientations in multidimensional signals , 2004, 2004 International Conference on Image Processing, 2004. ICIP '04..

[3]  Otmar Scherzer,et al.  Factors influencing the ill-posedness of nonlinear problems , 1994 .

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

[5]  H. Voss,et al.  Parameter estimation in nonlinear delayed feedback systems from noisy data , 2002 .

[6]  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.

[7]  W. Panzer,et al.  An adaptive algorithm for the detection of microcalcifications in simulated low-dose mammography. , 2003, Physics in medicine and biology.

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

[9]  A. Bakushinskii.  The problem of the convergence of the iteratively regularized Gauss-Newton method , 1992 .

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

[11]  M. Hanke.  Regularizing properties of a truncated Newton-CG algorithm for nonlinear inverse problems , 1997 .

[12]  A. Neubauer,et al.  On convergence rates for the Iteratively regularized Gauss-Newton method , 1997 .

[13]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

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

[15]  Thorsten Hohage,et al.  Logarithmic convergence rates of the iteratively regularized Gauss - Newton method for an inverse potential and an inverse scattering problem , 1997 .

[16]  U. Clarenz,et al.  Towards fast non-rigid registration , 2003 .

[17]  Otmar Scherzer,et al.  The use of Morozov's discrepancy principle for Tikhonov regularization for solving nonlinear ill-posed problems , 1993, Computing.

[18]  H. Voss,et al.  Phase synchronization from noisy univariate signals. , 2004, Physical review letters.

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

[20]  M. Rumpf,et al.  Morphological image sequence processing , 2004 .

[21]  Til Aach,et al.  Estimation of multiple local orientations in image signals , 2004, 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[22]  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).

[23]  Frank Natterer,et al.  Inversion of the attenuated Radon transform , 2001 .

[24]  H. Kantz,et al.  FREQUENCY DOMAIN CONDITIONS FOR THE EXISTENCE OF ALMOST PERIODIC SOLUTIONS IN EVOLUTIONARY VARIATIONAL INEQUALITIES , 2004 .

[25]  D. Lorenz.  Variational Denoising in Besov Spaces and Interpolation of Hard and Soft Wavelet Shrinkage , 2003 .

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

[27]  R. Ramlau.  A modified landweber method for inverse problems , 1999 .

[28]  M. Hanke,et al.  A convergence analysis of the Landweber iteration for nonlinear ill-posed problems , 1995 .

[29]  Otmar Scherzer,et al.  A convergence analysis of iterative methods for the solution of nonlinear ill-posed problems under affinely invariant conditions , 1998 .

[30]  Karsten Keller,et al.  Symbolic Analysis of High-Dimensional Time Series , 2003, Int. J. Bifurc. Chaos.

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

[32]  B. Kaltenbacher.  Some Newton-type methods for the regularization of nonlinear ill-posed problems , 1997 .

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

[34]  M. Hanke.  A regularizing Levenberg - Marquardt scheme, with applications to inverse groundwater filtration problems , 1997 .

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