Primal-dual interior point optimization for a regularized reconstruction of NMR relaxation time distributions

This paper deals with the reconstruction of relaxation time distributions in Nuclear Magnetic Resonance (NMR) spectroscopy. This large scale and ill-posed inverse problem is solved by the iterative minimization of a regularized objective function allowing to encode some prior assumptions on the sought distribution. The numerical optimization of the criterion is performed using a primal-dual interior point algorithm allowing to handle the non-negativity constraint. The performances of the proposed approach are illustrated through the processing of real data from a two-dimensional NMR experiment.

[1]  S. Moussaoui,et al.  Majorize–minimize linesearch for inversion methods involving barrier function optimization , 2012 .

[2]  T. Tsuchiya,et al.  On the formulation and theory of the Newton interior-point method for nonlinear programming , 1996 .

[3]  P. Armand,et al.  A Feasible BFGS Interior Point Algorithm for Solving Strongly Convex Minimization Problems , 1999 .

[4]  E. Sternin.  Use of inverse theory algorithms in the analysis of biomembrane NMR data. , 2007, Methods in molecular biology.

[5]  R. Henkelman,et al.  Quantitative Two‐Dimensional time Correlation Relaxometry , 1991, Magnetic resonance in medicine.

[6]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[7]  P. Toint,et al.  RAL-TR-96-096 A primal-dual algorithm for minimizinga non-convex function subject to boundand linear equality constraintsAndrew , 1996 .

[8]  S. Nash.  A survey of truncated-Newton methods , 2000 .

[9]  Émilie Chouzenoux,et al.  Efficient Maximum Entropy Reconstruction of Nuclear Magnetic Resonance T1-T2 Spectra , 2010, IEEE Transactions on Signal Processing.

[10]  James P. Butler,et al.  Estimating Solutions of First Kind Integral Equations with Nonnegative Constraints and Optimal Smoothing , 1981 .

[11]  Paul Armand,et al.  A Feasible BFGS Interior Point Algorithm for Solving Convex Minimization Problems , 2000, SIAM J. Optim..

[12]  W. Marsden.  I and J , 2012 .

[13]  Philippe Marchal,et al.  Chapter 10 Continuous relaxation time distribution decomposition by mem , 1996 .

[14]  P. J. Hore,et al.  Signal treatment and signal analysis in NMR , 1997 .

[15]  Lalitha Venkataramanan,et al.  Solving Fredholm integrals of the first kind with tensor product structure in 2 and 2.5 dimensions , 2002, IEEE Trans. Signal Process..

[16]  G. Bodenhausen,et al.  Principles of nuclear magnetic resonance in one and two dimensions , 1987 .