BlochSolver: A GPU-optimized fast 3D MRI simulator for experimentally compatible pulse sequences

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.

Brain MR image simulation for deep learning based medical image analysis networks

BACKGROUND AND OBJECTIVE

As large sets of annotated MRI data are needed for training and validating deep learning based medical image analysis algorithms, the lack of sufficient annotated data is a critical problem. A possible solution is the generation of artificial data by means of physics-based simulations. Existing brain simulation data is limited in terms of anatomical models, tissue classes, fixed tissue characteristics, MR sequences and overall realism.

METHODS

We propose a realistic simulation framework by incorporating patient-specific phantoms and Bloch equations-based analytical solutions for fast and accurate MRI simulations. A large number of labels are derived from open-source high-resolution T1w MRI data using a fully automated brain classification tool. The brain labels are taken as ground truth (GT) on which MR images are simulated using our framework. Moreover, we demonstrate that the T1w MR images generated from our framework along with GT annotations can be utilized directly to train a 3D brain segmentation network. To evaluate our model further on larger set of real multi-source MRI data without GT, we compared our model to existing brain segmentation tools, FSL-FAST and SynthSeg.

RESULTS

Our framework generates 3D brain MRI for variable anatomy, sequence, contrast, SNR and resolution. The brain segmentation network for WM/GM/CSF trained only on T1w simulated data shows promising results on real MRI data from MRBrainS18 challenge dataset with a Dice scores of 0.818/0.832/0.828. On OASIS data, our model exhibits a close performance to FSL, both qualitatively and quantitatively with a Dice scores of 0.901/0.939/0.937.

CONCLUSIONS

Our proposed simulation framework is the initial step towards achieving truly physics-based MRI image generation, providing flexibility to generate large sets of variable MRI data for desired anatomy, sequence, contrast, SNR, and resolution. Furthermore, the generated images can effectively train 3D brain segmentation networks, mitigating the reliance on real 3D annotated data.

Using Dictionary Matching to Improve the Accuracy of MOLLI Myocardial T1 Analysis and Measurements of Heart Rate Variability

We analyzed modified Look-Locker inversion recovery (MOLLI) T1 measurements by applying a dictionary matching strategy and aimed to acquire T1 measurements more accurately than those acquired by the conventional three-parameter matching analysis. We particularly clarified the robustness of this method for measuring heart rate (HR) variability. A phantom experiment using a 3T MRI system was performed for various HRs. The ideal MOLLI signal corresponding to the scan parameter in the MRI experiment was simulated over a wide range of T1 values according to the dictionary. The unknown T1 values were determined by finding the simulated signals in the dictionary corresponding to the measured signals using pattern matching. The measured T1 values showed that the proposed analysis improved the accuracy of T1 measurements compared to those acquired by traditional analysis by up to 10%. In addition, the variability of measurements at several HRs was reduced by up to 100 ms.

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

PURPOSE

To develop an open-source, high-performance, easy-to-use, extensible, cross-platform, and general MRI simulation framework (Koma).

METHODS

Koma was developed using the Julia programming language. Like other MRI simulators, it solves the Bloch equations with CPU and GPU parallelization. The inputs are the scanner parameters, the phantom, and the pulse sequence that is Pulseq-compatible. The raw data is stored in the ISMRMRD format. For the reconstruction, MRIReco.jl is used. A graphical user interface utilizing web technologies was also designed. Two types of experiments were performed: one to compare the quality of the results and the execution speed, and the second to compare its usability. Finally, the use of Koma in quantitative imaging was demonstrated by simulating Magnetic Resonance Fingerprinting (MRF) acquisitions.

RESULTS

Koma was compared to two well-known open-source MRI simulators, JEMRIS and MRiLab. Highly accurate results (with mean absolute differences below 0.1% compared to JEMRIS) and better GPU performance than MRiLab were demonstrated. In an experiment with students, Koma was proved to be easy to use, eight times faster on personal computers than JEMRIS, and 65% of test subjects recommended it. The potential for designing acquisition and reconstruction techniques was also shown through the simulation of MRF acquisitions, with conclusions that agree with the literature.

CONCLUSIONS

Koma's speed and flexibility have the potential to make simulations more accessible for education and research. Koma is expected to be used for designing and testing novel pulse sequences before implementing them in the scanner with Pulseq files, and for creating synthetic data to train machine learning models.

A fast and practical computation method for magnetic resonance simulators

PURPOSE

This work aims to develop a fast and practical computation method for MR simulations. The computational cost of MR simulations is often high because magnetizations of many isochromats are updated using a small step size on the order of microseconds. There are two types of subsequences to be processed for the simulations: subsequences with and without RF pulses. While straightforward implementations spend most of their time calculating subsequences with RF pulses, there is a method which efficiently reuses the computation for repetitive RF pulses.

THEORY AND METHODS

A new method for efficiently processing subsequences with RF pulses is proposed. Rather than using an iterative update approach, the proposed method computes the combined transition which combines all transitions applied iteratively for each subsequence with RF pulses. The combined transition is used again when the same subsequence is used later. The combined transitions are cached and managed using a least recently used algorithm.

RESULTS

The proposed method was found to accelerate the simulation by ˜20 times when 3.9 million isochromats were simulated using spin-echo sequences. Even on a laptop computer, the proposed method was able to simulate these sequences in ˜3.5 min.

CONCLUSION

An efficient method for simulating pulse sequences is proposed. The proposed method computes and manages combined transitions, making MR simulation practical on a wide range of computers, including laptops.

Magnetic Resonance Simulation in Education: Quantitative Evaluation of an Actual Classroom Experience

Magnetic resonance is an imaging modality that implies a high complexity for radiographers. Despite some simulators having been developed for training purposes, we are not aware of any attempt to quantitatively measure their educational performance. The present study gives an answer to the question: Does an MRI simulator built on specific functional and non-functional requirements help radiographers learn MRI theoretical and practical concepts better than traditional educational method based on lectures? Our study was carried out in a single day by a total of 60 students of a main hospital in Madrid, Spain. The experiment followed a randomized pre-test post-test design with a control group that used a traditional educational method, and an experimental group that used our simulator. Knowledge level was assessed by means of an instrument with evidence of validity in its format and content, while its reliability was analyzed after the experiment. Statistical differences between both groups were measured. Significant statistical differences were found in favor of the participants who used the simulator for both the post-test score and the gain (difference between post-test and pre-test scores). The effect size turned out to be significant as well. In this work we evaluated a magnetic resonance simulation paradigm as a tool to help in the training of radiographers. The study shows that a simulator built on specific design requirements is a valuable complement to traditional education procedures, backed up with significant quantitative results.

Bloch Simulation of a Three-point Dixon Experiment Using a Four-dimensional Numerical Phantom

A 4D numerical phantom, which is defined in the 3D spatial axes and the resonance frequency axis, is indispensable for Bloch simulations of biological tissues with complex distribution of materials. In this study, a 4D numerical phantom was created using MR image datasets of a biological sample containing water and fat, and the Bloch simulations were performed using the 4D numerical phantom. As a result, 3D images of the sample containing water and fat were successfully reproduced, which demonstrated the usefulness of the concept of the 4D numerical phantom.

Magnetic resonance imaging simulation with spin-lock preparations to detect tiny oscillatory magnetic fields

Spin-lock preparation was studied to detect tiny oscillatory magnetic fields such as a neural magnetic field without the blood oxygen level-dependent effect. This approach is a direct measurement and independent of static magnetic field strength. Accordingly, it is anticipated as a feasible functional magnetic resonance imaging (fMRI) in low and ultra-low-field MRI. Several reports have been published on spin-lock preparation but reports on imaging simulation are rare. Research in this area can assist in investigating magnetic resonance signal changes and, accordingly, can help to develop new spin-lock methods. In this study, we propose an imaging simulation method with an analytical solution using the Bloch equation. To demonstrate the feasibility of our proposed method, we compared simulated images with experimental results in which the number of sub-voxels and the amplitude and phase of the target oscillatory magnetic fields varied. In addition, we also applied graphics processing unit parallel computing and investigated the feasibility of avoiding an impracticable calculation time by doing so.

Development of a method for the Bloch image simulation of biological tissues

PURPOSE

The purpose of this study is to develop a method for the Bloch image simulation of biological tissues including various chemical components and T₂* distribution.

METHODS

The nuclear spins in the object material were modeled as a spectral intensity function Sr→ω defined by superposition of Lorentz functions with various central precession frequencies and the half width of 1/(πT₂'), where 1/T₂' is a relaxation rate attributable to microscopic field inhomogeneity in a voxel. Four-dimensional numerical phantoms were created to simulate Sr→ω and used for MRI simulations of the phantoms containing water and fat protons. Single slice multiple (16) gradient-echo sequences (ΔTE = 2.2 and 1.384 ms) were used for experiments at 1.5 T and 3 T and MRI simulations to evaluate the validity of the approach.

RESULTS

Experimentally measured image intensities of the multiple gradient-echo imaging sequences were well reproduced by the MRI simulations. The correlation coefficients between the experimentally measured image intensities and those numerically simulated were 0.9895 to 0.9992 for the 4-component phantom at 1.5 T and 0.9580 to 0.9996 for the 7-component phantom at 3 T.

CONCLUSION

T₂* and chemical shift effects were successfully implemented in the MRI simulator (BlochSolver). Because this approach can be applied to other MRI simulators, the method developed in this study is useful for MRI simulation of biological tissues containing water and fat protons.

Integration of an intelligent tutoring system in a magnetic resonance simulator for education: Technical feasibility and user experience

BACKGROUND AND OBJECTIVE

In this paper we propose to include an intelligent tutoring system (ITS) within a magnetic resonance (MR) simulator that has been developed in house. With this, we intend to measure the impact, in terms of user experience, of including an ITS in our simulator.

METHODS

We thoroughly describe the integration procedure and we have tested the benefits of this integration by means of two actual educational experiences, with one of them using the simulator as a standalone tool, and the other with the joint use of simulator+ITS. The experiences have consisted of two online courses with a number of students around 180 in both of them, where measurements of usability, perceived utility and likelihood to recommend were collected.

RESULTS

We have observed that the three measurements improved noticeably in the second course with respect to the first one; specifically, overall usability improved by 22.3%, perceived utility by an average of 55.1% and likelihood to recommend by 13.7%. In addition, quantitative measurements are complemented with comments in free text format directly provided by the students. Results show evidence on the benefits of integrating an ITS in terms of quantitative user experience, as well as qualitative comparative comments directly by students of both courses.

CONCLUSIONS

This is the first time that an ITS is used within the scope of MR simulation for training purposes. Benefits of integrating an ITS within an MR simulator have been evaluated in terms of user experience, with satisfactory comparative results.

An Accurate Dictionary Creation Method for MR Fingerprinting Using a Fast Bloch Simulator

This study proposes an accurate method for creating a dictionary for magnetic resonance fingerprinting (MRF) using a fast Bloch image simulator. An MRF sequence based on a fast imaging with steady precession sequence and a numerical phantom were used for dictionary generation. Cartesian and spiral readout gradients were used for the Bloch image simulation. The validity and usefulness of the method for accurate dictionary creation were demonstrated by MRF parameter maps obtained by pattern matching with the dictionaries generated by the proposed method.

Virtual clinical trials in medical imaging: a review

The accelerating complexity and variety of medical imaging devices and methods have outpaced the ability to evaluate and optimize their design and clinical use. This is a significant and increasing challenge for both scientific investigations and clinical applications. Evaluations would ideally be done using clinical imaging trials. These experiments, however, are often not practical due to ethical limitations, expense, time requirements, or lack of ground truth. Virtual clinical trials (VCTs) (also known as in silico imaging trials or virtual imaging trials) offer an alternative means to efficiently evaluate medical imaging technologies virtually. They do so by simulating the patients, imaging systems, and interpreters. The field of VCTs has been constantly advanced over the past decades in multiple areas. We summarize the major developments and current status of the field of VCTs in medical imaging. We review the core components of a VCT: computational phantoms, simulators of different imaging modalities, and interpretation models. We also highlight some of the applications of VCTs across various imaging modalities.

High-resolution in vivo MR-STAT using a matrix-free and parallelized reconstruction algorithm

MR-STAT is a recently proposed framework that allows the reconstruction of multiple quantitative parameter maps from a single short scan by performing spatial localisation and parameter estimation on the time-domain data simultaneously, without relying on the fast Fourier transform (FFT). To do this at high resolution, specialized algorithms are required to solve the underlying large-scale nonlinear optimisation problem. We propose a matrix-free and parallelized inexact Gauss-Newton based reconstruction algorithm for this purpose. The proposed algorithm is implemented on a high-performance computing cluster and is demonstrated to be able to generate high-resolution (1 mm  ×  1 mm in-plane resolution) quantitative parameter maps in simulation, phantom, and in vivo brain experiments. Reconstructed T1 and T2 values for the gel phantoms are in agreement with results from gold standard measurements and, for the in vivo experiments, the quantitative values show good agreement with literature values. In all experiments, short pulse sequences with robust Cartesian sampling are used, for which MR fingerprinting reconstructions are shown to fail.

A Web-Based Educational Magnetic Resonance Simulator: Design, Implementation and Testing

A new web-based education-oriented magnetic resonance (MR) simulator is presented. We have identified the main requirements that this simulator should comply with, so that trainees can face useful practical tasks such as setting the exact slice position and its properties, selecting the correct protocol or fitting the parameters to acquire an image. The tool follows the client-server model. The client contains the interface that mimics the console of a real machine and several of its features. The server stores anatomical models and executes the bulk of the simulation. This cross-platform simulator has been used in two real educational scenarios. The acceptance of the tool has been measured using two criteria, namely, the System Usability Scale and the Likelihood to Recommend, both with satisfactory results. Therefore, we conclude that given the potential of the tool, it may play a relevant role for the training of MRI operators and other involved personnel.

coreMRI: A high-performance, publicly available MR simulation platform on the cloud

Introduction

A Cloud-ORiented Engine for advanced MRI simulations (coreMRI) is presented in this study. The aim was to develop the first advanced MR simulation platform delivered as a web service through an on-demand, scalable cloud-based and GPU-based infrastructure. We hypothesized that such an online MR simulation platform could be utilized as a virtual MRI scanner but also as a cloud-based, high-performance engine for advanced MR simulations in simulation-based quantitative MR (qMR) methods.

Methods and results

The simulation framework of coreMRI was based on the solution of the Bloch equations and utilized a ground-up-approach design based on the principles already published in the literature. The development of a front-end environment allowed the connection of the end-users to the GPU-equipped instances on the cloud. The coreMRI simulation platform was based on a modular design where individual modules (such as the Gadgetron reconstruction framework and a newly developed Pulse Sequence Designer) could be inserted in the main simulation framework. Different types and sources of pulse sequences and anatomical models were utilized in this study revealing the flexibility that the coreMRI simulation platform offers to the users. The performance and scalability of coreMRI were also examined on multi-GPU configurations on the cloud, showing that a multi-GPU computer on the cloud equipped with a newer generation of GPU cards could significantly mitigate the prolonged execution times that accompany more realistic MRI and qMR simulations.

Conclusions

coreMRI is available to the entire MR community, whereas its high performance and scalability allow its users to configure advanced MRI experiments without the constraints imposed by experimentation in a true MRI scanner (such as time constraint and limited availability of MR scanners), without upfront investment for purchasing advanced computer systems and without any user expertise on computer programming or MR physics. coreMRI is available to the users through the webpage https://www.coreMRI.org.

A Fast GPU-optimized 3D MRI Simulator for Arbitrary k-space Sampling

Purpose:

To develop a fast 3D MRI simulator for arbitrary k-space sampling using a graphical processing unit (GPU) and demonstrate its performance by comparing simulation and experimental results in a real MRI system.

Materials and Methods:

A fast 3D MRI simulator using a GeForce GTX 1080 GPU (NVIDIA Corporation, Santa Clara, CA, USA) was developed using C++ and the CUDA 8.0 platform (NVIDIA Corporation). The unique advantage of this simulator was that it could use the same pulse sequence as used in the experiment. The performance of the MRI simulator was measured using two GTX 1080 GPUs and 3D Cones sequences. The MRI simulation results for 3D non-Cartesian sampling trajectories like 3D Cones sequences using a numerical 3D phantom were compared with the experimental results obtained with a real MRI system and a real 3D phantom.

Results:

The performance of the MRI simulator was about 3800–4900 gigaflops for 128- to 4-shot 3D Cones sequences with 256³ voxels, which was about 60% of the performance of the previous MRI simulator optimized for Cartesian sampling calculated for a Cartesian sampling gradient-echo sequence with 256³ voxels. The effects of the static magnetic field inhomogeneity, radio-frequency field inhomogeneity, gradient field nonlinearity, and fast repetition times on the MR images were reproduced in the simulated images as observed in the experimental images.

Conclusion:

The 3D MRI simulator developed for arbitrary k-space sampling optimized using GPUs is a powerful tool for the development and evaluation of advanced imaging sequences including both Cartesian and non-Cartesian k-space sampling.