Intravascular ultrasound image segmentation: a three-dimensional fast-marching method based on gray level distributions

Scale Mutualized Perception for Vessel Border Detection in Intravascular Ultrasound Images

Vessel border detection in IVUS images is essential for coronary disease diagnosis. It helps to obtain the clinical indices on the inner vessel morphology to indicate the stenosis. However, the existing methods suffer the challenge of scale-dependent interference. Early methods usually rely on the hand-crafted features, thus not robust to this interference. The existing deep learning methods are also ineffective to solve this challenge, because these methods aggregate multi-scale features in the top-down way. This aggregation may bring in interference from the non-adjacent scale. Besides, they only combine the features in all scales, and thus may weaken their complementary information. We propose the scale mutualized perception to solve this challenge by considering the adjacent scales mutually to preserve their complementary information. First, the adjacent small scales contain certain semantics to locate different vessel tissues. Then, they can also perceive the global context to assist the representation of the local context in the adjacent large scale, and vice versa. It helps to distinguish the objects with similar local features. Second, the adjacent large scales provide detailed information to refine the vessel boundaries. The experiments show the effectiveness of our method in 153 IVUS sequences, and its superiority to ten state-of-the-art methods.

3D pyramidal densely connected network with cross-frame uncertainty guidance for intravascular ultrasound sequence segmentation

Objective. Automatic extraction of external elastic membrane border (EEM) and lumen-intima border (LIB) in intravascular ultrasound (IVUS) sequences aids atherosclerosis diagnosis. Existing IVUS segmentation networks ignored longitudinal relations among sequential images and neglected that IVUS images of different vascular conditions vary largely in intricacy and informativeness. As a result, they suffered from performance degradation in complicated parts in IVUS sequences.Approach. In this paper, we develop a 3D Pyramidal Densely-connected Network (PDN) with Adaptive learning and post-Correction guided by a novel cross-frame uncertainty (CFU). The proposed method is named PDN-AC. Specifically, the PDN enables the longitudinal information exploitation and the effective perception of size-varied vessel regions in IVUS samples, by pyramidally connecting multi-scale 3D dilated convolutions. Additionally, the CFU enhances the robustness of the method to complicated pathology from the frame-level (f-CFU) and pixel-level (p-CFU) via exploiting cross-frame knowledge in IVUS sequences. The f-CFU weighs the complexity of IVUS frames and steers an adaptive sampling during the PDN training. The p-CFU visualizes uncertain pixels probably misclassified by the PDN and guides an active contour-based post-correction.Main results. Human and animal experiments were conducted on IVUS datasets acquired from atherosclerosis patients and pigs. Results showed that the f-CFU weighted adaptive sampling reduced the Hausdorff distance (HD) by 10.53%/7.69% in EEM/LIB detection. Improvements achieved by the p-CFU guided post-correction were 2.94%/5.56%.Significance. The PDN-AC attained mean Jaccard values of 0.90/0.87 and HD values of 0.33/0.34 mm in EEM/LIB detection, preferable to state-of-the-art IVUS segmentation methods.

AwCPM-Net: A Collaborative Constraint GAN for 3D Coronary Artery Reconstruction in Intravascular Ultrasound Sequences

3D coronary artery reconstruction (3D-CAR) in intravascular ultrasound (IVUS) sequences allows quantitative analyses of vessel properties. Existing methods treat two main tasks of the 3D-CAR separately, including the cardiac phase retrieval (CPR) and the membrane border extraction (MBE). They ignore the CPR-MBE connection that could achieve mutual promotions to both tasks. In this paper, we pioneer to achieve one-step 3D-CAR via a collaborative constraint generative adversarial network (GAN) named the AwCPM-Net. The AwCPM-Net consists of a dual-task collaborative generator and a dual-task constraint discriminator. The generator combines a self-supervised CPR branch with a semi-supervised MBE branch via a warming-up connection. The discriminator promotes dual-branch predictions simultaneously. The CPR branch requires no annotations and outputs inter-frame deformation fields used for identifying cardiac phases. Deformation fields are additionally constrained by the MBE branch and the discriminator. The MBE branch predicts membrane boundaries for each frame. Two aspects assist the semi-supervised segmentation: annotation augmentation by deformation fields of the CPR branch; information exploitation on unlabeled images enabled by GAN design. Trained and tested on an IVUS dataset acquired from atherosclerosis patients, the AwCPM-Net is effective in both CPR and MBE tasks, superior to state-of-the-art IVUS CPR or MBE methods. Hence, the AwCPM-Net reconstructs reliable 3D artery anatomy in the IVUS modality.

Extraction of media adventitia and luminal intima borders by reconstructing intravascular ultrasound image sequences with vascular structural continuity

PURPOSE

Most published methods directly achieve vessel membrane border detection on cross-sectional intravascular ultrasound (IVUS) images. The vascular structural continuity that exists in entire IVUS image sequences has been overlooked. However, this continuity can have a helpful role in the delineation of vessel membrane contours. To achieve the vessel membrane segmentation more effectively through employing this continuity, a strategy, referred to as multiangle reconstruction, segmentation, and recovery (RSR), is proposed in this paper.

METHODS

Four main steps are contained in the multiangle-RSR: first, a combination of sampling and interpolation is employed to reconstruct long-axis-model IVUS frames, in which continuity information becomes available. Second, a clustering algorithm is conducted on long-axis-model IVUS frames to roughly extract the media-adventitia (MA) and lumen-intima (LI) boundaries. Third, the segmentation results of cross-sectional IVUS frames are recovered based on the rough results of long-axis-model IVUS frames, and an optimization process that combines downsampling, fitting and smoothing is designed to reduce the interference of bifurcation and side vessels.

RESULTS

Multiangle-RSR is tested on a public dataset, and the Hausdorff distance (HD), Jaccard measure (JM), and percentage of area difference (PAD) are utilized as quantitative evaluation metrics. Mean HDs of 0.34 and 0.29 mm are obtained for MA border detection and LI border detection, respectively, which decrease by 43.3% and 9.4%, respectively, compared with their counterparts in previously published approaches. Furthermore, the mean JM is 0.87 for both MA border detection and LI border detection. The mean PADs of the MA contour extraction and the LI contour extraction are 0.10 and 0.11, respectively.

CONCLUSION

The results indicate that the proposed strategy effectively introduces vascular structural continuity by reconstructing long-axis-model IVUS frames and achieves more precise extraction of MA and LI borders.

Automatic Detection of Atherosclerotic Plaque and Calcification From Intravascular Ultrasound Images by Using Deep Convolutional Neural Networks

Atherosclerosis is the major cause of cardiovascular diseases (CVDs). Intravascular ultrasound (IVUS) is a common imaging modality for diagnosing CVDs. However, an efficient analyzer for IVUS image segmentation is required for assisting cardiologists. In this study, an end-to-end deep-learning convolutional neural network was developed for automatically detecting media-adventitia borders, luminal regions, and calcified plaque in IVUS images. A total of 713 grayscale IVUS images from 18 patients were used as training data for the proposed deep-learning model. The model is constructed using the three modified U-Nets and combined with the concept of cascaded networks to prevent errors in the detection of calcification owing to the interference of pixels outside the plaque regions. Three loss functions (Dice, Tversky, and focal loss) with various characteristics were tested to determine the best setting for the proposed model. The efficacy of the deep-learning model was evaluated by analyzing precision-recall curve. The average precision (AP), Dice score coefficient, precision, sensitivity, and specificity of the predicted and ground truth results were then compared. All training processes were validated using leave-one-subject-out cross-validation. The experimental results showed that the proposed deep-learning model exhibits high performance in segmenting the media-adventitia layers and luminal regions for all loss functions, with all tested metrics being higher than 0.90. For locating calcified tissues, the best result was obtained when the focal loss function was applied to the proposed model, with an AP of 0.73; however, the prediction efficacy was affected by the proportion of calcified tissues within the plaque region when the focal loss function was employed. Compared with commercial software, the proposed method exhibited high accuracy in segmenting IVUS images in some special cases, such as when shadow artifacts or side vessels surrounded the target vessel.

Supervised machine learning for coronary artery lumen segmentation in intravascular ultrasound images

Intravascular ultrasound (IVUS) has been widely used to capture cross sectional lumen frames of inner wall of coronary arteries. This kind of medical imaging modalities is capable of providing detailed and significant information of lumen contour shape, which is very important for clinical diagnosis and analysis of cardiovascular diseases. Numerous learning based techniques have recently become very popular for coronary artery segmentation due to their impressive results. In this work, a supervised machine learning method for coronary artery lumen segmentation with high accuracy and minimal user interaction is designed. The fully discriminative lumen segmentation method jointly learning a classifier the weak learners rely on and the features of the classifier is developed. Additionally, the theoretical supports of the Gradient Boosting framework used in this work and its quadratic approximation are presented. The proposed algorithm is tested on the public datasets of boundary detection of lumen in IVUS challenge held in MICCAI 2011 and achieves a higher average Jaccard similarity of 96.8% and a lower mean error distance of 0.55 (in Cartesian coordinates), which shows higher accuracy compared to the existing learning based methods. Moreover, three real patient IVUS datasets are used to evaluate the performance of the proposed coronary artery lumen segmentation algorithm, which is shown to achieve lower percent error of lumen area of 1.861% ± 0.965%, 1.968% ± 0.864%, and 1.671% ± 0.584%, respectively, compared to the manually measured lumen area (ground truth). The proposed lumen segmentation method is found to be superior to the latest learning based segmentation techniques. Given the efficiency and robustness, our method has great potential in IVUS images processing and coronary artery segmentation and quantification. NOVELTY STATEMENT: The main contributions are summarized in the following aspects: A detailed review of related work about learning based coronary artery lumen segmentation in intravascular ultrasound images is presented. A fully discriminative lumen segmentation method jointly learning a classifier our weak learners rely on and the features of the classifier is developed. The theoretical supports of the Gradient Boosting framework and its quadratic approximation used in this work are presented.

Vessel membrane segmentation and calcification location in intravascular ultrasound images using a region detector and an effective selection strategy

BACKGROUND AND OBJECTIVE

Segmenting vessel membranes and locating the calcific region in intravascular ultrasound (IVUS) images aid physicians in the diagnosis of atherosclerosis. However, the manual extraction of the media adventitia (MA)/lumen border and calcification location are cumbersome due to the excessive number of IVUS frames. Moreover, most existing (semi-)automatic detection methods cannot achieve both vessel membrane extraction and calcification location simultaneously, and they are unable to detect vessel membranes in IVUS frames from different acquisition systems.

METHOD

A fully automatic approach is proposed based on extremal regions and a flexible selection strategy to extract vessel membranes in different IVUS frames and locate the calcific region in high-frequency ones. Three main steps are included in the algorithm. First, a region detector is employed to extract extremal regions from an IVUS image. Then, according to the selection strategy, a part of the extracted regions is selected. At the same time, the calcification is located according to its special acoustic properties. Next, approximate MA and lumen border segmentation is achieved based on the selected extremal regions and the located calcification in polar coordinates. Finally, the final segmentation results are obtained by smoothing the approximate values.

RESULT

To demonstrate the feasibility of the method, it was evaluated based on a standard public dataset. Furthermore, to quantitatively evaluate the segmentation performance, the Hausdorff distance (HD), Jaccard measure (JM) and percentage of area difference (PAD) were used. The results show that a mean HD of 1.13/1.21 mm, a mean JM of 0.83/0.77 and a mean PAD of 0.11/0.23 are achieved for MA/lumen border detection in 77 40-MHz IVUS images. For MA/lumen border extraction in 435 20-MHz IVUS frames, the average HD, JM and PAD values are 0.47/0.28 mm, 0.84/0.89 and 0.13/0.10, respectively. In addition, the approach successfully achieves calcification location in 40-MHz IVUS frames. In comparison with other published methods, the method proposed in this study is competitive.

CONCLUSION

According to these results, our strategy can extract MA/lumen borders in different IVUS frames and effectively locate calcification in high-frequency IVUS frames.

Multi-Grid Phase Field Skin Tumor Segmentation in 3D Ultrasound Images

The aim of this paper is to present a new method for skin tumor segmentation in the 3D ultrasound images. We consider a variational formulation, the energy of which combines a diffuse interface phase field model (regularization term) and a log-likelihood computed using nonparametric estimates (data attachment term). We propose a multi-grid implementation with the exact solutions which has the advantage to avoid space discretization and numerical instabilities. The resulting algorithm is simple and easy to implement in multi-dimensions. Concerning applications, we focus on skin tumor segmentation. The clinical dataset used for the experiments is composed of 12 images with the ground truth given by a dermatologist. Comparisons with the reference methods show that the proposed method is more robust to the choice of the volume initialization. Moreover, thanks to the flexibility introduced by the diffuse interface, the sensitivity increases by 12% if the initialization is inside the lesion, and the Dice index increases by 59%, if the initialization covers the entire lesion. These results show that this new method is well designed to tackle the problem of underestimation of tumor volumes.

Robust segmentation of arterial walls in intravascular ultrasound images using Dual Path U-Net

A Fully Convolutional Network (FCN) based deep architecture called Dual Path U-Net (DPU-Net) is proposed for automatic segmentation of the lumen and media-adventitia in IntraVascular UltraSound (IVUS) frames, which is crucial for diagnosis of many cardiovascular diseases and also for facilitating 3D reconstructions of human arteries. One of the most prevalent problems in medical image analysis is the lack of training data. To overcome this limitation, we propose a twofold solution. First, we introduce a deep architecture that is able to learn using a small number of training images and still achieves a high degree of generalization ability. Second, we strengthen the proposed DPU-Net by having a real-time augmentor control the image augmentation process. Our real-time augmentor contains specially-designed operations that simulate three types of IVUS artifacts and integrate them into the training images. We exhaustively assessed our twofold contribution over Balocco's standard publicly available IVUS 20 MHz and 40 MHz B-mode dataset, which contain 109 training image, 326 test images and 19 training images, 59 test images, respectively. Models are trained from scratch with the training images provided and evaluated with two commonly used metrics in the IVUS segmentation literature, namely Jaccard Measure (JM) and Hausdorff Distance (HD). Experimental results show that DPU-Net achieves 0.87 JM, 0.82 mm HD and 0.86 JM, 1.07 mm HD over 40 MHz dataset for segmenting the lumen and the media, respectively. Also, DPU-Net achieves 0.90 JM, 0.25 mm HD and 0.92 JM, 0.30 mm HD over 20 MHz images for segmenting the lumen and the media, respectively. In addition, DPU-Net outperforms existing methods by 8-15% in terms of HD distance. DPU-Net also shows a strong generalization property for predicting images in the test sets that contain a significant amount of major artifacts such as bifurcations, shadows, and side branches that are not common in the training set. Furthermore, DPU-Net runs within 0.03 s to segment each frame with a single modern GPU (Nvidia GTX 1080). The proposed work leverages modern deep learning-based method for segmentation of lumen and the media vessel walls in both 20 MHz and 40 MHz IVUS B-mode images and achieves state-of-the-art results without any manual intervention. The code is available online at https://github.com/Kulbear/IVUS-Ultrasonic.

Detecting the Media-adventitia Border in Intravascular Ultrasound Images through a Classification-based Approach

The detection of the media-adventitia (MA) border in intravascular ultrasound (IVUS) images is essential for vessel assessment and disease diagnosis. However, it remains a challenging task, considering the existence of plaque, calcification, and various artifacts. In this article, an effective method based on classification is proposed to extract the MA border in IVUS images. First, a novel morphologic feature describing the relative position of each structure relative to the MA border, called RPES for short, is proposed. Then, the RPES feature and other features are employed in a multiclass extreme learning machine (ELM) to classify IVUS images into nine classes including the MA border and other structures. At last, a modified snake model is employed to effectively detect the MA border in the rectangular domain, in which a modified external force field is constructed on the basis of local border appearances and classification results. The proposed method is evaluated on a public dataset with 77 IVUS images by three indicators in eight situations, such as calcification and a guide wire artifact. With the proposed RPES feature, detection performances are improved by more than 39 percent, which shows an apparent advantage in comparative experiments. Furthermore, compared with two other existing methods used on the same dataset, the proposed method achieves 18 of the best indicators among 24, demonstrating its higher capability in detecting the MA border.

Optimal surface segmentation with convex priors in irregularly sampled space

Optimal surface segmentation is a state-of-the-art method used for segmentation of multiple globally optimal surfaces in volumetric datasets. The method is widely used in numerous medical image segmentation applications. However, nodes in the graph based optimal surface segmentation method typically encode uniformly distributed orthogonal voxels of the volume. Thus the segmentation cannot attain an accuracy greater than a single unit voxel, i.e. the distance between two adjoining nodes in graph space. Segmentation accuracy higher than a unit voxel is achievable by exploiting partial volume information in the voxels which shall result in non-equidistant spacing between adjoining graph nodes. This paper reports a generalized graph based multiple surface segmentation method with convex priors which can optimally segment the target surfaces in an irregularly sampled space. The proposed method allows non-equidistant spacing between the adjoining graph nodes to achieve subvoxel segmentation accuracy by utilizing the partial volume information in the voxels. The partial volume information in the voxels is exploited by computing a displacement field from the original volume data to identify the subvoxel-accurate centers within each voxel resulting in non-equidistant spacing between the adjoining graph nodes. The smoothness of each surface modeled as a convex constraint governs the connectivity and regularity of the surface. We employ an edge-based graph representation to incorporate the necessary constraints and the globally optimal solution is obtained by computing a minimum s-t cut. The proposed method was validated on 10 intravascular multi-frame ultrasound image datasets for subvoxel segmentation accuracy. In all cases, the approach yielded highly accurate results. Our approach can be readily extended to higher-dimensional segmentations.

Fractional-order Darwinian PSO-based feature selection for media-adventitia border detection in intravascular ultrasound images

Media-adventitia (MA) border delineates the outer appearance of arterial wall in intravascular ultrasound (IVUS) image. The detection of MA border is a challenging topic due to many difficulties such as complicated intravascular structures, intrinsic artifacts and image noises. We propose a classification-based MA border detection method with an embedded feature selection technique. The feature selection technique is based on Fractional-order Darwinian particle swarm optimization (FODPSO) algorithm. By employing feature selection, 293-dimension features including multi-scale features, gray-scale features and morphological feature are reducing to 37-dimension. The border detection method with feature selection is tested on a public dataset extracted from in-vivo pullbacks of human coronary arteries, which contains 77 IVUS images. Three indicators, Jaccard (JACC), Hausdorff Distance (HD) and Percentage of Area Difference (PAD), are measured for quantitative evaluation. Detection with 293-dimension features obtains JACC 0.79, HD 1.41 and PAD 0.16, while detection with 37-dimension features obtains JACC 0.83, HD 1.27 and PAD 0.12, indicating that the FODPSO-based feature selection method improves MA border detection by JACC 0.04, HD 0.14 and PAD 0.04. Furthermore, the proposed border detection method acquires better performances compared with two other automatic methods conducted on the same dataset available in literature.

Anatomical Structure Segmentation in Ultrasound Volumes Using Cross Frame Belief Propagating Iterative Random Walks

Ultrasound (US) is widely used as a low-cost alternative to computed tomography or magnetic resonance and primarily for preliminary imaging. Since speckle intensity in US images is inherently stochastic, readers are often challenged in their ability to identify the pathological regions in a volume of a large number of images. This paper introduces a generalized approach for volumetric segmentation of structures in US images and volumes. We employ an iterative random walks (IRW) solver, a random forest learning model, and a gradient vector flow (GVF) based interframe belief propagation technique for achieving cross-frame volumetric segmentation. At the start, a weak estimate of the tissue structure is obtained using estimates of parameters of a statistical mechanics model of US tissue interaction. Ensemble learning of these parameters further using a random forest is used to initialize the segmentation pipeline. IRW is used for correcting the contour in various steps of the algorithm. Subsequently, a GVF-based interframe belief propagation is applied to adjacent frames based on the initialization of contour using information in the current frame to segment the complete volume by frame-wise processing. We have experimentally evaluated our approach using two different datasets. Intravascular ultrasound (IVUS) segmentation was evaluated using 10 pullbacks acquired at 20 MHz and thyroid US segmentation is evaluated on 16 volumes acquired at [Formula: see text] MHz. Our approach obtains a Jaccard score of [Formula: see text] for IVUS segmentation and [Formula: see text] for thyroid segmentation while processing each frame in [Formula: see text] for the IVUS and in [Formula: see text] for thyroid segmentation without the need of any computing accelerators such as GPUs.

Segmentation of arterial walls in intravascular ultrasound cross-sectional images using extremal region selection

Intravascular Ultrasound (IVUS) is an intra-operative imaging modality that facilitates observing and appraising the vessel wall structure of the human coronary arteries. Segmentation of arterial wall boundaries from the IVUS images is not only crucial for quantitative analysis of the vessel walls and plaque characteristics, but is also necessary for generating 3D reconstructed models of the artery. The aim of this study is twofold. Firstly, we investigate the feasibility of using a recently proposed region detector, namely Extremal Region of Extremum Level (EREL) to delineate the luminal and media-adventitia borders in IVUS frames acquired by 20 MHz probes. Secondly, we propose a region selection strategy to label two ERELs as lumen and media based on the stability of their textural information. We extensively evaluated our selection strategy on the test set of a standard publicly available dataset containing 326 IVUS B-mode images. We showed that in the best case, the average Hausdorff Distances (HD) between the extracted ERELs and the actual lumen and media were 0.22  mm and 0.45 mm, respectively. The results of our experiments revealed that our selection strategy was able to segment the lumen with ⩽0.3 mm HD to the gold standard even though the images contained major artifacts such as bifurcations, shadows, and side branches. Moreover, when there was no artifact, our proposed method was able to delineate media-adventitia boundaries with 0.31 mm HD to the gold standard. Furthermore, our proposed segmentation method runs in time that is linear in the number of pixels in each frame. Based on the results of this work, by using a 20 MHz IVUS probe with controlled pullback, not only can we now analyze the internal structure of human arteries more accurately, but also segment each frame during the pullback procedure because of the low run time of our proposed segmentation method.

Automated localization and segmentation techniques for B-mode ultrasound images: A review

B-mode ultrasound imaging is used extensively in medicine. Hence, there is a need to have efficient segmentation tools to aid in computer-aided diagnosis, image-guided interventions, and therapy. This paper presents a comprehensive review on automated localization and segmentation techniques for B-mode ultrasound images. The paper first describes the general characteristics of B-mode ultrasound images. Then insight on the localization and segmentation of tissues is provided, both in the case in which the organ/tissue localization provides the final segmentation and in the case in which a two-step segmentation process is needed, due to the desired boundaries being too fine to locate from within the entire ultrasound frame. Subsequenly, examples of some main techniques found in literature are shown, including but not limited to shape priors, superpixel and classification, local pixel statistics, active contours, edge-tracking, dynamic programming, and data mining. Ten selected applications (abdomen/kidney, breast, cardiology, thyroid, liver, vascular, musculoskeletal, obstetrics, gynecology, prostate) are then investigated in depth, and the performances of a few specific applications are compared. In conclusion, future perspectives for B-mode based segmentation, such as the integration of RF information, the employment of higher frequency probes when possible, the focus on completely automatic algorithms, and the increase in available data are discussed.

Visualization of plaque distribution in a curved artery: three-dimensional intravascular ultrasound imaging

Intravascular ultrasound (IVUS) imaging provides an excellent tool for evaluation of the type, morphology, extent, and severity of an atheromatous plaque. 3 D IVUS imaging offers additive information pertaining to morphology of the arterial structures and volumetric plaque distributions. A new 3 D IVUS visualization technique was developed to provide 3 D structural information of a curved artery. A virtual 3 D curved arterial phantom consisting of varying cross-sectional shapes, wall thicknesses, and acoustic intensity information was utilized to validate the nonlinear interpolation technique to create intermediary 2 D IVUS images. IVUS imaging was performed for the iliofemoral arterial segment of an atherosclerotic Yucatan miniswine model. These in-vivo IVUS data were utilized for intermediary IVUS image generation and volumetric 3 D IVUS visualization. Smooth transitional changes of cross-sectional shape, wall thickness and grayscale intensity were found between the intermediary images and the original arterial phantom slices. The 3 D IVUS imaging of the unfolded curved iliofemoral artery provided realistic 3 D luminal surface images of the arteries with physiologic grayscale intensity information. This unique 3 D IVUS imaging technique may help with assessment of 3 D plaque distribution across the curved arterial structure, and improve 3 D visualization of atheromatous components.

Automatic segmentation of the lumen region in intravascular images of the coronary artery

Image assessment of the arterial system plays an important role in the diagnosis of cardiovascular diseases. The segmentation of the lumen and media-adventitia in intravascular (IVUS) images of the coronary artery is the first step towards the evaluation of the morphology of the vessel under analysis and the identification of possible atherosclerotic lesions. In this study, a fully automatic method for the segmentation of the lumen in IVUS images of the coronary artery is presented. The proposed method relies on the K-means algorithm and the mean roundness to identify the region corresponding to the potential lumen. An approach to identify and eliminate side branches on bifurcations is also proposed to delimit the area with the potential lumen regions. Additionally, an active contour model is applied to refine the contour of the lumen region. In order to evaluate the segmentation accuracy, the results of the proposed method were compared against manual delineations made by two experts in 326 IVUS images of the coronary artery. The average values of the Jaccard measure, Hausdorff distance, percentage of area difference and Dice coefficient were 0.88 ± 0.06, 0.29 ± 0.17  mm, 0.09 ± 0.07 and 0.94 ± 0.04, respectively, in 324 IVUS images successfully segmented. Additionally, a comparison with the studies found in the literature showed that the proposed method is slight better than the majority of the related methods that have been proposed. Hence, the new automatic segmentation method is shown to be effective in detecting the lumen in IVUS images without using complex solutions and user interaction.

Segmentation of Skin Tumors in High-Frequency 3-D Ultrasound Images

High-frequency 3-D ultrasound imaging is an informative tool for diagnosis, surgery planning and skin lesion examination. The purpose of this article was to describe a semi-automated segmentation tool providing easy access to the extent, shape and volume of a lesion. We propose an adaptive log-likelihood level-set segmentation procedure using non-parametric estimates of the intensity distribution. The algorithm has a single parameter to control the smoothness of the contour, and we describe how a fixed value yields satisfactory segmentation results with an average Dice coefficient of D = 0.76. The algorithm is implemented on a grid, which increases the speed by a factor of 100 compared with a standard pixelwise segmentation. We compare the method with parametric methods making the hypothesis of Rayleigh or Nakagami distributed signals, and illustrate that our method has greater robustness with similar computational speed. Benchmarks are made on realistic synthetic ultrasound images and a data set of nine clinical 3-D images acquired with a 50-MHz imaging system. The proposed algorithm is suitable for use in a clinical context as a post-processing tool.

Multi-scale texture-based level-set segmentation of breast B-mode images

Automatic segmentation of ultrasonographic breast lesions is very challenging, due to the lesions' spiculated nature and the variance in shape and texture of the B-mode ultrasound images. Many studies have tried to answer this challenge by applying a variety of computational methods including: Markov random field, artificial neural networks, and active contours and level-set techniques. These studies focused on creating an automatic contour, with maximal resemblance to a manual contour, delineated by a trained radiologist. In this study, we have developed an algorithm, designed to capture the spiculated boundary of the lesion by using the properties from the corresponding ultrasonic image. This is primarily achieved through a unique multi-scale texture identifier (inspired by visual system models) integrated in a level-set framework. The algorithm׳s performance has been evaluated quantitatively via contour-based and region-based error metrics. We compared the algorithm-generated contour to a manual contour delineated by an expert radiologist. In addition, we suggest here a new method for performance evaluation where corrections made by the radiologist replace the algorithm-generated (original) result in the correction zones. The resulting corrected contour is then compared to the original version. The evaluation showed: (1) Mean absolute error of 0.5 pixels between the original and the corrected contour; (2) Overlapping area of 99.2% between the lesion regions, obtained by the algorithm and the corrected contour. These results are significantly better than those previously reported. In addition, we have examined the potential of our segmentation results to contribute to the discrimination between malignant and benign lesions.

Reliable and Accurate Calcium Volume Measurement in Coronary Artery Using Intravascular Ultrasound Videos

Quantitative assessment of calcified atherosclerotic volume within the coronary artery wall is vital for cardiac interventional procedures. The goal of this study is to automatically measure the calcium volume, given the borders of coronary vessel wall for all the frames of the intravascular ultrasound (IVUS) video. Three soft computing fuzzy classification techniques were adapted namely Fuzzy c-Means (FCM), K-means, and Hidden Markov Random Field (HMRF) for automated segmentation of calcium regions and volume computation. These methods were benchmarked against previously developed threshold-based method. IVUS image data sets (around 30,600 IVUS frames) from 15 patients were collected using 40 MHz IVUS catheter (Atlantis® SR Pro, Boston Scientific®, pullback speed of 0.5 mm/s). Calcium mean volume for FCM, K-means, HMRF and threshold-based method were 37.84 ± 17.38 mm(3), 27.79 ± 10.94 mm(3), 46.44 ± 19.13 mm(3) and 35.92 ± 16.44 mm(3) respectively. Cross-correlation, Jaccard Index and Dice Similarity were highest between FCM and threshold-based method: 0.99, 0.92 ± 0.02 and 0.95 + 0.02 respectively. Student's t-test, z-test and Wilcoxon-test are also performed to demonstrate consistency, reliability and accuracy of the results. Given the vessel wall region, the system reliably and automatically measures the calcium volume in IVUS videos. Further, we validated our system against a trained expert using scoring: K-means showed the best performance with an accuracy of 92.80%. Out procedure and protocol is along the line with method previously published clinically.

Methods for 2-D and 3-D Endobronchial Ultrasound Image Segmentation

Endobronchial ultrasound (EBUS) is now commonly used for cancer-staging bronchoscopy. Unfortunately, EBUS is challenging to use and interpreting EBUS video sequences is difficult. Other ultrasound imaging domains, hampered by related difficulties, have benefited from computer-based image-segmentation methods. Yet, so far, no such methods have been proposed for EBUS. We propose image-segmentation methods for 2-D EBUS frames and 3-D EBUS sequences. Our 2-D method adapts the fast-marching level-set process, anisotropic diffusion, and region growing to the problem of segmenting 2-D EBUS frames. Our 3-D method builds upon the 2-D method while also incorporating the geodesic level-set process for segmenting EBUS sequences. Tests with lung-cancer patient data showed that the methods ran fully automatically for nearly 80% of test cases. For the remaining cases, the only user-interaction required was the selection of a seed point. When compared to ground-truth segmentations, the 2-D method achieved an overall Dice index = 90.0% ±4.9%, while the 3-D method achieved an overall Dice index = 83.9 ± 6.0%. In addition, the computation time (2-D, 0.070 s/frame; 3-D, 0.088 s/frame) was two orders of magnitude faster than interactive contour definition. Finally, we demonstrate the potential of the methods for EBUS localization in a multimodal image-guided bronchoscopy system.

Automated Framework for Detecting Lumen and Media-Adventitia Borders in Intravascular Ultrasound Images

An automated framework for detecting lumen and media-adventitia borders in intravascular ultrasound images was developed on the basis of an adaptive region-growing method and an unsupervised clustering method. To demonstrate the capability of the framework, linear regression, Bland-Altman analysis and distance analysis were used to quantitatively investigate the correlation, agreement and spatial distance, respectively, between our detected borders and manually traced borders in 337 intravascular ultrasound images in vivo acquired from six patients. The results of these investigations revealed good correlation (r = 0.99), good agreement (>96.82% of results within the 95% confidence interval) and small average distance errors (lumen border: 0.08 mm, media-adventitia border: 0.10 mm) between the borders generated by the automated framework and the manual tracing method. The proposed framework was found to be effective in detecting lumen and media-adventitia borders in intravascular ultrasound images, indicating its potential for use in routine studies of vascular disease.

Computer Vision Techniques for Transcatheter Intervention

Minimally invasive transcatheter technologies have demonstrated substantial promise for the diagnosis and the treatment of cardiovascular diseases. For example, transcatheter aortic valve implantation is an alternative to aortic valve replacement for the treatment of severe aortic stenosis, and transcatheter atrial fibrillation ablation is widely used for the treatment and the cure of atrial fibrillation. In addition, catheter-based intravascular ultrasound and optical coherence tomography imaging of coronary arteries provides important information about the coronary lumen, wall, and plaque characteristics. Qualitative and quantitative analysis of these cross-sectional image data will be beneficial to the evaluation and the treatment of coronary artery diseases such as atherosclerosis. In all the phases (preoperative, intraoperative, and postoperative) during the transcatheter intervention procedure, computer vision techniques (e.g., image segmentation and motion tracking) have been largely applied in the field to accomplish tasks like annulus measurement, valve selection, catheter placement control, and vessel centerline extraction. This provides beneficial guidance for the clinicians in surgical planning, disease diagnosis, and treatment assessment. In this paper, we present a systematical review on these state-of-the-art methods. We aim to give a comprehensive overview for researchers in the area of computer vision on the subject of transcatheter intervention. Research in medical computing is multi-disciplinary due to its nature, and hence, it is important to understand the application domain, clinical background, and imaging modality, so that methods and quantitative measurements derived from analyzing the imaging data are appropriate and meaningful. We thus provide an overview on the background information of the transcatheter intervention procedures, as well as a review of the computer vision techniques and methodologies applied in this area.

Value of C-Arm Computed Tomography to Evaluate Stent Deployment During Femoro-Popliteal Revascularization

PURPOSE

To compare the accuracy of C-arm computed tomography (CT) and digital subtraction angiography (DSA) in detecting incomplete stent expansion (ISE) after superficial femoral artery (SFA) stenting using intravascular ultrasound (IVUS) as a gold standard.

MATERIALS

Fifty patients with symptomatic SFA occlusive disease requiring angioplasty were prospectively included. Once technical success (<30 % residual stenosis) was obtained on post-procedural DSA, C-arm CT and IVUS were acquired. DSA and C-arm CT examinations were reviewed by 2 investigators and correlated with IVUS. C-arm CT image quality was rated on a four-point scale. Doppler ultrasound was performed at 1-year follow-up.

RESULTS

The ankle-brachial index was 0.69 ± 0.10 and 0.99 ± 0.40, respectively, pre- and post-procedure. C-arm CT imaging quality was rated as good or excellent in 80%. In-stent minimal luminal diameter (MLD) was evaluated at 4.71 ± 0.7 mm on DSA, 3.39 ± 0.6 mm on IVUS, and 3.12 ± 0.9 mm on C-arm CT. Compared to IVUS, DSA demonstrated an overestimation of MLD (p = 0.0001), an underestimation of ISE (DSA = 18.8% ± 7.6; IVUS = 29.8% ± 9) (p < 0.0001), and a poor inter-technique intra-class correlation coefficient (ICC = 0.24). No difference was observed between IVUS and C-arm CT in ISE as calculated by diameter (29.8 ± 9 vs. 28.2 ± 12.5%, p = 0.5) and area (30.2 ± 8.4 vs. 33.3 ± 9.5%, p = 0.2). Inter-technique ICC between C-arm CT and IVUS was 0.72 [95%CI 0.49; 0.85] for MLA measurements. The inter-observer ICC for MLD and MLA measurements on C-arm CT were, respectively, estimated at 0.75 [95% CI 0.40, 0.89] and 0.77 [95% CI 0.43, 0.90)].

CONCLUSIONS

C-arm CT presents a better correlation with IVUS than DSA to determine lumen diameter and ISE immediately after percutaneous revascularization.

Segmentation of tumor ultrasound image in HIFU therapy based on texture and boundary encoding

It is crucial in high intensity focused ultrasound (HIFU) therapy to detect the tumor precisely with less manual intervention for enhancing the therapy efficiency. Ultrasound image segmentation becomes a difficult task due to signal attenuation, speckle effect and shadows. This paper presents an unsupervised approach based on texture and boundary encoding customized for ultrasound image segmentation in HIFU therapy. The approach oversegments the ultrasound image into some small regions, which are merged by using the principle of minimum description length (MDL) afterwards. Small regions belonging to the same tumor are clustered as they preserve similar texture features. The mergence is completed by obtaining the shortest coding length from encoding textures and boundaries of these regions in the clustering process. The tumor region is finally selected from merged regions by a proposed algorithm without manual interaction. The performance of the method is tested on 50 uterine fibroid ultrasound images from HIFU guiding transducers. The segmentations are compared with manual delineations to verify its feasibility. The quantitative evaluation with HIFU images shows that the mean true positive of the approach is 93.53%, the mean false positive is 4.06%, the mean similarity is 89.92%, the mean norm Hausdorff distance is 3.62% and the mean norm maximum average distance is 0.57%. The experiments validate that the proposed method can achieve favorable segmentation without manual initialization and effectively handle the poor quality of the ultrasound guidance image in HIFU therapy, which indicates that the approach is applicable in HIFU therapy.

Automated detection framework of the calcified plaque with acoustic shadowing in IVUS images

Intravascular Ultrasound (IVUS) is one ultrasonic imaging technology to acquire vascular cross-sectional images for the visualization of the inner vessel structure. This technique has been widely used for the diagnosis and treatment of coronary artery diseases. The detection of the calcified plaque with acoustic shadowing in IVUS images plays a vital role in the quantitative analysis of atheromatous plaques. The conventional method of the calcium detection is manual drawing by the doctors. However, it is very time-consuming, and with high inter-observer and intra-observer variability between different doctors. Therefore, the computer-aided detection of the calcified plaque is highly desired. In this paper, an automated method is proposed to detect the calcified plaque with acoustic shadowing in IVUS images by the Rayleigh mixture model, the Markov random field, the graph searching method and the prior knowledge about the calcified plaque. The performance of our method was evaluated over 996 in-vivo IVUS images acquired from eight patients, and the detected calcified plaques are compared with manually detected calcified plaques by one cardiology doctor. The experimental results are quantitatively analyzed separately by three evaluation methods, the test of the sensitivity and specificity, the linear regression and the Bland-Altman analysis. The first method is used to evaluate the ability to distinguish between IVUS images with and without the calcified plaque, and the latter two methods can respectively measure the correlation and the agreement between our results and manual drawing results for locating the calcified plaque in the IVUS image. High sensitivity (94.68%) and specificity (95.82%), good correlation and agreement (>96.82% results fall within the 95% confidence interval in the Student t-test) demonstrate the effectiveness of the proposed method in the detection of the calcified plaque with acoustic shadowing in IVUS images.

Endovascular shear strain elastography for the detection and characterization of the severity of atherosclerotic plaques: in vitro validation and in vivo evaluation

This work explores the potential of shear strain elastograms to identify vulnerable atherosclerotic plaques. The Lagrangian speckle model estimator (LSME) elasticity imaging method was further developed to estimate shear strain elasticity (SSE). Three polyvinyl alcohol cryogel vessel phantoms were imaged with an intravascular ultrasound (IVUS) scanner. The estimated SSE maps were validated against finite-element results. Atherosclerosis was induced in carotid arteries of eight Sinclair mini-pigs using a combination of surgical techniques, diabetes and a high-fat diet. IVUS images were acquired in vivo in 14 plaques before euthanasia and histology. All plaques were characterized by high magnitudes in SSE maps that correlated with American Heart Association atherosclerosis stage classifications (r = 0.97, p < 0.001): the worse the plaque condition the higher was the absolute value of SSE, i.e. |SSE| (e.g., mean |SSE| was 3.70 ± 0.40% in Type V plaques, whereas it was reduced to 0.11 ± 0.01% in normal walls). This study indicates the feasibility of using SSE to highlight atherosclerotic plaque vulnerability characteristics.

A 3-D ultrasound imaging robotic system to detect and quantify lower limb arterial stenoses: in vivo feasibility

The degree of stenosis is the most common criterion used to assess the severity of lower limb peripheral arterial disease. Two-dimensional ultrasound (US) imaging is the first-line diagnostic method for investigating lesions, but it cannot render a 3-D map of the entire lower limb vascular tree required for therapy planning. We propose a prototype 3-D US imaging robotic system that can potentially reconstruct arteries from the iliac in the lower abdomen down to the popliteal behind the knee. A realistic multi-modal vascular phantom was first conceptualized to evaluate the system's performance. Geometric accuracies were assessed in surface reconstruction and cross-sectional area in comparison to computed tomography angiography (CTA). A mean surface map error of 0.55 mm was recorded for 3-D US vessel representations, and cross-sectional lumen areas were congruent with CTA geometry. In the phantom study, stenotic lesions were properly localized and severe stenoses up to 98.3% were evaluated with -3.6 to 11.8% errors. The feasibility of the in vivo system in reconstructing the normal femoral artery segment of a volunteer and detecting stenoses on a femoral segment of a patient was also investigated and compared with that of CTA. Together, these results encourage future developments to increase the robot's potential to adequately represent lower limb vessels and clinically evaluate stenotic lesions for therapy planning and recurrent non-invasive and non-ionizing follow-up examinations.

Standardized evaluation methodology and reference database for evaluating IVUS image segmentation

This paper describes an evaluation framework that allows a standardized and quantitative comparison of IVUS lumen and media segmentation algorithms. This framework has been introduced at the MICCAI 2011 Computing and Visualization for (Intra)Vascular Imaging (CVII) workshop, comparing the results of eight teams that participated. We describe the available data-base comprising of multi-center, multi-vendor and multi-frequency IVUS datasets, their acquisition, the creation of the reference standard and the evaluation measures. The approaches address segmentation of the lumen, the media, or both borders; semi- or fully-automatic operation; and 2-D vs. 3-D methodology. Three performance measures for quantitative analysis have been proposed. The results of the evaluation indicate that segmentation of the vessel lumen and media is possible with an accuracy that is comparable to manual annotation when semi-automatic methods are used, as well as encouraging results can be obtained also in case of fully-automatic segmentation. The analysis performed in this paper also highlights the challenges in IVUS segmentation that remains to be solved.

Graph-based IVUS segmentation with efficient computer-aided refinement

A new graph-based approach for segmentation of luminal and external elastic lamina (EEL) surface of coronary vessels in gated 20 MHz intravascular ultrasound (IVUS) image sequences (volumes) is presented. The approach consists of a fully automated segmentation stage ("new automated" or NA) and a user-guided computer-aided refinement ("new refinement" or NR) stage. Both approaches are based on the LOGISMOS approach for simultaneous dual-surface graph-based segmentation. This combination allows the user to efficiently combine general information about IVUS image appearance and case-specific IVUS morphology and therefore deal with frequently occurring issues like calcified plaque-causing signal shadowing-and imaging artifacts. The automated segmentation stage starts with pre-segmenting the lumen to automatically define the lumen centerline, which is used to transform the segmentation task into a LOGISMOS-family graph optimization problem. Following the automated segmentation, the user can inspect the result and correct local or regional segmentation inaccuracies by (iteratively) providing approximate clues regarding the location of the desired surface locations. This expert information is utilized to modify the previously calculated cost functions, locally re-optimizing the underlying modified graph without a need to start the new optimization from scratch. Validation of our method was performed on 41 gated 20 MHz IVUS data sets for which an expert-defined independent standard was available. Resulting from the automated stage of the approach (NA), the mean and standard deviation of the root mean square area errors for the luminal and external elastic lamina surfaces were 1.12 ±0.67 mm (2) and 2.35 ±1.61 mm (2) , respectively. Following the refinement stage (NR), the root mean square area errors significantly decreased to 0.82 ±0.44 mm (2) and 1.17 ±0.65 mm (2) for the same surfaces, respectively ( for both surfaces). The approach is delivering a previously unachievable speed of obtaining clinically relevant segmentations compared to the current approaches of automated segmentation followed by manual editing.

A four-criterion selection procedure for atherosclerotic plaque elasticity reconstruction based on in vivo coronary intravascular ultrasound radial strain sequences

Plaque elasticity (i.e., modulogram) and morphology are good predictors of plaque vulnerability. Recently, our group developed an intravascular ultrasound (IVUS) elasticity reconstruction method which was successfully implemented in vitro using vessel phantoms. In vivo IVUS modulography, however, remains a major challenge as the motion of the heart prevents accurate strain field estimation. We therefore designed a technique to extract accurate strain fields and modulograms from recorded IVUS sequences. We identified a set of four criteria based on tissue overlapping, RF-correlation coefficient between two successive frames, performance of the elasticity reconstruction method to recover the measured radial strain, and reproducibility of the computed modulograms over the cardiac cycle. This four-criterion selection procedure (4-CSP) was successfully tested on IVUS sequences obtained in twelve patients referred for a directional coronary atherectomy intervention. This study demonstrates the potential of the IVUS modulography technique based on the proposed 4-CSP to detect vulnerable plaques in vivo.

Improving quantification of intravascular fluorescence imaging using structural information

Intravascular near-infrared fluorescence (iNIRF) imaging can enable the in vivo visualization of biomarkers of vascular pathology, including high-risk plaques. The technique resolves the bio-distribution of systemically administered fluorescent probes with molecular specificity in the vessel wall. However, the geometrical variations that may occur in the distance between fibre-tip and vessel wall can lead to signal intensity variations and challenge quantification. Herein we examined whether the use of anatomical information of the cross-section vessel morphology, obtained from co-registered intravascular ultrasound (IVUS), can lead to quantification improvements when fibre-tip and vessel wall distance variations are present. The algorithm developed employs a photon propagation model derived from phantom experiments that is used to calculate the relative attenuation of fluorescence signals as they are collected over 360° along the vessel wall, and utilizes it to restore accurate fluorescence readings. The findings herein point to quantification improvements when employing hybrid iNIRF, with possible implications to the clinical detection of high-risk plaques or blood vessel theranostics.

Segmentation of skin lesions in 2-D and 3-D ultrasound images using a spatially coherent generalized Rayleigh mixture model

This paper addresses the problem of jointly estimating the statistical distribution and segmenting lesions in multiple-tissue high-frequency skin ultrasound images. The distribution of multiple-tissue images is modeled as a spatially coherent finite mixture of heavy-tailed Rayleigh distributions. Spatial coherence inherent to biological tissues is modeled by enforcing local dependence between the mixture components. An original Bayesian algorithm combined with a Markov chain Monte Carlo method is then proposed to jointly estimate the mixture parameters and a label-vector associating each voxel to a tissue. More precisely, a hybrid Metropolis-within-Gibbs sampler is used to draw samples that are asymptotically distributed according to the posterior distribution of the Bayesian model. The Bayesian estimators of the model parameters are then computed from the generated samples. Simulation results are conducted on synthetic data to illustrate the performance of the proposed estimation strategy. The method is then successfully applied to the segmentation of in vivo skin tumors in high-frequency 2-D and 3-D ultrasound images.

Automated tracing of the adventitial contour of aortoiliac and peripheral arterial walls in CT angiography (CTA) to allow calculation of non-calcified plaque burden

Aortoiliac and lower extremity arterial atherosclerotic plaque burden is a risk factor for the development of visceral and peripheral ischemic and aneurismal vascular disease. While prior research allows automated quantification of calcified plaque in these body regions using CT angiograms, no automated method exists to quantify soft plaque. We developed an automatic algorithm that defines the outer wall contour and wall thickness of vessels to quantify non-calcified plaque in CT angiograms of the chest, abdomen, pelvis, and lower extremities. The algorithm encodes the search space as a constrained graph and calculates the outer wall contour by deriving a minimum cost path through the graph, following the visible outer wall contour while minimizing path tortuosity. Our algorithm was statistically equivalent to a reference standard made by two reviewers. Absolute error was 1.9 ± 2.3% compared to the inter-observer variability of 3.9 ± 3.6%. Wall thickness in vessels with atherosclerosis was 3.4 ± 1.6 mm compared to 1.2 ± 0.4 mm in normal vessels. The algorithm shows promise as a tool for quantification of non-calcified plaque in CT angiography. When combined with previous research, our method has the potential to quantify both non-calcified and calcified plaque in all clinically significant systemic arteries, from the thoracic aorta to the arteries of the calf, over a wide range of diameters. This algorithm has the potential to enable risk stratification of patients and facilitate investigations into the relationships between asymptomatic atherosclerosis and a variety of behavioral, physiologic, pathologic, and genotypic conditions.

A state-of-the-art review on segmentation algorithms in intravascular ultrasound (IVUS) images

Over the past two decades, intravascular ultrasound (IVUS) image segmentation has remained a challenge for researchers while the use of this imaging modality is rapidly growing in catheterization procedures and in research studies. IVUS provides cross-sectional grayscale images of the arterial wall and the extent of atherosclerotic plaques with high spatial resolution in real time. In this paper, we review recently developed image processing methods for the detection of media-adventitia and luminal borders in IVUS images acquired with different transducers operating at frequencies ranging from 20 to 45 MHz. We discuss methodological challenges, lack of diversity in reported datasets, and weaknesses of quantification metrics that make IVUS segmentation still an open problem despite all efforts. In conclusion, we call for a common reference database, validation metrics, and ground-truth definition with which new and existing algorithms could be benchmarked.

A Parallel Method for Segmenting Intravascular Ultrasound Image Sequence

Intravascular Ultrasound (IVUS) is one of interventional imaging modalities widely used in clinical diagnosis of vascular diseases, especially coronary artery diseases. Segmentation of IVUS images to extract vessel wall boundaries is of importance for quantitative analysis and 3D vessel reconstruction. A 3D parallel method for segmenting IVUS image sequence is proposed in this paper. Firstly, original images are preprocessed to reduce possible noises and eliminate ring-down artifacts. Then, several longitudinal cuts are obtained and intima-lumen and media-adventitia boundaries are detected. Once these boundaries are mapped onto each cross-sectional slice, initial plan of vessel wall boundaries in each frame is obtained. Finally, these initial contours evolve continuously until stop at target contours. Consequently, segmentation of each IVUS tomographic frame is implemented simultaneously and the efficiency is greatly raised compared with 2D sequential approaches.