Semiautomatic carotid lumen segmentation for quantification of lumen geometry in multispectral MRI

A Dual-Stream Centerline-Guided Network for Segmentation of the Common and Internal Carotid Arteries From 3D Ultrasound Images

Segmentation of the carotid section encompassing the common carotid artery (CCA), the bifurcation and the internal carotid artery (ICA) from three-dimensional ultrasound (3DUS) is required to measure the vessel wall volume (VWV) and localized vessel-wall-plus-plaque thickness (VWT), shown to be sensitive to treatment effect. We proposed an approach to combine a centerline extraction network (CHG-Net) and a dual-stream centerline-guided network (DSCG-Net) to segment the lumen-intima (LIB) and media-adventitia boundaries (MAB) from 3DUS images. Correct arterial location is essential for successful segmentation of the carotid section encompassing the bifurcation. We addressed this challenge by using the arterial centerline to enhance the localization accuracy of the segmentation network. The CHG-Net was developed to generate a heatmap indicating high probability regions for the centerline location, which was then integrated with the 3DUS image by the DSCG-Net to generate the MAB and LIB. The DSCG-Net includes a scale-based and a spatial attention mechanism to fuse multi-level features extracted by the encoder, and a centerline heatmap reconstruction side-branch connected to the end of the encoder to increase the generalization ability of the network. Experiments involving 224 3DUS volumes produce a Dice similarity coefficient (DSC) of 95.8±1.9% and 92.3±5.4% for CCA MAB and LIB, respectively, and 93.2±4.4% and 89.0±10.0% for ICA MAB and LIB, respectively. Our approach outperformed four state-of-the-art 3D CNN models, even after their performances were boosted by centerline guidance. The efficiency afforded by the framework would allow it to be incorporated into the clinical workflow for improved quantification of plaque change.

VWI-APP: Vessel wall imaging-dedicated automated processing pipeline for intracranial atherosclerotic plaque quantification

BACKGROUND

Quantitative plaque assessment based on 3D magnetic resonance (MR) vessel wall imaging (VWI) has been shown to provide valuable numerical markers of the burden and risk of intracranial atherosclerotic disease (ICAD). However, plaque quantification is currently time-consuming and observer-dependent due to the demand for heavy manual effort. A VWI-dedicated automated processing pipeline (VWI-APP) is desirable.

PURPOSE

To develop and evaluate a VWI-APP for end-to-end quantitative analysis of intracranial atherosclerotic plaque.

METHODS

We retrospectively enrolled 91 subjects with ICAD (80 for pipeline development, 10 for an end-to-end pipeline evaluation, and 1 for demonstrating longitudinal plaque assessment) who had undergone VWI and MR angiography. In an end-to-end evaluation, diameter stenosis (DS), normalized wall index (NWI), remodeling ratio (RR), plaque wall contrast ratio (CR), and total plaque volume (TPV) were quantified at each culprit lesion using the developed VWI-APP and a computer-aided manual approach by a neuroradiologist, respectively. The time consumed in each quantification approach was recorded. Two-sided paired t-tests and intraclass correlation coefficient (ICC) were used to determine the difference and agreement in each plaque metric between VWI-APP and manual quantification approaches.

RESULTS

There was no significant difference between VWI-APP and manual quantification in each plaque metric. The ICC was 0.890, 0.813, 0.827, 0.891, and 0.991 for DS, NWI, RR, CR, and TPV, respectively, suggesting good to excellent accuracy of the pipeline method in plaque quantification. Quantitative analysis of each culprit lesion on average took 675.7 s using the manual approach but shortened to 238.3 s with the aid of VWI-APP.

CONCLUSIONS

VWI-APP is an accurate and efficient approach to intracranial atherosclerotic plaque quantification. Further clinical assessment of this automated tool is warranted to establish its utility in the risk assessment of ICAD lesions.

3D carotid artery segmentation using shape-constrained active contours

Reconstruction of the carotid artery is demanded in the detection and characterization of atherosclerosis. This study proposes a shape-constrained active contour model for segmenting the carotid artery from MR images, which embeds the output of the deep learning network into the active contour. First the centerline of the carotid artery is localized and then modified active contour initialized from the centerline is used to extract the vessel lumen, finally the probability atlas generated by the deep learning network in polar representation domain is integrated into the active contour as a prior information to detect the outer wall. The results showed that the proposed active contour model was efficient and comparable to manual segmentation.

Assessment of fully automatic segmentation of pulmonary artery and aorta on noncontrast CT with optimal surface graph cuts

Purpose

Accurate segmentation of the pulmonary arteries and aorta is important due to the association of the diameter and the shape of these vessels with several cardiovascular diseases and with the risk of exacerbations and death in patients with chronic obstructive pulmonary disease. We propose a fully automatic method based on an optimal surface graph‐cut algorithm to quantify the full 3D shape and the diameters of the pulmonary arteries and aorta in noncontrast computed tomography (CT) scans.

Methods

The proposed algorithm first extracts seed points in the right and left pulmonary arteries, the pulmonary trunk, and the ascending and descending aorta by using multi‐atlas registration. Subsequently, the centerlines of the pulmonary arteries and aorta are extracted by a minimum cost path tracking between the extracted seed points, with a cost based on a combination of lumen intensity similarity and multiscale medialness in three planes. The centerlines are refined by applying the path tracking algorithm to curved multiplanar reformatted scans and are then smoothed and dilated nonuniformly according to the extracted local vessel radius from the medialness filter. The resulting coarse estimates of the vessels are used as initialization for a graph‐cut segmentation. Once the vessels are segmented, the diameters of the pulmonary artery (PA) and the ascending aorta (AA) and the PA:AA ratio are automatically calculated both in a single axial slice and in a 10 mm volume around the automatically extracted PA bifurcation level. The method is evaluated on noncontrast CT scans from the Danish Lung Cancer Screening Trial (DLCST). Segmentation accuracy is determined by comparing with manual annotations on 25 CT scans. Intraclass correlation (ICC) between manual and automatic diameters, both measured in axial slices at the PA bifurcation level, is computed on an additional 200 CT scans. Repeatability of the automated 3D volumetric diameter and PA:AA ratio calculations (perpendicular to the vessel axis) are evaluated on 118 scan–rescan pairs with an average in‐between time of 3 months.

Results

We obtained a Dice segmentation overlap of 0.94 ± 0.02 for pulmonary arteries and 0.96 ± 0.01 for the aorta, with a mean surface distance of 0.62 ± 0.33 mm and 0.43 ± 0.07 mm, respectively. ICC between manual and automatic in‐slice diameter measures was 0.92 for PA, 0.97 for AA, and 0.90 for the PA:AA ratio, and for automatic diameters in 3D volumes around the PA bifurcation level between scan and rescan was 0.89, 0.95, and 0.86, respectively.

Conclusion

The proposed automatic segmentation method can reliably extract diameters of the large arteries in non‐ECG‐gated noncontrast CT scans such as are acquired in lung cancer screening.

Automated segmentation of the individual branches of the carotid arteries in contrast-enhanced MR angiography using DeepMedic

Background

Non-invasive imaging is of interest for tracking the progression of atherosclerosis in the carotid bifurcation, and segmenting this region into its constituent branch arteries is necessary for analyses. The purpose of this study was to validate and demonstrate a method for segmenting the carotid bifurcation into the common, internal, and external carotid arteries (CCA, ICA, ECA) in contrast-enhanced MR angiography (CE-MRA) data.

Methods

A segmentation pipeline utilizing a convolutional neural network (DeepMedic) was tailored and trained for multi-class segmentation of the carotid arteries in CE-MRA data from the Swedish CardioPulmonsary bioImage Study (SCAPIS). Segmentation quality was quantitatively assessed using the Dice similarity coefficient (DSC), Matthews Correlation Coefficient (MCC), F₂, F_0.5, and True Positive Ratio (TPR). Segmentations were also assessed qualitatively, by three observers using visual inspection. Finally, geometric descriptions of the carotid bifurcations were generated for each subject to demonstrate the utility of the proposed segmentation method.

Results

Branch-level segmentations scored DSC = 0.80 ± 0.13, MCC = 0.80 ± 0.12, F₂ = 0.82 ± 0.14, F_0.5 = 0.78 ± 0.13, and TPR = 0.84 ± 0.16, on average in a testing cohort of 46 carotid bifurcations. Qualitatively, 61% of segmentations were judged to be usable for analyses without adjustments in a cohort of 336 carotid bifurcations without ground-truth. Carotid artery geometry showed wide variation within the whole cohort, with CCA diameter 8.6 ± 1.1 mm, ICA 7.5 ± 1.4 mm, ECA 5.7 ± 1.0 mm and bifurcation angle 41 ± 21°.

Conclusion

The proposed segmentation method automatically generates branch-level segmentations of the carotid arteries that are suitable for use in further analyses and help enable large-cohort investigations.

Cascaded residual U-net for fully automatic segmentation of 3D carotid artery in high-resolution multi-contrast MR images

Accurate and automatic carotid artery segmentation for magnetic resonance (MR) images is eagerly expected, which can greatly assist a comprehensive study of atherosclerosis and accelerate the translation. Although many efforts have been made, identification of the inner lumen and outer wall in diseased vessels is still a challenging task due to complex vascular deformation, blurred wall boundary, and confusing componential expression. In this paper, we introduce a novel fully automatic 3D framework for simultaneously segmenting the carotid artery from high-resolution multi-contrast MR sequences based on deep learning. First, an optimal channel fitting structure is designed for identity mapping, and a novel 3D residual U-net is used as a basic network. Second, high-resolution MR images are trained using both patch-level and global-level strategies, and the two pre-segmentation results are optimized based on structural characteristics. Third, the optimized pre-segmentation results are cascaded with the patch-cropped MR volume data and trained to segment the carotid lumen and wall. Extensive experiments demonstrate the proposed method outperforms the state-of-the-art 3D Unet-based segmentation models.

Towards quantitative evaluation of wall shear stress from 4D flow imaging

Wall shear stress (WSS) is a relevant hemodynamic indicator of the local stress applied on the endothelium surface. More specifically, its spatiotemporal distribution reveals crucial in the evolution of many pathologies such as aneurysm, stenosis, and atherosclerosis. This paper introduces a new solution, called PaLMA, to quantify the WSS from 4D Flow MRI data. It relies on a two-step local parametric model, to accurately describe the vessel wall and the velocity-vector field in the neighborhood of a given point of interest. Extensive validations have been performed on synthetic 4D Flow MRI data, including four datasets generated from patient specific computational fluid dynamics simulations on carotids. The validation tests are focused on the impact of the noise component, of the resolution level, and of the segmentation accuracy concerning the vessel position in the context of complex flow patterns. In simulated cases aimed to reproduce clinical acquisition conditions, the WSS quantification performance reached by PaLMA is significantly higher (with a gain in RMSE of 12 to 27%) than the reference one obtained using the smoothing B-spline method proposed by Potters et al. (2015) method, while the computation time is equivalent for both WSS quantification methods.

A deep learning oriented method for automated 3D reconstruction of carotid arterial trees from MR imaging

The scope of this paper is to present a new carotid vessel segmentation algorithm implementing the U-net based convolutional neural network architecture. With carotid atherosclerosis being the major cause of stroke in Europe, new methods that can provide more accurate image segmentation of the carotid arterial tree and plaque tissue can help improve early diagnosis, prevention and treatment of carotid disease. Herein, we present a novel methodology combining the U-net model and morphological active contours in an iterative framework that accurately segments the carotid lumen and outer wall. The method automatically produces a 3D meshed model of the carotid bifurcation and smaller branches, using multispectral MR image series obtained from two clinical centres of the TAXINOMISIS study. As indicated by a validation study, the algorithm succeeds high accuracy (99.1% for lumen area and 92.6% for the perimeter) for lumen segmentation. The proposed algorithm will be used in the TAXINOMISIS study to obtain more accurate 3D vessel models for improved computational fluid dynamics simulations and the development of models of atherosclerotic plaque progression.

Deep morphology aided diagnosis network for segmentation of carotid artery vessel wall and diagnosis of carotid atherosclerosis on black-blood vessel wall MRI

PURPOSE

Early detection of carotid atherosclerosis on the vessel wall (VW) magnetic resonance imaging (MRI) (VW-MRI) images can prevent the progression of cardiovascular disease. However, the manual inspection process of the VW-MRI images is cumbersome and has low reproducibility. Therefore in this paper, by using the convolutional neural networks (CNNs), we develop a deep morphology aided diagnosis (DeepMAD) network for automated segmentation of the VW of carotid artery and for automated diagnosis of the carotid atherosclerosis with the black-blood (BB) VW-MRI (i.e., the T1-weighted MRI) in a slice-by-slice manner.

METHODS

The proposed DeepMAD network consists of a segmentation subnetwork and a diagnosis subnetwork for performing the segmentation and diagnosis tasks on the BB-VW-MRI images, where the manual labeled lumen area, the manual labeled outer wall area and the manual labeled lesion Types based on the modified American Heart Association (AHA) criteria are used as the ground-truth. Specifically, a deep U-shape CNN with a weighted fusion layer is designed as the segmentation subnetwork, where the lumen area and the outer wall area can be simultaneously segmented under the supervision of the triple Dice loss to provide the vessel wall map as morphological information. Then, the image stream from the BB-VWMRI image and the morphology stream from the obtained vessel wall map are extracted from two deep CNNs and combined to obtain the diagnosis results of atherosclerosis in the diagnosis subnetwork. In addition, the triple input set is formed by three carotid regions of interest (ROIs) from three consecutive slices of the MRI sequence and input to the DeepMAD network, where the first and last slices used as additional adjacent slices to provide 2.5D spatial information along the carotid artery centerline for the intermediate slice, which is the target slice for segmentation and diagnosis in the study.

RESULTS

Compared to other existing methods, the DeepMAD network can achieve promising segmentation performances (0.9594 Dice for the lumen and 0.9657 Dice for the outer wall) and better diagnosis Accuracy of the carotid atherosclerosis (0.9503 AUC and 0.8916 Accuracy) in the test dataset (including invisible subjects) from same source as the training dataset. In addition, the trained DeepMAD model can be successfully transferred to another test dataset for segmentation and diagnosis tasks with remarkable performance (0.9475 Dice for the lumen and 0.9542 Dice for the outer wall, 0. 9227 AUC and 0.8679 Accuracy for diagnosis).

CONCLUSIONS

Even without the intervention of reviewers required for previous works, the proposed DeepMAD network automatically segments the lumen and the outer wall together and diagnoses the carotid atherosclerosis with high performances. The DeepMAD network can be used in clinical trials to help radiologists get rid of tedious reading tasks, such as screening review to separate the normal carotid from the atherosclerotic arteries and outlining the vessel wall contours.

Automated 3D segmentation and diameter measurement of the thoracic aorta on non-contrast enhanced CT

Objectives

To develop and evaluate a fully automatic method to measure diameters of the ascending and descending aorta on non-ECG-gated, non-contrast computed tomography (CT) scans.

Material and methods

The method combines multi-atlas registration to obtain seed points, aorta centerline extraction, and an optimal surface segmentation approach to extract the aorta surface around the centerline. From the extracted 3D aorta segmentation, the diameter of the ascending and descending aorta was calculated at cross-sectional slices perpendicular to the extracted centerline, at the level of the pulmonary artery bifurcation, and at 1-cm intervals up to 3 cm above and below this level. Agreement with manual annotations was evaluated by dice similarity coefficient (DSC) for segmentation overlap, mean surface distance (MSD), and intra-class correlation (ICC) of diameters on 100 CT scans from a lung cancer screening trial. Repeatability of the diameter measurements was evaluated on 617 baseline-one year follow-up CT scan pairs.

Results

The agreement between manual and automatic segmentations was good with 0.95 ± 0.01 DSC and 0.56 ± 0.08 mm MSD. ICC between the diameters derived from manual and from automatic segmentations was 0.97, with the per-level ICC ranging from 0.87 to 0.94. An ICC of 0.98 for all measurements and per-level ICC ranging from 0.91 to 0.96 were obtained for repeatability.

Conclusion

This fully automatic method can assess diameters in the thoracic aorta reliably even in non-ECG-gated, non-contrast CT scans. This could be a promising tool to assess aorta dilatation in screening and in clinical practice.

Key Points

• Fully automatic method to assess thoracic aorta diameters.

• High agreement between fully automatic method and manual segmentations.

• Method is suitable for non-ECG-gated CT and can therefore be used in screening.

Electronic supplementary material

The online version of this article (10.1007/s00330-018-5931-z) contains supplementary material, which is available to authorized users.

Cooperative carotid artery centerline extraction in MRI

Centerline extraction of the carotid artery in MRI is important to analyze the artery geometry and to provide input for further processing such as registration and segmentation. The centerline of the artery bifurcation is often extracted by means of two independent minimum cost paths ranging from the common to the internal and the external carotid artery. Often the cost is not well defined at the artery bifurcation, leading to centerline errors. To solve this problem, we developed a method to cooperatively extract both centerlines, where in the cost to extract each centerline, we integrate a constraint region derived from the estimated position of the neighbor centerline. This method avoids that both centerlines follow the same cheapest path after the bifurcation, which is a common error when the paths are extracted independently. We show that this method results in less error compared to extracting them independently: 10 failed centerlines Vs. 3 failures in a data set of 161 arteries with manual annotations. Additionally, we show that the new method improves the non-cooperative approach in 28 cases (p < 0.0001) in a data set of 3,904 arteries.

Blood vessel segmentation algorithms - Review of methods, datasets and evaluation metrics

BACKGROUND

Blood vessel segmentation is a topic of high interest in medical image analysis since the analysis of vessels is crucial for diagnosis, treatment planning and execution, and evaluation of clinical outcomes in different fields, including laryngology, neurosurgery and ophthalmology. Automatic or semi-automatic vessel segmentation can support clinicians in performing these tasks. Different medical imaging techniques are currently used in clinical practice and an appropriate choice of the segmentation algorithm is mandatory to deal with the adopted imaging technique characteristics (e.g. resolution, noise and vessel contrast).

OBJECTIVE

This paper aims at reviewing the most recent and innovative blood vessel segmentation algorithms. Among the algorithms and approaches considered, we deeply investigated the most novel blood vessel segmentation including machine learning, deformable model, and tracking-based approaches.

METHODS

This paper analyzes more than 100 articles focused on blood vessel segmentation methods. For each analyzed approach, summary tables are presented reporting imaging technique used, anatomical region and performance measures employed. Benefits and disadvantages of each method are highlighted.

DISCUSSION

Despite the constant progress and efforts addressed in the field, several issues still need to be overcome. A relevant limitation consists in the segmentation of pathological vessels. Unfortunately, not consistent research effort has been addressed to this issue yet. Research is needed since some of the main assumptions made for healthy vessels (such as linearity and circular cross-section) do not hold in pathological tissues, which on the other hand require new vessel model formulations. Moreover, image intensity drops, noise and low contrast still represent an important obstacle for the achievement of a high-quality enhancement. This is particularly true for optical imaging, where the image quality is usually lower in terms of noise and contrast with respect to magnetic resonance and computer tomography angiography.

CONCLUSION

No single segmentation approach is suitable for all the different anatomical region or imaging modalities, thus the primary goal of this review was to provide an up to date source of information about the state of the art of the vessel segmentation algorithms so that the most suitable methods can be chosen according to the specific task.

Automated Registration of Freehand B-Mode Ultrasound and Magnetic Resonance Imaging of the Carotid Arteries Based on Geometric Features

An automated method for registering B-mode ultrasound (US) and magnetic resonance imaging (MRI) of the carotid arteries is proposed. The registration uses geometric features, namely, lumen centerlines and lumen segmentations, which are extracted fully automatically from the images after manual annotation of three seed points in US and MRI. The registration procedure starts with alignment of the lumen centerlines using a point-based registration algorithm. The resulting rigid transformation is used to initialize a rigid and subsequent non-rigid registration procedure that jointly aligns centerlines and segmentations by minimizing a weighted sum of the Euclidean distance between centerlines and the dissimilarity between segmentations. The method was evaluated in 28 carotid arteries from eight patients and six healthy volunteers. First, the automated US lumen segmentation method was validated and optimized in a cross-validation experiment. Next, the effect of the weighting parameter of the proposed registration dissimilarity metric and the control point spacing in the non-rigid registration was evaluated. Finally, the proposed registration method was evaluated in comparison to an existing intensity-and-point-based method, a registration using only the centerlines and a registration using manual US lumen segmentations. Registration accuracy was measured in terms of the mean surface distance between manual US segmentations and the registered MRI segmentations. The average mean surface distance was 0.78 ± 0.34 mm for all subjects, 0.65 ± 0.09 mm for healthy volunteers and 0.87 ± 0.42 mm for patients. The results for the complete set were significantly better (Wilcoxon test, p < 0.01) than the results for the intensity-and-point-based method and the centerline-based registration method. We conclude that the proposed method can robustly and accurately register US and MR images of the carotid artery, allowing multimodal analysis of the carotid plaque to improve plaque assessment.

Carotid Artery Wall Segmentation in Multispectral MRI by Coupled Optimal Surface Graph Cuts

We present a new three-dimensional coupled optimal surface graph-cut algorithm to segment the wall of the carotid artery bifurcation from Magnetic Resonance (MR) images. The method combines the search for both inner and outer borders into a single graph cut and uses cost functions that integrate information from multiple sequences. Our approach requires manual localization of only three seed points indicating the start and end points of the segmentation in the internal, external, and common carotid artery. We performed a quantitative validation using images of 57 carotid arteries. Dice overlap of 0.86 ± 0.06 for the complete vessel and 0.89 ± 0.05 for the lumen compared to manual annotation were obtained. Reproducibility tests were performed in 60 scans acquired with an interval of 15 ± 9 days, showing good agreement between baseline and follow-up segmentations with intraclass correlations of 0.96 and 0.74 for the lumen and complete vessel volumes respectively.

A Robust and Accurate Two-Step Auto-Labeling Conditional Iterative Closest Points (TACICP) Algorithm for Three-Dimensional Multi-Modal Carotid Image Registration

Atherosclerosis is among the leading causes of death and disability. Combining information from multi-modal vascular images is an effective and efficient way to diagnose and monitor atherosclerosis, in which image registration is a key technique. In this paper a feature-based registration algorithm, Two-step Auto-labeling Conditional Iterative Closed Points (TACICP) algorithm, is proposed to align three-dimensional carotid image datasets from ultrasound (US) and magnetic resonance (MR). Based on 2D segmented contours, a coarse-to-fine strategy is employed with two steps: rigid initialization step and non-rigid refinement step. Conditional Iterative Closest Points (CICP) algorithm is given in rigid initialization step to obtain the robust rigid transformation and label configurations. Then the labels and CICP algorithm with non-rigid thin-plate-spline (TPS) transformation model is introduced to solve non-rigid carotid deformation between different body positions. The results demonstrate that proposed TACICP algorithm has achieved an average registration error of less than 0.2mm with no failure case, which is superior to the state-of-the-art feature-based methods.

IT Infrastructure to support the secondary use of routinely acquired clinical imaging data for research

We propose an infrastructure for the automated anonymization, extraction and processing of image data stored in clinical data repositories to make routinely acquired imaging data available for research purposes. The automated system, which was tested in the context of analyzing routinely acquired MR brain imaging data, consists of four modules: subject selection using PACS query, anonymization of privacy sensitive information and removal of facial features, quality assurance on DICOM header and image information, and quantitative imaging biomarker extraction. In total, 1,616 examinations were selected based on the following MRI scanning protocols: dementia protocol (246), multiple sclerosis protocol (446) and open question protocol (924). We evaluated the effectiveness of the infrastructure in accessing and successfully extracting biomarkers from routinely acquired clinical imaging data. To examine the validity, we compared brain volumes between patient groups with positive and negative diagnosis, according to the patient reports. Overall, success rates of image data retrieval and automatic processing were 82.5 %, 82.3 % and 66.2 % for the three protocol groups respectively, indicating that a large percentage of routinely acquired clinical imaging data can be used for brain volumetry research, despite image heterogeneity. In line with the literature, brain volumes were found to be significantly smaller (p-value <0.001) in patients with a positive diagnosis of dementia (915 ml) compared to patients with a negative diagnosis (939 ml). This study demonstrates that quantitative image biomarkers such as intracranial and brain volume can be extracted from routinely acquired clinical imaging data. This enables secondary use of clinical images for research into quantitative biomarkers at a hitherto unprecedented scale.

Joint segmentation of lumen and outer wall from femoral artery MR images: Towards 3D imaging measurements of peripheral arterial disease

Three-dimensional (3D) measurements of peripheral arterial disease (PAD) plaque burden extracted from fast black-blood magnetic resonance (MR) images have shown to be more predictive of clinical outcomes than PAD stenosis measurements. To this end, accurate segmentation of the femoral artery lumen and outer wall is required for generating volumetric measurements of PAD plaque burden. Here, we propose a semi-automated algorithm to jointly segment the femoral artery lumen and outer wall surfaces from 3D black-blood MR images, which are reoriented and reconstructed along the medial axis of the femoral artery to obtain improved spatial coherence between slices of the long, thin femoral artery and to reduce computation time. The developed segmentation algorithm enforces two priors in a global optimization manner: the spatial consistency between the adjacent 2D slices and the anatomical region order between the femoral artery lumen and outer wall surfaces. The formulated combinatorial optimization problem for segmentation is solved globally and exactly by means of convex relaxation using a coupled continuous max-flow (CCMF) model, which is a dual formulation to the convex relaxed optimization problem. In addition, the CCMF model directly derives an efficient duality-based algorithm based on the modern multiplier augmented optimization scheme, which has been implemented on a GPU for fast computation. The computed segmentations from the developed algorithm were compared to manual delineations from experts using 20 black-blood MR images. The developed algorithm yielded both high accuracy (Dice similarity coefficients ≥ 87% for both the lumen and outer wall surfaces) and high reproducibility (intra-class correlation coefficient of 0.95 for generating vessel wall area), while outperforming the state-of-the-art method in terms of computational time by a factor of ≈ 20.

Topology adaptive vessel network skeleton extraction with novel medialness measuring function

Vessel tree skeleton extraction is widely applied in vascular structure segmentation, however, conventional approaches often suffer from the adjacent interferences and poor topological adaptability. To avoid these problems, a robust, topology adaptive tree-like structure skeleton extraction framework is proposed in this paper. Specifically, to avoid the adjacent interferences, a local message passing procedure called Gaussian affinity voting (GAV) is proposed to realize adaptive scale-growing of vessel voxels. Then the medialness measuring function (MMF) based on GAV, namely GAV-MMF, is constructed to extract medialness patterns robustly. In order to improve topological adaptability, a level-set graph embedded with GAV-MMF is employed to build initial curve skeletons without any user interaction. Furthermore, the GAV-MMF is embedded in stretching open active contours (SOAC) to drive the initial curves to the expected location, maintaining smoothness and continuity. In addition, to provide an accurate and smooth final skeleton tree topology, topological checks and skeleton network reconfiguration is proposed. The continuity and scalability of this method is validated experimentally on synthetic and clinical images for multi-scale vessels. Experimental results show that the proposed method achieves acceptable topological adaptability for skeleton extraction of vessel trees.

Joint intensity-and-point based registration of free-hand B-mode ultrasound and MRI of the carotid artery

PURPOSE

To introduce a semiautomatic algorithm to perform the registration of free-hand B-Mode ultrasound (US) and magnetic resonance imaging (MRI) of the carotid artery.

METHODS

The authors' approach combines geometrical features and intensity information. The only user interaction consists of placing three seed points in US and MRI. First, the lumen centerlines are used as landmarks for point based registration. Subsequently, in a joint optimization the distance between centerlines and the dissimilarity of the image intensities is minimized. Evaluation is performed in left and right carotids from six healthy volunteers and five patients with atherosclerosis. For the validation, the authors measure the Dice similarity coefficient (DSC) and the mean surface distance (MSD) between carotid lumen segmentations in US and MRI after registration. The effect of several design parameters on the registration accuracy is investigated by an exhaustive search on a training set of five volunteers and three patients. The optimum configuration is validated on the remaining images of one volunteer and two patients.

RESULTS

On the training set, the authors achieve an average DSC of 0.74 and a MSD of 0.66 mm on volunteer data. For the patient data, the authors obtain a DSC of 0.77 and a MSD of 0.69 mm. In the independent set composed of patient and volunteer data, the DSC is 0.69 and the MSD is 0.87 mm. The experiments with different design parameters show that nonrigid registration outperforms rigid registration, and that the combination of intensity and point information is superior to approaches that use intensity or points only.

CONCLUSIONS

The proposed method achieves an accurate registration of US and MRI, and may thus enable multimodal analysis of the carotid plaque.

Lumen segmentation and stenosis quantification of atherosclerotic carotid arteries in CTA utilizing a centerline intensity prior

PURPOSE

The degree of stenosis is an important biomarker in assessing the severity of cardiovascular disease. The purpose of our work is to develop and evaluate a semiautomatic method for carotid lumen segmentation and subsequent carotid artery stenosis quantification in CTA images.

METHODS

The authors present a semiautomatic stenosis detection and quantification method following lumen segmentation. The lumen of the carotid arteries is segmented in three steps. First, centerlines of the internal and external carotid arteries are extracted with an iterative minimum cost path approach in which the costs are based on a measure of medialness and intensity similarity to lumen. Second, the lumen boundary is delineated using a level set procedure which is steered by gradient information, regional intensity information, and spatial information. Special effort is made in adding terms based on local centerline intensity prior so as to exclude all possible plaque tissues from the segmentation. Third, side branches in the segmented lumen are removed by applying a shape constraint to the envelope of the maximum inscribed spheres of the segmentation. From the segmented lumen, the authors detect and quantify the cross-sectional area-based and cross-sectional diameter-based stenosis degrees according to the North American Symptomatic Carotid En-darterectomy Trial criterion.

RESULTS

The method is trained and tested on a publicly available database from the cls2009 challenge. For the segmentation, the authors obtain a Dice similarity coefficient of 90.2% and a mean absolute surface distance of 0.34 mm. For the stenosis quantification, the authors obtain an average error of 15.7% for cross-sectional diameter-based stenosis and 19.2% for cross-sectional area-based stenosis quantification.

CONCLUSIONS

With these results, the method ranks second in terms of carotid lumen segmentation accuracy, and first in terms of carotid artery stenosis quantification.

Model-guided directional minimal path for fully automatic extraction of coronary centerlines from cardiac CTA

Extracting centerlines of coronary arteries is a challenging but important task in clinical applications of cardiac CTA. In this paper, we propose a model-guided approach, the directional minimal path, for the centerline extraction. The proposed method is based on the minimal path algorithm and a prior coronary model is used. The model is first registered to the unseen image. Then, the start point and end point for the minimal path algorithm are provided by the model to automate the centerline extraction process. Also, the direction information of the coronary model is used to guide the path tracking of the minimal path procedure. This directional tracking improves the robustness and accuracy of the centerline extraction. Finally, the proposed method can automatically recognize the branches of the extracted coronary artery using the prior information in the model. We validated the proposed method by extracting the three main coronary branches. The mean accuracy of the 56 cases was 1.32+/-0.81 mm and the detection ratio was 88.7%.