Sebastien Ourselin

Fast free-form deformation using graphics processing units

A large number of algorithms have been developed to perform non-rigid registration and it is a tool commonly used in medical image analysis. The free-form deformation algorithm is a well-established technique, but is extremely time consuming. In this paper we present a parallel-friendly formulation of the algorithm suitable for graphics processing unit execution. Using our approach we perform registration of T1-weighted MR images in less than 1 min and show the same level of accuracy as a classical serial implementation when performing segmentation propagation. This technology could be of significant utility in time-critical applications such as image-guided interventions, or in the processing of large data sets.

Reconstructing a 3D structure from serial histological sections

Abstract We consider the problem of aligning histological sections for 3D reconstruction and analysis. The method we propose is based on a block-matching strategy that allows us to compute local displacements between the sections. We then collect these local measures to estimate a rigid transformation. Our emphasis is on the necessity to use a robust approach for this estimation step. The process is integrated within a multi-scale scheme to improve both accuracy and computation time. We prove experimentally that we can reach sub-pixel accuracy and we show some results of aligning histological sections from a rats brain and a rhesus monkeys brain.

Biopsy Site Re-localisation Based on the Computation of Epipolar Lines from Two Previous Endoscopic Images

Tracking biopsy sites in endoscopic images can be useful to provide a visual aid for the guidance of surgical tools, for example when endoscopic guided biopsy is required. A new method for re-localisation of these sites is presented in this paper. It makes use of epipolar geometry properties between three images of the same site observed from different viewpoints with an endoscope. Two epipolar lines are derived from the two first images in the third image where the site needs to be re-localised. Their intersection corresponds to the location of the biopsy site. This method was tested with gastroscopic data from 2 patients with 9 series of three images of the oesophagus. The re-localisation error was estimated at less than 1.5 millimetres by a clinical endoscopist, which is sufficient for most clinical endoscopic applications.

Evaluation of Registration Methods on Thoracic CT: The EMPIRE10 Challenge

EMPIRE10 (Evaluation of Methods for Pulmonary Image REgistration 2010) is a public platform for fair and meaningful comparison of registration algorithms which are applied to a database of intra patient thoracic CT image pairs. Evaluation of nonrigid registration techniques is a nontrivial task. This is compounded by the fact that researchers typically test only on their own data, which varies widely. For this reason, reliable assessment and comparison of different registration algorithms has been virtually impossible in the past. In this work we present the results of the launch phase of EMPIRE10, which comprised the comprehensive evaluation and comparison of 20 individual algorithms from leading academic and industrial research groups. All algorithms are applied to the same set of 30 thoracic CT pairs. Algorithm settings and parameters are chosen by researchers expert in the con figuration of their own method and the evaluation is independent, using the same criteria for all participants. All results are published on the EMPIRE10 website (http://empire10.isi.uu.nl). The challenge remains ongoing and open to new participants. Full results from 24 algorithms have been published at the time of writing. This paper details the organization of the challenge, the data and evaluation methods and the outcome of the initial launch with 20 algorithms. The gain in knowledge and future work are discussed.

A three-dimensional, histological and deformable atlas of the human basal ganglia. I. Atlas construction based on immunohistochemical and MRI data.

This paper describes the construction of an atlas of the human basal ganglia. The successive steps of the construction were as follows. First a postmortem specimen was subjected to a MRI acquisition prior to extraction of the brain from the skull. The brain was then cryosectioned (70 microm thickness). One section out of ten (80 sections) was Nissl-stained with cresyl violet, another series of 80 sections was immunostained for the calcium binding protein calbindin. Contours of basal ganglia nuclei including their calbindin-stained functional subdivisions, fiber bundles and ventricles (n=80 structures) were traced from histological sections and digitized. A novelty of this atlas is the MRI acquisition, which represents the core data element of the study. MRI was used for the coregistration of the atlas data and permitted, through multimodal (Nissl, calbindin, images of cryosectioning, T1 and T2 MRI) and 3D optimization, the production of anatomically and geometrically consistent 3D surfaces, which can be sliced through any desired orientation. The atlas MRI is also used for its deformation to provide accurate conformation to the MRI of living patients, thus adding information at the histological level to the patients MRI volume. This latter aspect will be presented in a forthcoming paper.

Block Matching: A General Framework to Improve Robustness of Rigid Registration of Medical Images

In order to improve the robustness of rigid registration algorithms in various medical imaging problems, we propose in this article a general framework built on block matching strategies. This framework combines two stages in a multi-scale hierarchy. The first stage consists in finding for each block (or subregion) of the first image, the most similar subregion in the other image, using a similarity criterion which depends on the nature of the images. The second stage consists in finding the global rigid transformation which best explains most of these local correspondances. This is done with a robust procedure which allows up to 50% of false matches. We show that this approach, besides its simplicity, provides a robust and efficient way to rigidly register images in various situations. This includes for instance the alignment of 2D histological sections for the 3D reconstructions of trimmed organs and tissues, the automatic computation of the mid-sagittal plane in multimodal 3D images of the brain, and the multimodal registration of 3D CT and MR images of the brain. A quantitative evaluation of the results is provided for this last example, as well as a comparison with the classical approaches involving the minimization of a global measure of similarity based on Mutual Information or the Correlation Ratio. This shows a significant improvement of the robustness, for a comparable final accuracy. Although slightly more expensive in terms of computational requirements, the proposed approach can easily be implemented on a parallel architecture, which opens potentialities for real time applications using a large number of processors.

Head size, age and gender adjustment in MRI studies: a necessary nuisance?

Imaging studies of cerebral volumes often adjust for factors such as age that may confound between-subject comparisons. However the use of nuisance covariates in imaging studies is inconsistent, which can make interpreting results across studies difficult. Using magnetic resonance images of 78 healthy controls we assessed the effects of age, gender, head size and scanner upgrade on region of interest (ROI) volumetry, cortical thickness and voxel-based morphometric (VBM) measures. We found numerous significant associations between these variables and volumetric measures: cerebral volumes and cortical thicknesses decreased with increasing age, men had larger volumes and smaller thicknesses than women, and increasing head size was associated with larger volumes. The relationships between most ROIs and head size volumes were non-linear. With age, gender, head size and upgrade in one model we found that volumes and thicknesses decreased with increasing age, women had larger volumes than men (VBM, whole-brain and white matter volumes), increasing head size was associated with larger volumes but not cortical thickness, and scanner upgrade had an effect on thickness and some volume measures. The effects of gender on cortical thickness when adjusting for head size, age and upgrade showed some non-significant effect (women>men), whereas the independent effect of head size showed little pattern. We conclude that age and head size should be considered in ROI volume studies, age, gender and upgrade should be considered for cortical thickness studies and all variables require consideration for VBM analyses. Division of all volumes by head size is unlikely to be adequate owing to their non-proportional relationship.

Automated cross-sectional and longitudinal hippocampal volume measurement in mild cognitive impairment and Alzheimer's disease

Volume and change in volume of the hippocampus are both important markers of Alzheimers disease (AD). Delineation of the structure on MRI is time-consuming and therefore reliable automated methods are required. We describe an improvement (multiple-atlas propagation and segmentation (MAPS)) to our template library-based segmentation technique. The improved technique uses non-linear registration of the best-matched templates from our manually segmented library to generate multiple segmentations and combines them using the simultaneous truth and performance level estimation (STAPLE) algorithm. Change in volume over 12months (MAPS-HBSI) was measured by applying the boundary shift integral using MAPS regions. Methods were developed and validated against manual measures using subsets from Alzheimers Disease Neuroimaging Initiative (ADNI). The best method was applied to 682 ADNI subjects, at baseline and 12-month follow-up, enabling assessment of volumes and atrophy rates in control, mild cognitive impairment (MCI) and AD groups, and within MCI subgroups classified by subsequent clinical outcome. We compared our measures with those generated by Surgical Navigation Technologies (SNT) available from ADNI. The accuracy of our volumes was one of the highest reported (mean(SD) Jaccard Index 0.80(0.04) (N=30)). Both MAPS baseline volume and MAPS-HBSI atrophy rate distinguished between control, MCI and AD groups. Comparing MCI subgroups (reverters, stable and converters): volumes were lower and rates higher in converters compared with stable and reverter groups (p< or =0.03). MAPS-HBSI required the lowest sample sizes (78 subjects) for a hypothetical trial. In conclusion, the MAPS and MAPS-HBSI methods give accurate and reliable volumes and atrophy rates across the clinical spectrum from healthy aging to AD.

β-Amyloid burden in the temporal neocortex is related to hippocampal atrophy in elderly subjects without dementia

Objective: To investigate whether global and regional β-amyloid (Aβ) burden as measured with 11C Pittsburgh compound B (PIB) PET is associated with hippocampal atrophy characterized using MRI in healthy controls and patients with amnestic mild cognitive impairment (aMCI) or Alzheimer disease (AD). Methods: Ninety-two elderly healthy controls, 32 subjects with aMCI, and 35 patients with AD were imaged using 11C-PIB PET and MRI. Hippocampal volume was measured and PIB standardized uptake value ratio was extracted after partial volume correction within 41 regions of interest. Global, regional, and voxel-based correlations between PIB and hippocampal volume were computed for each group. Results: In healthy control participants with elevated neocortex PIB retention, significant correlation was found between PIB retention in the inferior temporal region and hippocampal volume using both region-based and voxel-based approaches. No correlation was found in any other group. Conclusions: The strong correlation between hippocampal atrophy and β-amyloid (Aβ) burden in the Pittsburgh compound B–positive healthy control group suggests that Aβ deposition in the inferior temporal neocortex is related to hippocampal synaptic and neuronal degeneration.

Clinical and neuroanatomical signatures of tissue pathology in frontotemporal lobar degeneration

Relating clinical symptoms to neuroanatomical profiles of brain damage and ultimately to tissue pathology is a key challenge in the field of neurodegenerative disease and particularly relevant to the heterogeneous disorders that comprise the frontotemporal lobar degeneration spectrum. Here we present a retrospective analysis of clinical, neuropsychological and neuroimaging (volumetric and voxel-based morphometric) features in a pathologically ascertained cohort of 95 cases of frontotemporal lobar degeneration classified according to contemporary neuropathological criteria. Forty-eight cases (51%) had TDP-43 pathology, 42 (44%) had tau pathology and five (5%) had fused-in-sarcoma pathology. Certain relatively specific clinicopathological associations were identified. Semantic dementia was predominantly associated with TDP-43 type C pathology; frontotemporal dementia and motoneuron disease with TDP-43 type B pathology; young-onset behavioural variant frontotemporal dementia with FUS pathology; and the progressive supranuclear palsy syndrome with progressive supranuclear palsy pathology. Progressive non-fluent aphasia was most commonly associated with tau pathology. However, the most common clinical syndrome (behavioural variant frontotemporal dementia) was pathologically heterogeneous; while pathologically proven Picks disease and corticobasal degeneration were clinically heterogeneous, and TDP-43 type A pathology was associated with similar clinical features in cases with and without progranulin mutations. Volumetric magnetic resonance imaging, voxel-based morphometry and cluster analyses of the pathological groups here suggested a neuroanatomical framework underpinning this clinical and pathological diversity. Frontotemporal lobar degeneration-associated pathologies segregated based on their cerebral atrophy profiles, according to the following scheme: asymmetric, relatively localized (predominantly temporal lobe) atrophy (TDP-43 type C); relatively symmetric, relatively localized (predominantly temporal lobe) atrophy (microtubule-associated protein tau mutations); strongly asymmetric, distributed atrophy (Picks disease); relatively symmetric, predominantly extratemporal atrophy (corticobasal degeneration, fused-in-sarcoma pathology). TDP-43 type A pathology was associated with substantial individual variation; however, within this group progranulin mutations were associated with strongly asymmetric, distributed hemispheric atrophy. We interpret the findings in terms of emerging network models of neurodegenerative disease: the neuroanatomical specificity of particular frontotemporal lobar degeneration pathologies may depend on an interaction of disease-specific and network-specific factors.