Development and characterization of techniques for neuro-imaging registration

Three automated techniques were developed for the alignment of Neuro-Images acquired during distinct scanning periods and their performance were characterized. The techniques are based on the assumption that the human brain is a rigid body and will assume different positions during different scanning periods. One technique uses three fiducial markers, while the other two uses eigenvectors of the inertia matrix of the Neuro-Image, to compute the three angles (pitch, yaw and roll) needed to register the test Neuro-Image to the reference Neuro-Image. A rigid body transformation is computed and applied to the test Neuro-Image such that it results aligned to the reference Neuro-Image. These techniques were tested by applying known rigid body transformations to given Neuro-Images. The transformations were retrieved automatically on the basis of unit vectors or eigenvectors. The results show that the precision of two techniques is dependent on the axial resolution of the Neuro-Images and for one of them also on the imaging modality, while the precision of one technique is also dependent on the interpolation. Such methods can be applied to any Neuro-Imaging modality and have been tested for both fMRI and MRI

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

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

Techniques for Analysis and Motion Correction of Arterial Spin Labelling (ASL) Data from Dementia Group Studies

This investigation examines how Arterial Spin Labelling (ASL) Magnetic Resonance Imaging can be optimised to assist in the early diagnosis of diseases which cause dementia, by considering group study analysis and control of motion artefacts. ASL can produce quantitative cerebral blood flow maps noninvasively - without a radioactive or paramagnetic contrast agent being injected. ASL studies have already shown perfusion changes which correlate with the metabolic changes measured by Positron Emission Tomography in the early stages of dementia, before structural changes are evident. But the clinical use of ASL for dementia diagnosis is not yet widespread, due to a combination of a lack of protocol consistency, lack of accepted biomarkers, and sensitivity to motion artefacts. Applying ASL to improve early diagnosis of dementia may allow emerging treatments to be administered earlier, thus with greater effect. In this project, ASL data acquired from two separate patient cohorts ( (i) Young Onset Alzheimer’s Disease (YOAD) study, acquired at Queen Square; and (ii) Incidence and RISk of dementia (IRIS) study, acquired in Rotterdam) were analysed using a pipeline optimised for each acquisition protocol, with several statistical approaches considered including support-vector machine learning. Machine learning was also applied to improve the compatibility of the two studies, and to demonstrate a novel method to disentangle perfusion changes measured by ASL from grey matter atrophy. Also in this project, retrospective motion correction techniques for specific ASL sequences were developed, based on autofocusing and exploiting parallel imaging algorithms. These were tested using a specially developed simulation of the 3D GRASE ASL protocol, which is capable of modelling motion. The parallel imaging based approach was verified by performing a specifically designed MRI experiment involving deliberate motion, then applying the algorithm to demonstrably reduce motion artefacts retrospectively

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

Quality data assessment and improvement in pre-processing pipeline to minimize impact of spurious signals in functional magnetic imaging (fMRI)

In the recent years, the field of quality data assessment and signal denoising in functional magnetic resonance imaging (fMRI) is rapidly evolving and the identification and reduction of spurious signal with pre-processing pipeline is one of the most discussed topic. In particular, subject motion or physiological signals, such as respiratory or/and cardiac pulsatility, were showed to introduce false-positive activations in subsequent statistical analyses. Different measures for the evaluation of the impact of motion related artefacts, such as frame-wise displacement and root mean square of movement parameters, and the reduction of these artefacts with different approaches, such as linear regression of nuisance signals and scrubbing or censoring procedure, were introduced. However, we identify two main drawbacks: i) the different measures used for the evaluation of motion artefacts were based on user-dependent thresholds, and ii) each study described and applied their own pre-processing pipeline. Few studies analysed the effect of these different pipelines on subsequent analyses methods in task-based fMRI.The first aim of the study is to obtain a tool for motion fMRI data assessment, based on auto-calibrated procedures, to detect outlier subjects and outliers volumes, targeted on each investigated sample to ensure homogeneity of data for motion. The second aim is to compare the impact of different pre-processing pipelines on task-based fMRI using GLM based on recent advances in resting state fMRI preprocessing pipelines. Different output measures based on signal variability and task strength were used for the assessment

Processing strategies for functional magnetic resonance imaging data sets

Thesis (Ph.D.)--Harvard--Massachusetts Institute of Technology Division of Health Sciences and Technology, 1999.Includes bibliographical references (leaves 108-118).by Luis Carlos Maas, III.Ph.D

Real-time motion correction for Magnetic Resonance Imaging of the human brain at 7 Tesla

Magnetic resonance imaging (MRI) at the magnetic field strength of 7 Tesla (7T) enhances the quality of images available for research and clinical use. The improvements are however accompanied by novel challenges that are specific to ultra high-field MRI, which includes field strengths of 7T and above. Transmit B +1 field inhomogeneity is also higher, causing uneven signal intensity and linking to an uneven SAR distribution, which is also higher than at lower field strengths. The potential for higher spatial resolution imaging can also result in more pronounced motion artefacts. To address these issues in routine clinical use, motion correction strategies are required. This thesis describe the implementation of real-time, image-based Multislice Prospective Acquisition Correction (MS-PACE) technique for 7T MRI. Firstly, developmental work was done to establish a 7T-specific MS-PACE implementation. Pulse sequence and image reconstruction pipeline work was implemented using the Siemens Integrated Development Environment for Applications (IDEA) and Image Calculation Environment (ICE) framework. The technique was then validated in a task-based functional MRI study with healthy subjects. It was also integrated with parallel transmit imaging using slice-by-slice B+1 shimming. Validation experiments were performed in vivo using the Siemens MAGNETOM Terra 7T MRI scanner (Siemens Healthineers, Erlangen, Germany) at the Imaging Centre of Excellence (ICE)

Using an Internal Auditory Stimulus to Activate the Developing Primary Auditory Cortex: A Fetal fMRI Study

Insight into the rapidly developing brain in utero is scarce. Fetal functional magnetic resonance imaging (fMRI) is a technique used to gain awareness into the developmental process. Previous auditory task-based fMRI studies employed an external sound stimulus directly on the maternal abdomen. However, there has since been recommendation to cease doing so. We sought to investigate a reliable paradigm to study the development of fetal brain networks and postulate that by using an internal stimulus, such as the mother singing, it would result in activation of the fetal primary auditory cortex. Volunteers carrying singleton fetuses with a gestational age of 33-38 weeks underwent two stimulus-based block design BOLD fMRI series. All of the nine fetal subjects analyzed had activation in the right Heschl’s gyrus, and seven out of the nine fetal subjects had activation in the left Heschl’s gyrus when exposed to the internal acoustic stimulus. Ultimately, this internal auditory stimulus can be used to analyze the developing fetal brain

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