Search

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Down Syndrome is a chromosomal disorder that affects the development of cerebellar cortical lobules. Impaired neurogenesis in the cerebellum varies among different types of neuronal cells and neuronal layers. In this study, we developed an imaging analysis framework that utilizes gadolinium-enhanced ex vivo mouse brain MRI. We extracted the middle Purkinje layer of the mouse cerebellar cortex, enabling the estimation of the volume, thickness, and surface area of the entire cerebellar cortex, the internal granular layer, and the molecular layer in the Tc1 mouse model of Down Syndrome. The morphometric analysis of our method revealed that a larger proportion of the cerebellar thinning in this model of Down Syndrome resided in the inner granule cell layer, while a larger proportion of the surface area shrinkage was in the molecular layer., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Spinal cord atrophy measurements obtained from structural magnetic resonance imaging (MRI) are associated with disability in many neurological diseases and serve as in vivo biomarkers of neurodegeneration. Longitudinal spinal cord atrophy rate is commonly determined from the numerical difference between two volumes (based on 3D surface fitting) or two cross-sectional areas (CSA, based on 2D edge detection) obtained at different time-points. Being an indirect measure, atrophy rates are susceptible to variable segmentation errors at the edge of the spinal cord. To overcome those limitations, we developed a new registration-based pipeline that measures atrophy rates directly. We based our approach on the generalised boundary shift integral (GBSI) method, which registers 2 scans and uses a probabilistic XOR mask over the edge of the spinal cord, thereby measuring atrophy more accurately than segmentation-based techniques. Using a large cohort of longitudinal spinal cord images (610 subjects with multiple sclerosis from a multi-centre trial and 52 healthy controls), we demonstrated that GBSI is a sensitive, quantitative and objective measure of longitudinal spinal cord volume change. The GBSI pipeline is repeatable, reproducible, and provides more precise measurements of longitudinal spinal cord atrophy than segmentation-based methods in longitudinal spinal cord atrophy studies., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

In this paper we present a hierarchical generative model of medical image data, which can capture simultaneously the variability of both signal intensity and anatomical shapes across large populations. Such a model has a direct application for learning average-shaped probabilistic tissue templates in a fully automated manner. While in principle the generality of the proposed Bayesian approach makes it suitable to address a wide range of medical image computing problems, our work focuses primarily on neuroimaging applications. In particular we validate the proposed method on both real and synthetic brain MR scans including the cervical cord and demonstrate that it yields accurate alignment of brain and spinal cord structures, as compared to state-of-the-art tools for medical image registration. At the same time we illustrate how the resulting tissue probability maps can readily be used to segment, bias correct and spatially normalise unseen data, which are all crucial pre-processing steps for MR imaging studies., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Hospitals often hold historical MR image data printed on films without being able to make it accessible to modern image processing techniques. Having the possibility to recover geometrically consistent, volumetric images from scans acquired decades ago will enable more comprehensive, longitudinal studies to understand disease progressions. In this paper, we propose a consistent framework to reconstruct a volumetric representation from printed films holding thick single-slice brain MR image acquisitions dating back to the 1980's. We introduce a flexible framework based on semi-automatic slice extraction, followed by automated slice-to-volume registration with inter-slice transformation regularisation and slice intensity correction. Our algorithm is robust against numerous detrimental effects being present in archaic films. A subsequent, isotropic total variation deconvolution technique revitalises the visual appearance of the obtained volumes. We assess the accuracy and perform the validation of our reconstruction framework on a uniquely long-term MRI dataset where a ground-truth is available. This method will be used to facilitate a robust longitudinal analysis spanning 30 years of MRI scans., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Aim: To accurately quantify the radioactivity concentration measured by PET, emission data need to be corrected for photon attenuation; however, the MRI signal cannot easily be converted into attenuation values, making attenuation correction (AC) in PET/MRI challenging. In order to further improve the current vendor-implemented MR-AC methods for absolute quantification, a number of prototype methods have been proposed in the literature. These can be categorized into three types: template/atlas-based, segmentation-based, and reconstruction-based. These proposed methods in general demonstrated improvements compared to vendor-implemented AC, and many studies report deviations in PET uptake after AC of only a few percent from a gold standard CT-AC. Using a unified quantitative evaluation with identical metrics, subject cohort, and common CT-based reference, the aims of this study were to evaluate a selection of novel methods proposed in the literature, and identify the ones suitable for clinical use., Methods: In total, 11 AC methods were evaluated: two vendor-implemented (MR-AC DIXON and MR-AC UTE ), five based on template/atlas information (MR-AC SEGBONE (Koesters et al., 2016), MR-AC ONTARIO (Anazodo et al., 2014), MR-AC BOSTON (Izquierdo-Garcia et al., 2014), MR-AC UCL (Burgos et al., 2014), and MR-AC MAXPROB (Merida et al., 2015)), one based on simultaneous reconstruction of attenuation and emission (MR-AC MLAA (Benoit et al., 2015)), and three based on image-segmentation (MR-AC MUNICH (Cabello et al., 2015), MR-AC CAR-RiDR (Juttukonda et al., 2015), and MR-AC RESOLUTE (Ladefoged et al., 2015)). We selected 359 subjects who were scanned using one of the following radiotracers: [ 18 F]FDG (210), [ 11 C]PiB (51), and [ 18 F]florbetapir (98). The comparison to AC with a gold standard CT was performed both globally and regionally, with a special focus on robustness and outlier analysis., Results: The average performance in PET tracer uptake was within ±5% of CT for all of the proposed methods, with the average±SD global percentage bias in PET FDG uptake for each method being: MR-AC DIXON (-11.3±3.5)%, MR-AC UTE (-5.7±2.0)%, MR-AC ONTARIO (-4.3±3.6)%, MR-AC MUNICH (3.7±2.1)%, MR-AC MLAA (-1.9±2.6)%, MR-AC SEGBONE (-1.7±3.6)%, MR-AC UCL (0.8±1.2)%, MR-AC CAR-RiDR (-0.4±1.9)%, MR-AC MAXPROB (-0.4±1.6)%, MR-AC BOSTON (-0.3±1.8)%, and MR-AC RESOLUTE (0.3±1.7)%, ordered by average bias. The overall best performing methods (MR-AC BOSTON , MR-AC MAXPROB , MR-AC RESOLUTE and MR-AC UCL , ordered alphabetically) showed regional average errors within ±3% of PET with CT-AC in all regions of the brain with FDG, and the same four methods, as well as MR-AC CAR-RiDR , showed that for 95% of the patients, 95% of brain voxels had an uptake that deviated by less than 15% from the reference. Comparable performance was obtained with PiB and florbetapir., Conclusions: All of the proposed novel methods have an average global performance within likely acceptable limits (±5% of CT-based reference), and the main difference among the methods was found in the robustness, outlier analysis, and clinical feasibility. Overall, the best performing methods were MR-ACBOSTON, MR-ACMAXPROB, MR-ACRESOLUTE and MR-ACUCL, ordered alphabetically. These methods all minimized the number of outliers, standard deviation, and average global and local error. The methods MR-ACMUNICH and MR-ACCAR-RiDR were both within acceptable quantitative limits, so these methods should be considered if processing time is a factor. The method MR-ACSEGBONE also demonstrates promising results, and performs well within the likely acceptable quantitative limits. For clinical routine scans where processing time can be a key factor, this vendor-provided solution currently outperforms most methods. With the performance of the methods presented here, it may be concluded that the challenge of improving the accuracy of MR-AC in adult brains with normal anatomy has been solved to a quantitatively acceptable degree, which is smaller than the quantification reproducibility in PET imaging., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Multiple sclerosis lesions influence the process of image analysis, leading to tissue segmentation problems and biased morphometric estimates. Existing techniques try to reduce this bias by filling all lesions as normal-appearing white matter on T1-weighted images, considering each time-point separately. However, due to lesion segmentation errors and the presence of structures adjacent to the lesions, such as the ventricles and deep grey matter nuclei, filling all lesions with white matter-like intensities introduces errors and artefacts. In this paper, we present a novel lesion filling strategy inspired by in-painting techniques used in computer graphics applications for image completion. The proposed technique uses a five-dimensional (5D), patch-based (multi-modality and multi-time-point), Non-Local Means algorithm that fills lesions with the most plausible texture. We demonstrate that this strategy introduces less bias, fewer artefacts and spurious edges than the current, publicly available techniques. The proposed method is modality-agnostic and can be applied to multiple time-points simultaneously. In addition, it preserves anatomical structures and signal-to-noise characteristics even when the lesions are neighbouring grey matter or cerebrospinal fluid, and avoids excess of blurring or rasterisation due to the choice of the segmentation plane, shape of the lesions, and their size and/or location., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Structural MRI is widely used for investigating brain atrophy in many neurodegenerative disorders, with several research groups developing and publishing techniques to provide quantitative assessments of this longitudinal change. Often techniques are compared through computation of required sample size estimates for future clinical trials. However interpretation of such comparisons is rendered complex because, despite using the same publicly available cohorts, the various techniques have been assessed with different data exclusions and different statistical analysis models. We created the MIRIAD atrophy challenge in order to test various capabilities of atrophy measurement techniques. The data consisted of 69 subjects (46 Alzheimer's disease, 23 control) who were scanned multiple (up to twelve) times at nine visits over a follow-up period of one to two years, resulting in 708 total image sets. Nine participating groups from 6 countries completed the challenge by providing volumetric measurements of key structures (whole brain, lateral ventricle, left and right hippocampi) for each dataset and atrophy measurements of these structures for each time point pair (both forward and backward) of a given subject. From these results, we formally compared techniques using exactly the same dataset. First, we assessed the repeatability of each technique using rates obtained from short intervals where no measurable atrophy is expected. For those measures that provided direct measures of atrophy between pairs of images, we also assessed symmetry and transitivity. Then, we performed a statistical analysis in a consistent manner using linear mixed effect models. The models, one for repeated measures of volume made at multiple time-points and a second for repeated "direct" measures of change in brain volume, appropriately allowed for the correlation between measures made on the same subject and were shown to fit the data well. From these models, we obtained estimates of the distribution of atrophy rates in the Alzheimer's disease (AD) and control groups and of required sample sizes to detect a 25% treatment effect, in relation to healthy ageing, with 95% significance and 80% power over follow-up periods of 6, 12, and 24months. Uncertainty in these estimates, and head-to-head comparisons between techniques, were carried out using the bootstrap. The lateral ventricles provided the most stable measurements, followed by the brain. The hippocampi had much more variability across participants, likely because of differences in segmentation protocol and less distinct boundaries. Most methods showed no indication of bias based on the short-term interval results, and direct measures provided good consistency in terms of symmetry and transitivity. The resulting annualized rates of change derived from the model ranged from, for whole brain: -1.4% to -2.2% (AD) and -0.35% to -0.67% (control), for ventricles: 4.6% to 10.2% (AD) and 1.2% to 3.4% (control), and for hippocampi: -1.5% to -7.0% (AD) and -0.4% to -1.4% (control). There were large and statistically significant differences in the sample size requirements between many of the techniques. The lowest sample sizes for each of these structures, for a trial with a 12month follow-up period, were 242 (95% CI: 154 to 422) for whole brain, 168 (95% CI: 112 to 282) for ventricles, 190 (95% CI: 146 to 268) for left hippocampi, and 158 (95% CI: 116 to 228) for right hippocampi. This analysis represents one of the most extensive statistical comparisons of a large number of different atrophy measurement techniques from around the globe. The challenge data will remain online and publicly available so that other groups can assess their methods., (Copyright © 2015. Published by Elsevier Inc.)

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., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

As the number of people diagnosed with Alzheimer's disease (AD) reaches epidemic proportions, there is an urgent need to develop effective treatment strategies to tackle the social and economic costs of this fatal condition. Dozens of candidate therapeutics are currently being tested in clinical trials, and compounds targeting the aberrant accumulation of tau proteins into neurofibrillary tangles (NFTs) are the focus of substantial current interest. Reliable, translatable biomarkers sensitive to both tau pathology and its modulation by treatment along with animal models that faithfully reflect aspects of the human disease are urgently required. Magnetic resonance imaging (MRI) is well established as a valuable tool for monitoring the structural brain changes that accompany AD progression. However the descent into dementia is not defined by macroscopic brain matter loss alone: non-invasive imaging measurements sensitive to protein accumulation, white matter integrity and cerebral haemodynamics probe distinct aspects of AD pathophysiology and may serve as superior biomarkers for assessing drug efficacy. Here we employ a multi-parametric array of five translatable MRI techniques to characterise the in vivo pathophysiological phenotype of the rTg4510 mouse model of tauopathy (structural imaging, diffusion tensor imaging (DTI), arterial spin labelling (ASL), chemical exchange saturation transfer (CEST) and glucose CEST). Tau-induced pathological changes included grey matter atrophy, increased radial diffusivity in the white matter, decreased amide proton transfer and hyperperfusion. We demonstrate that the above markers unambiguously discriminate between the transgenic group and age-matched controls and provide a comprehensive profile of the multifaceted neuropathological processes underlying the rTg4510 model. Furthermore, we show that ASL and DTI techniques offer heightened sensitivity to processes believed to precede detectable structural changes and, as such, provides a platform for the study of disease mechanisms and therapeutic intervention., (Copyright © 2015 Elsevier Inc. All rights reserved.)

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., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

We introduce a novel image-processing framework for tracking longitudinal changes in white matter microstructure using diffusion tensor imaging (DTI). Charting the trajectory of such temporal changes offers new insight into disease progression but to do so accurately faces a number of challenges. Recent developments have highlighted the importance of processing each subject's data at multiple time points in an unbiased way. In this paper, we aim to highlight a different challenge critical to the processing of longitudinal DTI data, namely the approach to image alignment. Standard approaches in the literature align DTI data by registering the corresponding scalar-valued fractional anisotropy (FA) maps. We propose instead a DTI registration algorithm that leverages full tensor information to drive improved alignment. This proposed pipeline is evaluated against the standard FA-based approach using a DTI dataset from an ongoing study of Alzheimer's disease (AD). The dataset consists of subjects scanned at two time points and at each time point the DTI acquisition consists of two back-to-back repeats in the same scanning session. The repeated scans allow us to evaluate the specificity of each pipeline, using a test-retest design, and assess precision, using bootstrap-based method. The results show that the tensor-based pipeline achieves both higher specificity and precision than the standard FA-based approach. Tensor-based registration for longitudinal processing of DTI data in clinical studies may be of particular value in studies assessing disease progression., (Copyright © 2013 Elsevier Inc. All rights reserved.)

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., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Cortical thickness estimation performed in-vivo via magnetic resonance imaging is an important technique for the diagnosis and understanding of the progression of neurodegenerative diseases. Currently, two different computational paradigms exist, with methods generally classified as either surface or voxel-based. This paper provides a much needed comparison of the surface-based method FreeSurfer and two voxel-based methods using clinical data. We test the effects of computing regional statistics using two different atlases and demonstrate that this makes a significant difference to the cortical thickness results. We assess reproducibility, and show that FreeSurfer has a regional standard deviation of thickness difference on same day scans that is significantly lower than either a Laplacian or Registration based method and discuss the trade off between reproducibility and segmentation accuracy caused by bending energy constraints. We demonstrate that voxel-based methods can detect similar patterns of group-wise differences as well as FreeSurfer in typical applications such as producing group-wise maps of statistically significant thickness change, but that regional statistics can vary between methods. We use a Support Vector Machine to classify patients against controls and did not find statistically significantly different results with voxel based methods compared to FreeSurfer. Finally we assessed longitudinal performance and concluded that currently FreeSurfer provides the most plausible measure of change over time, with further work required for voxel based methods., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Thickness measurements of the cerebral cortex can aid diagnosis and provide valuable information about the temporal evolution of diseases such as Alzheimer's, Huntington's, and schizophrenia. Methods that measure the thickness of the cerebral cortex from in-vivo magnetic resonance (MR) images rely on an accurate segmentation of the MR data. However, segmenting the cortex in a robust and accurate way still poses a challenge due to the presence of noise, intensity non-uniformity, partial volume effects, the limited resolution of MRI and the highly convoluted shape of the cortical folds. Beginning with a well-established probabilistic segmentation model with anatomical tissue priors, we propose three post-processing refinements: a novel modification of the prior information to reduce segmentation bias; introduction of explicit partial volume classes; and a locally varying MRF-based model for enhancement of sulci and gyri. Experiments performed on a new digital phantom, on BrainWeb data and on data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) show statistically significant improvements in Dice scores and PV estimation (p<10(-3)) and also increased thickness estimation accuracy when compared to three well established techniques., (Copyright © 2011 Elsevier Inc. All rights reserved.)