A set of command line tools and a MATLAB (R) toolbox/library to process (diffusion) magnetic resonance imaging (MRI) data. The idea behind MRIToolkit is to serve as a stable distribution platform for the methods I develop as well as to other state-of-the-art technologies in diffusion MRI and beyond.
- The MATLAB (R) toolbox is available here on Github!
- The command line tools will be uploaded soon!
Please, see this guide
Examples of some functionalities can be found in Demos
Main functionalities:
- Complete diffusion MRI pre-processing (signal drift correction, Gibbs ringing correction, motion and eddy currents correction, B-matrix rotation, EPI correction, MPPCA denoising)
- Diffusion Tensor Imaging (DTI) and Diffusion Kurtosis Imaging (DKI) fit, including the MK-curve method;
- Fiber tractography, both with DTI and Constrained Spherical Deconvolution (CSD);
- Spherical deconvolution using the damped Richardson Lucy, the Generalized Richardson Lucy and mFOD methods;
- T1 / T2 quantification
I have coded most functions to accept Python-like name/value argument couples. To know which arguments to specify, just try the MATLAB help as, for instance, "help EDTI.LoadNifti"
- 'ExploreDTIInterface': I am pleased to announce that MRIToolkit now contains, distributes and develops many functions originally developed as part of ExploreDTI. They are here available as a consolidated library and are planned to also become command line tools. A big thank to Alexander Leemans and Ben Jeurissen for this!
- 'SphericalDeconvolution': Methods used for two novel deconvolution methods we developed, namely the Generalized Richardson Lucy and mFOD. Some of the functions here included have been taken from Hardi Tools of Erick Canales-Rodriguez.
-
- 'LesionEditor': a graphical user interface to visualise and edit segmentations of 3D MRI images, originally designed for delineation of multiple-sclerosis lesions on fluid attenuated inversion recovery (FLAIR) images. Requires MATLAB R2018a or newer. Documentation coming soon
- 'NiftiIO_basic': Basic Nifti input/output, including code originally written by Jimmy Shen
- 'DW_basic': Utilities to load / manipulate / save dMRI data
- 'OptimizationMethods': Classes and functions for numeric optimisation methods
- 'Relaxometry': Classes for T1/T2 quantification using inversion-recovery / spin-echo multi-echo data
- 'ThirdParty': Utilities from third parties used in other scripts. Includes EPG code from Brian Hargreaves
- 'ImageRegistrations': Image registration utils based on Elastix
- 'Demos': Small examples showcasing functionalities of MRIToolkit.
- 'Diffusion_basic': Class for (basic) dMRI quantification
- 'DW_IVIMDTDK_I': Diffusion MRI fit utilities - IVIM, DT, DKI - used in De Luca et al. 2017
- 'Dicom_utils': Tools for handling unconventional or buggy DICOMs
- 'DW_Laplacian_NNLS': Tools for spectral multi-compartment fit (NNLS/L2NNLS/PL2NNLS) used in De Luca et al. 2018
- 'MixedCodeUtils': 'Useful general purpose functions
- 'MRIfoundation': Classes for MRI sequences abstraction
- 'EPG_simulator': Classes for EPG simulations
- 'dfMRI', Diffusion fMRI utilities used in De Luca et al. 2019
- 'EPGFitT2Muscle.m': in Relaxometry, this function fits bi-exponential T2 (water/fat) with EPG - previously used in Arrigoni et al. 2018
- 'DESPOT12Fit.m': an example script to fit DESPOT1 and DESPOT2 as in Deoni et al. 2005
- 'MultiEchoTFit.m': an example script to fit mono-exponential T2 with multi-echo data. More advanced fit approaches (bi-exponential, EPG) will come soon
- 'InversionRecoveryT1Fit.m': an example script to fit T1 with inversion recovery data.
- 'CSD_Tractography_Script.m': an example script showcasing how to use the ExploreDTIInterface toolbox.
- More examples will come in the next days.
- The file naming convention is to always indicate Niftis as .nii, even when they are actually compressed in .nii.gz. The code takes care of that, but expects only .nii as arguments in function calls!
- Not everything has been checked yet - expect many bug fixes and new releases in the next time.
- MRIToolkit relies on a couple of third party dependencies:
- Elastix: 1) Either compile your own version or grab the executables for your platform here. 2) Copy the file "TemplateMRIToolkitDefineLocalVars.m" to your MATLAB default folder (user/MATLAB or Documents/MATLAB), rename the file as "MRIToolkitDefineLocalVars.m". 3) Edit the script, adjusting the variable MRIToolkit.Elastix.Location as needed.
- NODDI toolbox: if you would like to try the mFOD method, you will need to add the NODDI toolbox to the MATLAB path.
- ExploreDTI: While MRIToolkit is entirely self-sufficient (e.g. all needed ExploreDTI functions are bundled and adapted), the visualization of fiber tractograhy and other results will need ExploreDTI. Get it for free from Alexander Leemans.
- This code is a work in progress. It will be updated without notice to ensure bug-fixes and the inclusion of best available methods
- Most code is poorly commented and not general, but will be improved over releases.
- This software is distributed under the LGPLv3 license (https://opensource.org/licenses/lgpl-3.0.html).
- Magnetic Resonance Imaging (MRI)
- Image segmentation
- T1 quantification, Inversion Recovery
- T2 quantification, spin echo multi echo
- Extended Phase Graphs
- Diffusion MRI (dMRI) - Diffusion Tensor Imaging (DTI) - Diffusion Kurtosis Imaging (DKI)
- dMRI preprocessing - motion correction - eddy currents correction - EPI distortions correction
- Fiber tractography - Constrained Spherical Deconvolution (CSD) - Generalized Richardson Lucy (GRL) - mFOD
Alberto De Luca - First published in 2019