Andrew Melbourne

Registration of dynamic contrast-enhanced MRI using a progressive principal component registration (PPCR)

Registration of dynamic contrast-enhanced magnetic resonance images (DCE-MRI) of soft tissue is difficult. Conventional registration cost functions that depend on information content are compromised by the changing intensity profile, leading to misregistration. We present a new data-driven model of uptake patterns formed from a principal components analysis (PCA) of time-series data, avoiding the need for a physiological model. We term this process progressive principal component registration (PPCR). Registration is performed repeatedly to an artificial time series of target images generated using the principal components of the current best-registered time-series data. The aim is to produce a dataset that has had random motion artefacts removed but long-term contrast enhancement implicitly preserved. The procedure is tested on 22 DCE-MRI datasets of the liver. Preliminary assessment of the images is by expert observer comparison with registration to the first image in the sequence. The PPCR is preferred in all cases where a preference exists. The method requires neither segmentation nor a pharmacokinetic uptake model and can allow successful registration in the presence of contrast enhancement.

Geodesic Information Flows: Spatially-Variant Graphs and Their Application to Segmentation and Fusion

Clinical annotations, such as voxel-wise binary or probabilistic tissue segmentations, structural parcellations, pathological regions-of-interest and anatomical landmarks are key to many clinical studies. However, due to the time consuming nature of manually generating these annotations, they tend to be scarce and limited to small subsets of data. This work explores a novel framework to propagate voxel-wise annotations between morphologically dissimilar images by diffusing and mapping the available examples through intermediate steps. A spatially-variant graph structure connecting morphologically similar subjects is introduced over a database of images, enabling the gradual diffusion of information to all the subjects, even in the presence of large-scale morphological variability. We illustrate the utility of the proposed framework on two example applications: brain parcellation using categorical labels and tissue segmentation using probabilistic features. The application of the proposed method to categorical label fusion showed highly statistically significant improvements when compared to state-of-the-art methodologies. Significant improvements were also observed when applying the proposed framework to probabilistic tissue segmentation of both synthetic and real data, mainly in the presence of large morphological variability.

Respiratory motion correction in dynamic MRI using robust data decomposition registration - application to DCE-MRI.

Motion correction in Dynamic Contrast Enhanced (DCE-) MRI is challenging because rapid intensity changes can compromise common (intensity based) registration algorithms. In this study we introduce a novel registration technique based on robust principal component analysis (RPCA) to decompose a given time-series into a low rank and a sparse component. This allows robust separation of motion components that can be registered, from intensity variations that are left unchanged. This Robust Data Decomposition Registration (RDDR) is demonstrated on both simulated and a wide range of clinical data. Robustness to different types of motion and breathing choices during acquisition is demonstrated for a variety of imaged organs including liver, small bowel and prostate. The analysis of clinically relevant regions of interest showed both a decrease of error (15-62% reduction following registration) in tissue time-intensity curves and improved areas under the curve (AUC60) at early enhancement.

AdaPT: An adaptive preterm segmentation algorithm for neonatal brain MRI

Advances in neonatal care have improved the survival of infants born prematurely although these infants remain at increased risk of adverse neurodevelopmental outcome. The measurement of white matter structure and features of the cortical surface can help define biomarkers that predict this risk. The measurement of these structures relies upon accurate automated segmentation routines, but these are often confounded by neonatal-specific imaging difficulties including poor contrast, low resolution, partial volume effects and the presence of significant natural and pathological anatomical variability. In this work we develop and evaluate an adaptive preterm multi-modal maximum a posteriori expectation-maximisation segmentation algorithm (AdaPT) incorporating an iterative relaxation strategy that adapts the tissue proportion priors toward the subject data. Also incorporated are intensity non-uniformity correction, a spatial homogeneity term in the form of a Markov random field and furthermore, the proposed method explicitly models the partial volume effect specifically mitigating the neonatal specific grey and white matter contrast inversion. Spatial priors are iteratively relaxed, enabling the segmentation of images with high anatomical disparity from a normal population. Experiments performed on a clinical cohort of 92 infants are validated against manual segmentation of normal and pathological cortical grey matter, cerebellum and ventricular volumes. Dice overlap scores increase significantly when compared to a widely-used maximum likelihood expectation maximisation algorithm for pathological cortical grey matter, cerebellum and ventricular volumes. Adaptive maximum a posteriori expectation maximisation is shown to be a useful tool for accurate and robust neonatal brain segmentation.

Evaluation of automatic neonatal brain segmentation algorithms:the NeoBrainS12 challenge

A number of algorithms for brain segmentation in preterm born infants have been published, but a reliable comparison of their performance is lacking. The NeoBrainS12 study (http://neobrains12.isi.uu.nl), providing three different image sets of preterm born infants, was set up to provide such a comparison. These sets are (i) axial scans acquired at 40 weeks corrected age, (ii) coronal scans acquired at 30 weeks corrected age and (iii) coronal scans acquired at 40 weeks corrected age. Each of these three sets consists of three T1- and T2-weighted MR images of the brain acquired with a 3T MRI scanner. The task was to segment cortical grey matter, non-myelinated and myelinated white matter, brainstem, basal ganglia and thalami, cerebellum, and cerebrospinal fluid in the ventricles and in the extracerebral space separately. Any team could upload the results and all segmentations were evaluated in the same way. This paper presents the results of eight participating teams. The results demonstrate that the participating methods were able to segment all tissue classes well, except myelinated white matter.

Preterm birth affects the developmental synergy between cortical folding and cortical connectivity observed on multimodal MRI.

The survival rates of infants born prematurely have improved as a result of advances in neonatal care, although there remains an increased risk of subsequent disability. Accurate measurement of the shape and appearance of the very preterm brain at term-equivalent age may guide the development of predictive biomarkers of neurological outcome. We demonstrate in 92 preterm infants (born at an average gestational age of 27.0±2.7weeks) scanned at term equivalent age (scanned at 40.4±1.74weeks) that the cortical sulcation ratio varies spatially over the cortical surface at term equivalent age and correlates significantly with gestational age at birth (r=0.49,p<0.0001). In the underlying white matter, fractional anisotropy of local white matter regions correlated significantly with gestational age at birth at term equivalent age (for the genu of the corpus callosum r=0.26,p=0.02 and for the splenium r=0.52,p<0.001) and in addition the fractional anisotropy in these local regions varies according to location. Finally, we demonstrate that connectivity measurements from tractography correlate significantly and specifically with the sulcation ratio of the overlying cortical surface at term equivalent age in a subgroup of 20 infants (r={0.67,0.61,0.86}, p={0.004,0.01,0.00002}) for tract systems emanating from the left and right corticospinal tracts and the corpus callosum respectively). Combined, these results suggest a close relationship between the cortical surface phenotype and underlying white matter structure assessed by diffusion weighted MRI. The spatial surface pattern may allow inference on the connectivity and developmental trajectory of the underlying white matter complementary to diffusion imaging and this result may guide the development of biomarkers of functional outcome.

Longitudinal measurement of the developing grey matter in preterm subjects using multi-modal MRI.

Preterm birth is a major public health concern, with the severity and occurrence of adverse outcome increasing with earlier delivery. Being born preterm disrupts a time of rapid brain development: in addition to volumetric growth, the cortex folds, myelination is occurring and there are changes on the cellular level. These neurological events have been imaged non-invasively using diffusion-weighted (DW) MRI. In this population, there has been a focus on examining diffusion in the white matter, but the grey matter is also critically important for neurological health. We acquired multi-shell high-resolution diffusion data on 12 infants born at ≤ 28 weeks of gestational age at two time-points: once when stable after birth, and again at term-equivalent age. We used the Neurite Orientation Dispersion and Density Imaging model (NODDI) (Zhang et al., 2012) to analyse the changes in the cerebral cortex and the thalamus, both grey matter regions. We showed region-dependent changes in NODDI parameters over the preterm period, highlighting underlying changes specific to the microstructure. This work is the first time that NODDI parameters have been evaluated in both the cortical and the thalamic grey matter as a function of age in preterm infants, offering a unique insight into neuro-development in this at-risk population.

Motion artifact correction in free‐breathing abdominal MRI using overlapping partial samples to recover image deformations

This article presents a method to reconstruct liver MRI data acquired continuously during free breathing, without any external sensor or navigator measurements. When the deformations associated with k‐space data are known, generalized matrix inversion reconstruction has been shown to be effective in reducing the ghosting and blurring artifacts of motion. This article describes a novel method to obtain these nonrigid deformations. A breathing model is built from a fast dynamic series: low spatial resolution images are registered and their deformations parameterized by overall superior–inferior displacement. The correct deformation for each subset of the subsequent imaging data is then found by comparing a few lines of k‐space with the equivalent lines from a deformed reference image while varying the deformation over the model parameter. This procedure is known as image deformation recovery using overlapping partial samples (iDROPS). Simulations using 10 rapid dynamic studies from volunteers showed the average error in iDROPS‐derived deformations within the liver to be 1.43 mm. A further four volunteers were imaged at higher spatial resolution. The complete reconstruction process using data from throughout several breathing cycles was shown to reduce blurring and ghosting in the liver. Retrospective respiratory gating was also demonstrated using the iDROPS parameterization. Magn Reson Med, 2009.

The effect of motion correction on pharmacokinetic parameter estimation in dynamic-contrast-enhanced MRI

A dynamic-contrast-enhanced magnetic resonance imaging (DCE-MRI) dataset consists of many imaging frames, often acquired both before and after contrast injection. Due to the length of time spent acquiring images, patient motion is likely and image re-alignment or registration is required before further analysis such as pharmacokinetic model fitting. Non-rigid image registration procedures may be used to correct motion artefacts; however, a careful choice of registration strategy is required to reduce misregistration artefacts associated with enhancing features. This work investigates the effect of registration on the results of model-fitting algorithms for 52 DCE-MR mammography cases for 14 patients. Results are divided into two sections: a comparison of registration strategies in which a DCE-MRI-specific algorithm is preferred in 50% of cases, followed by an investigation of parameter changes with known applied deformations, inspecting the effect of magnitude and timing of motion artefacts. Increased motion magnitude correlates with increased model-fit residual and is seen to have a strong influence on the visibility of strongly enhancing features. Motion artefacts in images close to the contrast agent arrival have a disproportionate effect on discrepancies in parameter estimation. The choice of algorithm, magnitude of motion and timing of the motion are each shown to influence estimated pharmacokinetic parameters even when motion magnitude is small.

Adaptive neonate brain segmentation

Babies born prematurely are at increased risk of adverse neurodevelopmental outcomes. Recent advances suggest that measurement of brain volumes can help in defining biomarkers for neurodevelopmental outcome. These techniques rely on an accurate segmentation of the MRI data. However, due to lack of contrast, partial volume (PV) effect, the existence of both hypo- and hyper-intensities and significant natural and pathological anatomical variability, the segmentation of neonatal brain MRI is challenging. We propose a pipeline for image segmentation that uses a novel multi-model Maximum a posteriori Expectation Maximisation (MAP-EM) segmentation algorithm with a prior over both intensities and the tissue proportions, a B0 inhomogeneity correction, and a spatial homogeneity term through the use of a Markov Random Field. This robust and adaptive technique enables the segmentation of images with high anatomical disparity from a normal population. Furthermore, the proposed method implicitly models Partial Volume, mitigating the problem of neonatal white/grey matter intensity inversion. Experiments performed on a clinical cohort show expected statistically significant correlations with gestational age at birth and birthweight. Furthermore, the proposed method obtains statistically significant improvements in Dice scores when compared to the a Maximum Likelihood EM algorithm.