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

Improved fMRI Time-Series Registration Using Joint Probability Density Priors

Functional MRI (fMRI) time-series studies are plagued by varying degrees of subject head motion. Faithful head motion correction is essential to accurately detect brain activation using statistical analyses of these time-series. Mutual information (MI) based slice-to-volume (SV) registration is used for motion estimation when the rate of change of head position is large. SV registration accounts for head motion between slice acquisitions by estimating an independent rigid transformation for each slice in the time-series. Consequently each MI optimization uses intensity counts from a single time-series slice, making the algorithm susceptible to noise for low complexity endslices (i.e., slices near the top of the head scans). This work focuses on improving the accuracy of MI-based SV registration of end-slices by using joint probability density priors derived from registered high complexity centerslices (i.e., slices near the middle of the head scans). Results show that the use of such priors can significantly improve SV registration accuracy.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85928/1/Fessler236.pd

Motion correction in fMRI via registration of individual slices into an anatomical volume

An automated retrospective image registration based on mutual information is adapted to a multislice functional magnetic resonance imaging (fMRI) acquisition protocol to provide accurate motion correction. Motion correction is performed by mapping each slice to an anatomic volume data set acquired in the same fMRI session to accommodate inter-slice head motion. Accuracy of the registration parameters was assessed by registration of simulated MR data of the known truth. The widely used rigid body volume registration approach based on stacked slices from the time series data may hinder statistical accuracy by introducing inaccurate assumptions of no motion between slices for multislice fMRI data. Improved sensitivity and specificity of the fMRI signal from mapping-each-slice-to-volume method is demonstrated in comparison with a stacked-slice correction method by examining functional data from two normal volunteers. The data presented in a standard anatomical coordinate system suggest the reliability of the mapping-each-slice-to-volume method to detect the activation signals consistent between the two subjects. Magn Reson Med 41:964–972, 1999. © 1999 Wiley-Liss, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/34927/1/16_ftp.pd

Neural Representations of Visual Motion Processing in the Human Brain Using Laminar Imaging at 9.4 Tesla

During natural behavior, much of the motion signal falling into our eyes is due to our own movements. Therefore, in order to correctly perceive motion in our environment, it is important to parse visual motion signals into those caused by self-motion such as eye- or head-movements and those caused by external motion. Neural mechanisms underlying this task, which are also required to allow for a stable perception of the world during pursuit eye movements, are not fully understood. Both, perceptual stability as well as perception of real-world (i.e. objective) motion are the product of integration between motion signals on the retina and efference copies of eye movements. The central aim of this thesis is to examine whether different levels of cortical depth or distinct columnar structures of visual motion regions are differentially involved in disentangling signals related to self-motion, objective, or object motion. Based on previous studies reporting segregated populations of voxels in high level visual areas such as V3A, V6, and MST responding predominantly to either retinal or extra- retinal (‘real’) motion, we speculated such voxels to reside within laminar or columnar functional units. We used ultra-high field (9.4T) fMRI along with an experimental paradigm that independently manipulated retinal and extra-retinal motion signals (smooth pursuit) while controlling for effects of eye-movements, to investigate whether processing of real world motion in human V5/MT, putative MST (pMST), and V1 is associated to differential laminar signal intensities. We also examined motion integration across cortical depths in human motion areas V3A and V6 that have strong objective motion responses. We found a unique, condition specific laminar profile in human area V6, showing reduced mid-layer responses for retinal motion only, suggestive of an inhibitory retinal contribution to motion integration in mid layers or alternatively an excitatory contribution in deep and superficial layers. We also found evidence indicating that in V5/MT and pMST, processing related to retinal, objective, and pursuit motion are either integrated or colocalized at the scale of our resolution. In contrast, in V1, independent functional processes seem to be driving the response to retinal and objective motion on the one hand, and to pursuit signals on the other. The lack of differential signals across depth in these regions suggests either that a columnar rather than laminar segregation governs these functions in these areas, or that the methods used were unable to detect differential neural laminar processing. Furthermore, the thesis provides a thorough analysis of the relevant technical modalities used for data acquisition and data analysis at ultra-high field in the context of laminar fMRI. Relying on our technical implementations we were able to conduct two high-resolution fMRI experiments that helped us to further investigate the laminar organization of self-induced and externally induced motion cues in human high-level visual areas and to form speculations about the site and the mechanisms of their integration

Development and characterization of methodology and technology for the alignment of fMRI time series

This dissertation has developed, implemented and tested a novel computer based system (AUTOALIGN) that incorporates an algorithm for the alignment of functional Magnetic Resonance Image (fMRI) time series. The algorithm assumes the human brain to be a rigid body and computes a head coordinate system on the basis of three reference points that lie on the directions correspondent to two of the eigenvectors of inertia of the volume, at the intersections with the head boundary. The eigenvectors are found weighting the inertia components with the voxel\u27s intensity values assumed as mass. The three reference points are found in the same position, relative to the origin of the head coordinate system, in both test and reference brain images. Intensity correction is performed at sub-voxel accuracy by tri-linear interpolation. A test fMR brain volume in which controlled simulations of rigid-body transformations have been introduced has preliminarily assessed system performance. Further experimentation has been conducted with real fMRI time series. Rigid-body transformations have been retrieved automatically and the values of the motion parameters compared to those obtained by the Statistical Parametric Mapping (SPM99), and the Automatic Image Registration (AIR 3.08). Results indicated that AUTOALIGN offers subvoxel accuracy in correcting both misalignment and intensity among time points in fMR images time series, and also that its performance is comparable to that of SPM99 and AIR3.08

A blind deconvolution approach to recover effective connectivity brain networks from resting state fMRI data

A great improvement to the insight on brain function that we can get from fMRI data can come from effective connectivity analysis, in which the flow of information between even remote brain regions is inferred by the parameters of a predictive dynamical model. As opposed to biologically inspired models, some techniques as Granger causality (GC) are purely data-driven and rely on statistical prediction and temporal precedence. While powerful and widely applicable, this approach could suffer from two main limitations when applied to BOLD fMRI data: confounding effect of hemodynamic response function (HRF) and conditioning to a large number of variables in presence of short time series. For task-related fMRI, neural population dynamics can be captured by modeling signal dynamics with explicit exogenous inputs; for resting-state fMRI on the other hand, the absence of explicit inputs makes this task more difficult, unless relying on some specific prior physiological hypothesis. In order to overcome these issues and to allow a more general approach, here we present a simple and novel blind-deconvolution technique for BOLD-fMRI signal. Coming to the second limitation, a fully multivariate conditioning with short and noisy data leads to computational problems due to overfitting. Furthermore, conceptual issues arise in presence of redundancy. We thus apply partial conditioning to a limited subset of variables in the framework of information theory, as recently proposed. Mixing these two improvements we compare the differences between BOLD and deconvolved BOLD level effective networks and draw some conclusions

Synthesizing pseudo-T2w images to recapture missing data in neonatal neuroimaging with applications in rs-fMRI

T1- and T2-weighted (T1w and T2w) images are essential for tissue classification and anatomical localization in Magnetic Resonance Imaging (MRI) analyses. However, these anatomical data can be challenging to acquire in non-sedated neonatal cohorts, which are prone to high amplitude movement and display lower tissue contrast than adults. As a result, one of these modalities may be missing or of such poor quality that they cannot be used for accurate image processing, resulting in subject loss. While recent literature attempts to overcome these issues in adult populations using synthetic imaging approaches, evaluation of the efficacy of these methods in pediatric populations and the impact of these techniques in conventional MR analyses has not been performed. In this work, we present two novel methods to generate pseudo-T2w images: the first is based in deep learning and expands upon previous models to 3D imaging without the requirement of paired data, the second is based in nonlinear multi-atlas registration providing a computationally lightweight alternative. We demonstrate the anatomical accuracy of pseudo-T2w images and their efficacy in existing MR processing pipelines in two independent neonatal cohorts. Critically, we show that implementing these pseudo-T2w methods in resting-state functional MRI analyses produces virtually identical functional connectivity results when compared to those resulting from T2w images, confirming their utility in infant MRI studies for salvaging otherwise lost subject data

Analysis and Strategies to Enhance Intensity-Base Image Registration.

The availability of numerous complementary imaging modalities allows us to obtain a detailed picture of the body and its functioning. To aid diagnostics and surgical planning, all available information can be presented by visually aligning images from different modalities using image registration. This dissertation investigates strategies to improve the performance of image registration algorithms that use intensity-based similarity metrics. Nonrigid warp estimation using intensity-based registration can be very time consuming. We develop a novel framework based on importance sampling and stochastic approximation techniques to accelerate nonrigid registration methods while preserving their accuracy. Registration results for simulated brain MRI data and human lung CT data demonstrate the efficacy of the proposed framework. Functional MRI (fMRI) is used to non-invasively detect brain-activation by acquiring a series of brain images, called a time-series, while the subject performs tasks designed to stimulate parts of the brain. Consequently, these studies are plagued by subject head motion. Mutual information (MI) based slice-to-volume (SV) registration algorithms used to estimate time-series motion are less accurate for end-slices (i.e., slices near the top of the head scans), where a loss in image complexity yields noisy MI estimates. We present a strategy, dubbed SV-JP, to improve SV registration accuracy for time-series end-slices by using joint pdf priors derived from successfully registered high complexity slices near the middle of the head scans to bolster noisy MI estimates. Although fMRI time-series registration can estimate head motion, this motion also spawns extraneous intensity fluctuations called spin saturation artifacts. These artifacts hamper brain-activation detection. We describe spin saturation using mathematical expressions and develop a weighted-average spin saturation (WASS) correction scheme. An algorithm to identify time-series voxels affected by spin saturation and to implement WASS correction is outlined. The performance of registration methods is dependant on the tuning parameters used to implement their similarity metrics. To facilitate finding optimal tuning parameters, we develop a computationally efficient linear approximation of the (co)variance of MI-based registration estimates. However, empirically, our approximation was satisfactory only for a simple mono-modality registration example and broke down for realistic multi-modality registration where the MI metric becomes strongly nonlinear.Ph.D.Electrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/61552/1/rbhagali_1.pd

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