KomaMRI.jl: An open-source framework for general MRI simulations with GPU acceleration

Phase distribution graphs for fast, differentiable, and spatially encoded Bloch simulations of arbitrary MRI sequences

PURPOSE

An analytical approach to Bloch simulations for MRI sequences is introduced that enables time efficient calculations of signals free of Monte-Carlo noise, while providing full flexibility and differentiability in RF flip angles, RF phases, magnetic field gradients and time, as well as insights into image formation.

THEORY AND METHODS

We present an implementation of the extended phase graph (EPG) concept implemented in PyTorch, including an efficient state selection algorithm. This simulation is compared with an isochromat-based Bloch simulation with random isochromat distribution as well as with in vivo measurements using the Pulseq standard. Additionally, different sequences are tested and analyzed using this novel simulation approach.

RESULTS

Our simulation outperforms isochromat-based simulations in terms of computation time as well as signal quality, without exhibiting any kind of Monte-Carlo noise. The novel approach allows extracting useful information about the magnetization evolution not available otherwise. Our approach extends the common state-tensor-based EPG simulation approach for the contribution of dephased states including spatial encoding andT 2 ' $$ {T}_2^{\prime } $$ effects, and arbitrary timing. This allows calculation of echo shapes in addition to echo amplitudes only. Our implementation provides full differentiability in all input parameters allowing gradient descent optimization. Simulation of non-instantaneous pulses via hard-pulse approximation is left for future work, as the performance and accuracy characteristics are not yet analyzed.

CONCLUSIONS

Phase distribution graphs provide fast, differentiable, and spatially encoded Bloch simulations for most MRI sequences. It allows efficient simulation and optimization of arbitrary MRI sequences, which was previously only possible via high isochromat counts.

GPU-accelerated Bloch simulations and MR-STAT reconstructions using the Julia programming language

PURPOSE

MR-STAT is a relatively new multiparametric quantitative MRI technique in which quantitative paramater maps are obtained by solving a large-scale nonlinear optimization problem. Managing reconstruction times is one of the main challenges of MR-STAT. In this work we leverage GPU hardware to reduce MR-STAT reconstruction times. A highly optimized, GPU-compatible Bloch simulation toolbox is developed as part of this work that can be utilized for other quantitative MRI techniques as well.

METHODS

The Julia programming language was used to develop a flexible yet highly performant and GPU-compatible Bloch simulation toolbox called BlochSimulators.jl. The runtime performance of the toolbox is benchmarked against other Bloch simulation toolboxes. Furthermore, a (partially matrix-free) modification of a previously presented (matrix-free) MR-STAT reconstruction algorithm is proposed and implemented using the Julia language on GPU hardware. The proposed algorithm is combined with BlochSimulators.jl and the resulting MR-STAT reconstruction times on GPU hardware are compared to previously presented MR-STAT reconstruction times.

RESULTS

The BlochSimulators.jl package demonstrates superior runtime performance on both CPU and GPU hardware when compared to other existing Bloch simulation toolboxes. The GPU-accelerated partially matrix-free MR-STAT reconstruction algorithm, which relies on BlochSimulators.jl, allows for reconstructions of 68 seconds per two-dimensional (2D slice).

CONCLUSION

By combining the proposed Bloch simulation toolbox and the partially matrix-free reconstruction algorithm, 2D MR-STAT reconstructions can be performed in the order of one minute on a modern GPU card. The Bloch simulation toolbox can be utilized for other quantitative MRI techniques as well, for example for online dictionary generation for MR Fingerprinting.

PhoenixMR: A GPU-based MRI simulation framework with runtime-dynamic code execution

BACKGROUND

Simulations of physical processes and behavior can provide unique insights and understanding of real-world problems. Magnetic Resonance Imaging (MRI) is an imaging technique with several components of complexity. Several of these components have been characterized and simulated in the past. However, several computational challenges prevent simulations from being simultaneously fast, flexible, and accurate.

PURPOSE

The simulation of MRI experiments is underutilized by medical physicists and researchers using currently available simulators due to reasons including speed, accuracy, and extensibility constraints. This paper introduces an innovative MRI simulation engine and framework that aims to overcome these issues making available realistic and fast MRI simulation.

METHODS

Using the CUDA C/C++ programing language, an MRI simulation engine (PhoenixMR), incorporating a Turing-complete virtual machine (VM) to simulate abstract spatiotemporal complexities, was developed. This engine solves a set of time-discrete Bloch equations using the symmetric operator splitting technique. An extensible front-end framework package (written in Python) aids the use of PhoenixMR to simplify simulation development.

RESULTS

The PhoenixMR library and front-end codes have been developed and tested. A set of example simulations were performed to demonstrate the ease of use and flexibility of simulation components such as geometrical setup, pulse sequence design, phantom design, and so forth. Initial validation of PhoenixMR is performed by comparing its accuracy and performance against a widely used MRI simulator using identical simulation parameters. Validation results show PhoenixMR simulations are three orders of magnitude faster. There is also strong agreement between models.

CONCLUSIONS

A novel MRI simulation platform called PhoenixMR has been introduced. This research tool is designed to be usable by physicists and engineers interested in performing MRI simulations. Examples are shown demonstrating the accuracy, flexibility, and usability of PhoenixMR in several key areas of MRI simulation.

Open-source magnetic resonance imaging: Improving access, science, and education through global collaboration

Open-source practices and resources in magnetic resonance imaging (MRI) have increased substantially in recent years. This trend started with software and data being published open-source and, more recently, open-source hardware designs have become increasingly available. These developments towards a culture of sharing and establishing nonexclusive global collaborations have already improved the reproducibility and reusability of code and designs, while providing a more inclusive approach, especially for low-income settings. Community-driven standardization and documentation efforts are further strengthening and expanding these milestones. The future of open-source MRI is bright and we have just started to discover its full collaborative potential. In this review we will give an overview of open-source software and open-source hardware projects in human MRI research.

CMRsim-A python package for cardiovascular MR simulations incorporating complex motion and flow

PURPOSE

To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction.

METHODS

CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided.

RESULTS

The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain.

CONCLUSION

CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research.