Efficient Maximum Entropy Reconstruction of Nuclear Magnetic Resonance T1-T2 Spectra
- 2 September 2010
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Signal Processing
- Vol. 58 (12) , 6040-6051
Abstract
This paper deals with the reconstruction of T1-T2 correlation spectra in nuclear magnetic resonance relaxometry. The ill-posed character and the large size of this inverse problem are the main difficulties to tackle. While maximum entropy is retained as an adequate regularization approach, the choice of an efficient optimization algorithm remains a challenging task. Our proposal is to apply a truncated Newton algorithm with two original features. First, a theoretically sound line search strategy suitable for the entropy function is applied to ensure the convergence of the algorithm. Second, an appropriate preconditioning structure based on a singular value decomposition of the forward model matrix is used to speed up the algorithm convergence. Furthermore, we exploit the specific structures of the observation model and the Hessian of the criterion to reduce the computation cost of the algorithm. The performances of the proposed strategy are illustrated by means of synthetic and real data processing.Keywords
All Related Versions
This publication has 33 references indexed in Scilit:
- Convergence of Conjugate Gradient Methods with a Closed-Form Stepsize FormulaJournal of Optimization Theory and Applications, 2007
- An Expanded Theoretical Treatment of Iteration-Dependent Majorize-Minimize AlgorithmsIEEE Transactions on Image Processing, 2007
- Convergence of line search methods for unconstrained optimizationApplied Mathematics and Computation, 2004
- A Tutorial on MM AlgorithmsThe American Statistician, 2004
- Iterative Methods for Sparse Linear Systems2003
- Global Convergence of Conjugate Gradient Methods without Line SearchAnnals of Operations Research, 2001
- Rank-Deficient and Discrete Ill-Posed Problems1998
- Line search algorithms with guaranteed sufficient decreaseACM Transactions on Mathematical Software, 1994
- Maximum entropy image reconstruction: general algorithmMonthly Notices of the Royal Astronomical Society, 1984
- Estimating Solutions of First Kind Integral Equations with Nonnegative Constraints and Optimal SmoothingSIAM Journal on Numerical Analysis, 1981