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

AM-FM Analysis of Structural and Functional Magnetic Resonance Images

This thesis proposes the application of multi-dimensional Amplitude-Modulation Frequency-Modulation (AM-FM) methods to magnetic resonance images (MRI). The basic goal is to provide a framework for exploring non-stationary characteristics of structural and functional MRI (sMRI and fMRI). First, we provide a comparison framework for the most popular AM-FM methods using different filterbank configurations that includes Gabor, Equirriple and multi-scale directional designs. We compare the performance and robustness to Gaussian noise using synthetic FM image examples. We show that the multi-dimensional quasi-local method (QLM) with an equiripple filterbank gave the best results in terms of instantaneous frequency (IF) estimation. We then apply the best performing AM-FM method to sMRI to compute the 3D IF features. We use a t-test on the IF magnitude for each voxel to find evidence of significant differences between healthy controls and patients diagnosed with schizophrenia (n=353) can be found in the IF. We also propose the use of the instantaneous phase (IP) as a new feature for analyzing fMRI images. Using principal component analysis and independent component analysis on the instantaneous phase from fMRI, we built spatial maps and identified brain regions that are biologically coherent with the task performed by the subject. This thesis provides the first application of AM-FM models to fMRI and sMRI

Improving the accuracy of brain activation maps in the group-level analysis of fMRI data utilizing spatiotemporal Gaussian process model

OBJECTIVE: Accuracy and precision of the statistical analysis methods used for brain activation maps are essential. Adjusting models to consider spatiotemporal correlation embedded in fMRI data may increase their accuracy, but it also introduces a high computational cost. The present study aimed to apply and assess the spatiotemporal Gaussian process (STGP) model to improve accuracy and reduce cost. METHODS: We applied the spatiotemporal Gaussian process (STGP) model for both simulated and experimental memory tfMRI data and compared the findings with fast, fully Bayesian, and General Linear Models (GLM). To assess their accuracy and precision, the models were fitted to the simulated data (1000 voxels,100 times point for 50 people), and an average of accuracy indexes of 100 repetitions was computed. Functional and activation maps for all models were calculated in experimental data analysis. RESULTS: STGP model resulted in a higher Z-score in the whole brain, in the 1000 most activated voxels, and in the frontal lobe as the approved memory area. Based on the simulated data, the STGP model showed more accuracy and precision than the other two models. However, its computational time was more than the GLM, as the price of model correction, but much less than that of the fast, fully Bayesian model. CONCLUSION: Spatiotemporal correlation further improved the accuracy of the STGP compared to the GLM and fast, fully Bayesian model. This can result in more accurate activation maps. Moreover, the STGP model’s computational speed appears to be reasonable for model application

Estimating the effective degrees of freedom in univariate multiple regression analysis

The general linear model provides the most widely applied statistical framework for analyzing functional MRI (fMRI) data. With the increasing temporal resolution of recent scanning protocols, and more elaborate data preprocessing schemes, data independency is no longer a valid assumption. In this paper, we revise the statistical background of the general linear model in the presence of temporal autocorrelations. First, when detecting the activation signal, we explicitly account for the temporal autocorrelation structure, which yields a generalized F-test and the associated corrected (or effective) degrees of freedom (DOF). The proposed approach is data driven and thus independent of any specific preprocessing method. Then, for event-related protocols, we propose a new model for the temporal autocorrelations (“damped oscillator” model) and compare this model to another, previously used in the field (first-order autoregressive model, or AR(1) model). In the case of long fMRI time series, an efficient approximation for the number of effective DOF is provided for both models. Finally, the validity of our approach is assessed using simulated and real fMRI data and is compared with more conventional methods

Fusion analysis of functional MRI data for classification of individuals based on patterns of activation.

Classification of individuals based on patterns of brain activity observed in functional MRI contrasts may be helpful for diagnosis of neurological disorders. Prior work for classification based on these patterns have primarily focused on using a single contrast, which does not take advantage of complementary information that may be available in multiple contrasts. Where multiple contrasts are used, the objective has been only to identify the joint, distinct brain activity patterns that differ between groups of subjects; not to use the information to classify individuals. Here, we use joint Independent Component Analysis (jICA) within a Support Vector Machine (SVM) classification method, and take advantage of the relative contribution of activation patterns generated from multiple fMRI contrasts to improve classification accuracy. Young (age: 19-26) and older (age: 57-73) adults (16 each) were scanned while listening to noise alone and to speech degraded with noise, half of which contained meaningful context that could be used to enhance intelligibility. Functional contrasts based on these conditions (and a silent baseline condition) were used within jICA to generate spatially independent joint activation sources and their corresponding modulation profiles. Modulation profiles were used within a non-linear SVM framework to classify individuals as young or older. Results demonstrate that a combination of activation maps across the multiple contrasts yielded an area under ROC curve of 0.86, superior to classification resulting from individual contrasts. Moreover, class separability, measured by a divergence criterion, was substantially higher when using the combination of activation maps