Search

Multi-Atlas based Segmentation (MAS) algorithms have been successfully applied to many medical image segmentation tasks, but their success relies on a large number of atlases and good image registration performance. Choosing well-registered atlases for label fusion is vital for an accurate segmentation. This choice becomes even more crucial when the segmentation involves organs characterized by a high anatomical and pathological variability. In this paper, we propose a new genetic atlas selection strategy (GAS) that automatically chooses the best subset of atlases to be used for segmenting the target image, on the basis of both image similarity and segmentation overlap. More precisely, the key idea of GAS is that if two images are similar, the performances of an atlas for segmenting each image are similar. Since the ground truth of each atlas is known, GAS first selects a predefined number of similar images to the target, then, for each one of them, finds a near-optimal subset of atlases by means of a genetic algorithm. All these near-optimal subsets are then combined and used to segment the target image. GAS was tested on single-label and multi-label segmentation problems. In the first case, we considered the segmentation of both the whole prostate and of the left ventricle of the heart from magnetic resonance images. Regarding multi-label problems, the zonal segmentation of the prostate into peripheral and transition zone was considered. The results showed that the performance of MAS algorithms statistically improved when GAS is used., (Copyright © 2018 Elsevier B.V. All rights reserved.)

This paper introduces a novel method for inferring spatially varying regularisation in non-linear registration. This is achieved through full Bayesian inference on a probabilistic registration model, where the prior on the transformation parameters is parameterised as a weighted mixture of spatially localised components. Such an approach has the advantage of allowing the registration to be more flexibly driven by the data than a traditional globally defined regularisation penalty, such as bending energy. The proposed method adaptively determines the influence of the prior in a local region. The strength of the prior may be reduced in areas where the data better support deformations, or can enforce a stronger constraint in less informative areas. Consequently, the use of such a spatially adaptive prior may reduce unwanted impacts of regularisation on the inferred transformation. This is especially important for applications where the deformation field itself is of interest, such as tensor based morphometry. The proposed approach is demonstrated using synthetic images, and with application to tensor based morphometry analysis of subjects with Alzheimer's disease and healthy controls. The results indicate that using the proposed spatially adaptive prior leads to sparser deformations, which provide better localisation of regional volume change. Additionally, the proposed regularisation model leads to more data driven and localised maps of registration uncertainty. This paper also demonstrates for the first time the use of Bayesian model comparison for selecting different types of regularisation., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

We present a framework for intrinsic comparison of surface metric structures and curvatures. This work parallels the work of Kurtek et al. on parameterization-invariant comparison of genus zero shapes. Here, instead of comparing the embedding of spherically parameterized surfaces in space, we focus on the first fundamental form. To ensure that the distance on spherical metric tensor fields is invariant to parameterization, we apply the conjugation-invariant metric arising from the L2 norm on symmetric positive definite matrices. As a reparameterization changes the metric tensor by a congruent Jacobian transform, this metric perfectly suits our purpose. The result is an intrinsic comparison of shape metric structure that does not depend on the specifics of a spherical mapping. Further, when restricted to tensors of fixed volume form, the manifold of metric tensor fields and its quotient of the group of unitary diffeomorphisms becomes a proper metric manifold that is geodesically complete. Exploiting this fact, and augmenting the metric with analogous metrics on curvatures, we derive a complete Riemannian framework for shape comparison and reconstruction. A by-product of our framework is a near-isometric and curvature-preserving mapping between surfaces. The correspondence is optimized using the fast spherical fluid algorithm. We validate our framework using several subcortical boundary surface models from the ADNI dataset.

The advent of large of multi-modal imaging databases opens up the opportunity to learn how local intensity patterns covariate between multiple modalities. These models can then be used to describe expected intensities in an unseen image modalities given one or multiple observations, or to detect deviations (e.g. pathology) from the expected intensity patterns. In this work, we propose a template-based multi-modal generative mixture-model of imaging data and apply it to the problems of inlier/outlier pattern classification and image synthesis. Results on synthetic and patient data demonstrate that the proposed method is able to synthesise unseen data and accurately localise pathological regions, even in the presence of large abnormalities. It also demonstrates that the proposed model can provide accurate and uncertainty-aware intensity estimates of expected imaging patterns.

Magnetic Resonance Imaging (MRI), a reference examination for cardiac morphology and function in humans, allows to image the cardiac right ventricle (RV) with high spatial resolution. The segmentation of the RV is a difficult task due to the variable shape of the RV and its ill-defined borders in these images. The aim of this paper is to evaluate several RV segmentation algorithms on common data. More precisely, we report here the results of the Right Ventricle Segmentation Challenge (RVSC), concretized during the MICCAI'12 Conference with an on-site competition. Seven automated and semi-automated methods have been considered, along them three atlas-based methods, two prior based methods, and two prior-free, image-driven methods that make use of cardiac motion. The obtained contours were compared against a manual tracing by an expert cardiac radiologist, taken as a reference, using Dice metric and Hausdorff distance. We herein describe the cardiac data composed of 48 patients, the evaluation protocol and the results. Best results show that an average 80% Dice accuracy and a 1cm Hausdorff distance can be expected from semi-automated algorithms for this challenging task on the datasets, and that an automated algorithm can reach similar performance, at the expense of a high computational burden. Data are now publicly available and the website remains open for new submissions (http://www.litislab.eu/rvsc/)., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Multiple Sclerosis lesions influence the process of image analysis, leading to tissue segmentation problems and biased morphometric estimates. With the aim of reducing this bias, existing techniques fill segmented lesions as normal appearing white matter. However, due to lesion segmentation errors or the presence of neighbouring structures, such as the ventricles and deep grey matter structures, filling all lesions as white matter like intensities is prone to introduce errors and artefacts. In this paper, we present a novel lesion filling strategy based on in-painting techniques for image completion. This technique makes use of a patch-based Non-Local Means algorithm that fills the lesions with the most plausible texture, rather than normal appearing white matter. We demonstrate that this strategy introduces less bias and fewer artefacts and spurious edges than previous techniques. The advantages of the proposed methodology are that it preserves both anatomical structure and signal-to-noise characteristics even when the lesions are neighbouring grey matter and cerebrospinal fluid, and avoids excess blurring or rasterisation due to the choice of segmentation plane, and lesion shape, size and/or position.

StereoEEG implantation is performed in patients with epilepsy to determine the site of the seizure onset zone. Intracranial haemorrhage is the most common complication associated to implantation carrying a risk that ranges from 0.6 to 2.7%, with significant associated morbidity. SEEG planning is done pre-operatively to identify avascular trajectories for the electrodes. In current practice neurosurgeons have no assistance in the planning of the electrode trajectories. There is great interest in developing computer assisted planning systems that can optimize the safety profile of electrode trajectories, maximizing the distance to critical brain structures. In this work, we address the problem of blood vessel extraction for SEEG trajectory planning. The proposed method exploits the availability of multi-modal images within a trajectory planning system to formulate a vessel extraction framework that combines the scale and the neighbouring structure of an object. We validated the proposed method in twelve multi-modal patient image sets. The mean Dice similarity coefficient (DSC) was 0.88 ± 0.03, representing a statistically significantly improvement when compared to the semi-automated single rater, single modality segmentation protocol used in current practice (DSC = 0.78 ± 0.02).

Preterm birth is a significant public health concern. For infants born very preterm (≤ 32 weeks completed gestation), there is a high instance of developmental disability. Due to the heterogeneity of patient outcomes, it is important to investigate early markers of future ability to provide effective and targeted intervention. As a neuronal relay centre, the thalamus is critical for effective cognitive function and, thus, development of white matter connections between the thalamus and cortex is vital. By non-invasively examining the state of the thalamus we can monitor development in the preterm period. To track the development we develop a novel registration technique to combine data from multiple modalities, in order to derive the transformation from a preterm scan, to a scan of the same infant at term-equivalent age. By measuring the changes in diffusion parameters over this period on a per-voxel basis, we hope to provide unique insight into neurodevelopment.

The detection of abnormal intensities in brain images caused by the presence of pathologies is currently under great scrutiny. Selecting appropriate models for pathological data is of critical importance for an unbiased and biologically plausible model fit, which in itself enables a better understanding of the underlying data and biological processes. Besides, it impacts on one's ability to extract pathologically meaningful imaging biomarkers. With this aim in mind, this work proposes a fully unsupervised hierarchical model selection framework for neuroimaging data which permits the stratification of different types of abnormal image atterns without prior knowledge about the subject's pathological status.

Multi-atlas segmentation has been widely used to segment various anatomical structures. The success of this technique partly relies on the selection of atlases that are best mapped to a new target image after registration. Recently, manifold learning has been proposed as a method for atlas selection. Each manifold learning technique seeks to optimize a unique objective function. Therefore, different techniques produce different embeddings even when applied to the same data set. Previous studies used a single technique in their method and gave no reason for the choice of the manifold learning technique employed nor the theoretical grounds for the choice of the manifold parameters. In this study, we compare side-by-side the results given by 3 manifold learning techniques (Isomap, Laplacian Eigenmaps and Locally Linear Embedding) on the same data set. We assess the ability of those 3 different techniques to select the best atlases to combine in the framework of multi-atlas segmentation. First, a leave-one-out experiment is used to optimize our method on a set of 110 manually segmented atlases of hippocampi and find the manifold learning technique and associated manifold parameters that give the best segmentation accuracy. Then, the optimal parameters are used to automatically segment 30 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI). For our dataset, the selection of atlases with Locally Linear Embedding gives the best results. Our findings show that selection of atlases with manifold learning leads to segmentation accuracy close to or significantly higher than the state-of-the-art method and that accuracy can be increased by fine tuning the manifold learning process.

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.)

This paper introduces a novel method for inferring spatially varying regularisation in non-rigid registration. This is achieved through full Bayesian inference on a probabilistic registration model, where the prior on transformations is parametrised as a weighted mixture of spatially localised components. Such an approach has the advantage of allowing the registration to be more flexibly driven by the data than a more traditional global regularisation scheme, such as bending energy. The proposed method adaptively determines the influence of the prior in a local region. The importance of the prior may be reduced in areas where the data better supports deformations, or can enforce a stronger constraint in less informative areas. Consequently, the use of such a spatially adaptive prior may reduce the unwanted impact of regularisation on the inferred deformation field. This is especially important for applications such as tensor based morphometry, where the features of interest are directly derived from the deformation field. The proposed approach is demonstrated with application to tensor based morphometry analysis of subjects with Alzheimer's disease and healthy controls. The results show that using the proposed spatially adaptive prior leads to deformation fields that have a substantially lower average complexity, but which also provide more accurate localisation of statistical group differences.

Quantifying local changes to the airway wall surfaces from computed tomography images is important in the study of diseases such as chronic obstructive pulmonary disease. Current approaches segment the airways in the individual time point images and subsequently aggregate per airway generation or perform branch matching to assess regional changes. In contrast, we propose an integrated approach analysing the time points simultaneously using a subject-specific groupwise space and 4D optimal surface segmentation. The method combines information from all time points and measurements are matched locally at any position on the resulting surfaces. Visual inspection of the scans of 10 subjects showed increased tree length compared to the state of the art with little change in the amount of false positives. A large scale analysis of the airways of 374 subjects including a total of 1870 images showed significant correlation with lung function and high reproducibility of the measurements.

The combination of functional and anatomical imaging technologies such as Positron Emission Tomography (PET) and Computed Tomography (CT) has shown its value in the preclinical and clinical fields. In PET/CT hybrid acquisition systems, CT-derived attenuation maps enable a more accurate PET reconstruction. However, CT provides only very limited soft-tissue contrast and exposes the patient to an additional radiation dose. In comparison, Magnetic Resonance Imaging (MRI) provides good soft-tissue contrast and the ability to study functional activation and tissue microstructures, but does not directly provide patient-specific electron density maps for PET reconstruction. The aim of the proposed work is to improve PET/MR reconstruction by generating synthetic CTs and attenuation-maps. The synthetic images are generated through a multi-atlas information propagation scheme, locally matching the MRI-derived patient's morphology to a database of pre-acquired MRI/CT pairs. Results show improvements in CT synthesis and PET reconstruction accuracy when compared to a segmentation method using an Ultrashort-Echo-Time MRI sequence.

X-ray mammography is routinely used in national screening programmes and as a clinical diagnostic tool. Magnetic Resonance Imaging (MRI) is commonly used as a complementary modality, providing functional information about the breast and a 3D image that can overcome ambiguities caused by the superimposition of fibro-glandular structures associated with X-ray imaging. Relating findings between these modalities is a challenging task however, due to the different imaging processes involved and the large deformation that the breast undergoes. In this work we present a registration method to determine spatial correspondence between pairs of MR and X-ray images of the breast, that is targeted for clinical use. We propose a generic registration framework which incorporates a volume-preserving affine transformation model and validate its performance using routinely acquired clinical data. Experiments on simulated mammograms from 8 volunteers produced a mean registration error of 3.8±1.6mm for a mean of 12 manually identified landmarks per volume. When validated using 57 lesions identified on routine clinical CC and MLO mammograms (n=113 registration tasks) from 49 subjects the median registration error was 13.1mm. When applied to the registration of an MR image to CC and MLO mammograms of a patient with a localisation clip, the mean error was 8.9mm. The results indicate that an intensity based registration algorithm, using a relatively simple transformation model, can provide radiologists with a clinically useful tool for breast cancer diagnosis., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Cortical thinning is a widely used and powerful biomarker for measuring disease progression in Alzheimer's disease (AD). However, there has been little work on the effect of atrophy on the cortical folding patterns. In this study, we examined whether the cortical folding could be used as a biomarker of AD. Cortical folding metrics were computed on 678 patients from the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. For each subject, the boundary between grey matter and white matter was extracted using a level set technique. At each point on the boundary two metrics characterising folding, curvedness and shape index, were generated. Joint histograms using these metrics were calculated for five regions of interest (ROIs): frontal, temporal, occipital, and parietal lobes as well as the cingulum. Pixelwise statistical maps were generated from the joint histograms using permutations tests. In each ROI, a significant reduction was observed between controls and AD in areas associated with the sulcal folds, suggesting a sulcal opening associated with neurodegeneration. When comparing to MCI patients, the regions of significance were smaller but overlapping with those regions found comparing controls to AD. It indicates that the differences in cortical folding are progressive and can be detected before formal diagnosis of AD. Our preliminary analysis showed a viable signal in the cortical folding patterns for Alzheimer's disease that should be explored further.

Automatic pectoral muscle removal on medio-lateral oblique view of mammogram is an essential step for many mammographic processing algorithms. However, the wide variability in the position of the muscle contour, together with the similarity between in muscle and breast tissues makes the detection a difficult task. In this paper, we propose a two step procedure to detect the muscle contour. In a first step, the endpoints of the contour are predicted with a pair of support vector regression models; one model is trained to predict the intersection point of the contour with the top row while the other is designed for the prediction of the endpoint of the contour on the left column. Next, the muscle contour is computed as the shortest path between the two endpoints. A comprehensive comparison with manually-drawn contours reveals the strength of the proposed method.

The automatic detection and segmentation of the pectoral muscle in the medio-lateral oblique view of mammograms is essential for further analysis of breast anormalies. However, it is still a very difficult task since the sizes, shapes and intensity contrasts of pectoral muscles change greatly from image to image. In this paper, an algorithm based on the shortest path on a graph is proposed to automatically detect the pectoral muscle contour. To overcome the difficulties of searching for the path between a lateral and the top margins of the image, this is first transformed, using polar coordinates. In the transformed image, the muscle boundary in amongst the shortest paths between the top and the bottom rows. A comprehensive comparison with manually-drawn contours reveals the strength of the proposed method.

Thickness measurements of the cerebral cortex can aid diagnosis and provide valuable information about the temporal evolution of several diseases such as Alzheimer's, Huntington's, Schizophrenia, as well as normal ageing. The presence of deep sulci and 'collapsed gyri' (caused by the loss of tissue in patients with neurodegenerative diseases) complicates the tissue segmentation due to partial volume (PV) effects and limited resolution of MRI. We extend existing work to improve the segmentation and thickness estimation in a single framework. We model the PV effect using a maximum a posteriori approach with novel iterative modification of the prior information to enhance deep sulci and gyri delineation. We use a voxel based approach to estimate thickness using the Laplace equation within a Lagrangian-Eulerian framework leading to sub-voxel accuracy. Experiments performed on a new digital phantom and on clinical Alzheimer's disease MR images show improvements in both accuracy and robustness of the thickness measurements, as well as a reduction of errors in deep sulci and collapsed gyri.