Pathological lung segmentation in chest CT images based on improved random walker.

BACKGROUND AND OBJECTIVE: Pathological lung segmentation as a pretreatment step in the diagnosis of lung diseases has been widely explored. Because of the complexity of pathological lung structures and gray blur of the border, accurate lung segmentation in clinical 3D computed tomography images is a challenging task. In view of the current situation, the work proposes a fast and accurate pathological lung segmentation method. The following contributions have been made: First, the edge weights introduce spatial information and clustering information, so that walkers can use more image information during walking. Second, a Gaussian Distribution of seed point set is established to further expand the possibility of selection between fake seed points and real seed points. Finally, the pre-parameter is calculated using original seed points, and the final results are fitted with new seed points. METHODS: This study proposes a segmentation method based on an improved random walker algorithm. The proposed method consists of the following steps: First, a gray value is used as the sample distribution. Gaussian mixture model is used to obtain the clustering probability of an image. Thus, the spatial distance and clustering result are added as new weights, and the new edge weights are used to construct a random walker map. Second, a large number of marked points are automatically selected, and the intermediate results are obtained from the newly constructed map and retained only as pre-parameters. When new seed points are introduced, the probability value of the walker is quickly calculated from the new parameters and pre-parameters, and the final segmentation result can be obtained. RESULTS: The proposed method was tested on 65 sets of CT cases. Quantitative evaluation with different methods confirms the high accuracy on our dataset (98.55%) and LOLA11 dataset (97.41%). Similarly, the average segmentation time (10.5s) is faster than random walker (1,332.5s). CONCLUSIONS: The comparison of the experimental results show that the proposed method can accurately and quickly obtain pathological lung processing results. Therefore, it has potential clinical applications.

Segmentation of the placenta and its vascular tree in Doppler ultrasound for fetal surgery planning.

PURPOSE: Twin-to-twin transfusion syndrome (TTTS) is a serious condition that occurs in about 10-15% of monochorionic twin pregnancies. In most instances, the blood flow is unevenly distributed throughout the placenta anastomoses leading to the death of both fetuses if no surgical procedure is performed. Fetoscopic laser coagulation is the optimal therapy to considerably improve co-twin prognosis by clogging the abnormal anastomoses. Notwithstanding progress in recent years, TTTS surgery is highly risky. Computer-assisted planning of the intervention can thus improve the outcome. METHODS: In this work, we implement a GPU-accelerated random walker (RW) algorithm to detect the placenta, both umbilical cord insertions and the placental vasculature from Doppler ultrasound (US). Placenta and background seeds are manually initialized in 10-20 slices (out of 245). Vessels are automatically initialized in the same slices by means of Otsu thresholding. The RW finds the boundaries of the placenta and reconstructs the vasculature. RESULTS: We evaluate our semiautomatic method in 5 monochorionic and 24 singleton pregnancies. Although satisfactory performance is achieved on placenta segmentation (Dice ≥ 84.0%), some vascular connections are still neglected due to the presence of US reverberation artifacts (Dice ≥ 56.9%). We also compared inter-user variability and obtained Dice coefficients of ≥ 76.8% and ≥ 97.42% for placenta and vasculature, respectively. After a 3-min manual initialization, our GPU approach speeds the computation 10.6 times compared to the CPU. CONCLUSIONS: Our semiautomatic method provides a near real-time user experience and requires short training without compromising the segmentation accuracy. A powerful approach is thus presented to rapidly plan the fetoscope insertion point ahead of TTTS surgery.

A self-supervised strategy for fully automatic segmentation of renal dynamic contrast-enhanced magnetic resonance images.

PURPOSE: An automated accurate segmentation for dynamic contrast-enhanced magnetic resonance (DCE-MR) image sequences is essential for quantification of renal function. A self-supervised strategy is proposed for fully automatic segmentation of the renal DCE-MR images without using manually labeled data. METHODS: The proposed strategy employed both temporal and spatial information of the DCE-MR image sequences. First, the kidney area, the seed regions of the cortex, the medulla, and the pelvis were automatically detected in the spatial domain. Subsequently, all the pixels in the kidney were automatically labeled as the cortex, the medulla, or the pelvis based on their time-intensity signal and spatial position using a supervised classifier. The feasibility of the proposed strategy was verified on a dataset of renal DCE-MR images of 14 subjects without history of kidney diseases. Furthermore, the self-supervised strategy and the commonly used traditional unsupervised method were quantitatively compared with a reference manual segmentation by an experienced radiologist, using similarity indexes. RESULTS: The average Dice coefficient (ADC) for the segmentations of the proposed self-supervised method is 0.92 using a ransom walker model as the classifier or 0.86 using a K-nearest neighbor model as the classifier. The ADC of the Kmeans-based unsupervised methods with three and six clusters were 0.65 and 0.79, respectively. The Dice coefficients of the self-supervised method were remarkably higher than that of the unsupervised method (one-tailed paired-sample t-test, P-values <10-3 ). CONCLUSIONS: The results indicate that the proposed self-supervised approach yields a satisfactory similarity with the reference manual segmentation. Compared with the traditional unsupervised clustering method, the new strategy does not require manual intervention during the segmentation process and achieves better results for the segmentation of renal DCE-MR images.

A fully automatic approach for multimodal PET and MR image segmentation in gamma knife treatment planning.

BACKGROUND AND OBJECTIVES: Nowadays, clinical practice in Gamma Knife treatments is generally based on MRI anatomical information alone. However, the joint use of MRI and PET images can be useful for considering both anatomical and metabolic information about the lesion to be treated. In this paper we present a co-segmentation method to integrate the segmented Biological Target Volume (BTV), using [11C]-Methionine-PET (MET-PET) images, and the segmented Gross Target Volume (GTV), on the respective co-registered MR images. The resulting volume gives enhanced brain tumor information to be used in stereotactic neuro-radiosurgery treatment planning. GTV often does not match entirely with BTV, which provides metabolic information about brain lesions. For this reason, PET imaging is valuable and it could be used to provide complementary information useful for treatment planning. In this way, BTV can be used to modify GTV, enhancing Clinical Target Volume (CTV) delineation. METHODS: A novel fully automatic multimodal PET/MRI segmentation method for Leksell Gamma Knife® treatments is proposed. This approach improves and combines two computer-assisted and operator-independent single modality methods, previously developed and validated, to segment BTV and GTV from PET and MR images, respectively. In addition, the GTV is utilized to combine the superior contrast of PET images with the higher spatial resolution of MRI, obtaining a new BTV, called BTVMRI. A total of 19 brain metastatic tumors, undergone stereotactic neuro-radiosurgery, were retrospectively analyzed. A framework for the evaluation of multimodal PET/MRI segmentation is also presented. Overlap-based and spatial distance-based metrics were considered to quantify similarity concerning PET and MRI segmentation approaches. Statistics was also included to measure correlation among the different segmentation processes. Since it is not possible to define a gold-standard CTV according to both MRI and PET images without treatment response assessment, the feasibility and the clinical value of BTV integration in Gamma Knife treatment planning were considered. Therefore, a qualitative evaluation was carried out by three experienced clinicians. RESULTS: The achieved experimental results showed that GTV and BTV segmentations are statistically correlated (Spearman's rank correlation coefficient: 0.898) but they have low similarity degree (average Dice Similarity Coefficient: 61.87 ± 14.64). Therefore, volume measurements as well as evaluation metrics values demonstrated that MRI and PET convey different but complementary imaging information. GTV and BTV could be combined to enhance treatment planning. In more than 50% of cases the CTV was strongly or moderately conditioned by metabolic imaging. Especially, BTVMRI enhanced the CTV more accurately than BTV in 25% of cases. CONCLUSIONS: The proposed fully automatic multimodal PET/MRI segmentation method is a valid operator-independent methodology helping the clinicians to define a CTV that includes both metabolic and morphologic information. BTVMRI and GTV should be considered for a comprehensive treatment planning.

A graph-based approach for spatio-temporal segmentation of coronary arteries in X-ray angiographic sequences.

The segmentation and tracking of coronary arteries (CAs) are critical steps for the computation of biophysical measurements in pediatric interventional cardiology. In the literature, most methods are focused on either segmenting the vessel lumen or on tracking the vessel centerline. However, they do not simultaneously combine the segmentation and tracking of a specific CA. This paper introduces a novel algorithm for CA segmentation and tracking from 2D X-ray angiography sequences. The proposed algorithm is based on the Temporal Vessel Walker (TVW) segmentation method, which combines graph-based formulation and temporal priors. Moreover, superpixel groups are used by TVW as image primitives to ensure a better extraction of the CA. The proposed algorithm, TVW with superpixels (SP-TVW), returns an accurate result to segment and track the artery along the angiogram. Quantitative results over 12 sequences of young patients show the accuracy of the proposed framework. The results return a mean recall of 84% in the dataset. In addition, the proposed method returned a Dice index of 70% in segmenting and tracking right coronary arteries and circumflex arteries. The performance of the proposed method surpasses the existing polyline method in tracking the centerline of CA with a more precise localization of the centerline, resulting in a smaller distance error of 0.23mm compared to 0.94mm.

Automatic liver tumor segmentation on computed tomography for patient treatment planning and monitoring.

Segmentation of liver tumors from Computed Tomography (CT) and tumor burden analysis play an important role in the choice of therapeutic strategies for liver diseases and treatment monitoring. In this paper, a new segmentation method for liver tumors from contrast-enhanced CT imaging is proposed. As manual segmentation of tumors for liver treatment planning is both labor intensive and time-consuming, a highly accurate automatic tumor segmentation is desired. The proposed framework is fully automatic requiring no user interaction. The proposed segmentation evaluated on real-world clinical data from patients is based on a hybrid method integrating cuckoo optimization and fuzzy c-means algorithm with random walkers algorithm. The accuracy of the proposed method was validated using a clinical liver dataset containing one of the highest numbers of tumors utilized for liver tumor segmentation containing 127 tumors in total with further validation of the results by a consultant radiologist. The proposed method was able to achieve one of the highest accuracies reported in the literature for liver tumor segmentation compared to other segmentation methods with a mean overlap error of 22.78 % and dice similarity coefficient of 0.75 in 3Dircadb dataset and a mean overlap error of 15.61 % and dice similarity coefficient of 0.81 in MIDAS dataset. The proposed method was able to outperform most other tumor segmentation methods reported in the literature while representing an overlap error improvement of 6 % compared to one of the best performing automatic methods in the literature. The proposed framework was able to provide consistently accurate results considering the number of tumors and the variations in tumor contrast enhancements and tumor appearances while the tumor burden was estimated with a mean error of 0.84 % in 3Dircadb dataset.

Automatic liver segmentation on Computed Tomography using random walkers for treatment planning.

Segmentation of the liver from Computed Tomography (CT) volumes plays an important role during the choice of treatment strategies for liver diseases. Despite lots of attention, liver segmentation remains a challenging task due to the lack of visible edges on most boundaries of the liver coupled with high variability of both intensity patterns and anatomical appearances with all these difficulties becoming more prominent in pathological livers. To achieve a more accurate segmentation, a random walker based framework is proposed that can segment contrast-enhanced livers CT images with great accuracy and speed. Based on the location of the right lung lobe, the liver dome is automatically detected thus eliminating the need for manual initialization. The computational requirements are further minimized utilizing rib-caged area segmentation, the liver is then extracted by utilizing random walker method. The proposed method was able to achieve one of the highest accuracies reported in the literature against a mixed healthy and pathological liver dataset compared to other segmentation methods with an overlap error of 4.47 % and dice similarity coefficient of 0.94 while it showed exceptional accuracy on segmenting the pathological livers with an overlap error of 5.95 % and dice similarity coefficient of 0.91.

3-D Ultrasound Segmentation of the Placenta Using the Random Walker Algorithm: Reliability and Agreement.

Volumetric segmentation of the placenta using 3-D ultrasound is currently performed clinically to investigate correlation between organ volume and fetal outcome or pathology. Previously, interpolative or semi-automatic contour-based methodologies were used to provide volumetric results. We describe the validation of an original random walker (RW)-based algorithm against manual segmentation and an existing semi-automated method, virtual organ computer-aided analysis (VOCAL), using initialization time, inter- and intra-observer variability of volumetric measurements and quantification accuracy (with respect to manual segmentation) as metrics of success. Both semi-automatic methods require initialization. Therefore, the first experiment compared initialization times. Initialization was timed by one observer using 20 subjects. This revealed significant differences (p < 0.001) in time taken to initialize the VOCAL method compared with the RW method. In the second experiment, 10 subjects were used to analyze intra-/inter-observer variability between two observers. Bland-Altman plots were used to analyze variability combined with intra- and inter-observer variability measured by intra-class correlation coefficients, which were reported for all three methods. Intra-class correlation coefficient values for intra-observer variability were higher for the RW method than for VOCAL, and both were similar to manual segmentation. Inter-observer variability was 0.94 (0.88, 0.97), 0.91 (0.81, 0.95) and 0.80 (0.61, 0.90) for manual, RW and VOCAL, respectively. Finally, a third observer with no prior ultrasound experience was introduced and volumetric differences from manual segmentation were reported. Dice similarity coefficients for observers 1, 2 and 3 were respectively 0.84 ± 0.12, 0.94 ± 0.08 and 0.84 ± 0.11, and the mean was 0.87 ± 0.13. The RW algorithm was found to provide results concordant with those for manual segmentation and to outperform VOCAL in aspects of observer reliability. The training of an additional untrained observer was investigated, and results revealed that with the appropriate initialization protocol, results for observers with varying levels of experience were concordant. We found that with appropriate training, the RW method can be used for fast, repeatable 3-D measurement of placental volume.

Semi-automatic segmentation for 3D motion analysis of the tongue with dynamic MRI.

Dynamic MRI has been widely used to track the motion of the tongue and measure its internal deformation during speech and swallowing. Accurate segmentation of the tongue is a prerequisite step to define the target boundary and constrain the tracking to tissue points within the tongue. Segmentation of 2D slices or 3D volumes is challenging because of the large number of slices and time frames involved in the segmentation, as well as the incorporation of numerous local deformations that occur throughout the tongue during motion. In this paper, we propose a semi-automatic approach to segment 3D dynamic MRI of the tongue. The algorithm steps include seeding a few slices at one time frame, propagating seeds to the same slices at different time frames using deformable registration, and random walker segmentation based on these seed positions. This method was validated on the tongue of five normal subjects carrying out the same speech task with multi-slice 2D dynamic cine-MR images obtained at three orthogonal orientations and 26 time frames. The resulting semi-automatic segmentations of a total of 130 volumes showed an average dice similarity coefficient (DSC) score of 0.92 with less segmented volume variability between time frames than in manual segmentations.

Endocardial surface delineation in 3-D transesophageal echocardiography.

We describe and compare several methods for recovering endocardial walls from 3-D transesophageal echocardiography (3-D TEE), which can help with diagnostics or providing input into biomechanical models. We employ a segmentation method based on 3-D level sets that maximizes enclosed volume while minimizing surface area and uses a growth inhibition function that includes 3-D gradient magnitude (to locate the endocardial walls) and a thin tissue detector (for the mitral valve leaflets). We also study delineation using a graph cut method that performs automated seeding by leveraging a fast radial symmetry transform to determine a central axis along which the 3-D volume is warped into a cylindrical coordinate space. Finally, a random walker approach is also used for automated delineation. The methods are used to estimate clinically relevant cardiovascular volumetric parameters such as stroke volume and left ventricular ejection fraction. Experiments are performed on clinical data collected from patients undergoing cardiothoracic surgery. Performance evaluation includes comparisons of the automated delineations against expert-defined ground truth using a number of error metrics, as well as errors between automatically computed and expert-derived physiologic parameters.

SEMI-AUTOMATIC SEGMENTATION OF THE TONGUE FOR 3D MOTION ANALYSIS WITH DYNAMIC MRI.

Accurate segmentation is an important preprocessing step for measuring the internal deformation of the tongue during speech and swallowing using 3D dynamic MRI. In an MRI stack, manual segmentation of every 2D slice and time frame is time-consuming due to the large number of volumes captured over the entire task cycle. In this paper, we propose a semi-automatic segmentation workflow for processing 3D dynamic MRI of the tongue. The steps comprise seeding a few slices, seed propagation by deformable registration, random walker segmentation of the temporal stack of images and 3D super-resolution volumes. This method was validated on the tongue of two subjects carrying out the same speech task with multi-slice 2D dynamic cine-MR images obtained at three orthogonal orientations and 26 time frames. The resulting semi-automatic segmentations of 52 volumes showed an average dice similarity coefficient (DSC) score of 0.9 with reduced segmented volume variability compared to manual segmentations.