Interleaved EPI based fMRI improved by multiplexed sensitivity encoding (MUSE) and simultaneous multi-band imaging

© 2014 Chang et al. Functional magnetic resonance imaging (fMRI) is a non-invasive and powerful imaging tool for detecting brain activities. The majority of fMRI studies are performed with single-shot echo-planar imaging (EPI) due to its high temporal resolution. Recent studies have demonstrated that, by increasing the spatial-resolution of fMRI, previously unidentified neuronal networks can be measured. However, it is challenging to improve the spatial resolution of conventional single-shot EPI based fMRI. Although multi-shot interleaved EPI is superior to single-shot EPI in terms of the improved spatial-resolution, reduced geometric distortions, and sharper point spread function (PSF), interleaved EPI based fMRI has two main limitations: 1) the imaging throughput is lower in interleaved EPI; 2) the magnitude and phase signal variations among EPI segments (due to physiological noise, subject motion, and B0 drift) are translated to significant in-plane aliasing artifact across the field of view (FOV). Here we report a method that integrates multiple approaches to address the technical limitations of interleaved EPI-based fMRI. Firstly, the multiplexed sensitivity-encoding (MUSE) post-processing algorithm is used to suppress in-plane aliasing artifacts resulting from time-domain signal instabilities during dynamic scans. Secondly, a simultaneous multi-band interleaved EPI pulse sequence, with a controlled aliasing scheme incorporated, is implemented to increase the imaging throughput. Thirdly, the MUSE algorithm is then generalized to accommodate fMRI data obtained with our multi-band interleaved EPI pulse sequence, suppressing both in-plane and through-plane aliasing artifacts. The blood-oxygenation-level-dependent (BOLD) signal detectability and the scan throughput can be significantly improved for interleaved EPI-based fMRI. Our human fMRI data obtained from 3 Tesla systems demonstrate the effectiveness of the developed methods. It is expected that future fMRI studies requiring high spatial-resolvability and fidelity will largely benefit from the reported techniques.published_or_final_versio

Applications of the golden angle in cardiovascular MRI

The use of radial trajectories has been seen as a potential solution to highly efficient cardiovascular magnetic resonance imaging (MRI). By acquiring a broad range of spatial frequencies per repetition time, the acquisition is time-efficient and robust against motion. Of particular interest is the golden angle profile order, which promises a near-uniform k-space coverage for an arbitrary number of readouts, enabling flexible data resorting, which is critical for efficient cardiovascular MRI. In Study I the use of 2D golden angle profile ordering is explored for imaging pulmonary embolisms. The insensitivity to motion and flow is used to reduce the artifacts that otherwise degrade images of the pulmonary vasculature when imaging with thin slices. It was found that the proposed technique could improve the image quality. Another source of artifacts arises when gradients are rapidly switched, and local induction of eddy currents may perturb spin equilibrium. In Study II, we propose a generalized golden angle profile orderings in 3D which reduces eddy-current artifacts. We demonstrate the efficacy of our generalization through numerical simulations, phantom imaging and imaging of a healthy volunteer. In Study III an improved 2D golden angle profile ordering was explored which resulted in a higher degree of k-space uniformity after physiological binning. This novel profile ordering was used in combination with a phase-contrast readout to enable quantification of myocardial tissue velocity and transmitral blood flow velocity, which are essential parameters for diastolic function assessment. When compared to echocardiography, it was found that MRI could accurately quantify myocardial tissue velocity, whereas transmitral blood flow velocity was underestimated. Study IV explored a further development of Study III by proposing a 3D version of the improved profile ordering. This novel ordering was used to acquire whole-heart functional images during free-breathing in less than one minute. Together, these results indicate that golden-angle-based imaging has the potential to improve cardiovascular MRI in several areas

Optimal Correction of The Slice Timing Problem and Subject Motion Artifacts in fMRI

Functional magnetic resonance imaging (fMRI) is an extremely popular investigative and clinical imaging tool that allows safe and noninvasive study of the functional living brain. Fundamentally, fMRI measures a physiological signal as it changes over time. The manner in which this spatio-temporal signal is acquired can create technical challenges during image reconstruction that must be corrected for if any meaningful information is to be extracted from the data. Two particular challenges that are fundamentally intertwined with each other are temporal misalignment and spatial misalignment. Temporal misalignment is due to the nature of fMRI acquisition protocols themselves: a 3D volume is created by sampling and stacking multiple 2D slices. However, these slices are not acquired simultaneously or sequentially, and therefore will always be temporally misaligned with each other. Spatial misalignment arises when subject motion is present during the scan, resulting in individual volumes being spatially misaligned with each other. Spatial and temporal misalignment are not independent from each other, and their interaction can cause additional artifacts and reconstruction challenges if not addressed properly. The purpose of this thesis is to critically examine the problem of both spatial and temporal misalignment from a signal processing perspective, while considering the physical nature and origin of the signal itself, and develop optimal correction routines for spatial and temporal misalignment and their associated artifacts. One of the most immediate problems associated with temporal misalignment is that the order in which the slices are acquired must be known in order for correction to be possible. Surprisingly, this information is rarely provided with old or shared data, meaning that this critical preprocessing step must be skipped, significantly lowering the value of the data. We use the spatio-temporal properties of the fMRI signal to develop a robust and accurate algorithm to infer the slice acquisition order retrospectively from any fMRI scan. The ability to extract the interleave parameter from any data set allows us to perform slice timing correction even if this information had been lost, or was not provided with the scan. In the next section of this work, we develop a new optimal method of slice timing correction (Filter-Shift) based on the fundamental properties of sampling theory in digital signal processing. By examining the properties of the signal of interest (The blood oxygen level depended signal: BOLD signal), we are able to design and implement an effective FIR filter to simultaneously remove noise and reconstruct the signal of interest at any shifted offset, without the need for sub-optimal interpolation. In the final section, we investigate the effects of different motion types on the MR signal based on the Bloch equation, in order to develop a theoretical foundation from which we can create an optimal correction method. We devise a novel method to remove these artifacts: Discrete reconstruction of irregular fMRI trajectory (DRIFT). Our method calculates the exact displacement of the k-space samples due to motion at each dwell time and retrospectively corrects each slice of the fMRI volume using an inverse nonuniform Fourier transform. We conclude that a hybrid approach with both prospective and retrospective components are essentially required for optimal removal of motion artifacts from the fMRI data. The combined work of this thesis provides two theoretically sound and extremely effective correction routines, that both remove artifacts and restore the underlying sampled signal. Motion correction and slice timing correction are typically the first two preprocessing steps to be applied to any fMRI data, and thus provide the foundation for any further analysis. While many other preprocessing steps can be omitted or included depending on the analysis, motion correction and slice timing correction are unequivocally beneficial and necessary for accurate and reliable results. This work provides a theoretical and quantitative framework that describes the optimal removal of artifacts associated with motion and slice timing

Methods for cleaning the BOLD fMRI signal

Available online 9 December 2016 http://www.sciencedirect.com/science/article/pii/S1053811916307418?via%3Dihubhttp://www.sciencedirect.com/science/article/pii/S1053811916307418?via%3DihubBlood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.This work was supported by the Spanish Ministry of Economy and Competitiveness [Grant PSI 2013–42343 Neuroimagen Multimodal], the Severo Ochoa Programme for Centres/Units of Excellence in R & D [SEV-2015-490], and the research and writing of the paper were supported by the NIMH and NINDS Intramural Research Programs (ZICMH002888) of the NIH/HHS

The current state-of-the-art of spinal cord imaging: methods.

A first-ever spinal cord imaging meeting was sponsored by the International Spinal Research Trust and the Wings for Life Foundation with the aim of identifying the current state-of-the-art of spinal cord imaging, the current greatest challenges, and greatest needs for future development. This meeting was attended by a small group of invited experts spanning all aspects of spinal cord imaging from basic research to clinical practice. The greatest current challenges for spinal cord imaging were identified as arising from the imaging environment itself; difficult imaging environment created by the bone surrounding the spinal canal, physiological motion of the cord and adjacent tissues, and small cross-sectional dimensions of the spinal cord, exacerbated by metallic implants often present in injured patients. Challenges were also identified as a result of a lack of "critical mass" of researchers taking on the development of spinal cord imaging, affecting both the rate of progress in the field, and the demand for equipment and software to manufacturers to produce the necessary tools. Here we define the current state-of-the-art of spinal cord imaging, discuss the underlying theory and challenges, and present the evidence for the current and potential power of these methods. In two review papers (part I and part II), we propose that the challenges can be overcome with advances in methods, improving availability and effectiveness of methods, and linking existing researchers to create the necessary scientific and clinical network to advance the rate of progress and impact of the research

Semi-PROPELLER Compressed Sensing Image Reconstruction with Enhanced Resolution in MRI

Magnetic Resonance Imaging (MRI) reconstruction algorithm using semi-PROPELLER compressed sensing is pre- sented in this paper. It is exhibited that introduced algorithm for estimating data shifts is feasible when super- resolution is applied. The offered approach utilizes compressively sensed MRI PROPELLER sequences and improves MR images spatial resolution in circumstances when highly undersampled k-space trajectories are applied. Compressed Sensing (CS) aims at signal and images reconstructing from significantly fewer measurements than were conventionally assumed necessary. Compressed sensing (CS) aims at signal and images reconstructing from significantly fewer measurements than were traditionally thought necessary. It is shown that the presented approach improves MR spatial resolution in cases when Compressed Sensing (CS) sequences are used. The application of CS in medical modalities has the potential for significant scan time reductions, with visible benefits for patients and health care economics. These methods emphasize on maximizing image sparsity on known sparse transform do- main and minimizing fidelity. This diagnostic modality struggles with an inherently slow data acquisition process. The use of CS to MRI leads to substantial scan time reductions and visible benefits for patients and economic factors. In this report the objective is to combine Super-Resolution image enhancement algorithm with both PROPELLER sequence and CS framework. All the techniques emphasize on maximizing image sparsity on known sparse transform domain and minimizing fidelity. The motion estimation algorithm being a part of super resolution reconstruction (SRR) estimates shifts for all blades jointly, emphasizing blade-pair correlations that are both strong and more robust to noise.

Doctor of Philosophy

dissertationThis dissertation presents original research that improves the ability of magnetic resonance imaging (MRI) to measure temperature in aqueous tissue using the proton resonance frequency (PRF) shift and T1 measurements in fat tissue in order to monitor focused ultrasound (FUS) treatments. The inherent errors involved in measuring the longitudinal relaxation time T1 using the variable flip angle method with a two-dimensional (2D) acquisition are presented. The edges of the slice profile can contribute a significant amount of signal for large flip angles at steady state, which causes significant errors in the T1 estimate. Only a narrow range of flip angle combinations provided accurate T1 estimates. Respiration motion causes phase artifacts, which lead to errors when measuring temperature changes using the PRF method. A respiration correction method for 3D imaging temperature of the breast is presented. Free induction decay (FID) navigators were used to measure and correct phase offsets induced by respiration. The precision of PRF temperature measurements within the breast was improved by an average factor of 2.1 with final temperature precision of approximately 1 °C. Locating the position of the ultrasound focus in MR coordinates of an ultrasound transducer with multiple degrees of freedom can be difficult. A rapid method for predicting the position using 3 tracker coils with a special MRI pulse iv sequence is presented. The Euclidean transformation of the coil's current positions to their calibration positions was used to predict the current focus position. The focus position was predicted to within approximately 2.1 mm in less than 1 s. MRI typically has tradeoffs between imaging field of view and spatial and temporal resolution. A method for acquiring a large field of view with high spatial and temporal resolution is presented. This method used a multiecho pseudo-golden angle stack of stars imaging sequence to acquire the large field of view with high spatial resolution and k-space weighted image contrast (KWIC) to increase the temporal resolution. The pseudo-golden angle allowed for removal of artifacts introduced by the KWIC reconstruction algorithm. The multiple echoes allowed for high readout bandwidth to reduce blurring due to off resonance and chemical shift as well as provide separate water/fat images, estimates of the initial signal magnitude M(0), T2 * time constant, and combination of echo phases. The combined echo phases provided significant improvement to the PRF temperature precision, and ranged from ~0.3-1.0 °C within human breast. M(0) and T2 * values can possibly be used as a measure of temperature in fat