DEcomposition and Component Analysis of Exponential Signals (DECAES)

Table of contents

Introduction

DECAES provides tools for decomposing multi-exponential signals which arise from multi spin-echo magnetic resonance imaging (MRI) scans into exponential components. The main decomposition method used is an inverse Laplace transform-based technique which involves solving the regularized nonnegative least squares (NNLS) inverse problem

\[X_{\mu} = \underset{x \ge 0}{\operatorname{argmin}}\; ||Ax - b||_2^2 + \mu^2 ||x||_2^2\]

where $b$ is the signal magnitude data, $A$ is a matrix of exponential decay bases, and $\mu$ is a regularization parameter. $A$ is constructed using the extended phase graph algorithm with stimulated echo correction. The columns of $A$ are exponential bases with differing characteristic $T_2$ decay times $T_{2, j}$.

The output $X_{\mu}$ is the spectrum of (nonnegative) exponential decay amplitudes. Amplitude $X_{\mu, j}$ of the spectrum $X_{\mu}$ is therefore interpreted physically as the amount of the signal $b$ which decays with time constant $T_{2, j}$. For this reason, the spectrum $X_{\mu}$ is commonly referred to as the $T_2$ distribution. DECAES provides methods for computing $T_2$-distributions.

Installation

Using Julia v1.9 or later you can install DECAES as follows:

$ julia --project=@decaes -e 'import Pkg; Pkg.add("DECAES"); Pkg.build("DECAES")'

This will add DECAES.jl to a named Julia project environment separate from your global environment, and build the decaes executable at ~/.julia/bin for running DECAES from the command line.

Updating DECAES

DECAES can similarly be updated to the latest version as follows:

$ julia --project=@decaes -e 'import Pkg; Pkg.update("DECAES"); Pkg.build("DECAES")'

Myelin water imaging

Myelin water imaging (MWI) is an MRI technique used to visualize the myelin water contained within the sheaths of myelinated axons within the body, such as within the brain's white matter.

Through analysing $T_2$-distributions computed from multi spin-echo MRI scans, one can separate the contribution due to the myelin water from the intra- and extra-cellular water and compute the myelin water fraction (MWF). The MWF describes the fraction of water trapped between myelin lipid bilayers relative to the total water in the region. DECAES provides methods for computing the MWF.

MWI was pioneered at the University of British Columbia by Alex MacKay and Ken Whittal.

Basics of myelin water imaging:

Validation of myelin water imaging:

Some applications of myelin water imaging:

Acknowledgements

  • Porting to Julia was done by Jonathan Doucette (email: jdoucette@physics.ubc.ca) in November 2019. This work was funded by NSERC (016-05371) and CIHR (RN382474-418628) under PI Alexander Rauscher at the University of British Columbia
  • Christian Kames (email: ckames@physics.ubc.ca) contributed to optimizing the Julia port for both speed and memory efficiency, as well as writing the PAR/XML/REC file reader used internally
  • Original MATLAB code was written by Thomas Prasloski (email: tprasloski@gmail.com). Modifications to the MATLAB code were made by Vanessa Wiggermann to enable processing on various MATLAB versions in February 2019. The Julia port is based on this modified version

Citing this work

Z Med Phys

If you use DECAES in your research, please cite the following:

@article{DECAES.jl-2020,
  title = {{{DECAES}} - {{DEcomposition}} and {{Component Analysis}} of {{Exponential Signals}}},
  author = {Doucette, Jonathan and Kames, Christian and Rauscher, Alexander},
  year = {2020},
  month = may,
  issn = {1876-4436},
  doi = {10.1016/j.zemedi.2020.04.001},
  journal = {Zeitschrift Fur Medizinische Physik},
  keywords = {Brain,Luminal Water Imaging,MRI,Myelin Water Imaging,Prostate},
  language = {eng},
  pmid = {32451148}
}