Parsimonious discretization for characterizing multi‐exponential decay in magnetic resonance

We address the problem of analyzing noise‐corrupted magnetic resonance transverse decay signals as a superposition of underlying independently decaying monoexponentials of positive amplitude. First, we indicate the manner in which this is an ill‐conditioned inverse problem, rendering the analysis un...

Full description

Saved in:
Bibliographic Details
Published in:NMR in biomedicine 2020-12, Vol.33 (12), p.e4366-n/a
Main Authors: Bonny, Jean‐Marie, Traore, Amidou, Bouhrara, Mustapha, Spencer, Richard G., Pages, Guilhem
Format: Article
Language:English
Subjects:
Biological products
cumulative distribution function
Decay
Discretization
inverse Laplace transform
Inverse problems
Life Sciences
Magnetic resonance
Myelin
Neuroimaging
Noise
Plateaus
Regularization
Resonance
simulation
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We address the problem of analyzing noise‐corrupted magnetic resonance transverse decay signals as a superposition of underlying independently decaying monoexponentials of positive amplitude. First, we indicate the manner in which this is an ill‐conditioned inverse problem, rendering the analysis unstable with respect to noise. Second, we define an approach to this analysis, stabilized solely by the nonnegativity constraint without regularization. This is made possible by appropriate discretization, which is coarser than that often used in practice. Thirdly, we indicate further stabilization by inspecting the plateaus of cumulative distributions. We demonstrate our approach through analysis of simulated myelin water fraction measurements, and compare the accuracy with more conventional approaches. Finally, we apply our method to brain imaging data obtained from a human subject, showing that our approach leads to maps of the myelin water fraction which are much more stable with respect to increasing noise than those obtained with conventional approaches. A discretization coarser than that often used in practice and a positivity constraint suffice for stabilizing the estimation of multi‐exponential amplitudes. The amplitudes are better estimated by the plateaus of the cumulative density function derived from the outputs of the nonnegative least squares algorithm. This parsimonious approach leads to improved myelin water fraction estimation.
ISSN:0952-3480
1099-1492
DOI:10.1002/nbm.4366