Numerical Phantoms
Gary Zhang1

1Computer Science, University College London

Synopsis

Numerical phantoms have played and will continue to play an important role in the development and validation of advanced diffusion MRI techniques. They complement biological phantoms (in vitro, ex vivo and in vivo) with their controllability and physical phantoms with their flexibility. This talk will review the aspects of diffusion MRI techniques that have benefited from validation with numerical phantoms and the range of numerical phantoms currently available. Examples of using numerical phantoms for validating the mapping of tissue microstructure and structural connectivity in the brain will be presented.

Target Audience

Scientists and clinicians interested in learning what has been done in validation of commonly used MRI diffusion methods (tractography and tissue microstructure) and applying those methods in their own scientific inquiries and the use of these measurements in both basic research and clinical applications.

Educational Objectives

This talk aims to provide its audience with the following:

1. Overview of the aspects of diffusion MRI techniques that can be validated with numerical phantoms.

2. Overview of the range of numerical phantoms that are currently available.

3. Examples of using numerical phantoms for validating diffusion MRI techniques.

Overview

Over the past two decades, diffusion MRI has firmly established itself as the modality of choice for the non-invasive mapping of tissue microstructure and structural connectivity in the brain. Its success is a testament to the relentless efforts within the diffusion MRI research community in advancing the state-of-the-art to mitigate the limitations of earlier developments. A common denominator among all these advancements has been the adoption of numerical phantoms to provide rigorous validation.

Numerical phantoms play a unique role in validation that is complementary of other forms of phantoms (physical, in vitro, ex vivo and in vivo). Although they necessarily represent a model of the real world based on our current understanding, they provide a framework for validation that is much more controlled (with known ground truth) than in vitro, ex vivo and in vivo phantoms, and is much more flexible than physical phantoms.

This talk will give an overview of the numerical phantoms that are currently available for validating diffusion MRI techniques. It will give examples of using numerical phantoms to validate

1. microstructure imaging,

2. tractography/structural connectivity,

3. pre-processing and group comparison.

Validation of Microstructure Imaging

This is an area in which numerical phantoms have played a major role. They allow recent model-based microstructure imaging techniques to be evaluated in a controlled and flexible manner that is not possible otherwise.

Model-based microstructure imaging techniques (Stanisz et al 1997, Assaf et al 2004, Jespersen 2007, Alexander et al 2010, Fieremans et al 2011, Zhang et al 2012, Scherrer et al In Press) work by prescribing a model that uses a set of salient microstructural features (e.g. axon density and diameter) to predict the measured diffusion MRI signals, and when presented with a set of measurements, the prescribed model is used to estimate the salient features from the measurements. To be of practical use, the models chosen generally need to be sufficiently simple, thus representing a crude approximation to the reality. To understand the limitations of such approximations, numerical phantoms have been developed. They implement much more accurate models that are computationally too expensive for practical use but provide the ground truth mapping between microstructural features and the diffusion signal.

Numerical phantoms for this purpose can be divided into several types based on the models implemented. The most common class uses Monte Carlo simulation of diffusion (Szafer et al 1995, Ford et al 1997, Liu et al 2004, Peled 2007, Balls et al 2009, Hall et al 2009, Nilsson et al 2009, Fieremans et al 2010, Landman et al 2010, Yeh et al 2013), with the system implemented in the Open Source Camino software being most widely used (Hall et al 2009). This class of phantoms has been used to validate the estimation of axon density and diameter (Alexander et al 2010), crossing-fibre resolution (Ramirez-Manzanares et al 2011), exchange (Fieremans et al 2010, Nilsson et al 2010), neurite beading (Budde et al 2010) and axonal undulation (Nilsson et al 2011), time-dependent extra-axonal diffusion (Burcaw et al 2015, Lam et al 2015, De Santis et al 2016), and g-ratio mapping (Stikov et al 2011).

An alternative approach is to directly solve the diffusion equations using either finite differences (Chin et al 2002, Hwang et al 2003, Russell et al 2012, Li et al 2014) or finite element method (Hagslätt et al 2003, Moroney et al 2013, Nguyen et al 2014, Beltrachini et al 2015).

Validation of Tractography/Structural Connectivity

The tractography evaluation is dominated by physical and biological phantoms but numerical phantoms still have a role to play with its flexibility (Leemans et al 2005, Close et al 2009, Neher et al 2014). In particular, Fiberfox (Neher et al 2014), combined with Tractometer (Côté et al 2013), has recently been chosen for the Tractography Challenge 2015 hosted by the ISMRM Diffusion Study Group.

Validation of Pre-processing and Group Comparison

Numerical phantoms have also been developed to validate common pre-processing and group comparison tools. Graham et al (2016) have recently developed a tool for simulating realistic imaging artefacts in diffusion MRI and demonstrated its application to the evaluation of eddy-current distortion and motion correction techniques. Van Hecke et al (2007) have used synthetic diffusion MRI data (Leemans et al 2005) to evaluate a diffusion tensor image registration framework. The authors have also developed a framework based on simulated diffusion MRI data to assess various aspect of population-based group comparison using diffusion MRI (Van Hecke et al 2008-2011).

Acknowledgements

No acknowledgement found.

References

Alexander et al, Orientationally invariant indices of axon diameter and density from diffusion MRI, NeuroImage (2010) 52:1374-89.

Assaf et al, New Modeling and Experimental Framework to Characterize Hindered and Restricted Water Diffusion in Brain White Matter, MRM (2004) 52:965-78.

Balls et al, A Simulation Environment for Diffusion Weighted MR Experiments in Complex Media, MRM (2009) 62:771-8.

Beltrachini et al, A parametric finite element solution of the generalised Bloch–Torrey equation for arbitrary domains, JMR (2015) 259:126-34.

Budde et al, Neurite beading is sufficient to decrease the apparent diffusion coefficient after ischemic stroke, PNAS (2010) 107:14472-7.

Burcaw et al, Mesoscopic structure of neuronal tracts from time-dependent diffusion, NeuroImage (2015) 114:18-37.

Chin et al, Biexponential Diffusion Attenuation in the Rat Spinal Cord: Computer Simulations Based on Anatomic Images of Axonal Architecture, MRM (2002) 47:455-60.

Close et al, A software tool to generate simulated white matter structures for the assessment of fibre-tracking algorithms, NeuroImage (2009) 47:1288-1300.

De Santis et al, Including diffusion time dependence in the extra-axonal space improves in vivo estimates of axonal diameter and density in human white matter, NeuroImage (2016) 130:91-103.

Fieremans et al, Monte Carlo study of a two-compartment exchange model of diffusion, NMR Biomed (2010) 23:711-24.

Fieremans et al, White matter characterization with diffusional kurtosis imaging, NeuroImage (2011) 58:177-88.

Ford and Hackney, Numerical Model for Calculation of Apparent Diffusion Coefficients (ADC) in Permeable Cylinders-Comparison with Measured ADC in Spinal Cord White Matter, MRM (1997) 37:387-94.

Graham et al, Realistic simulation of artefacts in diffusion MRI for validating post-processing correction techniques, NeuroImage (2016) 125:1079-94.

Hagslätt et al, Predictions of pulsed field gradient NMR echo-decays for molecules diffusing in various restrictive geometries. Simulations of diffusion propagators based on a finite element method, JMR (2003) 161:138-47.

Hall and Alexander, Convergence and Parameter Choice for Monte-Carlo Simulations of Diffusion MRI, IEEE TMI (2009) 28:1354-64.

Hwang et al, An Image-Based Finite Difference Model for Simulating Restricted Diffusion, MRM (2003) 50:373-82.

Jespersen et al, Modeling dendrite density from magnetic resonance diffusion measurements, NeuroImage (2007) 34:1473-86.

Lam et al, A Model for Extra-Axonal Diffusion Spectra with Frequency-Dependent Restriction, MRM (2015) 73:2306-20.

Landman et al, Complex geometric models of diffusion and relaxation in healthy and damaged white matter, NMR Biomed (2010) 23:152-62.

Leemans et al, Mathematical Framework for Simulating Diffusion Tensor MR Neural Fiber Bundles, MRM (2005) 53:944-53.

Li et al, Numerical simulation of diffusion MRI signals using an adaptive time-stepping method, Phys Med Biol (2014) 59:441-54.

Liu et al, Characterizing Non-Gaussian Diffusion by Using Generalized Diffusion Tensors, MRM (2004) 51:924-37.

Moroney et al, Numerical analysis of NMR diffusion measurements in the short gradient pulse limit, JMR (2013) 234:165-75.

Neher et al, Fiberfox: Facilitating the Creation of Realistic White Matter Software Phantoms, MRM (2014) 72:1460-70.

Nguyen et al, A finite elements method to solve the Bloch–Torrey equation applied to diffusion magnetic resonance imaging, J Comp Phys (2014) 263:283-302.

Nilsson et al, On the effects of a varied diffusion time in vivo: is the diffusion in white matter restricted? MRI (2009) 27:176-87.

Nilsson et al, Evaluating the accuracy and precision of a two-compartment Kärger model using Monte Carlo simulations, JMR (2010) 206:59-67.

Nilsson et al, The importance of axonal undulation in diffusion MR measurements: a Monte Carlo simulation study, NMR Biomed (2012) 25:795-805.

Peled, New Perspectives on the Sources of White Matter DTI Signal, IEEE TMI (2007) 26:1448-55.

Ramirez-Manzanares et al, Resolving axon fiber crossings at clinical b-values: An evaluation study, Medical Physics (2011) 38:5239-53.

Russell et al, A finite difference method with periodic boundary conditions for simulations of diffusion-weighted magnetic resonance experiments in tissue, Phys Med Biol (2012) 57:N35-46.

Scherrer et al, Characterizing brain tissue by assessment of the distribution of anisotropic microstructural environments in diffusion-compartment imaging (DIAMOND), MRM (In Press).

Stanisz et al, An Analytical Model of Restricted Diffusion in Bovine Optic Nerve, MRM (1997) 37:103-11.

Szafer et al, Theoretical Model for Water Diffusion in Tissues, MRM (1995) 33:697-712.

Van Hecke et al, Nonrigid Coregistration of Diffusion Tensor Images Using a Viscous Fluid Model and Mutual Information, IEEE TMI (2007) 26:1598-612.

Van Hecke et al, On the construction of an inter-subject diffusion tensor magnetic resonance atlas of the healthy human brain, NeuroImage (2008) 43:69-80.

Van Hecke et al, On the construction of a ground truth framework for evaluating voxel-based diffusion tensor MRI analysis methods, NeuroImage (2009) 46:692-707.

Van Hecke et al, Comparing Isotropic and Anisotropic Smoothing for Voxel-Based DTI Analyses: A Simulation Study, HBM (2010) 31:98-114.

Van Hecke et al, The effect of template selection on diffusion tensor voxel-based analysis results, NeuroImage (2011) 55:566-73.

Xu et al, Numerical study of water diffusion in biological tissues using an improved finite difference method, Phys Med Biol (2007) 52:N111-26.

Yeh et al, Diffusion Microscopist Simulator: A General Monte Carlo Simulation System for Diffusion Magnetic Resonance Imaging, PLoS One (2013) e76626.

Zhang et al, NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain, NeuroImage (2012) 61:1000-16.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)