Search

Quadrantanopia caused by inadvertent severing of Meyer's Loop of the optic radiation is a well-recognised complication of temporal lobectomy for conditions such as epilepsy. Dissection studies indicate that the anterior extent of Meyer's Loop varies considerably between individuals. Quantifying this for individual patients is thus an important step to improve the safety profile of temporal lobectomies. Previous attempts to delineate Meyer's Loop using diffusion MRI tractography have had difficulty estimating its full anterior extent, required manual ROI placement, and/or relied on advanced diffusion sequences that cannot be acquired routinely in most clinics. Here we present CONSULT: a pipeline that can delineate the optic radiation from raw DICOM data in a completely automated way via a combination of robust pre-processing, segmentation, and alignment stages, plus simple improvements that bolster the efficiency and reliability of standard tractography. We tested CONSULT on 696 scans of predominantly healthy participants (539 unique brains), including both advanced acquisitions and simpler acquisitions that could be acquired in clinically acceptable timeframes. Delineations completed without error in 99.4% of the scans. The distance between Meyer's Loop and the temporal pole closely matched both averages and ranges reported in dissection studies for all tested sequences. Median scan-rescan error of this distance was 1 mm. When tested on two participants with considerable pathology, delineations were successful and realistic. Through this, we demonstrate not only how to identify Meyer's Loop with clinically feasible sequences, but also that this can be achieved without fundamental changes to tractography algorithms or complex post-processing methods., (© 2021 Commonwealth of Australia. Human Brain Mapping published by Wiley Periodicals LLC.)

Background: Centiloid was introduced to harmonise β-Amyloid (Aβ) PET quantification across different tracers, scanners and analysis techniques. Unfortunately, Centiloid still suffers from some quantification disparities in longitudinal analysis when normalising data from different tracers or scanners. In this work, we aim to reduce this variability using a different analysis technique applied to the existing calibration data., Method: All PET images from the Centiloid calibration dataset, along with 3762 PET images from the AIBL study were analysed using the recommended SPM pipeline. The PET images were SUVR normalised using the whole cerebellum. All SUVR normalised PiB images from the calibration dataset were decomposed using non-negative matrix factorisation (NMF). The NMF coefficients related to the first component were strongly correlated with global SUVR and were subsequently used as a surrogate for Aβ retention. For each tracer of the calibration dataset, the components of the NMF were computed in a way such that the coefficients of the first component would match those of the corresponding PiB. Given the strong correlations between the SUVR and the NMF coefficients on the calibration dataset, all PET images from AIBL were subsequently decomposed using the computed NMF, and their coefficients transformed into Centiloids., Results: Using the AIBL data, the correlation between the standard Centiloid and the novel NMF-based Centiloid was high in each tracer. The NMF-based Centiloids showed a reduction of outliers, and improved longitudinal consistency. Furthermore, it removed the effects of switching tracers from the longitudinal variance of the Centiloid measure, when assessed using a linear mixed effects model., Conclusion: We here propose a novel image driven method to perform the Centiloid quantification. The methods is highly correlated with standard Centiloids while improving the longitudinal reliability when switching tracers. Implementation of this method across multiple studies may lend to more robust and comparable data for future research., (Copyright © 2020. Published by Elsevier Inc.)

Accurate quantification of white matter hyperintensities (WMH) from Magnetic Resonance Imaging (MRI) is a valuable tool for the analysis of normal brain ageing or neurodegeneration. Reliable automatic extraction of WMH lesions is challenging due to their heterogeneous spatial occurrence, their small size and their diffuse nature. In this paper, we present an automatic method to segment these lesions based on an ensemble of overcomplete patch-based neural networks. The proposed method successfully provides accurate and regular segmentations due to its overcomplete nature while minimizing the segmentation error by using a boosted ensemble of neural networks. The proposed method compared favourably to state of the art techniques using two different neurodegenerative datasets., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Purpose: To validate a fully automated scheme to extract biochemical information from the hip joint cartilages using MR T2 mapping images incorporating segmentation of co-registered three-dimensional Fast-Spin-Echo (3D-SPACE) images., Methods: Manual analyses of unilateral hip (3 Tesla) MR images of 24 asymptomatic volunteers were used to validate a 3D deformable model method for automated cartilage segmentation of SPACE scans, partitioning of the individual femoral and acetabular cartilage plates into clinically defined sub-regions and propagating these results to T2 maps to calculate region-wise T2 value statistics. Analyses were completed on a desktop computer (∼ 10 min per case)., Results: The mean voxel overlap between automated A and manual M segmentations of the cartilage volumes in the (clinically based) SPACE images was 73% (100 × 2|A∩M|/[|A|+|M|]). The automated and manual analyses demonstrated a relative difference error <10% in the median "T2 average signal" for each cartilage plate. The automated and manual analyses showed consistent patterns between significant differences in T2 data across the hip cartilage sub-regions., Conclusion: The good agreement between the manual and automatic analyses of T2 values indicates the use of structural 3D-SPACE MR images with the proposed method provides a promising approach for automated quantitative T2 assessment of hip joint cartilages., (© 2015 Wiley Periodicals, Inc.)

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

Many medical image processing techniques rely on accurate shape modeling of anatomical features. The presence of shape abnormalities challenges traditional processing algorithms based on strong morphological priors. In this work, a sparse shape reconstruction from a statistical shape model is presented. It combines the advantages of traditional statistical shape models (defining a 'normal' shape space) and previously presented sparse shape composition (providing localized descriptors of anomalies). The algorithm was incorporated into our image segmentation and classification software. Evaluation was performed on simulated and clinical MRI data from 22 sciatica patients with intervertebral disc herniation, containing 35 herniated and 97 normal discs. Moderate to high correlation (R=0.73) was achieved between simulated and detected herniations. The sparse reconstruction provided novel quantitative features describing the herniation morphology and MRI signal appearance in three dimensions (3D). The proposed descriptors of local disc morphology resulted to the 3D segmentation accuracy of 1.07±1.00mm (mean absolute vertex-to-vertex mesh distance over the posterior disc region), and improved the intervertebral disc classification from 0.888 to 0.931 (area under receiver operating curve). The results show that the sparse shape reconstruction may improve computer-aided diagnosis of pathological conditions presenting local morphological alterations, as seen in intervertebral disc herniation., (Crown Copyright © 2015. Published by Elsevier Ltd. All rights reserved.)

Accurate segmentation of hip joint cartilage from magnetic resonance (MR) images offers opportunities for quantitative investigations of pathoanatomical conditions such as osteoarthritis. In this paper, we present a fully automatic scheme for the segmentation of the individual femoral and acetabular cartilage plates in the human hip joint from high-resolution 3D MR images. The developed scheme uses an improved optimal multi-object multi-surface graph search framework with an arc-weighted graph representation that incorporates prior morphological knowledge as a basis for segmentation of the individual femoral and acetabular cartilage plates despite weak or incomplete boundary interfaces. This automated scheme was validated against manual segmentations from 3D true fast imaging with steady-state precession (TrueFISP) MR examinations of the right hip joints in 52 asymptomatic volunteers. Compared with expert manual segmentations of the combined, femoral and acetabular cartilage volumes, the automatic scheme obtained mean (± standard deviation) Dice's similarity coefficients of 0.81 (± 0.03), 0.79 (± 0.03) and 0.72 (± 0.05). The corresponding mean absolute volume difference errors were 8.44% (± 6.36), 9.44% (± 7.19) and 9.05% (± 8.02). The mean absolute differences between manual and automated measures of cartilage thickness for femoral and acetabular cartilage plates were 0.13 mm (± 0.12) and 0.11 mm (± 0.11), respectively.

Deformable models incorporating shape priors have proved to be a successful approach in segmenting anatomical regions and specific structures in medical images. This paper introduces weighted shape priors for deformable models in the context of 3D magnetic resonance (MR) image segmentation of the bony elements of the human hip joint. The fully automated approach allows the focusing of the shape model energy to a priori selected anatomical structures or regions of clinical interest by preferentially ordering the shape representation (or eigen-modes) within this type of model to the highly weighted areas. This focused shape model improves accuracy of the shape constraints in those regions compared to standard approaches. The proposed method achieved femoral head and acetabular bone segmentation mean absolute surface distance errors of 0.55±0.18mm and 0.75±0.20mm respectively in 35 3D unilateral MR datasets from 25 subjects acquired at 3T with different limited field of views for individual bony components of the hip joint., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Objective: As large-scale medical imaging studies are becoming more common, there is an increasing reliance on automated software to extract quantitative information from these images. As the size of the cohorts keeps increasing with large studies, there is a also a need for tools that allow results from automated image processing and analysis to be presented in a way that enables fast and efficient quality checking, tagging and reporting on cases in which automatic processing failed or was problematic., Materials and Methods: MilxXplore is an open source visualization platform, which provides an interface to navigate and explore imaging data in a web browser, giving the end user the opportunity to perform quality control and reporting in a user friendly, collaborative and efficient way., Discussion: Compared to existing software solutions that often provide an overview of the results at the subject's level, MilxXplore pools the results of individual subjects and time points together, allowing easy and efficient navigation and browsing through the different acquisitions of a subject over time, and comparing the results against the rest of the population., Conclusions: MilxXplore is fast, flexible and allows remote quality checks of processed imaging data, facilitating data sharing and collaboration across multiple locations, and can be easily integrated into a cloud computing pipeline. With the growing trend of open data and open science, such a tool will become increasingly important to share and publish results of imaging analysis.

Due to physical limitations inherent in magnetic resonance imaging scanners, three dimensional volumetric scans are often acquired with anisotropic voxel resolution. We investigate several interpolation approaches to reduce the anisotropy and present a novel approach - constrained reverse diffusion for thick slice interpolation. This technique was compared to common methods: linear and cubic B-Spline interpolation and a technique based on non-rigid registration of neighboring slices. The methods were evaluated on artificial MR phantoms and real MR scans of human brain. The constrained reverse diffusion approach delivered promising results and provides an alternative for thick slice interpolation, especially for higher anisotropy factors., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Identifying mild cognitive impairment (MCI) subjects who will convert to clinical Alzheimer's disease (AD) is important for therapeutic decisions, patient counselling and clinical trials. Hippocampal volume and rate of atrophy predict clinical decline at the MCI stage and progression to AD. In this paper, we create p-maps from the differences in the shape of the hippocampus between 60 normal controls and 60 AD subjects using statistical shape models, and generate different regions of interest (ROI) by thresholding the p-maps at different significance levels. We demonstrate increased statistical power to classify 86 MCI converters and 128 MCI stable subjects using the hippocampal atrophy rates calculated by the boundary shift integral within these ROIs.

Accurate cortical thickness estimation is important for the study of many neurodegenerative diseases. Many approaches have been previously proposed, which can be broadly categorised as mesh-based and voxel-based. While the mesh-based approaches can potentially achieve subvoxel resolution, they usually lack the computational efficiency needed for clinical applications and large database studies. In contrast, voxel-based approaches, are computationally efficient, but lack accuracy. The aim of this paper is to propose a novel voxel-based method based upon the Laplacian definition of thickness that is both accurate and computationally efficient. A framework was developed to estimate and integrate the partial volume information within the thickness estimation process. Firstly, in a Lagrangian step, the boundaries are initialized using the partial volume information. Subsequently, in an Eulerian step, a pair of partial differential equations are solved on the remaining voxels to finally compute the thickness. Using partial volume information significantly improved the accuracy of the thickness estimation on synthetic phantoms, and improved reproducibility on real data. Significant differences in the hippocampus and temporal lobe between healthy controls (NC), mild cognitive impaired (MCI) and Alzheimer's disease (AD) patients were found on clinical data from the ADNI database. We compared our method in terms of precision, computational speed and statistical power against the Eulerian approach. With a slight increase in computation time, accuracy and precision were greatly improved. Power analysis demonstrated the ability of our method to yield statistically significant results when comparing AD and NC. Overall, with our method the number of samples is reduced by 25% to find significant differences between the two groups.

This study presents a novel method for the automatic segmentation of the quadratus lumborum (QL) muscle from axial magnetic resonance (MR) images using a hybrid scheme incorporating the use of non-rigid registration with probabilistic atlases (PAs) and geodesic active contours (GACs). The scheme was evaluated on an MR database of 7mm axial images of the lumbar spine from 20 subjects (fast bowlers and athletic controls). This scheme involved several steps, including (i) image pre-processing, (ii) generation of PAs for the QL, psoas (PS) and erector spinae+multifidus (ES+MT) muscles and (iii) segmentation, using 3D GACs initialized and constrained by the propagation of the PAs using non-rigid registration. Pre-processing of the images involved bias field correction based on local entropy minimization with a bicubic spline model and a reverse diffusion interpolation algorithm to increase the slice resolution to 0.98 x 0.98 x 1.75mm. The processed images were then registered (affine and non-rigid) and used to generate an average atlas. The PAs for the QL, PS and ES+MT were then generated by propagation of manual segmentations. These atlases were further analysed with specialised filtering to constrain the QL segmentation from adjacent non-muscle tissues (kidney, fat). This information was then used in 3D GACs to obtain the final segmentation of the QL. The automatic segmentation results were compared with the manual segmentations using the Dice similarity metric (DSC), with a median DSC for the right and left QL muscles of 0.78 (mean = 0.77, sd=0.07) and 0.75 (mean =0.74, sd=0.07), respectively.

Accurate cortical thickness estimation in-vivo is important for the study of many neurodegenerative diseases. When using magnetic resonance images (MRI), accuracy may be hampered by artifacts such as partial volume (PV) as the cortex spans only a few voxels. In zones of opposed sulcal banks (tight sulci) the measurement can be even more difficult. The aim of this work is to propose a voxel-based cortical thickness estimation method from MR by integrating a mechanism for correcting sulci delineation after an improved partial volume classification. First, an efficient and accurate framework was developed to enhance partial volume classification with structural information. Then, the correction of sulci delineation is performed after a homotopic thinning of a cost function image. Integrated to our voxel-based cortical thickness estimation pipeline, the overall method showed a better estimate of thickness and a high reproducibility on real data (R2 > 0.9). A quantitative analysis on clinical data from an Alzheimer's disease study showed significant differences between normal controls and Alzheimer's disease patients.

Beta-amyloid (Abeta) plaques are one of the neuropathological hallmarks of Alzheimer's disease (AD) and can be quantified using the marker 11C PiB. As l1C PiB PET images have limited anatomical information, an Magnetic Resonance Image (MRI) is usually acquired to perform the spatial normalization needed for population analysis. We designed and evaluated a high dimensional spatial normalization approach that only uses the 11C PiB PET image. The non-rigid registration (NRR) is based on free form deformation (FFD) modelled using B-splines. To compensate for the limited anatomical information, the FFD is constrained to an allowable transform space using a model trained from MR registrations. Abeta deposition is dependent on disease staging, so a spatially normalized 11C PiB PET appearance model selects and refines the atlas. The approach was compared with MR NRR using data from healthy elderly, mild cognitive impaired and Alzheimer disease participants. Using segmentation propagation, an average Dice similarity coefficient of 0.64 and 0.73 was obtained for white and gray matter. The R-squared correlation between the uptake obtained in the frontal, parietal, occipital and temporal was 0.789, 0.843, 0.871 and 0.964. These are very promising results, considering the low resolution of 11C PiB PET images.

To perform cartilage quantitative analysis requires the accurate segmentation of each individual cartilage. In this paper we present a model based scheme that can automatically and accurately segment each individual cartilage in healthy knees from a clinical MR sequence (fat suppressed spoiled gradient recall). This scheme consists of three stages; the automatic segmentation of the bones, the extraction of the bone-cartilage interfaces (BCI) and segmentation of the cartilages. The bone segmentation is performed using three-dimensional active shape models. The BCI is extracted using image information and prior knowledge about the likelihood of each point belonging to the interface. A cartilage thickness model then provides constraints and regularizes the cartilage segmentation performed from the BCI. The accuracy and robustness of the approach was experimentally validated, with (patellar, tibial and femoral) cartilage segmentations having a median DSC of (0.870, 0.855, 0.870), performing significantly better than non-rigid registration (0.787, 0.814, 0.795). The total cartilage segmentation had an average DSC of (0.891), close to the (0.896) obtained using a semi-automatic watershed algorithm. The error in quantitative volume and thickness measures was (8.29, 4.94, 5.56)% and (0.19, 0.33, 0.10) mm respectively.

This paper considers the problem of automatic classification of textured tissues in 3D MRI. More specifically, it aims at validating the use of features extracted from the phase of the MR signal to improve texture discrimination in bone segmentation. This extra information provides better segmentation, compared to using magnitude only features. We also present a novel multiscale scheme to improve the speed of pixel based classification algorithm, such as support vector machines. This algorithm dramatically increases the speed of the segmentation process by an order of magnitude through a reduction of the number of pixels that needs to be classified in the image.

This paper considers the problem of tissue classification in 3D MRI. More specifically, a new set of texture features, based on phase information, is used to perform the segmentation of the bones of the knee. The phase information provides a very good discrimination between the bone and the surrounding tissues, but is usually not used due to phase unwrapping problems. We present a method to extract textural information from the phase that does not require phase unwrapping. The textural information extracted from the magnitude and the phase can be combined to perform tissue classification, and used to initialise an active shape model, leading to a more precise segmentation.