Spatio-temporal wavelet regularization for parallel MRI reconstruction: application to functional MRI

Parallel MRI is a fast imaging technique that enables the acquisition of highly resolved images in space or/and in time. The performance of parallel imaging strongly depends on the reconstruction algorithm, which can proceed either in the original k-space (GRAPPA, SMASH) or in the image domain (SENSE-like methods). To improve the performance of the widely used SENSE algorithm, 2D- or slice-specific regularization in the wavelet domain has been deeply investigated. In this paper, we extend this approach using 3D-wavelet representations in order to handle all slices together and address reconstruction artifacts which propagate across adjacent slices. The gain induced by such extension (3D-Unconstrained Wavelet Regularized -SENSE: 3D-UWR-SENSE) is validated on anatomical image reconstruction where no temporal acquisition is considered. Another important extension accounts for temporal correlations that exist between successive scans in functional MRI (fMRI). In addition to the case of 2D+t acquisition schemes addressed by some other methods like kt-FOCUSS, our approach allows us to deal with 3D+t acquisition schemes which are widely used in neuroimaging. The resulting 3D-UWR-SENSE and 4D-UWR-SENSE reconstruction schemes are fully unsupervised in the sense that all regularization parameters are estimated in the maximum likelihood sense on a reference scan. The gain induced by such extensions is illustrated on both anatomical and functional image reconstruction, and also measured in terms of statistical sensitivity for the 4D-UWR-SENSE approach during a fast event-related fMRI protocol. Our 4D-UWR-SENSE algorithm outperforms the SENSE reconstruction at the subject and group levels (15 subjects) for different contrasts of interest (eg, motor or computation tasks) and using different parallel acceleration factors (R=2 and R=4) on 2x2x3mm3 EPI images.Comment: arXiv admin note: substantial text overlap with arXiv:1103.353

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

Functional MRI Data Analysis Techniques and Strategies to Map the Olfactory System of a Rat Brain.

Understanding mysteries of a brain represents one of the great challenges for modern science. Functional magnetic resonance imaging (fMRI) has two features that make it unique amongst other imaging modalities used in behavioral neuroscience. First, it can be entirely non-invasive and second, fMRI has the spatial and temporal resolution to resolve patterns of neuronal activity across the entire brain in less than a minute. fMRI indirectly detects neural activity in different parts of the brain by comparing contrast in MR signal intensity prior to and following stimulation. Areas of the brain with increased synaptic and neuronal activity require increased levels of oxygen to sustain this activity. Enhanced brain activity is accompanied by an increase in metabolism followed by increases in blood flow and blood volume. The enhanced blood flow usually exceeds the metabolic demand exposing the active brain area to high level of oxygenated hemoglobin. Oxygenated hemoglobin increases the MR signal intensity that can be detected in MR scanner. This relatively straight forward scenario is, unfortunately, oversimplified. The fMRI signal change to noise ratio is extremely small. In this work a quantitative analysis strategy to analyze fMRI data was successfully developed, implemented and optimized for the rat brain. Therein, each subject is registered or aligned to a complete volume-segmented rat atlas. The matrices that transformed the subject\u27s anatomy to the atlas space are used to embed each slice within the atlas. All transformed pixel locations of the anatomy images are tagged with the segmented atlas major and minor regions creating a fully segmented representation of each subject. This task required the development of a full 3D surface atlas based upon 2D non-uniformly spaced 2D slices from an existing atlas. A multiple materials marching cube (M3C) algorithm was used to generate these 1277 subvolumes. After this process, they were coalesced into a dozen major zones of the brain (amygdaloid complex, cerebrum, cerebellum, hypothalamus, etc.). Each major brain category was subdivided into approximately 10 sub-major zones. Many scientists are interested in behavior and reactions to pain, pleasure, smell, for example. Consequently, the 3D volume atlas was segmented into functional zones as well as the anatomical regions. A utility (program) called Tree Browser was developed to interactively display and choose different anatomical and/or functional areas. Statistical t-tests are performed to determine activation on each subject within their original coordinate system. Due to the multiple t-test analyses performed, a false-positive detection controlling mechanism was introduced. A statistical composite of five components was created for each group. The individual analyses were summed within groups. The strategy developed in this work is unique as it registers segments and analyzes multiple subjects and presents a composite response of the whole group. This strategy is robust, incredibly fast and statistically powerful. The power of this system was demonstrated by mapping the olfactory system of a rat brain. Synchronized changes in neuronal activity across multiple subjects and brain areas can be viewed as functional neuro-anatomical circuits coordinating the thoughts, memories and emotions for particular behaviors using this fMRI module

Application of bootstrap resampling in fMRI

This thesis demonstrates the use of the bootstrap resampling technique considering temporal dependency in the fMRI data to determine the reliability and confidence interval of fMRI parameters. Traditionally, the test-retest method has been used to reliably detect active voxels in the fMRI image of the brain, which is based on repetitive experimentation. The main concern with the test-retest method is the reproducibility of data over these multiple repetitions. Fatigue, habituation, motion artifacts, and repositioning errors are few of the factors, which can affect the reproducibility of data. The conventional bootstrap resampling technique is based on the assumption that the dataset is independent and identically distributed over time. However, studies have shown temporal dependency in the fMRI images of the brain acquired from subjects in the resting phase. This study demonstrates the use of the bootstrap resampling technique, incorporating the criterion of temporal dependency in the fMRI data set, to detect reliable active voxels in the fMRI images acquired during a task activated motor paradigm, where the subject is instructed to perform bilateral finger tapping. The results of the study showed that the active regions detected using the bootstrap resampling technique considering temporal dependency in the fMRI data were more reliable than the active regions detected using the bootstrap resampling technique without considering any temporal dependency in the fMRI data

Performance Comparison of 3D Sinc Interpolation for fMRI Motion Correction by Language of Implementation and Hardware Platform

Substantial effort is devoted to improving neuroimaging data processing; this effort however, is typically from the algorithmic perspective only. I demonstrate that substantive running time performance improvements to neuroscientific data processing algorithms can be realized by considering their implementation. Focusing specifically on 3D sinc interpolation, an algorithm used for processing functional magnetic resonance imaging (fMRI) data, I compare the performance of Python, C and OpenCL implementations of this algorithm across multiple hardware platforms. I also benchmark the performance of a novel implementation of 3D sinc interpolation on a field programmable gate array (FPGA). Together, these comparisons demonstrate that the performance of a neuroimaging data processing algorithm is significantly impacted by its implementation. I also present a case study demonstrating the practical benefits of improving a neuroscientific data processing algorithm\u27s implementation, then conclude by addressing threats to the validity of the study and discussing future directions

A Novel Synergistic Model Fusing Electroencephalography and Functional Magnetic Resonance Imaging for Modeling Brain Activities

Study of the human brain is an important and very active area of research. Unraveling the way the human brain works would allow us to better understand, predict and prevent brain related diseases that affect a significant part of the population. Studying the brain response to certain input stimuli can help us determine the involved brain areas and understand the mechanisms that characterize behavioral and psychological traits. In this research work two methods used for the monitoring of brain activities, Electroencephalography (EEG) and functional Magnetic Resonance (fMRI) have been studied for their fusion, in an attempt to bridge together the advantages of each one. In particular, this work has focused in the analysis of a specific type of EEG and fMRI recordings that are related to certain events and capture the brain response under specific experimental conditions. Using spatial features of the EEG we can describe the temporal evolution of the electrical field recorded in the scalp of the head. This work introduces the use of Hidden Markov Models (HMM) for modeling the EEG dynamics. This novel approach is applied for the discrimination of normal and progressive Mild Cognitive Impairment patients with significant results. EEG alone is not able to provide the spatial localization needed to uncover and understand the neural mechanisms and processes of the human brain. Functional Magnetic Resonance imaging (fMRI) provides the means of localizing functional activity, without though, providing the timing details of these activations. Although, at first glance it is apparent that the strengths of these two modalities, EEG and fMRI, complement each other, the fusion of information provided from each one is a challenging task. A novel methodology for fusing EEG spatiotemporal features and fMRI features, based on Canonical Partial Least Squares (CPLS) is presented in this work. A HMM modeling approach is used in order to derive a novel feature-based representation of the EEG signal that characterizes the topographic information of the EEG. We use the HMM model in order to project the EEG data in the Fisher score space and use the Fisher score to describe the dynamics of the EEG topography sequence. The correspondence between this new feature and the fMRI is studied using CPLS. This methodology is applied for extracting features for the classification of a visual task. The results indicate that the proposed methodology is able to capture task related activations that can be used for the classification of mental tasks. Extensions on the proposed models are examined along with future research directions and applications

Doctor of Philosophy

dissertationFunctional magnetic resonance imaging (fMRI) measures the change of oxygen consumption level in the blood vessels of the human brain, hence indirectly detecting the neuronal activity. Resting-state fMRI (rs-fMRI) is used to identify the intrinsic functional patterns of the brain when there is no external stimulus. Accurate estimation of intrinsic activity is important for understanding the functional organization and dynamics of the brain, as well as differences in the functional networks of patients with mental disorders. This dissertation aims to robustly estimate the functional connectivities and networks of the human brain using rs-fMRI data of multiple subjects. We use Markov random field (MRF), an undirected graphical model to represent the statistical dependency among the functional network variables. Graphical models describe multivariate probability distributions that can be factorized and represented by a graph. By defining the nodes and the edges along with their weights according to our assumptions, we build soft constraints into the graph structure as prior information. We explore various approximate optimization methods including variational Bayesian, graph cuts, and Markov chain Monte Carlo sampling (MCMC). We develop the random field models to solve three related problems. In the first problem, the goal is to detect the pairwise connectivity between gray matter voxels in a rs-fMRI dataset of the single subject. We define a six-dimensional graph to represent our prior information that two voxels are more likely to be connected if their spatial neighbors are connected. The posterior mean of the connectivity variables are estimated by variational inference, also known as mean field theory in statistical physics. The proposed method proves to outperform the standard spatial smoothing and is able to detect finer patterns of brain activity. Our second work aims to identify multiple functional systems. We define a Potts model, a special case of MRF, on the network label variables, and define von Mises-Fisher distribution on the normalized fMRI signal. The inference is significantly more difficult than the binary classification in the previous problem. We use MCMC to draw samples from the posterior distribution of network labels. In the third application, we extend the graphical model to the multiple subject scenario. By building a graph including the network labels of both a group map and the subject label maps, we define a hierarchical model that has richer structure than the flat single-subject model, and captures the shared patterns as well as the variation among the subjects. All three solutions are data-driven Bayesian methods, which estimate model parameters from the data. The experiments show that by the regularization of MRF, the functional network maps we estimate are more accurate and more consistent across multiple sessions

Magnetic Resonance Imaging of the Brain in Moving Subjects. Application of Fetal, Neonatal and Adult Brain Studies

Imaging in the presence of subject motion has been an ongoing challenge for magnetic resonance imaging (MRI). Motion makes MRI data inconsistent, causing artifacts in conventional anatomical imaging as well as invalidating diffusion tensor imaging (DTI) reconstruction. In this thesis some of the important issues regarding the acquisition and reconstruction of anatomical and DTI imaging of moving subjects are addressed; methods to achieve high resolution and high signalto- noise ratio (SNR) volume data are proposed. An approach has been developed that uses multiple overlapped dynamic single shot slice by slice imaging combined with retrospective alignment and data fusion to produce self consistent 3D volume images under subject motion. We term this method as snapshot MRI with volume reconstruction or SVR. The SVR method has been performed successfully for brain studies on subjects that cannot stay still, and in some cases were moving substantially during scanning. For example, awake neonates, deliberately moved adults and, especially, on fetuses, for which no conventional high resolution 3D method is currently available. Fine structure of the in-utero fetal brain is clearly revealed for the first time with substantially improved SNR. The SVR method has been extended to correct motion artifacts from conventional multi-slice sequences when the subject drifts in position during data acquisition. Besides anatomical imaging, the SVR method has also been further extended to DTI reconstruction when there is subject motion. This has been validated successfully from an adult who was deliberately moving and then applied to inutero fetal brain imaging, which no conventional high resolution 3D method is currently available. Excellent fetal brain 3D apparent diffusion coefficient (ADC) maps in high resolution have been achieved for the first time as well as promising fractional Anisotropy (FA) maps. Pilot clinical studies using SVR reconstructed data to study fetal brain development in-utero have been performed. Growth curves for the normally developing fetal brain have been devised by the quantification of cerebral and cerebellar volumes as well as some one dimensional measurements. A Verhulst model is proposed to describe these growth curves, and this approach has achieved a correlation over 0.99 between the fitted model and actual data