The MRI Coil-reconstruct Toolbox
Project description
The MRI Coil-reconstruct Toolbox, MCT, is a small toolbox for combining the channels of a multi-channel MRI acquisition. Where possible, this toolbox uses GPU accelerated routines to speed-up the processing. For example, the weights of the STARC (STAbility-weighted Rf-coil Combination) reconstruction model are fitted using the GPU or using multi-threaded CPU (utilizing OpenCL). At the moment MCT supports rSoS (root Sum Of Squares), rCovSoS (same as SoS but than with additional usage of a noise covariance matrix) and STARC.
Summary
features multiple coil combine / reconstruction methods
command line and python interface
GPU/multicore-CPU accelerated STARC
Free Open Source Software: LGPL v3 license
Python and OpenCL based
Full documentation: https://mct.readthedocs.io/
Project home: https://github.com/robbert-harms/MCT
Tags: MRI, coil-reconstruct, image reconstruction, opencl, python
Data reconstruction
This software contains various reconstruction methods that can be used to combine your channels into one (or more) volumes. Not all reconstruction methods may be applicable to your data, for example the STARC [1] method only works when working with fMRI data.
Console
To reconstruct your data using the command line, after installation you can use:
$ mct-reconstruct <method> {0..15}.nii
Where method at the moment is one of “rSoS”, “rCovSoS” or “STARC”. Some methods require more information to combine the channels, please see the full documentation for this.
If you only want to use certain volumes of your data, use the “–volumes” or “-v” switch on the command line:
$ mct-reconstruct <method> {0..15}.nii -v odd
To use (for example) only the odd volumes. Available options are “odd”, “even” or a list of indices, such as “0 2 4 5” (space separated).
Python
It is also possible to reconstruct your data using the Python API, for example:
from mct.reconstruction_methods import rSoS, rCovSoS, STARC
input_path = '/data/'
output_path = '/data/output/'
nmr_channels = 16
input_filenames = [input_path + str(ind)
for ind in range(nmr_channels)]
method = rSoS(input_filenames)
method.reconstruct(output_path, volumes='odd')
This would reconstruct your data using rSoS using only the odd volumes.
References:
Simple approach to improve time series fMRI stability: STAbility-weighted Rf-coil Combination (STARC), L. Huber et al. ISMRM 2017 abstract #0586.
Kashyap S, Fritz FJ, Harms RL, Huber L, Ivanov D, Roebroeck A, Poser BA, Uludag, K. Effect of optimised coil-combinations on high-resolution laminar fMRI at 9.4T, Proceedings of ISMRM 2018 Abstract #5442, (https://doi.org/10.7490/f1000research.1117449.1)
Quick installation guide
The basic requirements for MCT are:
Python 3.x (recommended) or Python 2.7
OpenCL 1.2 (or higher) support in GPU driver or CPU runtime
Linux
For Ubuntu >= 16 you can use:
sudo add-apt-repository ppa:robbert-harms/cbclab
sudo apt-get update
sudo apt-get install python3-mct
For Debian users and Ubuntu < 16 users, install MDT with:
sudo apt-get install python3 python3-pip python3-pyopencl python3-numpy python3-nibabel python3-pyqt5 python3-matplotlib python3-yaml python3-argcomplete libpng-dev libfreetype6-dev libxft-dev
sudo pip3 install mct
Note that python3-nibabel may need NeuroDebian to be available on your machine. An alternative is to use pip3 install nibabel instead.
Windows
The installation on Windows is a little bit more complex and the following is only a quick reference guide. To save duplication of information and since this package depends on MDT and MOT, the complete install instructions can be copied from the MDT documentation. After following that guide, installation of MCT is simply done using pip install mct. The quick overview is:
Install Anaconda Python 3.5
Install MOT using the guide at https://mot.readthedocs.io
Open an Anaconda shell and type: pip install mct
Mac
Install Anaconda Python 3.5
Open a terminal and type: pip install mct
Please note that Mac support is experimental due to the unstable nature of the OpenCL drivers in Mac, that is, users running MDT with the GPU as selected device may experience crashes. Running MDT in the CPU seems to work though.
For more information and full installation instructions please see the documentation of the MDT package https://maastrichtdiffusiontoolbox.readthedocs.org
Roadmap
Add a few more reconstruction methods like:
Roemer
GRAPPA
SENSE
Improve the data handling and memory usage.
Interruptable processing (saving intermediate data to HD and start from there)
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.