Development of a functional magnetic resonance imaging simulator for modeling realistic rigid-body motion artifacts

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.

Outliers in diffusion-weighted MRI: Exploring detection models and mitigation strategies

Diffusion-weighted MRI (dMRI) is a medical imaging method that can be used to investigate the brain microstructure and structural connections between different brain regions. The method, however, requires relatively complex data processing frameworks and analysis pipelines. Many of these approaches are vulnerable to signal dropout artefacts that can originate from subjects moving their head during the scan. To combat these artefacts and eliminate such outliers, researchers have proposed two approaches: to replace outliers or to downweight outliers during modelling and analysis. With the rising interest in dMRI for clinical research, these types of corrections are increasingly important. Therefore, we set out to investigate the differences between outlier replacement and weighting approaches to help the dMRI community to select the best tool for their data processing pipelines. We evaluated dMRI motion correction registration and single tensor model fit pipelines using Gaussian Process and Spherical Harmonic based replacement approaches and outlier downweighting using highly realistic whole-brain simulations. As a proof of concept, we applied these approaches to dMRI infant data sets that contained varying numbers of dropout artefacts. Based on our results, we concluded that the Gaussian Process based outlier replacement provided similar tensor fit results to Gaussian Process based outlier detection and downweighting. Therefore, if only the least-squares estimate of the single tensor model is of interest, our recommendation is to use outlier replacement. However, outlier downweighting can potentially provide a more accurate estimate of the model precision which could be relevant for applications such as probabilistic tractoraphy.

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.

Effect of subject-specific head morphometry on specific absorption rate estimates in parallel-transmit MRI at 7 T

PURPOSE

To assess the accuracy of morphing an established reference electromagnetic head model to a subject-specific morphometry for the estimation of specific absorption rate (SAR) in 7T parallel-transmit (pTx) MRI.

METHODS

Synthetic T1 -weighted MR images were created from three high-resolution open-source electromagnetic head voxel models. The accuracy of morphing a "reference" (multimodal image-based detailed anatomical [MIDA]) electromagnetic model into a different subject's native space (Duke and Ella) was compared. Both linear and nonlinear registration methods were evaluated. Maximum 10-g averaged SAR was estimated for circularly polarized mode and for 5000 random RF shim sets in an eight-channel transmit head coil, and comparison made between the morphed MIDA electromagnetic models and the native Duke and Ella electromagnetic models, respectively.

RESULTS

The averaged error in maximum 10-g averaged SAR estimation across pTx MRI shim sets between the MIDA and the Duke target model was reduced from 17.5% with only rigid-body registration, to 11.8% when affine linear registration was used, and further reduced to 10.7% when nonlinear registration was used. The corresponding figures for the Ella model were 16.7%, 11.2%, and 10.1%.

CONCLUSION

We found that morphometry accounts for up to half of the subject-specific differences in pTx SAR. Both linear and nonlinear morphing of an electromagnetic model into a target subject improved SAR agreement by better matching head size, morphometry, and position. However, differences remained, likely arising from details in tissue composition estimation. Thus, the uncertainty of the head morphometry and tissue composition may need to be considered separately to achieve personalized SAR estimation.

Deep Learning for Image Enhancement and Correction in Magnetic Resonance Imaging-State-of-the-Art and Challenges

Magnetic resonance imaging (MRI) provides excellent soft-tissue contrast for clinical diagnoses and research which underpin many recent breakthroughs in medicine and biology. The post-processing of reconstructed MR images is often automated for incorporation into MRI scanners by the manufacturers and increasingly plays a critical role in the final image quality for clinical reporting and interpretation. For image enhancement and correction, the post-processing steps include noise reduction, image artefact correction, and image resolution improvements. With the recent success of deep learning in many research fields, there is great potential to apply deep learning for MR image enhancement, and recent publications have demonstrated promising results. Motivated by the rapidly growing literature in this area, in this review paper, we provide a comprehensive overview of deep learning-based methods for post-processing MR images to enhance image quality and correct image artefacts. We aim to provide researchers in MRI or other research fields, including computer vision and image processing, a literature survey of deep learning approaches for MR image enhancement. We discuss the current limitations of the application of artificial intelligence in MRI and highlight possible directions for future developments. In the era of deep learning, we highlight the importance of a critical appraisal of the explanatory information provided and the generalizability of deep learning algorithms in medical imaging.

Estimation of undistorted images in brain echo-planar images with distortions using the conjugate gradient method with anatomical regularization

PURPOSE

Although echo-planar imaging (EPI) is widely used for diffusion magnetic resonance (MR) imaging, EPI images suffer from susceptibility-induced geometric distortions. We herein propose a new estimation method for undistorted EPI images using anatomical T₁ -weighted images (T₁ WIs) based on the physics of MR imaging.

METHODS

Our proposed method estimates the undistorted EPI image in the image domain while estimating the magnetic field inhomogeneity map using the conjugate gradient method with anatomical regularization. Our method synthesizes the distorted image to match the measured EPI image containing geometric distortions by alternately updating the undistorted EPI image and the magnetic field inhomogeneity map. We evaluated our proposed method and compared it with a nonrigid registration-based distortion correction method using simulated data and using real data. In the evaluation of the estimation of the magnetic field inhomogeneity map, we used the normalized root-mean-squared error (NRMSE) between the estimated results and the ground truth. In the evaluation of the estimation of undistorted images, we used mutual information (MI) between the undistorted EPI image and the anatomical T₁ WI.

RESULTS

Using the simulated data, the means and standard deviations of the NRMSE values in the nonrigid registration-based method and proposed method were 1.29 ± 0.63 and 0.64 ± 0.30, respectively. The MI values in the proposed method were larger than those in the nonrigid registration-based method in all evaluated slices. For the real data, the proposed method improved the distortion, and the MI values in the proposed method were larger than those in the nonrigid registration-based method. In the estimation of the magnetic field inhomogeneity map, the NRMSE values in our method were smaller than those in the nonrigid registration-based method.

CONCLUSIONS

We demonstrated that our proposed method can estimate the regions with compressed distortions that are not well represented by the nonrigid registration-based methods. The results suggest that the proposed method could be useful in analyses combining EPI images with T₁ WIs.

A Fetal Brain magnetic resonance Acquisition Numerical phantom (FaBiAN)

Accurate characterization of in utero human brain maturation is critical as it involves complex and interconnected structural and functional processes that may influence health later in life. Magnetic resonance imaging is a powerful tool to investigate equivocal neurological patterns during fetal development. However, the number of acquisitions of satisfactory quality available in this cohort of sensitive subjects remains scarce, thus hindering the validation of advanced image processing techniques. Numerical phantoms can mitigate these limitations by providing a controlled environment with a known ground truth. In this work, we present FaBiAN, an open-source Fetal Brain magnetic resonance Acquisition Numerical phantom that simulates clinical T2-weighted fast spin echo sequences of the fetal brain. This unique tool is based on a general, flexible and realistic setup that includes stochastic fetal movements, thus providing images of the fetal brain throughout maturation comparable to clinical acquisitions. We demonstrate its value to evaluate the robustness and optimize the accuracy of an algorithm for super-resolution fetal brain magnetic resonance imaging from simulated motion-corrupted 2D low-resolution series compared to a synthetic high-resolution reference volume. We also show that the images generated can complement clinical datasets to support data-intensive deep learning methods for fetal brain tissue segmentation.

Physical and digital phantoms for validating tractography and assessing artifacts

Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: "What are dMRI phantoms and what are they good for?", "What would the ideal phantom for validation fiber tractography look like?" and "What phantoms, phantom datasets and tools used for their creation are available to the research community?". We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take.

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.

Improving B 1 + parametric estimation in the brain from multispin-echo sequences using a fusion bootstrap moves solver

PURPOSE

To simultaneously estimate the B 1 + field (along with the T₂ ) in the brain with multispin-echo (MSE) sequences and dictionary matching.

METHODS

T₂ mapping provides clinically relevant information such as in the assessment of brain degenerative diseases. It is commonly obtained with MSE sequences, and accuracy can be further improved by matching the MSE signal to a precomputed dictionary of echo-modulation curves. For additional T₁ quantification, transmit B 1 + field knowledge is also required. Preliminary work has shown that although simultaneous brain B 1 + estimation along with T₂ is possible, it presents a bimodal distribution with the main peak coinciding with the true value. By taking advantage of this, the B 1 + maps are expected to be spatially smooth by applying an iterative method that takes into account each pixel neighborhood known as the fusion bootstrap moves solver (FBMS). The effect of the FBMS on B 1 + accuracy and piecewise smoothness is investigated and different spatial regularization levels are compared. Total variation regularization was used for both B 1 + and T₂ simultaneous estimation because of its simplicity as an initial proof-of-concept; future work could explore non edge-preserving regularization independently for B 1 + .

RESULTS

Improvements in B 1 + accuracy (up to 45.37% and 16.81% B 1 + error decrease) and recovery of spatially homogeneous maps are shown in simulations and in vivo 3.0T brain data, respectively.

CONCLUSION

Accurate B 1 + estimated values can be obtained from widely available MSE sequences while jointly estimating T₂ maps with the use of echo-modulation curve matching and FBMS at no further cost.

Attenuation of motion artifacts in fMRI using discrete reconstruction of irregular fMRI trajectories (DRIFT)

PURPOSE

Numerous studies report motion as the most detrimental source of noise and artifacts in fMRI. Current motion correction methods fail to completely address the motion problem. Retrospective techniques such as spatial realignment can correct for between-volume misalignment but fail to address within volume contamination and spin-history artifacts. Prospective motion correction can prevent spin-history artifacts but currently cannot update the gradients fast enough to remove k-space filling artifacts, calling for a hybrid approach to fully address these problems.

THEORY AND METHODS

Motion can be mathematically formulated into the MR signal equation to describe the motion artifacts at their origin in k-space. From these equations, it is demonstrated that different motions have different effects on the signal. A novel motion correction algorithm is designed from these equations to remove motion-induced artifacts directly in k-space, discrete reconstruction of irregular fMRI trajectory (DRIFT). This method is evaluated rigorously using fMRI simulations and data from a rotating phantom inside the scanner.

RESULTS

The results indicate that although some motion types have negligible effects on the MR signal, others produce catastrophic and lasting artifacts even after motion cessation. In simulation, DRIFT is able to remove motion artifacts in the absence of spin history. In a phantom scan, DRIFT significantly attenuates the motion artifacts in the fMRI data.

CONCLUSION

Neither prospective nor retrospective motion correction methods could completely remove the motion artifacts from the fMRI data. However, DRIFT, as a retrospective technique, when combined with prospective motion correction, can eliminate a significant portion of motion artifacts.

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.

R-fMRI reconstruction from k-t undersampled data using a subject-invariant dictionary model and VB-EM with nested minorization

Higher spatial resolution in resting-state functional magnetic resonance imaging (R-fMRI) can give reliable information about the functional networks in the cerebral cortex. Typical methods can achieve higher spatial or temporal resolution by speeding up scans using either (i) complex pulse-sequence designs or (ii) k-space undersampling coupled with priors on the signal. We propose to undersample the R-fMRI acquisition in k-space and time to speedup scans in order to improve spatial resolution. We propose a novel model-based R-fMRI reconstruction framework using a robust, subject-invariant, spatially regularized dictionary prior on the signal. Furthermore, we propose a novel inference framework based on variational Bayesian expectation maximization with nested minorization (VB-EM-NM). Our inference framework allows us to provide an estimate of uncertainty of the reconstruction, unlike typical reconstruction methods. Empirical evaluation of (i) simulated R-fMRI reconstruction and (ii) functional-network estimates from brain R-fMRI reconstructions demonstrate that our framework improves over the state of the art, and, additionally, enables significantly higher spatial resolution.

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.

Evaluation of Six Phase Encoding Based Susceptibility Distortion Correction Methods for Diffusion MRI

Purpose: Susceptibility distortions impact diffusion MRI data analysis and is typically corrected during preprocessing. Correction strategies involve three classes of methods: registration to a structural image, the use of a fieldmap, or the use of images acquired with opposing phase encoding directions. It has been demonstrated that phase encoding based methods outperform the other two classes, but unfortunately, the choice of which phase encoding based method to use is still an open question due to the absence of any systematic comparisons.

Methods: In this paper we quantitatively evaluated six popular phase encoding based methods for correcting susceptibility distortions in diffusion MRI data. We employed a framework that allows for the simulation of realistic diffusion MRI data with susceptibility distortions. We evaluated the ability for methods to correct distortions by comparing the corrected data with the ground truth. Four diffusion tensor metrics (FA, MD, eigenvalues and eigenvectors) were calculated from the corrected data and compared with the ground truth. We also validated two popular indirect metrics using both simulated data and real data. The two indirect metrics are the difference between the corrected LR and AP data, and the FA standard deviation over the corrected LR, RL, AP, and PA data.

Results: We found that DR-BUDDI and TOPUP offered the most accurate and robust correction compared to the other four methods using both direct and indirect evaluation metrics. EPIC and HySCO performed well in correcting b₀ images but produced poor corrections for diffusion weighted volumes, and also they produced large errors for the four diffusion tensor metrics. We also demonstrate that the indirect metric (the difference between corrected LR and AP data) gives a different ordering of correction quality than the direct metric.

Conclusion: We suggest researchers to use DR-BUDDI or TOPUP for susceptibility distortion correction. The two indirect metrics (the difference between corrected LR and AP data, and the FA standard deviation) should be interpreted together as a measure of distortion correction quality. The performance ranking of the various tools inferred from direct and indirect metrics differs slightly. However, across all tools, the results of direct and indirect metrics are highly correlated indicating that the analysis of indirect metrics may provide a good proxy of the performance of a correction tool if assessment using direct metrics is not feasible.

Comparison of quality control methods for automated diffusion tensor imaging analysis pipelines

The processing of brain diffusion tensor imaging (DTI) data for large cohort studies requires fully automatic pipelines to perform quality control (QC) and artifact/outlier removal procedures on the raw DTI data prior to calculation of diffusion parameters. In this study, three automatic DTI processing pipelines, each complying with the general ENIGMA framework, were designed by uniquely combining multiple image processing software tools. Different QC procedures based on the RESTORE algorithm, the DTIPrep protocol, and a combination of both methods were compared using simulated ground truth and artifact containing DTI datasets modeling eddy current induced distortions, various levels of motion artifacts, and thermal noise. Variability was also examined in 20 DTI datasets acquired in subjects with vascular cognitive impairment (VCI) from the multi-site Ontario Neurodegenerative Disease Research Initiative (ONDRI). The mean fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were calculated in global brain grey matter (GM) and white matter (WM) regions. For the simulated DTI datasets, the measure used to evaluate the performance of the pipelines was the normalized difference between the mean DTI metrics measured in GM and WM regions and the corresponding ground truth DTI value. The performance of the proposed pipelines was very similar, particularly in FA measurements. However, the pipeline based on the RESTORE algorithm was the most accurate when analyzing the artifact containing DTI datasets. The pipeline that combined the DTIPrep protocol and the RESTORE algorithm produced the lowest standard deviation in FA measurements in normal appearing WM across subjects. We concluded that this pipeline was the most robust and is preferred for automated analysis of multisite brain DTI data.

Agito ergo sum: Correlates of spatio-temporal motion characteristics during fMRI

The impact of in-scanner motion on functional magnetic resonance imaging (fMRI) data has a notorious reputation in the neuroimaging community. State-of-the-art guidelines advise to scrub out excessively corrupted frames as assessed by a composite framewise displacement (FD) score, to regress out models of nuisance variables, and to include average FD as a covariate in group-level analyses. Here, we studied individual motion time courses at time points typically retained in fMRI analyses. We observed that even in this set of putatively clean time points, motion exhibited a very clear spatio-temporal structure, so that we could distinguish subjects into separate groups of movers with varying characteristics. Then, we showed that this spatio-temporal motion cartography tightly relates to a broad array of anthropometric and cognitive factors. Convergent results were obtained from two different analytical perspectives: univariate assessment of behavioural differences across mover subgroups unraveled defining markers, while subsequent multivariate analysis broadened the range of involved factors and clarified that multiple motion/behaviour modes of covariance overlap in the data. Our results demonstrate that even the smaller episodes of motion typically retained in fMRI analyses carry structured, behaviourally relevant information. They call for further examinations of possible biases in current regression-based motion correction strategies.

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.

The Benefit of Slice Timing Correction in Common fMRI Preprocessing Pipelines

Due to the nature of fMRI acquisition protocols, slices cannot be acquired simultaneously, and as a result, are temporally misaligned from each other. To correct from this misalignment, preprocessing pipelines often incorporate slice timing correction (STC). However, evaluating the benefits of STC is challenging because it (1) is dependent on slice acquisition parameters, (2) interacts with head movement in a non-linear fashion, and (3) significantly changes with other preprocessing steps, fMRI experimental design, and fMRI acquisition parameters. Presently, the interaction of STC with various scan conditions has not been extensively examined. Here, we examine the effect of STC when it is applied with various other preprocessing steps such as motion correction (MC), motion parameter residualization (MPR), and spatial smoothing. Using 180 simulated and 30 real fMRI data, we quantitatively demonstrate that the optimal order in which STC should be applied depends on interleave parameters and motion level. We also demonstrate the benefit STC on sub-second-TR scans and for functional connectivity analysis. We conclude that STC is a critical part of the preprocessing pipeline that can be extremely beneficial for fMRI processing. However, its effectiveness interacts with other preprocessing steps and with other scan parameters and conditions which may obscure its significant importance in the fMRI processing pipeline.

A supervised learning approach for diffusion MRI quality control with minimal training data

Quality control (QC) is a fundamental component of any study. Diffusion MRI has unique challenges that make manual QC particularly difficult, including a greater number of artefacts than other MR modalities and a greater volume of data. The gold standard is manual inspection of the data, but this process is time-consuming and subjective. Recently supervised learning approaches based on convolutional neural networks have been shown to be competitive with manual inspection. A drawback of these approaches is they still require a manually labelled dataset for training, which is itself time-consuming to produce and still introduces an element of subjectivity. In this work we demonstrate the need for manual labelling can be greatly reduced by training on simulated data, and using a small amount of labelled data for a final calibration step. We demonstrate its potential for the detection of severe movement artefacts, and compare performance to a classifier trained on manually-labelled real data.

Flow MRI simulation in complex 3D geometries: Application to the cerebral venous network

PURPOSE

Develop and evaluate a complete tool to include 3D fluid flows in MRI simulation, leveraging from existing software. Simulation of MR spin flow motion is of high interest in the study of flow artifacts and angiography. However, at present, only a few simulators include this option and most are restricted to static tissue imaging.

THEORY AND METHODS

An extension of JEMRIS, one of the most advanced high performance open-source simulation platforms to date, was developed. The implementation of a Lagrangian description of the flow allows simulating any MR experiment, including both static tissues and complex flow data from computational fluid dynamics. Simulations of simple flow models are compared with real experiments on a physical flow phantom. A realistic simulation of 3D flow MRI on the cerebral venous network is also carried out.

RESULTS

Simulations and real experiments are in good agreement. The generality of the framework is illustrated in 2D and 3D with some common flow artifacts (misregistration and inflow enhancement) and with the three main angiographic techniques: phase contrast velocimetry (PC), time-of-flight, and contrast-enhanced imaging MRA.

CONCLUSION

The framework provides a versatile and reusable tool for the simulation of any MRI experiment including physiological fluids and arbitrarily complex flow motion.

Susceptibility-induced distortion that varies due to motion: Correction in diffusion MR without acquiring additional data

Because of their low bandwidth in the phase-encode (PE) direction, the susceptibility-induced off-resonance field causes distortions in echo planar imaging (EPI) images. It is therefore crucial to correct for susceptibility-induced distortions when performing diffusion studies using EPI. The susceptibility-induced field is caused by the object (head) disrupting the field and it is typically assumed that it remains constant within a framework defined by the object, (i.e. it follows the object as it moves in the scanner).

However, this is only approximately true. When a non-spherical object rotates around an axis other than that parallel with the magnetic flux (the z-axis) it changes the way it disrupts the field, leading to different distortions. Hence, if using a single field to correct for distortions there will be residual distortions in the volumes where the object orientation is substantially different to that when the field was measured.

In this paper we present a post-processing method for estimating the field as it changes with motion during the course of an experiment. It only requires a single measured field and knowledge of the orientation of the subject when that field was acquired. The volume-to-volume changes of the field as a consequence of subject movement are estimated directly from the diffusion data without the need for any additional or special acquisitions. It uses a generative model that predicts how each volume would look predicated on field change and inverts that model to yield an estimate of the field changes. It has been validated on both simulations and experimental data. The results show that we are able to track the field with high accuracy and that we are able to correct the data for the adverse effects of the changing field.

Quantitative assessment of the susceptibility artefact and its interaction with motion in diffusion MRI

In this paper we evaluate the three main methods for correcting the susceptibility-induced artefact in diffusion-weighted magnetic-resonance (DW-MR) data, and assess how correction is affected by the susceptibility field's interaction with motion. The susceptibility artefact adversely impacts analysis performed on the data and is typically corrected in post-processing. Correction strategies involve either registration to a structural image, the application of an acquired field-map or the use of additional images acquired with different phase-encoding. Unfortunately, the choice of which method to use is made difficult by the absence of any systematic comparisons of them. In this work we quantitatively evaluate these methods, by extending and employing a recently proposed framework that allows for the simulation of realistic DW-MR datasets with artefacts. Our analysis separately evaluates the ability for methods to correct for geometric distortions and to recover lost information in regions of signal compression. In terms of geometric distortions, we find that registration-based methods offer the poorest correction. Field-mapping techniques are better, but are influenced by noise and partial volume effects, whilst multiple phase-encode methods performed best. We use our simulations to validate a popular surrogate metric of correction quality, the comparison of corrected data acquired with AP and LR phase-encoding, and apply this surrogate to real datasets. Furthermore, we demonstrate that failing to account for the interaction of the susceptibility field with head movement leads to increased errors when analysing DW-MR data. None of the commonly used post-processing methods account for this interaction, and we suggest this may be a valuable area for future methods development.

Automated Detection of Motion Artefacts in MR Imaging Using Decision Forests

The acquisition of a Magnetic Resonance (MR) scan usually takes longer than subjects can remain still. Movement of the subject such as bulk patient motion or respiratory motion degrades the image quality and its diagnostic value by producing image artefacts like ghosting, blurring, and smearing. This work focuses on the effect of motion on the reconstructed slices and the detection of motion artefacts in the reconstruction by using a supervised learning approach based on random decision forests. Both the effects of bulk patient motion occurring at various time points in the acquisition on head scans and the effects of respiratory motion on cardiac scans are studied. Evaluation is performed on synthetic images where motion artefacts have been introduced by altering the k-space data according to a motion trajectory, using the three common k-space sampling patterns: Cartesian, radial, and spiral. The results suggest that a machine learning approach is well capable of learning the characteristics of motion artefacts and subsequently detecting motion artefacts with a confidence that depends on the sampling pattern.

Towards a comprehensive framework for movement and distortion correction of diffusion MR images: Within volume movement

Most motion correction methods work by aligning a set of volumes together, or to a volume that represents a reference location. These are based on an implicit assumption that the subject remains motionless during the several seconds it takes to acquire all slices in a volume, and that any movement occurs in the brief moment between acquiring the last slice of one volume and the first slice of the next. This is clearly an approximation that can be more or less good depending on how long it takes to acquire one volume and in how rapidly the subject moves. In this paper we present a method that increases the temporal resolution of the motion correction by modelling movement as a piecewise continous function over time. This intra-volume movement correction is implemented within a previously presented framework that simultaneously estimates distortions, movement and movement-induced signal dropout. We validate the method on highly realistic simulated data containing all of these effects. It is demonstrated that we can estimate the true movement with high accuracy, and that scalar parameters derived from the data, such as fractional anisotropy, are estimated with greater fidelity when data has been corrected for intra-volume movement. Importantly, we also show that the difference in fidelity between data affected by different amounts of movement is much reduced when taking intra-volume movement into account. Additional validation was performed on data from a healthy volunteer scanned when lying still and when performing deliberate movements. We show an increased correspondence between the “still” and the “movement” data when the latter is corrected for intra-volume movement. Finally we demonstrate a big reduction in the telltale signs of intra-volume movement in data acquired on elderly subjects.

•

We introduce a method to correct for intra-volume movement into an existing framework for movement and distortion correction.

•

It has been validated on realistic simulated data and on experimental data with deliberate movement.

•

The results indicate that one can reliably reverse the adverse effects of intra-volume movement.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

Realistic simulation of artefacts in diffusion MRI for validating post-processing correction techniques

In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. In the literature, improved alignment is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We apply this framework in three ways. Firstly we assess the most commonly used method for artefact correction, FSL's eddy_correct, and compare it to a recently proposed alternative, eddy. We demonstrate quantitatively that using eddy_correct leads to significant errors in the corrected data, whilst eddy is able to provide much improved correction. Secondly we investigate the datasets required to achieve good correction with eddy, by looking at the minimum number of directions required and comparing the recommended full-sphere acquisitions to equivalent half-sphere protocols. Finally, we investigate the impact of correction quality by examining the fits from microstructure models to real and simulated data.

A Simulation Framework for Quantitative Validation of Artefact Correction in Diffusion MRI

In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. The ultimate aim of these techniques is improved image alignment, yet in the literature this is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We demonstrate the application of this framework by testing one of the most commonly used methods for EC correction, registration of DWIs to b = 0, and reveal the systematic bias this introduces into corrected datasets.

Isolated assessment of translation or rotation severely underestimates the effects of subject motion in fMRI data

Subject motion has long since been known to be a major confound in functional MRI studies of the human brain. For resting-state functional MRI in particular, data corruption due to motion artefacts has been shown to be most relevant. However, despite 6 parameters (3 for translations and 3 for rotations) being required to fully describe the head's motion trajectory between timepoints, not all are routinely used to assess subject motion. Using structural (n = 964) as well as functional MRI (n = 200) data from public repositories, a series of experiments was performed to assess the impact of using a reduced parameter set (translation_only and rotation_only) versus using the complete parameter set. It could be shown that the usage of 65 mm as an indicator of the average cortical distance is a valid approximation in adults, although care must be taken when comparing children and adults using the same measure. The effect of using slightly smaller or larger values is minimal. Further, both translation_only and rotation_only severely underestimate the full extent of subject motion; consequently, both translation_only and rotation_only discard substantially fewer datapoints when used for quality control purposes (“motion scrubbing”). Finally, both translation_only and rotation_only severely underperform in predicting the full extent of the signal changes and the overall variance explained by motion in functional MRI data. These results suggest that a comprehensive measure, taking into account all available parameters, should be used to characterize subject motion in fMRI.

Dysregulated daily rhythmicity of neuronal resting-state networks in MCI patients

In young healthy participants, the degree of daily rhythmicity largely varies across different neuronal resting-state networks (RSNs), while it is to date unknown whether this temporal pattern of activity is conserved in healthy and pathological aging. Twelve healthy elderly (mean age=65.1±5.7 years) and 12 patients with amnestic mild cognitive impairment (aMCI; mean age=69.6±6.2 years) underwent four resting-state functional magnetic resonance imaging scans at fixed 2.5 h intervals throughout a day. Time courses of a RSN were extracted by a connectivity strength and a spatial extent approach performed individually for each participant. Highly rhythmic RSNs included a sensorimotor, a cerebellar and a visual network in healthy elderly; the least rhythmic RSNs in this group included a network associated with executive control and an orbitofrontal network. The degree of daily rhythmicity in aMCI patients was reduced and dysregulated. For healthy elderly, the findings are in accordance with results reported for young healthy participants suggesting a comparable distribution of daily rhythmicity across RSNs during healthy aging. In contrast, the reduction and dysregulation of daily rhythmicity observed in aMCI patients is presumably indicative of underlying neurodegenerative processes in this group.

Susceptibility artefact correction using dynamic graph cuts: application to neurosurgery

Echo Planar Imaging (EPI) is routinely used in diffusion and functional MR imaging due to its rapid acquisition time. However, the long readout period makes it prone to susceptibility artefacts which results in geometric and intensity distortions of the acquired image. The use of these distorted images for neuronavigation hampers the effectiveness of image-guided surgery systems as critical white matter tracts and functionally eloquent brain areas cannot be accurately localised. In this paper, we present a novel method for correction of distortions arising from susceptibility artefacts in EPI images. The proposed method combines fieldmap and image registration based correction techniques in a unified framework. A phase unwrapping algorithm is presented that can efficiently compute the B₀ magnetic field inhomogeneity map as well as the uncertainty associated with the estimated solution through the use of dynamic graph cuts. This information is fed to a subsequent image registration step to further refine the results in areas with high uncertainty. This work has been integrated into the surgical workflow at the National Hospital for Neurology and Neurosurgery and its effectiveness in correcting for geometric distortions due to susceptibility artefacts is demonstrated on EPI images acquired with an interventional MRI scanner during neurosurgery.

High performance MRI simulations of motion on multi-GPU systems

BACKGROUND

MRI physics simulators have been developed in the past for optimizing imaging protocols and for training purposes. However, these simulators have only addressed motion within a limited scope. The purpose of this study was the incorporation of realistic motion, such as cardiac motion, respiratory motion and flow, within MRI simulations in a high performance multi-GPU environment.

METHODS

Three different motion models were introduced in the Magnetic Resonance Imaging SIMULator (MRISIMUL) of this study: cardiac motion, respiratory motion and flow. Simulation of a simple Gradient Echo pulse sequence and a CINE pulse sequence on the corresponding anatomical model was performed. Myocardial tagging was also investigated. In pulse sequence design, software crushers were introduced to accommodate the long execution times in order to avoid spurious echoes formation. The displacement of the anatomical model isochromats was calculated within the Graphics Processing Unit (GPU) kernel for every timestep of the pulse sequence. Experiments that would allow simulation of custom anatomical and motion models were also performed. Last, simulations of motion with MRISIMUL on single-node and multi-node multi-GPU systems were examined.

RESULTS

Gradient Echo and CINE images of the three motion models were produced and motion-related artifacts were demonstrated. The temporal evolution of the contractility of the heart was presented through the application of myocardial tagging. Better simulation performance and image quality were presented through the introduction of software crushers without the need to further increase the computational load and GPU resources. Last, MRISIMUL demonstrated an almost linear scalable performance with the increasing number of available GPU cards, in both single-node and multi-node multi-GPU computer systems.

CONCLUSIONS

MRISIMUL is the first MR physics simulator to have implemented motion with a 3D large computational load on a single computer multi-GPU configuration. The incorporation of realistic motion models, such as cardiac motion, respiratory motion and flow may benefit the design and optimization of existing or new MR pulse sequences, protocols and algorithms, which examine motion related MR applications.

SimPACE: generating simulated motion corrupted BOLD data with synthetic-navigated acquisition for the development and evaluation of SLOMOCO: a new, highly effective slicewise motion correction

Head motion in functional MRI and resting-state MRI is a major problem. Existing methods do not robustly reflect the true level of motion artifact for in vivo fMRI data. The primary issue is that current methods assume that motion is synchronized to the volume acquisition and thus ignore intra-volume motion. This manuscript covers three sections in the use of gold-standard motion-corrupted data to pursue an intra-volume motion correction. First, we present a way to get motion corrupted data with accurately known motion at the slice acquisition level. This technique simulates important data acquisition-related motion artifacts while acquiring real BOLD MRI data. It is based on a novel motion-injection pulse sequence that introduces known motion independently for every slice: Simulated Prospective Acquisition CorrEction (SimPACE). Secondly, with data acquired using SimPACE, we evaluate several motion correction and characterization techniques, including several commonly used BOLD signal- and motion parameter-based metrics. Finally, we introduce and evaluate a novel, slice-based motion correction technique. Our novel method, SLice-Oriented MOtion COrrection (SLOMOCO) performs better than the volumetric methods and, moreover, accurately detects the motion of independent slices, in this case equivalent to the known injected motion. We demonstrate that SLOMOCO can model and correct for nearly all effects of motion in BOLD data. Also, none of the commonly used motion metrics was observed to robustly identify motion corrupted events, especially in the most realistic scenario of sudden head movement. For some popular metrics, performance was poor even when using the ideal known slice motion instead of volumetric parameters. This has negative implications for methods relying on these metrics, such as recently proposed motion correction methods such as data censoring and global signal regression.

MRISIMUL: a GPU-based parallel approach to MRI simulations

A new step-by-step comprehensive MR physics simulator (MRISIMUL) of the Bloch equations is presented. The aim was to develop a magnetic resonance imaging (MRI) simulator that makes no assumptions with respect to the underlying pulse sequence and also allows for complex large-scale analysis on a single computer without requiring simplifications of the MRI model. We hypothesized that such a simulation platform could be developed with parallel acceleration of the executable core within the graphic processing unit (GPU) environment. MRISIMUL integrates realistic aspects of the MRI experiment from signal generation to image formation and solves the entire complex problem for densely spaced isochromats and for a densely spaced time axis. The simulation platform was developed in MATLAB whereas the computationally demanding core services were developed in CUDA-C. The MRISIMUL simulator imaged three different computer models: a user-defined phantom, a human brain model and a human heart model. The high computational power of GPU-based simulations was compared against other computer configurations. A speedup of about 228 times was achieved when compared to serially executed C-code on the CPU whereas a speedup between 31 to 115 times was achieved when compared to the OpenMP parallel executed C-code on the CPU, depending on the number of threads used in multithreading (2-8 threads). The high performance of MRISIMUL allows its application in large-scale analysis and can bring the computational power of a supercomputer or a large computer cluster to a single GPU personal computer.

Bloch-based MRI system simulator considering realistic electromagnetic fields for calculation of signal, noise, and specific absorption rate

PURPOSE

To describe and introduce new software capable of accurately simulating MR signal, noise, and specific absorption rate (SAR) given arbitrary sample, sequence, static magnetic field distribution, and radiofrequency magnetic and electric field distributions for each transmit and receive coil.

THEORY AND METHODS

Using fundamental equations for nuclear precession and relaxation, signal reception, noise reception, and calculation of SAR, a versatile MR simulator was developed. The resulting simulator was tested with simulation of a variety of sequences demonstrating several common imaging contrast types and artifacts. The simulation of intravoxel dephasing and rephasing with both tracking of the first order derivatives of each magnetization vector and multiple magnetization vectors was examined to ensure adequate representation of the MR signal. A quantitative comparison of simulated and experimentally measured SNR was also performed.

RESULTS

The simulator showed good agreement with our expectations, theory, and experiment.

CONCLUSION

With careful design, an MR simulator producing realistic signal, noise, and SAR for arbitrary sample, sequence, and fields has been created. It is hoped that this tool will be valuable in a wide variety of applications.

A virtual imaging platform for multi-modality medical image simulation

This paper presents the Virtual Imaging Platform (VIP), a platform accessible at http://vip.creatis.insa-lyon.fr to facilitate the sharing of object models and medical image simulators, and to provide access to distributed computing and storage resources. A complete overview is presented, describing the ontologies designed to share models in a common repository, the workflow template used to integrate simulators, and the tools and strategies used to exploit computing and storage resources. Simulation results obtained in four image modalities and with different models show that VIP is versatile and robust enough to support large simulations. The platform currently has 200 registered users who consumed 33 years of CPU time in 2011.

Prospective active marker motion correction improves statistical power in BOLD fMRI

Group level statistical maps of blood oxygenation level dependent (BOLD) signals acquired using functional magnetic resonance imaging (fMRI) have become a basic measurement for much of systems, cognitive and social neuroscience. A challenge in making inferences from these statistical maps is the noise and potential confounds that arise from the head motion that occurs within and between acquisition volumes. This motion results in the scan plane being misaligned during acquisition, ultimately leading to reduced statistical power when maps are constructed at the group level. In most cases, an attempt is made to correct for this motion through the use of retrospective analysis methods. In this paper, we use a prospective active marker motion correction (PRAMMO) system that uses radio frequency markers for real-time tracking of motion, enabling on-line slice plane correction. We show that the statistical power of the activation maps is substantially increased using PRAMMO compared to conventional retrospective correction. Analysis of our results indicates that the PRAMMO acquisition reduces the variance without decreasing the signal component of the BOLD (beta). Using PRAMMO could thus improve the overall statistical power of fMRI based BOLD measurements, leading to stronger inferences of the nature of processing in the human brain.

How ignoring physiological noise can bias the conclusions from fMRI simulation results

Neuroimaging researchers use simulation studies to validate their statistical methods because it is acknowledged that this is the most feasible way to know the ground truth of the data. The noise model used in these studies typically varies from a simple Gaussian distribution to an estimate of the noise distribution from real data. However, although several studies point out the presence of physiological noise in fMRI data, this noise source is currently lacking in simulation studies. Therefore, we explored the impact of adding physiological noise to the simulated data. For several experimental designs, fMRI data were generated under different noise models while the signal-to-noise ratio was kept constant. The sensitivity and specificity of a standard statistical parametric mapping (SPM) analysis were determined by comparing the known activation with the detected activation. We show that by including physiological noise in the data generation process, the simulation results in terms of sensitivity and specificity drop dramatically. Additionally, we used the new proposed simulation model to compare a standard SPM analysis against the method proposed by Cabella et al. (2009). The results indicate that the analysis of data containing no physiological noise yields a better performance of the SPM analysis. However, if physiological noise is included in the data, the sensitivity and specificity of the Cabella method are higher compared to the SPM analysis. Based on these results, we argue that the results of current simulation studies are likely to be biased, especially when analysis methods are compared using ROC curves.

Classifying fMRI-derived resting-state connectivity patterns according to their daily rhythmicity

The vast majority of biological functions express rhythmic fluctuations across the 24-hour day. We investigated the degree of daily modulation across fMRI (functional Magnetic Resonance Imaging) derived resting-state data in 15 subjects by evaluating the time courses of 20 connectivity patterns over 8h (4 sessions). For each subject, we determined the chronotype, which describes the relationship between the individual circadian rhythm and the local time. We could therefore analyze the daily time course of the connectivity patterns controlling for internal time. Furthermore, as the participants' scan times were staggered as a function of their chronotype, we prevented sleep deprivation and kept time awake constant across subjects. Individual functional connectivity within each connectivity pattern was defined at each session as connectivity strength measured by a mean z-value and, in addition, as the spatial extent expressed by the number of activated voxels. Highly rhythmic connectivity patterns included two sub-systems of the Default-Mode Network (DMN) and a network extending over sensori-motor regions. The network characterized as the most stable across the day is mainly associated with processing of executive control. We conclude that the degree of daily modulation largely varies across fMRI derived resting-state connectivity patterns, ranging from highly rhythmic to stable. This finding should be considered when interpreting results from fMRI studies.

Simultaneous multislice inversion contrast imaging using power independent of the number of slices (PINS) and delays alternating with nutation for tailored excitation (DANTE) radio frequency pulses

A method for simultaneous multislice (SMS) inversion contrast imaging is presented using a combination of the delays alternating with nutation for tailored excitation (DANTE) and the power independent of the number of slices (PINS) techniques. In SMS imaging, simultaneously excited slices result in an aliased image that is disentangled using parallel imaging reconstruction techniques. At high-magnetic field strengths, the peak amplitude and specific absorption rate of conventional (summed) SMS radio frequency pulses can be prohibitively high. Using the PINS approach, specific absorption rate is independent of the number of slices allowing high SMS acceleration factors even at high fields. Using DANTE, adiabatic SMS radio frequency pulses can be created to be combined with PINS. This allows 2D imaging protocols that employ adiabatic pulses to also reap the benefits of low specific absorption rate SMS acceleration. As a proof-of-concept, simulations and measurements using hyperbolic secant inversion pulses are shown.

SimTB, a simulation toolbox for fMRI data under a model of spatiotemporal separability

We introduce SimTB, a MATLAB toolbox designed to simulate functional magnetic resonance imaging (fMRI) datasets under a model of spatiotemporal separability. The toolbox meets the increasing need of the fMRI community to more comprehensively understand the effects of complex processing strategies by providing a ground truth that estimation methods may be compared against. SimTB captures the fundamental structure of real data, but data generation is fully parameterized and fully controlled by the user, allowing for accurate and precise comparisons. The toolbox offers a wealth of options regarding the number and configuration of spatial sources, implementation of experimental paradigms, inclusion of tissue-specific properties, addition of noise and head movement, and much more. A straightforward data generation method and short computation time (3-10 seconds for each dataset) allow a practitioner to simulate and analyze many datasets to potentially understand a problem from many angles. Beginning MATLAB users can use the SimTB graphical user interface (GUI) to design and execute simulations while experienced users can write batch scripts to automate and customize this process. The toolbox is freely available at http://mialab.mrn.org/software together with sample scripts and tutorials.

On the numerically predicted spatial BOLD fMRI specificity for spin echo sequences

This work utilises general numerical magnetic resonance imaging MRI simulations to predict the spatial specificity of the blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signal. A Monte Carlo simulation approach was utilized on a microvascular structure consisting of randomly oriented cylinders representing blood vessels. This framework was employed to numerically investigate the spatial specificity, defined as ratio of pial vessel to microvascular signal, of the spin echo BOLD fMRI signal as a function of field strength, echo time and tissue types [grey matter (GM) and cerebrospinal fluid (CSF), respectively]. Spatial specificity of spin echo BOLD fMRI signal was determined to increase with field strength up to 16 T and with maximal specificity for echo time shorter than tissue T(2). In addition, it was found that, for large pial vessels, the extravascular signal decay could not be described using the oversimplified but nevertheless commonly employed mono-exponential signal decay approximation (MEA). Consequently, a recently proposed model relying on the MEA deviates substantially from our results on the spatial specificity. A refinement of this model is proposed based on an available, more detailed signal description. Finally, the effect of CSF on the spatial specificity was investigated. While a large spatial specificity of the spin echo BOLD fMRI signal is observed if a pial vessel is surrounded by grey matter, this is greatly reduced for a pial vessel situated on a GM/CSF interface, rendering the suppression of pial vessels on the cortex surface unlikely.

An iterative two-threshold analysis for single-subject functional MRI of the human brain

OBJECTIVES

Current thresholding strategies for the analysis of functional MRI (fMRI) datasets may suffer from specific limitations (e.g. with respect to the required smoothness) or lead to reduced performance for a low signal-to-noise ratio (SNR). Although a previously proposed two-threshold (TT) method offers a promising solution to these problems, the use of preset settings limits its performance. This work presents an optimised TT approach that estimates the required parameters in an iterative manner.

METHODS

The iterative TT (iTT) method is compared with the original TT method, as well as other established voxel-based and cluster-based thresholding approaches and spatial mixture modelling (SMM) for both simulated data and fMRI of a hometown walking task at different experimental settings (spatial resolution, filtering and SNR).

RESULTS

In general, the iTT method presents with remarkable sensitivity and good specificity that outperforms all conventional approaches tested except for SMM in a few cases. This also holds true for challenging conditions such as high spatial resolution, the absence of filtering, high noise level, or a low number of task repetitions.

CONCLUSION

Thus, iTT emerges as a good candidate for both scientific fMRI studies at high spatial resolution and more routine applications for clinical purposes.

An empirical investigation of motion effects in eMRI of interictal epileptiform spikes

We recently developed a functional neuroimaging technique called encephalographic magnetic resonance imaging (eMRI). Our method acquires rapid single-shot gradient-echo echo-planar MRI (repetition time=47 ms); it attempts to measure an MR signal more directly linked to neuronal electromagnetic activity than existing methods. To increase the likelihood of detecting such an MR signal, we recorded concurrent MRI and scalp electroencephalography (EEG) during fast (20-200 ms), localized, high-amplitude (>50 μV on EEG) cortical discharges in a cohort of focal epilepsy patients. Seen on EEG as interictal spikes, these discharges occur in between seizures and induced easily detectable MR magnitude and phase changes concurrent with the spikes with a lag of milliseconds to tens of milliseconds. Due to the time scale of the responses, localized changes in blood flow or hemoglobin oxygenation are unlikely to cause the MR signal changes that we observed. While the precise underlying mechanisms are unclear, in this study, we empirically investigate one potentially important confounding variable - motion. Head motion in the scanner affects both EEG and MR recording. It can produce brief "spike-like" artifacts on EEG and induce large MR signal changes similar to our interictal spike-related signal changes. In order to explore the possibility that interictal spikes were associated with head motions (although such an association had never been reported), we had previously tracked head position in epilepsy patients during interictal spikes and explicitly demonstrated a lack of associated head motion. However, that study was performed outside the MR scanner, and the root-mean-square error in the head position measurement was 0.7 mm. The large inaccuracy in this measurement therefore did not definitively rule out motion as a possible signal generator. In this study, we instructed healthy subjects to make deliberate brief (<500 ms) head motions inside the MR scanner and imaged these head motions with concurrent EEG and MRI. We compared these artifactual MR and EEG data to genuine interictal spikes. While per-voxel MR and per-electrode EEG time courses for the motion case can mimic the corresponding time courses associated with a genuine interictal spike, head motion can be unambiguously differentiated from interictal spikes via scalp EEG potential maps. Motion induces widespread changes in scalp potential, whereas interictal spikes are localized and have a regional fall-off in amplitude. These findings make bulk head motion an unlikely generator of the large spike-related MR signal changes that we had observed. Further work is required to precisely identify the underlying mechanisms.