SpinDoctor: A MATLAB toolbox for diffusion MRI simulation.
Apparent diffusion coefficient
Bloch-torrey equation
Diffusion magnetic resonance imaging
Finite elements
Simulation
Journal
NeuroImage
ISSN: 1095-9572
Titre abrégé: Neuroimage
Pays: United States
ID NLM: 9215515
Informations de publication
Date de publication:
15 11 2019
15 11 2019
Historique:
received:
26
01
2019
revised:
21
08
2019
accepted:
22
08
2019
pubmed:
31
8
2019
medline:
29
9
2020
entrez:
31
8
2019
Statut:
ppublish
Résumé
The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium can be modeled by the multiple compartment Bloch-Torrey partial differential equation. Under the assumption of negligible water exchange between compartments, the time-dependent apparent diffusion coefficient can be directly computed from the solution of a diffusion equation subject to a time-dependent Neumann boundary condition. This paper describes a publicly available MATLAB toolbox called SpinDoctor that can be used 1) to solve the Bloch-Torrey partial differential equation in order to simulate the diffusion magnetic resonance imaging signal; 2) to solve a diffusion partial differential equation to obtain directly the apparent diffusion coefficient; 3) to compare the simulated apparent diffusion coefficient with a short-time approximation formula. The partial differential equations are solved by P1 finite elements combined with built-in MATLAB routines for solving ordinary differential equations. The finite element mesh generation is performed using an external package called Tetgen. SpinDoctor provides built-in options of including 1) spherical cells with a nucleus; 2) cylindrical cells with a myelin layer; 3) an extra-cellular space enclosed either a) in a box or b) in a tight wrapping around the cells; 4) deformation of canonical cells by bending and twisting; 5) permeable membranes; Built-in diffusion-encoding pulse sequences include the Pulsed Gradient Spin Echo and the Oscillating Gradient Spin Echo. We describe in detail how to use the SpinDoctor toolbox. We validate SpinDoctor simulations using reference signals computed by the Matrix Formalism method. We compare the accuracy and computational time of SpinDoctor simulations with Monte-Carlo simulations and show significant speed-up of SpinDoctor over Monte-Carlo simulations in complex geometries. We also illustrate several extensions of SpinDoctor functionalities, including the incorporation of T
Identifiants
pubmed: 31470126
pii: S1053-8119(19)30711-6
doi: 10.1016/j.neuroimage.2019.116120
pii:
doi:
Types de publication
Journal Article
Research Support, Non-U.S. Gov't
Langues
eng
Sous-ensembles de citation
IM
Pagination
116120Informations de copyright
Copyright © 2019 Elsevier Inc. All rights reserved.