Simple MATLAB and Python scripts for multi‐exponential analysis

Multi‐exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi‐exponential analysis methods for 1D and 2D data and example test problems and experiments. Regularized least‐squ...

Full description

Saved in:
Bibliographic Details
Published in:Magnetic resonance in chemistry Vol. 62; no. 10; pp. 698 - 711
Main Authors: Afrough, Armin, Mokhtari, Rasoul, Feilberg, Karen L.
Format: Journal Article
Language:English
Published: England Wiley Subscription Services, Inc 01.10.2024
Subjects:
ISSN:0749-1581, 1097-458X, 1097-458X
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Multi‐exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi‐exponential analysis methods for 1D and 2D data and example test problems and experiments. Regularized least‐squares solutions provide production‐quality outputs with robust stopping rules in ~5 and ~20 lines of code for 1D and 2D inversions, respectively. The software provides an open‐architecture simple solution for transforming exponential decay data to the distribution of their decay lifetimes. Examples from magnetic resonance relaxation of a complex fluid, a Danish North Sea crude oil, and fluid mixtures in porous materials—brine/crude oil mixture in North Sea reservoir chalk—are presented. Developed codes may be incorporated in other software or directly used by other researchers, in magnetic resonance relaxation, diffusion, and imaging or other physical phenomena that require multi‐exponential analysis. Multi‐exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi‐exponential analysis methods for 1D and 2D data and example test problems and experiments. Regularized least‐squares solutions provide production‐quality outputs with robust stopping rules in ~5 and ~20 lines of code for 1D and 2D inversions, respectively. The software provides an open‐architecture simple solution for transforming exponential decay data to the distribution of their decay lifetimes.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
content type line 23
ISSN:0749-1581
1097-458X
1097-458X
DOI:10.1002/mrc.5453