Geometric Deformable Model Research Articles

Multi-object 3D segmentation of brain structures using a geometric deformable model with a priori knowledge

ABSTRACT Brain structure segmentation in 3D Magnetic Resonance Images is crucial for understanding neurodegenerative disorders. Manual segmentation is error-prone, necessitating robust automated techniques. In this paper, we introduce a novel and robust approach for the simultaneous segmentation of multiple brain structures in MRI images. Our method involves the concurrent evolution of 3D surfaces toward predefined anatomical targets, employing an efficient multi-object generalized fast marching method (MOGFMM) for simultaneous object detection. Additionally, we propose an effective evolution function that integrates prior knowledge from anatomical and probabilistic atlases, as well as spatial relationships among the segmented structures. Each deformable surface corresponds to a specific structure. To validate our approach, we conducted experiments on a dataset of real brain images (IBSR) and compared the results with several state-of-the-art methods. The obtained results were promising, demonstrating the effectiveness and superiority of our developed method.

Reconstruct Dynamic Soft-Tissue With Stereo Endoscope Based on a Single-Layer Network.

In dynamic minimally invasive surgery environments, 3D reconstruction of deformable soft-tissue surfaces with stereo endoscopic images is very challenging. A simple self-supervised stereo reconstruction framework is proposed to address this issue, which bridges the traditional geometric deformable models and the newly revived neural networks. The equivalence between the classical thin plate spline (TPS) model and a single-layer fully-connected or convolutional network is studied. By alternating training of two TPS equivalent networks within the self-supervised framework, disparity priors are learnt from the past stereo frames of target tissues to form an optimized disparity basis, on which disparity maps of subsequent frames can be estimated more accurately without sacrificing computational efficiency and robustness. The proposed method was verified on stereo-endoscopic videos recorded by the da Vinci® surgical robots.

Automatic brain extraction and hemisphere segmentation in rat brain MR images after stroke using deformable models.

Experimental ischemic stroke models play an essential role in understanding the mechanisms of cerebral ischemia and evaluating the development of pathological extent. An important precursor to the investigation of ischemic strokes associated with rodents is the brain extraction and hemisphere segmentation in rat brain diffusion-weighted imaging (DWI) and T2-weighted MRI (T2WI) images. Accurate and reliable image segmentation tools for extracting the rat brain and hemispheres in the MR images are critical in subsequent processes, such as lesion identification and injury analysis. This study is an attempt to investigate rat brain extraction and hemisphere segmentation algorithms that are practicable in both DWI and T2WI images. To automatically perform brain extraction, the proposed framework is based on an efficient geometric deformable model. By introducing an additional image force in response to the rat brain characteristics into the skull stripping model, we establish a unique rat brain extraction scheme in DWI and T2WI images. For the subsequent hemisphere segmentation, we develop an efficient brain feature detection algorithm to approximately separate the rat brain. A refinement process is enforced by constructing a gradient vector flow in the proximity of the midsurface, where a parametric active contour is attracted to achieve hemisphere segmentation. Extensive experiments with 55 DWI and T2WI subjects were executed in comparison with the state-of-the-art methods. Experimental results indicated that our rat brain extraction and hemisphere segmentation schemes outperformed the competitive methods and exhibited high performance both qualitatively and quantitatively. For rat brain extraction, the average Dice scores were 97.13% and 97.42% in DWI and T2WI image volumes, respectively. Rat hemisphere segmentation results based on the Hausdorff distance metric revealed average values of 0.17 and 0.15mm for DWI and T2WI subjects, respectively. We believe that the established frameworks are advantageous to facilitate preclinical stroke investigation and relevant neuroscience research that requires accurate brain extraction and hemisphere segmentation using rat DWI and T2WI images.

A topology constrained geometric deformable model for medical image segmentation

BackgroundThe analysis and segmentation of Magnetic Resonance images are usually performed by experts. This process is tedious and time consuming. While many automatic and semi-automatic methods have been proposed in the literature, most of them are not reliable or accurate enough for medical image segmentation. Geometric Deformable Models have been widely used as a method of medical image segmentation using only geometric measures. They have proved their efficiency, ease of implementation and ability to handle automatic topology changes. Since, in medical images, the topology of the target to be segmented is known beforehand, different concepts of topology preserving have been implemented along with Geometric Deformable Models as a term to constrain their topology while evolving. MethodIn this paper, we propose a new approach for image segmentation, called the Selective Binary Level Set function and a new variant of the Topology Preserving Selective Binary Level Set model. ResultsThe proposed approach successfully achieves either global or local segmentations. Different metrics were used to evaluate and compare our segmentation results of the heart ventricle to other methods. ConclusionResults show the effectiveness of our proposed method, which proved accurate as a primary step for medical image analysis.

A DEEP LEARNING-BASED CAD SYSTEM FOR RENAL ALLOGRAFT ASSESSMENT: DIFFUSION, BOLD, AND CLINICAL BIOMARKERS.

Recently, studies for non-invasive renal transplant evaluation have been explored to control allograft rejection. In this paper, a computer-aided diagnostic system has been developed to accommodate with an early-stage renal transplant status assessment, called RT-CAD. Our model of this system integrated multiple sources for a more accurate diagnosis: two image-based sources and two clinical-based sources. The image-based sources included apparent diffusion coefficients (ADCs) and the amount of deoxygenated hemoglobin (R2*). More specifically, these ADCs were extracted from 47 diffusion weighted magnetic resonance imaging (DW-MRI) scans at 11 different b-values (b0, b50, b100, …, b1000 s/mm2), while the R2* values were extracted from 30 blood oxygen level-dependent MRI (BOLD-MRI) scans at 5 different echo times (2ms, 7ms, 12ms, 17ms, and 22ms). The clinical sources included serum creatinine (SCr) and creatinine clearance (CrCl). First, the kidney was segmented through the RT-CAD system using a geometric deformable model called a level-set method. Second, both ADCs and R2* were estimated for common patients (N = 30) and then were integrated with the corresponding SCr and CrCl. Last, these integrated biomarkers were considered the discriminatory features to be used as trainers and testers for future deep learning-based classifiers such as stacked auto-encoders (SAEs). We used a k-fold cross-validation criteria to evaluate the RT-CAD system diagnostic performance, which achieved the following scores: 93.3%, 90.0%, and 95.0% in terms of accuracy, sensitivity, and specificity in differentiating between acute renal rejection (AR) and non-rejection (NR). The reliability and completeness of the RT-CAD system was further accepted by the area under the curve score of 0.92. The conclusions ensured that the presented RT-CAD system has a high reliability to diagnose the status of the renal transplant in a non-invasive way.

Segmentation of the lumen and media‐adventitial borders in intravascular ultrasound images using a geometric deformable model

This study presents a geometric deformable model-based segmentation approach to segmentation of the intima and media-adventitial (MA) borders in sequential intravascular ultrasound (IVUS) images. The initial estimation of the vessel borders was done manually only for the first frame of each sequence. After the border initialisation, pre-processing including edge preservation, noise reduction, and dead zone preservation was successively performed on each IVUS frame. To improve segmentation performance, the image masks were determined preliminarily by local binary pattern-based mask initialisation. Then, the inner and outer borders were approximated using a modified distance regularised level set evolution model. The results showed superior performance of the suggested approach for estimating intima and MA layers from the IVUS images. The corresponding correlation coefficients of area, vessel perimeter, maximum vessel diameter, and maximum lumen diameter were r = 0.782, r = 0.716, r = 0.956, and r = 0.874 for the 20 MHz images, respectively, and r = 0.990, r = 0.995, r = 0.989, and r = 0.996 for the 45 MHz images, respectively. In addition, linear regression analysis indicated that the manual segmentation had significantly high similarity at r > 0.967 and r > 0.993 for 20 and 45 MHz images, respectively.

3D kidney segmentation from abdominal diffusion MRI using an appearance-guided deformable boundary.

A new technique for more accurate automatic segmentation of the kidney from its surrounding abdominal structures in diffusion-weighted magnetic resonance imaging (DW-MRI) is presented. This approach combines a new 3D probabilistic shape model of the kidney with a first-order appearance model and fourth-order spatial model of the diffusion-weighted signal intensity to guide the evolution of a 3D geometric deformable model. The probabilistic shape model was built from labeled training datasets to produce a spatially variant, independent random field of region labels. A Markov-Gibbs random field spatial model with up to fourth-order interactions was adequate to capture the inhomogeneity of renal tissues in the DW-MRI signal. A new analytical approach estimated the Gibbs potentials directly from the DW-MRI data to be segmented, in order that the segmentation procedure would be fully automatic. Finally, to better distinguish the kidney object from the surrounding tissues, marginal gray level distributions inside and outside of the deformable boundary were modeled with adaptive linear combinations of discrete Gaussians (first-order appearance model). The approach was tested on a cohort of 64 DW-MRI datasets with b-values ranging from 50 to 1000 s/mm2. The performance of the presented approach was evaluated using leave-one-subject-out cross validation and compared against three other well-known segmentation methods applied to the same DW-MRI data using the following evaluation metrics: 1) the Dice similarity coefficient (DSC); 2) the 95-percentile modified Hausdorff distance (MHD); and 3) the percentage kidney volume difference (PKVD). High performance of the new approach was confirmed by the high DSC (0.95±0.01), low MHD (3.9±0.76) mm, and low PKVD (9.5±2.2)% relative to manual segmentation by an MR expert (a board certified radiologist).

Computer-Aided Diagnostic System for Early Detection of Acute Renal Transplant Rejection Using Diffusion-Weighted MRI

Early diagnosis of acute renal transplant rejection (ARTR) is critical for accurate treatment. Although the current gold standard, diagnostic technique is renal biopsy, it is not preferred due to its invasiveness, long recovery time (1-2 weeks), and potential for complications, e.g., bleeding and/or infection. This paper presents a computer-aided diagnostic (CAD) system for early ARTR detection using (3D + b-value) diffusion-weighted (DW) magnetic resonance imaging (MRI) data. The CAD process starts from kidney tissue segmentation with an evolving geometric (level-set-based) deformable model. The evolution is guided by a voxel-wise stochastic speed function, which follows from a joint kidney-background Markov-Gibbs random field model accounting for an adaptive kidney shape prior and on-going kidney-background visual appearances. A B-spline-based three-dimensional data alignment is employed to handle local deviations due to breathing and heart beating. Then, empirical cumulative distribution functions of apparent diffusion coefficients of the segmented DW-MRI at different b-values are collected as discriminatory transplant status features. Finally, a deep-learning-based classifier with stacked nonnegative constrained autoencoders is employed to distinguish between rejected and nonrejected renal transplants. In our initial "leave-one-subject-out" experiment on 100 subjects, [Formula: see text] of the subjects were correctly classified. The subsequent four-fold and ten-fold cross-validations gave the average accuracy of [Formula: see text] and [Formula: see text], respectively. These results demonstrate the promise of this new CAD system to reliably diagnose renal transplant rejection. The technology presented here can significantly impact the quality of care of renal transplant patients since it has the potential to replace the gold standard in kidney diagnosis, biopsy.

Hair detection and lesion segmentation in dermoscopic images using domain knowledge.

Automated segmentation and dermoscopic hair detection are one of the significant challenges in computer-aided diagnosis (CAD) of melanocytic lesions. Additionally, due to the presence of artifacts and variation in skin texture and smooth lesion boundaries, the accuracy of such methods gets hampered. The objective of this research is to develop an automated hair detection and lesion segmentation algorithm using lesion-specific properties to improve the accuracy. The aforementioned objective is achieved in two ways. Firstly, a novel hair detection algorithm is designed by considering the properties of dermoscopic hair. Second, a novel chroma-based geometric deformable model is used to effectively differentiate the lesion from the surrounding skin. The speed function incorporates the chrominance properties of the lesion to stop evolution at the lesion boundary. Automatic initialization of the initial contour and chrominance-based speed function aids in providing robust and flexible segmentation. The proposed approach is tested on 200 images from PH2 and 900 images from ISBI 2016 datasets. Average accuracy, sensitivity, specificity, and overlap scores of 93.4, 87.6, 95.3, and 11.52% respectively are obtained for the PH2 dataset. Similarly, the proposed method resulted in average accuracy, sensitivity, specificity, and overlap scores of 94.6, 82.4, 97.2, and 7.20% respectively for the ISBI 2016 dataset. Statistical and quantitative analyses prove the reliability of the algorithm for incorporation in CAD systems. Graphical Abstract Overview of proposed system.

Mass Segmentation in Automated 3-D Breast Ultrasound Using Adaptive Region Growing and Supervised Edge-Based Deformable Model.

Automated 3-D breast ultrasound has been proposed as a complementary modality to mammography for early detection of breast cancers. To facilitate the interpretation of these images, computer aided detection systems are being developed in which mass segmentation is an essential component for feature extraction and temporal comparisons. However, automated segmentation of masses is challenging because of the large variety in shape, size, and texture of these 3-D objects. In this paper, the authors aim to develop a computerized segmentation system, which uses a seed position as the only priori of the problem. A two-stage segmentation approach has been proposed incorporating shape information of training masses. At the first stage, a new adaptive region growing algorithm is used to give a rough estimation of the mass boundary. The similarity threshold of the proposed algorithm is determined using a Gaussian mixture model based on the volume and circularity of the training masses. In the second stage, a novel geometric edge-based deformable model is introduced using the result of the first stage as the initial contour. In a data set of 50 masses, including 38 malignant and 12 benign lesions, the proposed segmentation method achieved a mean Dice of 0.74 ± 0.19 which outperformed the adaptive region growing with a mean Dice of 0.65 ± 0.2 (p-value < 0.02). Moreover, the resulting mean Dice was significantly (p-value < 0.001) better than that of the distance regularized level set evolution method (0.52 ± 0.27). The supervised method presented in this paper achieved accurate mass segmentation results in terms of Dice measure. The suggested segmentation method can be utilized in two aspects: 1) to automatically measure the change in volume of breast lesions over time and 2) to extract features for a computer aided detection or diagnosis system.

Statistical analysis of ADCs and clinical biomarkers in detecting acute renal transplant rejection

The main goal of this study is to determine which parameters [e.g. clinical biomarkers, demographics and image-markers using 4D (3D + b-value) diffusion-weighted MRI (DW-MRI)] are more correlated with transplanted kidney status in patients who have undergone kidney transplantation, and can be used for early assessment of acute renal rejection. The study included 16 patients with stable graft function and 37 patients with acute rejection (AR), determined by renal biopsy post-transplantation. 3D DW-MRI of each allograft had been acquired using a series of b-values 50 and 100-1000 in steps of 100 smm-2. The kidney was automatically segmented and co-aligned across series for motion correction using geometric deformable models. Volume-averaged apparent diffusion coefficients (ADCs) at each b-value were calculated. All possible subsets of ADC were used, along with patient age, sex, serum plasma creatinine (SPCr) and creatinine clearance (CrCl), as predictors in 211 logistic regression models where AR was the outcome variable. Predictive value of ADC at each b-value was assessed using its Akaike weight. ANOVA of the saturated model found that odds of AR depended significantly on SPCr, CrCl and ADC at b = 500, 600, 700 and 900 smm-2. The model incorporating ADC at b = 100 and700 smm-2 had the lowest value of the Akaike information criterion; the same two b-values also had the greatest Akaike weights. For comparison, the top 10 submodels and the full model were reported. Preliminary findings suggest that ADC provides improved detection of AR than lab values alone. At least two non-zero gradient strengths should be used for optimal results. Advances in knowledge: This paper investigated possible correlations between image-based and clinical biomarkers, and the fusion of both with respect to biopsy diagnosis of AR.

A comprehensive non-invasive framework for diagnosing prostate cancer

Early detection of prostate cancer increases chances of patients' survival. Our automated non-invasive system for computer-aided diagnosis (CAD) of prostate cancer segments the prostate on diffusion-weighted magnetic resonance images (DW-MRI) acquired at different b-values, estimates its apparent diffusion coefficients (ADC), and classifies their descriptors – empirical cumulative distribution functions (CDF) – with a trained deep learning network. To segment the prostate, an evolving geometric (level-set-based) deformable model is guided by a speed function depending on intensity attributes extracted from the DW-MRI with nonnegative matrix factorization (NMF). For a more robust evolution, the attributes are fused with a probabilistic shape prior and estimated spatial dependencies between prostate voxels. To preserve continuity, the ADCs of the segmented prostate volume at different b-values are normalized and refined using a generalized Gauss-Markov random field image model. The CDFs of the refined ADCs at different b-values are considered global water diffusion features and used to distinguish between benign and malignant prostates. A deep learning network of stacked non-negativity-constrained auto-encoders (SNCAE) is trained to classify the benign or malignant prostates on the basis of the constructed CDFs. Our experiments on 53 clinical DW-MRI data sets resulted in 92.3% accuracy, 83.3% sensitivity, and 100% specificity, indicating that the proposed CAD system could be used as a reliable non-invasive diagnostic tool.

Age-related changes of thoracic aorta geometry used to predict the risk for acute type B dissection

AimsRisk models that use a single aortic diameter threshold have failed to successfully predict acute type B aortic dissection (TBAD). We sought to identify meaningful age-indexed anatomical variables to predict TBAD risk. Methods and resultsA geometric deformable model, consisting of virtual elastic balloons that inflate inside a vessel lumen, was developed to quantify thoracic aorta geometry. In the presence of TBAD, true and total artery lumen morphology were assessed. A stepwise logistic model was built to predict TBAD risk. Initial covariates included age, gender, body mass index and all anatomic variables not directly related to the dissected segment. Patients with acute TBAD (n=34, 62±12years old, 57% male gender) were compared with subjects with symptoms of dissection, but with a subsequent negative diagnosis (n=51, 62±12years old, 76% male gender). Patient risk factors did not differ between groups. Most aortic anatomical variables were age-dependent. Aortic size was larger in every segment of the dissected with respect to non-dissected aortas (p<0.001). Variables entering the TBAD risk prediction model were aortic arch diameter, thoracic aorta length and age (predictability=0.9764, r=0.85), confirmed by a bootstrap internal validation. In dissected aortas, the true lumen volume was correlated to age (r=0.72). ConclusionsTBAD probability increases with a larger aortic arch diameter and a longer thoracic aorta, whereas threshold values increase with age. The aortic morphology was age-dependent. After dissection, true lumen volume correlated to age. The use of threshold values indexed to age should be encouraged to better prevent and eventually treat TBAD.

Multiple-object geometric deformable model for segmentation of macular OCT: errata

Multiple-object geometric deformable model for segmentation of macular OCT: errata

Segmentation of the Cerebellar Peduncles Using a Random Forest Classifier and a Multi-object Geometric Deformable Model: Application to Spinocerebellar Ataxia Type 6.

The cerebellar peduncles, comprising the superior cerebellar peduncles (SCPs), the middle cerebellar peduncle (MCP), and the inferior cerebellar peduncles (ICPs), are white matter tracts that connect the cerebellum to other parts of the central nervous system. Methods for automatic segmentation and quantification of the cerebellar peduncles are needed for objectively and efficiently studying their structure and function. Diffusion tensor imaging (DTI) provides key information to support this goal, but it remains challenging because the tensors change dramatically in the decussation of the SCPs (dSCP), the region where the SCPs cross. This paper presents an automatic method for segmenting the cerebellar peduncles, including the dSCP. The method uses volumetric segmentation concepts based on extracted DTI features. The dSCP and noncrossing portions of the peduncles are modeled as separate objects, and are initially classified using a random forest classifier together with the DTI features. To obtain geometrically correct results, a multi-object geometric deformable model is used to refine the random forest classification. The method was evaluated using a leave-one-out cross-validation on five control subjects and four patients with spinocerebellar ataxia type 6 (SCA6). It was then used to evaluate group differences in the peduncles in a population of 32 controls and 11 SCA6 patients. In the SCA6 group, we have observed significant decreases in the volumes of the dSCP and the ICPs and significant increases in the mean diffusivity in the noncrossing SCPs, the MCP, and the ICPs. These results are consistent with a degeneration of the cerebellar peduncles in SCA6 patients.

Water flow based geometric active deformable model for road network

A width and color based geometric active deformable model is proposed for road network extraction from remote sensing images with minimal human interception. Orientation and width of road are computed from a single manual seed point, from which the propagation starts both right and left hand directions of the starting point, which extracts the interconnected road network from the aerial or high spatial resolution satellite image automatically. Here the propagation (like water flow in canal with defined boundary) is restricted with color and width of the road. Road extraction is done for linear, curvilinear (U shape and S shape) roads first, irrespective of width and color. Then, this algorithm is improved to extract road with junctions in a shape of L, T and X along with center line. Roads with small break or disconnected roads are also extracts by a modified version of this same algorithm. This methodology is tested and evaluated with various remote sensing images. The experimental results show that the proposed method is efficient and extracting roads accurately with less computation time. However, in complex urban areas, the identification accuracy declines due to the various sizes of obstacles, over bridges, multilane etc.

Multiple-object geometric deformable model for segmentation of macular OCT.

Optical coherence tomography (OCT) is the de facto standard imaging modality for ophthalmological assessment of retinal eye disease, and is of increasing importance in the study of neurological disorders. Quantification of the thicknesses of various retinal layers within the macular cube provides unique diagnostic insights for many diseases, but the capability for automatic segmentation and quantification remains quite limited. While manual segmentation has been used for many scientific studies, it is extremely time consuming and is subject to intra- and inter-rater variation. This paper presents a new computational domain, referred to as flat space, and a segmentation method for specific retinal layers in the macular cube using a recently developed deformable model approach for multiple objects. The framework maintains object relationships and topology while preventing overlaps and gaps. The algorithm segments eight retinal layers over the whole macular cube, where each boundary is defined with subvoxel precision. Evaluation of the method on single-eye OCT scans from 37 subjects, each with manual ground truth, shows improvement over a state-of-the-art method.

Biomedical image segmentation using geometric deformable models and metaheuristics

This paper describes a hybrid level set approach for medical image segmentation. This new geometric deformable model combines region- and edge-based information with the prior shape knowledge introduced using deformable registration. Our proposal consists of two phases: training and test. The former implies the learning of the level set parameters by means of a Genetic Algorithm, while the latter is the proper segmentation, where another metaheuristic, in this case Scatter Search, derives the shape prior. In an experimental comparison, this approach has shown a better performance than a number of state-of-the-art methods when segmenting anatomical structures from different biomedical image modalities.

A comprehensive non‐invasive framework for automated evaluation of acute renal transplant rejection using DCE‐MRI

The objective was to develop a novel and automated comprehensive framework for the non-invasive identification and classification of kidney non-rejection and acute rejection transplants using 2D dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). The proposed approach consists of four steps. First, kidney objects are segmented from the surrounding structures with a geometric deformable model. Second, a non-rigid registration approach is employed to account for any local kidney deformation. In the third step, the cortex of the kidney is extracted in order to determine dynamic agent delivery, since it is the cortex that is primarily affected by the perfusion deficits that underlie the pathophysiology of acute rejection. Finally, we use an analytical function-based model to fit the dynamic contrast agent kinetic curves in order to determine possible rejection candidates. Five features that map the data from the original data space to the feature space are chosen with a k-nearest-neighbor (KNN) classifier to distinguish between acute rejection and non-rejection transplants. Our study includes 50 transplant patients divided into two groups: 27 patients with stable kidney function and the remainder with impaired kidney function. All of the patients underwent DCE-MRI, while the patients in the impaired group also underwent ultrasound-guided fine needle biopsy. We extracted the kidney objects and the renal cortex from DCE-MRI for accurate medical evaluation with an accuracy of 0.97 ± 0.02 and 0.90 ± 0.03, respectively, using the Dice similarity metric. In a cohort of 50 participants, our framework classified all cases correctly (100%) as rejection or non-rejection transplant candidates, which is comparable to the gold standard of biopsy but without the associated deleterious side-effects. Both the 95% confidence interval (CI) statistic and the receiver operating characteristic (ROC) analysis document the ability to separate rejection and non-rejection groups. The average plateau (AP) signal magnitude and the gamma-variate model functional parameter α have the best individual discriminating characteristics.

SEGMENTATION OF THE COMPLETE SUPERIOR CEREBELLAR PEDUNCLES USING A MULTI-OBJECT GEOMETRIC DEFORMABLE MODEL.

The superior cerebellar peduncles (SCPs) are white matter tracts that serve as the major efferent pathways from the cerebellum to the thalamus. With diffusion tensor images (DTI), tractography algorithms or volumetric segmentation methods have been able to reconstruct part of the SCPs. However, when the fibers cross, the primary eigenvector (PEV) no longer represents the primary diffusion direction. Therefore, at the crossing of the left and right SCP, known as the decussation of the SCPs (dSCP), fiber tracts propagate incorrectly. To our knowledge, previous methods have not been able to segment the SCPs correctly. In this work, we explore the diffusion properties and seek to volumetrically segment the complete SCPs. The non-crossing SCPs and dSCP are modeled as different objects. A multi-object geometric deformable model is employed to define the boundaries of each piece of the SCPs, with the forces derived from diffusion properties as well as the PEV. We tested our method on a software phantom and real subjects. Results indicate that our method is able to the resolve the crossing and segment the complete SCPs with repeatability.