MAFT-SO: A novel multi-atlas fusion template based on spatial overlap for ASD diagnosis.

Autism spectrum disorder (ASD) is a common neurological condition. Early diagnosis and treatment are essential for enhancing the life quality of individuals with ASD. However, most existing studies either focus solely on the brain networks of subjects within a single atlas or merely employ simple matrix concatenation to represent the fusion of multi-atlas. These approaches neglected the natural spatial overlap that exists between brain regions across multi-atlas and did not fully capture the comprehensive information of brain regions under different atlases. To tackle this weakness, in this paper, we propose a novel multi-atlas fusion template based on spatial overlap degree of brain regions, which aims to obtain a comprehensive representation of brain networks. Specifically, we formally define a measurement of the spatial overlap among brain regions across different atlases, named spatial overlap degree. Then, we fuse the multi-atlas to obtain brain networks of each subject based on spatial overlap. Finally, the GCN is used to perform the final classification. The experimental results on Autism Brain Imaging Data Exchange (ABIDE) demonstrate that our proposed method achieved an accuracy of 0.757. Overall, our method outperforms SOTA methods in ASD/TC classification.

Subject-specific atlas for automatic brain tissue segmentation of neonatal magnetic resonance images.

Developing advanced systems for 3D brain tissue segmentation from neonatal magnetic resonance (MR) images is vital for newborn structural analysis. However, automatic segmentation of neonatal brain tissues is challenging due to smaller head size and inverted T1/T2 tissue contrast compared to adults. In this work, a subject-specific atlas based technique is presented for segmentation of gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) from neonatal MR images. It involves atlas selection, subject-specific atlas creation using random forest (RF) classifier, and brain tissue segmentation using the expectation maximization-Markov random field (EM-MRF) method. To increase the segmentation accuracy, different tissue intensity- and gradient-based features were used. Evaluation on 40 neonatal MR images (gestational age of 37-44 weeks) demonstrated an overall accuracy of 94.3% and an average Dice similarity coefficient (DSC) of 0.945 (GM), 0.947 (WM), and 0.912 (CSF). Compared to multi-atlas segmentation methods like SEGMA and EM-MRF with multiple atlases, our method improved accuracy by up to 4%, particularly in complex tissue regions. Our proposed method allows accurate brain tissue segmentation, a crucial step in brain magnetic resonance imaging (MRI) applications including brain surface reconstruction and realistic head model creation in neonates.

Dense Multi-Scale Graph Convolutional Network for Knee Joint Cartilage Segmentation.

In this paper, we propose a dense multi-scale adaptive graph convolutional network (DMA-GCN) method for automatic segmentation of the knee joint cartilage from MR images. Under the multi-atlas setting, the suggested approach exhibits several novelties, as described in the following. First, our models integrate both local-level and global-level learning simultaneously. The local learning task aggregates spatial contextual information from aligned spatial neighborhoods of nodes, at multiple scales, while global learning explores pairwise affinities between nodes, located globally at different positions in the image. We propose two different structures of building models, whereby the local and global convolutional units are combined by following an alternating or a sequential manner. Secondly, based on the previous models, we develop the DMA-GCN network, by utilizing a densely connected architecture with residual skip connections. This is a deeper GCN structure, expanded over different block layers, thus being capable of providing more expressive node feature representations. Third, all units pertaining to the overall network are equipped with their individual adaptive graph learning mechanism, which allows the graph structures to be automatically learned during training. The proposed cartilage segmentation method is evaluated on the entire publicly available Osteoarthritis Initiative (OAI) cohort. To this end, we have devised a thorough experimental setup, with the goal of investigating the effect of several factors of our approach on the classification rates. Furthermore, we present exhaustive comparative results, considering traditional existing methods, six deep learning segmentation methods, and seven graph-based convolution methods, including the currently most representative models from this field. The obtained results demonstrate that the DMA-GCN outperforms all competing methods across all evaluation measures, providing DSC=95.71% and DSC=94.02% for the segmentation of femoral and tibial cartilage, respectively.

A collection of 157 individual neuromelanin-sensitive images accompanied by non-linear neuromelanin-sensitive atlas and a probabilistic locus coeruleus atlas.

The current dataset aims to support and enhance the research reliability of neuromelanin regions in the brainstem, such as locus coeruleus (LC), by offering raw neuromelanin-sensitive images. The dataset includes raw neuromelanin-sensitive images from 157 healthy individuals (8-64 years old). In addition, leveraging individual neuromelanin-sensitive images, a non-linear neuromelanin-sensitive atlas, generated through an iterative warping process, is included to tackle the common challenge of a limited field of view in neuromelanin-sensitive images. Finally, the dataset encompasses a probabilistic LC atlas generated through a majority voting approach with pre-existing multiple atlas-based segmentations. This process entails warping pre-existing atlases onto individual spaces and identifying voxels with a majority consensus of over 50 % across the atlases. This LC probabilistic atlas can minimize uncertainty variance associated with choosing a specific single atlas.

An open source auto-segmentation algorithm for delineating heart and substructures - Development and validation within a multicenter lung cancer cohort.

BACKGROUND AND PURPOSE: Irradiation of the heart in thoracic cancers raises toxicity concerns. For accurate dose estimation, automated heart and substructure segmentation is potentially useful. In this study, a hybrid automatic segmentation is developed. The accuracy of delineation and dose predictions were evaluated, testing the method's potential within heart toxicity studies. MATERIALS AND METHODS: The hybrid segmentation method delineated the heart, four chambers, three large vessels, and the coronary arteries. The method consisted of a nnU-net heart segmentation and partly atlas- and model-based segmentation of the substructures. The nnU-net training and atlas segmentation was based on lung cancer patients and was validated against a national consensus dataset of 12 patients with breast cancer. The accuracy of dose predictions between manual and auto-segmented heart and substructures was evaluated by transferring the dose distribution of 240 previously treated lung cancer patients to the consensus data set. RESULTS: The hybrid auto-segmentation method performed well with a heart dice similarity coefficient (DSC) of 0.95, with no statistically significant difference between the automatic and manual delineations. The DSC for the chambers varied from 0.78-0.86 for the automatic segmentation and was comparable with the inter-observer variability. Most importantly, the automatic segmentation was as precise as the clinical experts in predicting the dose distribution to the heart and all substructures. CONCLUSION: The hybrid segmentation method performed well in delineating the heart and substructures. The prediction of dose by the automatic segmentation was aligned with the manual delineations, enabling measurement of heart and substructure dose in large cohorts. The delineation algorithm will be available for download.

HGM-cNet: Integrating hippocampal gray matter probability map into a cascaded deep learning framework improves hippocampus segmentation.

A robust cascaded deep learning framework with integrated hippocampal gray matter (HGM) probability map was developed to improve the hippocampus segmentation (called HGM-cNet) due to its significance in various neuropsychiatric disorders such as Alzheimer's disease (AD). Particularly, the HGM-cNet cascaded two identical convolutional neural networks (CNN), where each CNN was devised by incorporating Attention Block, Residual Block, and DropBlock into the typical encoder-decoder architecture. The two CNNs were skip-connected between encoder components at each scale. The adoption of the cascaded deep learning framework was to conveniently incorporate the HGM probability map with the feature map generated by the first CNN. Experiments on 135T1-weighted MRI scans and manual hippocampal labels from publicly available ADNI-HarP dataset demonstrated that the proposed HGM-cNet outperformed seven multi-atlas-based hippocampus segmentation methods and six deep learning methods under comparison in most evaluation metrics. The Dice (average > 0.89 for both left and right hippocampus) was increased by around or more than 1% over other methods. The HGM-cNet also achieved a superior hippocampus segmentation performance in each group of cognitive normal, mild cognitive impairment, and AD. The stability, conveniences and generalizability of the cascaded deep learning framework with integrated HGM probability map in improving hippocampus segmentation was validated by replacing the proposed CNN with 3D-UNet, Atten-UNet, HippoDeep, QuickNet, DeepHarp, and TransBTS models. The integration of the HGM probability map in the cascaded deep learning framework was also demonstrated to facilitate capturing hippocampal atrophy more accurately than alternative methods in AD analysis. The codes are publicly available at https://github.com/Liu1436510768/HGM-cNet.git.

Segmentation of the aorta in systolic phase from 4D flow MRI: multi-atlas vs. deep learning.

OBJECTIVE: In the management of the aortic aneurysm, 4D flow magnetic resonance Imaging provides valuable information for the computation of new biomarkers using computational fluid dynamics (CFD). However, accurate segmentation of the aorta is required. Thus, our objective is to evaluate the performance of two automatic segmentation methods on the calculation of aortic wall pressure. METHODS: Automatic segmentation of the aorta was performed with methods based on deep learning and multi-atlas using the systolic phase in the 4D flow MRI magnitude image of 36 patients. Using mesh morphing, isotopological meshes were generated, and CFD was performed to calculate the aortic wall pressure. Node-to-node comparisons of the pressure results were made to identify the most robust automatic method respect to the pressures obtained with a manually segmented model. RESULTS: Deep learning approach presented the best segmentation performance with a mean Dice similarity coefficient and a mean Hausdorff distance (HD) equal to 0.92+/- 0.02 and 21.02+/- 24.20 mm, respectively. At the global level HD is affected by the performance in the abdominal aorta. Locally, this distance decreases to 9.41+/- 3.45 and 5.82+/- 6.23 for the ascending and descending thoracic aorta, respectively. Moreover, with respect to the pressures from the manual segmentations, the differences in the pressures computed from deep learning were lower than those computed from multi-atlas method. CONCLUSION: To reduce biases in the calculation of aortic wall pressure, accurate segmentation is needed, particularly in regions with high blood flow velocities. Thus, the deep learning segmen-tation method should be preferred.

Deep label fusion: A generalizable hybrid multi-atlas and deep convolutional neural network for medical image segmentation.

Deep convolutional neural networks (DCNN) achieve very high accuracy in segmenting various anatomical structures in medical images but often suffer from relatively poor generalizability. Multi-atlas segmentation (MAS), while less accurate than DCNN in many applications, tends to generalize well to unseen datasets with different characteristics from the training dataset. Several groups have attempted to integrate the power of DCNN to learn complex data representations and the robustness of MAS to changes in image characteristics. However, these studies primarily focused on replacing individual components of MAS with DCNN models and reported marginal improvements in accuracy. In this study we describe and evaluate a 3D end-to-end hybrid MAS and DCNN segmentation pipeline, called Deep Label Fusion (DLF). The DLF pipeline consists of two main components with learnable weights, including a weighted voting subnet that mimics the MAS algorithm and a fine-tuning subnet that corrects residual segmentation errors to improve final segmentation accuracy. We evaluate DLF on five datasets that represent a diversity of anatomical structures (medial temporal lobe subregions and lumbar vertebrae) and imaging modalities (multi-modality, multi-field-strength MRI and Computational Tomography). These experiments show that DLF achieves comparable segmentation accuracy to nnU-Net (Isensee et al., 2020), the state-of-the-art DCNN pipeline, when evaluated on a dataset with similar characteristics to the training datasets, while outperforming nnU-Net on tasks that involve generalization to datasets with different characteristics (different MRI field strength or different patient population). DLF is also shown to consistently improve upon conventional MAS methods. In addition, a modality augmentation strategy tailored for multimodal imaging is proposed and demonstrated to be beneficial in improving the segmentation accuracy of learning-based methods, including DLF and DCNN, in missing data scenarios in test time as well as increasing the interpretability of the contribution of each individual modality.

DeepHeartCT: A fully automatic artificial intelligence hybrid framework based on convolutional neural network and multi-atlas segmentation for multi-structure cardiac computed tomography angiography image segmentation.

Cardiac computed tomography angiography (CTA) is an emerging imaging modality for assessing coronary artery as well as various cardiovascular structures. Recently, deep learning (DL) methods have been successfully applied to many applications of medical image analysis including cardiac CTA structure segmentation. However, DL requires a large amounts of data and high-quality labels for training which can be burdensome to obtain due to its labor-intensive nature. In this study, we aim to develop a fully automatic artificial intelligence (AI) system, named DeepHeartCT, for accurate and rapid cardiac CTA segmentation based on DL. The proposed system was trained using a large clinical dataset with computer-generated labels to segment various cardiovascular structures including left and right ventricles (LV, RV), left and right atria (LA, RA), and LV myocardium (LVM). This new system was trained directly using high-quality computer labels generated from our previously developed multi-atlas based AI system. In addition, a reverse ranking strategy was proposed to assess the segmentation quality in the absence of manual reference labels. This strategy allowed the new framework to assemble optimal computer-generated labels from a large dataset for effective training of a deep convolutional neural network (CNN). A large clinical cardiac CTA studies (n = 1,064) were used to train and validate our framework. The trained model was then tested on another independent dataset with manual labels (n = 60). The Dice score, Hausdorff distance and mean surface distance were used to quantify the segmentation accuracy. The proposed DeepHeartCT framework yields a high median Dice score of 0.90 [interquartile range (IQR), 0.90-0.91], a low median Hausdorff distance of 7 mm (IQR, 4-15 mm) and a low mean surface distance of 0.80 mm (IQR, 0.57-1.29 mm) across all segmented structures. An additional experiment was conducted to evaluate the proposed DL-based AI framework trained with a small vs. large dataset. The results show our framework also performed well when trained on a small optimal training dataset (n = 110) with a significantly reduced training time. These results demonstrated that the proposed DeepHeartCT framework provides accurate and rapid cardiac CTA segmentation that can be readily generalized for handling large-scale medical imaging applications.

A novel multi-atlas segmentation approach under the semi-supervised learning framework: Application to knee cartilage segmentation.

BACKGROUND AND OBJECTIVE: Multi-atlas based segmentation techniques, which rely on an atlas library comprised of training images labeled by an expert, have proven their effectiveness in multiple automatic segmentation applications. However, the usage of exhaustive patch libraries combined with the voxel-wise labeling incur a large computational cost in terms of memory requirements and execution times. METHODS: To confront this shortcoming, we propose a novel two-stage multi-atlas approach designed under the Semi-Supervised Learning (SSL) framework. The main properties of our method are as follows: First, instead of the voxel-wise labeling approach, the labeling of target voxels is accomplished here by exploiting the spectral content of globally sampled datasets from the target image, along with their spatially correspondent data collected from the atlases. Following SSL, voxels classification is boosted by incorporating unlabeled data from the target image, in addition to the labeled ones from atlas library. Our scheme integrates constructively fruitful concepts, including sparse reconstructions of voxels from linear neighborhoods, HOG feature descriptors of patches/regions, and label propagation via sparse graph constructions. Segmentation of the target image is carried out in two stages: stage-1 focuses on the sampling and labeling of global data, while stage-2 undertakes the above tasks for the out-of-sample data. Finally, we propose different graph-based methods for the labeling of global data, while these methods are extended to deal with the out-of-sample voxels. RESULTS: A thorough experimental investigation is conducted on 76 subjects provided by the publicly accessible Osteoarthritis Initiative (OAI) repository. Comparative results and statistical analysis demonstrate that the suggested methodology exhibits superior segmentation performance compared to the existing patch-based methods, across all evaluation metrics (DSC:88.89%, Precision: 89.86%, Recall: 88.12%), while at the same time it requires a considerably reduced computational load (>70% reduction on average execution time with respect to other patch-based). In addition, our approach is favorably compared against other non patch-based and deep learning methods in terms of performance accuracy (on the 3-class problem). A final experimentation on a 5-class setting of the problems demonstrates that our approach is capable of achieving performance comparable to existing state-of-the-art knee cartilage segmentation methods (DSC:88.22% and DSC:85.84% for femoral and tibial cartilage respectively).

Multi-atlas segmentation and quantification of muscle, bone and subcutaneous adipose tissue in the lower leg using peripheral quantitative computed tomography.

Accurate and reproducible tissue identification is essential for understanding structural and functional changes that may occur naturally with aging, or because of a chronic disease, or in response to intervention therapies. Peripheral quantitative computed tomography (pQCT) is regularly employed for body composition studies, especially for the structural and material properties of the bone. Furthermore, pQCT acquisition requires low radiation dose and the scanner is compact and portable. However, pQCT scans have limited spatial resolution and moderate SNR. pQCT image quality is frequently degraded by involuntary subject movement during image acquisition. These limitations may often compromise the accuracy of tissue quantification, and emphasize the need for automated and robust quantification methods. We propose a tissue identification and quantification methodology that addresses image quality limitations and artifacts, with increased interest in subject movement. We introduce a multi-atlas image segmentation (MAIS) framework for semantic segmentation of hard and soft tissues in pQCT scans at multiple levels of the lower leg. We describe the stages of statistical atlas generation, deformable registration and multi-tissue classifier fusion. We evaluated the performance of our methodology using multiple deformable registration approaches against reference tissue masks. We also evaluated the performance of conventional model-based segmentation against the same reference data to facilitate comparisons. We studied the effect of subject movement on tissue segmentation quality. We also applied the top performing method to a larger out-of-sample dataset and report the quantification results. The results show that multi-atlas image segmentation with diffeomorphic deformation and probabilistic label fusion produces very good quality over all tissues, even for scans with significant quality degradation. The application of our technique to the larger dataset reveals trends of age-related body composition changes that are consistent with the literature. Because of its robustness to subject motion artifacts, our MAIS methodology enables analysis of larger number of scans than conventional state-of-the-art methods. Automated analysis of both soft and hard tissues in pQCT is another contribution of this work.

An Efficient Optimization Approach for Glioma Tumor Segmentation in Brain MRI.

Glioma is an aggressive type of cancer that develops in the brain or spinal cord. Due to many differences in its shape and appearance, accurate segmentation of glioma for identifying all parts of the tumor and its surrounding cancerous tissues is a challenging task. In recent researches, the combination of multi-atlas segmentation and machine learning methods provides robust and accurate results by learning from annotated atlas datasets. To overcome the side effects of limited existing information on atlas-based segmentation, and the long training phase of learning methods, we proposed a semi-supervised unified framework for multi-label segmentation that formulates this problem in terms of a Markov Random Field energy optimization on a parametric graph. To evaluate the proposed framework, we apply it to publicly available BRATS datasets, including low- and high-grade glioma tumors. Experimental results indicate competitive performance compared to the state-of-the-art methods. Compared with the top ranked methods, the proposed framework obtains the best dice score for segmenting of "whole tumor" (WT), "tumor core" (TC ) and "enhancing active tumor" (ET) regions. The achieved accuracy is 94[Formula: see text] characterized by the mean dice score. The motivation of using MRF graph is to map the segmentation problem to an optimization model in a graphical environment. Therefore, by defining perfect graph structure and optimum constraints and flows in the continuous max-flow model, the segmentation is performed precisely.

Automatic macaque brain segmentation based on 7T MRI.

BACKGROUND: In monkey neuroimaging, particularly magnetic resonance imaging (MRI) studies, quick and accurate automatic macaque brain segmentation is essential. However, there are few processing and analysis tools dedicated to automatic brain tissue segmentation and labeling of the macaque brain on a subject-specific basis. As a result, currently most adopted methods are through direct implementation of existing tools that have been designed for human brain. However, the operation steps combining different functional modules of a variety of processing and analysis tool software are inevitably complicated, cumbersome, time-consuming and labor-intensive. NEW METHOD: In this study, we proposed a novel quick and accurate automatic macaque brain segmentation method based on multi-atlas registration and majority-vote algorithm. First, the single-atlas method based on S-HAMMER is used to register each template image of the reference atlas set (including brain tissue labeled images and brain anatomical structure labeled images) to the preprocessed image to be segmented. Thus, we obtain the corresponding deformation field and spatially transform the labeled image, and then get multiple segmentation results by local weighted voting method, which perform label fusion to obtain the final labeled images of brain structures segmentation result. RESULTS: By collecting high SNR and high spatial resolution images of macaque brain images from our 7T human MRI scanner, we have constructed two brain templates for each individual macaque subject, and macaque brain tissues and brain anatomical structure by one-atlas method. However, segmentation result of single-atlas method is not much accurate in some brain tissue area. It takes about 2 h and need more manual correction for segmentation. Automatic segmentation of macaque brain structure based on multi-atlas method was reasonably successful, the accuracy of segmentation was greatly improved without manual correction. Also, the proposed method provided good tissue fitting to V1 with smooth and continuous boundary. The Dice similarity of multi-atlas method showing 3.24%, 4.24%, 2.55%, 2.85%, 3.05%, and 0.35% improvement in image slices of 63, 66, 70, 71, 99 and 100, respectively. The entire processing time for the construction of a single template map took ~40 min. CONCLUSIONS: This study proposed a novel automatic segmentation method of individual macaque brain structure based on multi-atlas registration method, which is concise and reliable. It may offer a valuable tool to applications in the field of brain and neuroscience research using the macaque as an experimental animal model.

Synthesizing pseudo-T2w images to recapture missing data in neonatal neuroimaging with applications in rs-fMRI.

T1- and T2-weighted (T1w and T2w) images are essential for tissue classification and anatomical localization in Magnetic Resonance Imaging (MRI) analyses. However, these anatomical data can be challenging to acquire in non-sedated neonatal cohorts, which are prone to high amplitude movement and display lower tissue contrast than adults. As a result, one of these modalities may be missing or of such poor quality that they cannot be used for accurate image processing, resulting in subject loss. While recent literature attempts to overcome these issues in adult populations using synthetic imaging approaches, evaluation of the efficacy of these methods in pediatric populations and the impact of these techniques in conventional MR analyses has not been performed. In this work, we present two novel methods to generate pseudo-T2w images: the first is based in deep learning and expands upon previous models to 3D imaging without the requirement of paired data, the second is based in nonlinear multi-atlas registration providing a computationally lightweight alternative. We demonstrate the anatomical accuracy of pseudo-T2w images and their efficacy in existing MR processing pipelines in two independent neonatal cohorts. Critically, we show that implementing these pseudo-T2w methods in resting-state functional MRI analyses produces virtually identical functional connectivity results when compared to those resulting from T2w images, confirming their utility in infant MRI studies for salvaging otherwise lost subject data.

Anatomy Recognition in CT Images of Head & Neck Region via Precision Atlases.

Multi-atlas segmentation methods will benefit from atlases covering the complete spectrum of population patterns, while the difficulties in generating such large enough datasets and the computation burden required in the segmentation procedure reduce its practicality in clinical application. In this work, we start from a viewpoint that different parts of the target object can be recognized by different atlases and propose a precision atlas selection strategy. By comparing regional similarity between target image and atlases, precision atlases are ranked and selected by the frequency of regional best match, which have no need to be globally similar to the target subject at either image-level or object-level, largely increasing the implicit patterns contained in the atlas set. In the proposed anatomy recognition method, atlas building is first achieved by all-to-template registration, where the minimum spanning tree (MST) strategy is used to select a registration template from a subset of radiologically near-normal images. Then, a two-stage recognition process is conducted: in rough recognition, sub-image level similarity is calculated between the test image and each image of the whole atlas set, and only the atlas with the highest similarity contributes to the recognition map regionally; in refined recognition, the atlases with the highest frequencies of best match are selected as the precision atlases and are utilized to further increase the accuracy of boundary matching. The proposed method is demonstrated on 298 computed tomography (CT) images and 9 organs in the Head & Neck (H&N) body region. Experimental results illustrate that our method is effective for organs with different segmentation challenge and samples with different image quality, where remarkable improvement in boundary interpretation is made by refined recognition and most objects achieve a localization error within 2 voxels.

Automated Motion Analysis of Bony Joint Structures from Dynamic Computer Tomography Images: A Multi-Atlas Approach.

Dynamic computer tomography (CT) is an emerging modality to analyze in-vivo joint kinematics at the bone level, but it requires manual bone segmentation and, in some instances, landmark identification. The objective of this study is to present an automated workflow for the assessment of three-dimensional in vivo joint kinematics from dynamic musculoskeletal CT images. The proposed method relies on a multi-atlas, multi-label segmentation and landmark propagation framework to extract bony structures and detect anatomical landmarks on the CT dataset. The segmented structures serve as regions of interest for the subsequent motion estimation across the dynamic sequence. The landmarks are propagated across the dynamic sequence for the construction of bone embedded reference frames from which kinematic parameters are estimated. We applied our workflow on dynamic CT images obtained from 15 healthy subjects on two different joints: thumb base (n = 5) and knee (n = 10). The proposed method resulted in segmentation accuracies of 0.90 ± 0.01 for the thumb dataset and 0.94 ± 0.02 for the knee as measured by the Dice score coefficient. In terms of motion estimation, mean differences in cardan angles between the automated algorithm and manual segmentation, and landmark identification performed by an expert were below 1°. Intraclass correlation (ICC) between cardan angles from the algorithm and results from expert manual landmarks ranged from 0.72 to 0.99 for all joints across all axes. The proposed automated method resulted in reproducible and reliable measurements, enabling the assessment of joint kinematics using 4DCT in clinical routine.

SOMA: Subject-, object-, and modality-adapted precision atlas approach for automatic anatomy recognition and delineation in medical images.

PURPOSE: In the multi-atlas segmentation (MAS) method, a large enough atlas set, which can cover the complete spectrum of the whole population pattern of the target object will benefit the segmentation quality. However, the difficulty in obtaining and generating such a large set of atlases and the computational burden required in the segmentation procedure make this approach impractical. In this paper, we propose a method called SOMA to select subject-, object-, and modality-adapted precision atlases for automatic anatomy recognition in medical images with pathology, following the idea that different regions of the target object in a novel image can be recognized by different atlases with regionally best similarity, so that effective atlases have no need to be globally similar to the target subject and also have no need to be overall similar to the target object. METHODS: The SOMA method consists of three main components: atlas building, object recognition, and object delineation. Considering the computational complexity, we utilize an all-to-template strategy to align all images to the same image space belonging to the root image determined by the minimum spanning tree (MST) strategy among a subset of radiologically near-normal images. The object recognition process is composed of two stages: rough recognition and refined recognition. In rough recognition, subimage matching is conducted between the test image and each image of the whole atlas set, and only the atlas corresponding to the best-matched subimage contributes to the recognition map regionally. The frequency of best match for each atlas is recorded by a counter, and the atlases with the highest frequencies are selected as the precision atlases. In refined recognition, only the precision atlases are examined, and the subimage matching is conducted in a nonlocal manner of searching to further increase the accuracy of boundary matching. Delineation is based on a U-net-based deep learning network, where the original gray scale image together with the fuzzy map from refined recognition compose a two-channel input to the network, and the output is a segmentation map of the target object. RESULTS: Experiments are conducted on computed tomography (CT) images with different qualities in two body regions - head and neck (H&N) and thorax, from 298 subjects with nine objects and 241 subjects with six objects, respectively. Most objects achieve a localization error within two voxels after refined recognition, with marked improvement in localization accuracy from rough to refined recognition of 0.6-3 mm in H&N and 0.8-4.9 mm in thorax, and also in delineation accuracy (Dice coefficient) from refined recognition to delineation of 0.01-0.11 in H&N and 0.01-0.18 in thorax. CONCLUSIONS: The SOMA method shows high accuracy and robustness in anatomy recognition and delineation. The improvements from rough to refined recognition and further to delineation, as well as immunity of recognition accuracy to varying image and object qualities, demonstrate the core principles of SOMA where segmentation accuracy increases with precision atlases and gradually refined object matching.

Morphological analysis of subcortical structures for assessment of cognitive dysfunction in Parkinson's disease using multi-atlas based segmentation.

Cognitive impairment in Parkinson's Disease (PD) is the most prevalent non-motor symptom that requires analysis of anatomical associations to cognitive decline in PD. The objective of this study is to analyse the morphological variations of the subcortical structures to assess cognitive dysfunction in PD. In this study, T1 MR images of 58 Healthy Control (HC) and 135 PD subjects categorised as 91 Cognitively normal PD (NC-PD), 25 PD with Mild Cognitive Impairment (PD-MCI) and 19 PD with Dementia (PD-D) subjects, based on cognitive scores are utilised. The 132 anatomical regions are segmented using spatially localized multi-atlas model and volumetric analysis is carried out. The morphological alterations through textural features are captured to differentiate among the HC and PD subjects under different cognitive domains. The volumetric differences in the segmented subcortical structures of accumbens, amygdala, caudate, putamen and thalamus are able to predict cognitive impairment in PD. The volumetric distribution of the subcortical structures in PD-MCI subjects exhibit an overlap with the HC group due to lack of spatial specificity in their atrophy levels. The 3D GLCM features extracted from the significant subcortical structures could discriminate HC, NC-PD, PD-MCI and PD-D subjects with better classification accuracies. The disease related atrophy levels of the subcortical structures captured through morphological analysis provide sensitive evaluation of cognitive impairment in PD.

Validation of separate multi-atlases for auto segmentation of cardiac substructures in CT-scans acquired in deep inspiration breath hold and free breathing.

BACKGROUND AND PURPOSE: Developing NTCP-models for cardiac complications after breast cancer (BC) radiotherapy requires cardiac dose-volume parameters for many patients. These can be obtained by using multi-atlas based automatic segmentation (MABAS) of cardiac structures in planning CT scans. We investigated the relevance of separate multi-atlases for deep inspiration breath hold (DIBH) and free breathing (FB) CT scans. MATERIALS AND METHODS: BC patients scanned in DIBH (n = 10) and in FB (n = 20) were selected to create separate multi-atlases consisting of expert panel delineations of the whole heart, atria and ventricles. The accuracy of atlas-generated contours was validated with expert delineations in independent datasets (n = 10 for DIBH and FB) and reported as Dice coefficients, contour distances and dose-volume differences in relation to interobserver variability of manual contours. Dependency of MABAS contouring accuracy on breathing technique was assessed by validation of a FB atlas in DIBH patients and vice versa (cross-validation). RESULTS: For all structures the FB and DIBH atlases resulted in Dice coefficients with their respective reference contours ≥ 0.8 and average contour distances ≤ 2 mm smaller than slice thickness of (CTs). No significant differences were found for dose-volume parameters in volumes receiving relevant dose levels (WH, LV and RV). Accuracy of the DIBH atlas was at least similar to, and for the ventricles better than, the interobserver variation in manual delineation. Cross-validation between breathing techniques showed a reduced MABAS performance. CONCLUSION: Multi-atlas accuracy was at least similar to interobserver delineation variation. Separate atlases for scans made in DIBH and FB could benefit atlas performance because accuracy depends on breathing technique.

Multi-slice low-rank tensor decomposition based multi-atlas segmentation: Application to automatic pathological liver CT segmentation.

Liver segmentation from abdominal CT images is an essential step for liver cancer computer-aided diagnosis and surgical planning. However, both the accuracy and robustness of existing liver segmentation methods cannot meet the requirements of clinical applications. In particular, for the common clinical cases where the liver tissue contains major pathology, current segmentation methods show poor performance. In this paper, we propose a novel low-rank tensor decomposition (LRTD) based multi-atlas segmentation (MAS) framework that achieves accurate and robust pathological liver segmentation of CT images. Firstly, we propose a multi-slice LRTD scheme to recover the underlying low-rank structure embedded in 3D medical images. It performs the LRTD on small image segments consisting of multiple consecutive image slices. Then, we present an LRTD-based atlas construction method to generate tumor-free liver atlases that mitigates the performance degradation of liver segmentation due to the presence of tumors. Finally, we introduce an LRTD-based MAS algorithm to derive patient-specific liver atlases for each test image, and to achieve accurate pairwise image registration and label propagation. Extensive experiments on three public databases of pathological liver cases validate the effectiveness of the proposed method. Both qualitative and quantitative results demonstrate that, in the presence of major pathology, the proposed method is more accurate and robust than state-of-the-art methods.