Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Phase Image Texture Analysis for Motion Detection in Diffusion MRI (PITA-MDD)

Purpose Pronounced spin phase artifacts appear in diffusion-weighted imaging (DWI) with only minor subject motion. While DWI data corruption is often identified as signal drop out in diffusion-weighted (DW) magnitude images, DW phase images may have higher sensitivity for detecting subtle subject motion. Methods This article describes a novel method to return a metric of subject motion, computed using an image texture analysis of the DW phase image. This Phase Image Texture Analysis for Motion Detection in dMRI (PITA-MDD) method is computationally fast and reliably detects subject motion from diffusion-weighted images. A threshold of the motion metric was identified to remove motion-corrupted slices, and the effect of removing corrupted slices was assessed on the reconstructed FA maps and fiber tracts. Results Using a motion-metric threshold to remove the motion-corrupted slices results in superior fiber tracts and fractional anisotropy maps. When further compared to a state-of-the-art magnitude-based motion correction method, PITA-MDD was able to detect comparable corrupted slices in a more computationally efficient manner. Conclusion In this study, we evaluated the use of DW phase images to detect motion corruption. The proposed method can be a robust and fast alternative for automatic motion detection in the brain with multiple applications to inform prospective motion correction or as real-time feedback for data quality control during scanning, as well as after data is already acquired

High-resolution diffusion-weighted brain MRI under motion

Magnetic resonance imaging is one of the fastest developing medical imaging techniques. It provides excellent soft tissue contrast and has been a leading tool for neuroradiology and neuroscience research over the last decades. One of the possible MR imaging contrasts is the ability to visualize diffusion processes. The method, referred to as diffusion-weighted imaging, is one of the most common clinical contrasts but is prone to artifacts and is challenging to acquire at high resolutions. This thesis aimed to improve the resolution of diffusion weighted imaging, both in a clinical and in a research context. While diffusion-weighted imaging traditionally has been considered a 2D technique the manuscripts and methods presented here explore 3D diffusion acquisitions with isotropic resolution. Acquiring multiple small 3D volumes, or slabs, which are combined into one full volume has been the method of choice in this work. The first paper presented explores a parallel imaging driven multi-echo EPI readout to enable high resolution with reduced geometric distortions. The work performed on diffusion phase correction lead to an understanding that was used for the subsequent multi-slab papers. The second and third papers introduce the diffusion-weighted 3D multi-slab echo-planar imaging technique and explore its advantages and performance. As the method requires a slightly increased acquisition time the need for prospective motion correction became apparent. The forth paper suggests a new motion navigator using the subcutaneous fat surrounding the skull for rigid body head motion estimation, dubbed FatNav. The spatially sparse representation of the fat signal allowed for high parallel imaging acceleration factors, short acquisition times, and reduced geometric distortions of the navigator. The fifth manuscript presents a combination of the high-resolution 3D multi-slab technique and a modified FatNav module. Unlike our first FatNav implementation, using a single sagittal slab, this modified navigator acquired orthogonal projections of the head using the fat signal alone. The combined use of both presented methods provides a promising start for a fully motion corrected high-resolution diffusion acquisition in a clinical setting

Fast volume reconstruction from motion corrupted stacks of 2D slices

Capturing an enclosing volume of moving subjects and organs using fast individual image slice acquisition has shown promise in dealing with motion artefacts. Motion between slice acquisitions results in spatial inconsistencies that can be resolved by slice-to-volume reconstruction (SVR) methods to provide high quality 3D image data. Existing algorithms are, however, typically very slow, specialised to specific applications and rely on approximations, which impedes their potential clinical use. In this paper, we present a fast multi-GPU accelerated framework for slice-to-volume reconstruction. It is based on optimised 2D/3D registration, super-resolution with automatic outlier rejection and an additional (optional) intensity bias correction. We introduce a novel and fully automatic procedure for selecting the image stack with least motion to serve as an initial registration target. We evaluate the proposed method using artificial motion corrupted phantom data as well as clinical data, including tracked freehand ultrasound of the liver and fetal Magnetic Resonance Imaging. We achieve speed-up factors greater than 30 compared to a single CPU system and greater than 10 compared to currently available state-of-the-art multi-core CPU methods. We ensure high reconstruction accuracy by exact computation of the point-spread function for every input data point, which has not previously been possible due to computational limitations. Our framework and its implementation is scalable for available computational infrastructures and tests show a speed-up factor of 1.70 for each additional GPU. This paves the way for the online application of image based reconstruction methods during clinical examinations. The source code for the proposed approach is publicly available

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