Coronary vessel trees from 3D imagery: a topological approach

Automatic Coronary Artery Segmentation Using Active Search for Branches and Seemingly Disconnected Vessel Segments from Coronary CT Angiography

We propose a Bayesian tracking and segmentation method of coronary arteries on coronary computed tomographic angiography (CCTA). The geometry of coronary arteries including lumen boundary is estimated in Maximum A Posteriori (MAP) framework. Three consecutive sphere based filtering is combined with a stochastic process that is based on the similarity of the consecutive local neighborhood voxels and the geometric constraint of a vessel. It is also founded on the prior knowledge that an artery can be seen locally disconnected and consist of branches which may be seemingly disconnected due to plaque build up. For such problem, an active search method is proposed to find branches and seemingly disconnected but actually connected vessel segments. Several new measures have been developed for branch detection, disconnection check and planar vesselness measure. Using public domain Rotterdam CT dataset, the accuracy of extracted centerline is demonstrated and automatic reconstruction of coronary artery mesh is shown.

Digital Topology and Geometry in Medical Imaging: A Survey

Digital topology and geometry refers to the use of topologic and geometric properties and features for images defined in digital grids. Such methods have been widely used in many medical imaging applications, including image segmentation, visualization, manipulation, interpolation, registration, surface-tracking, object representation, correction, quantitative morphometry etc. Digital topology and geometry play important roles in medical imaging research by enriching the scope of target outcomes and by adding strong theoretical foundations with enhanced stability, fidelity, and efficiency. This paper presents a comprehensive yet compact survey on results, principles, and insights of methods related to digital topology and geometry with strong emphasis on understanding their roles in various medical imaging applications. Specifically, this paper reviews methods related to distance analysis and path propagation, connectivity, surface-tracking, image segmentation, boundary and centerline detection, topology preservation and local topological properties, skeletonization, and object representation, correction, and quantitative morphometry. A common thread among the topics reviewed in this paper is that their theory and algorithms use the principle of digital path connectivity, path propagation, and neighborhood analysis.

Automated generation of directed graphs from vascular segmentations

Automated feature extraction from medical images is an important task in imaging informatics. We describe a graph-based technique for automatically identifying vascular substructures within a vascular tree segmentation. We illustrate our technique using vascular segmentations from computed tomography pulmonary angiography images. The segmentations were acquired in a semi-automated fashion using existing segmentation tools. A 3D parallel thinning algorithm was used to generate the vascular skeleton and then graph-based techniques were used to transform the skeleton to a directed graph with bifurcations and endpoints as nodes in the graph. Machine-learning classifiers were used to automatically prune false vascular structures from the directed graph. Semantic labeling of portions of the graph with pulmonary anatomy (pulmonary trunk and left and right pulmonary arteries) was achieved with high accuracy (percent correct⩾0.97). Least-squares cubic splines of the centerline paths between nodes were computed and were used to extract morphological features of the vascular tree. The graphs were used to automatically obtain diameter measurements that had high correlation (r⩾0.77) with manual measurements made from the same arteries.

A fast seed detection using local geometrical feature for automatic tracking of coronary arteries in CTA

We propose a fast seed detection for automatic tracking of coronary arteries in coronary computed tomographic angiography (CCTA). To detect vessel regions, Hessian-based filtering is combined with a new local geometric feature that is based on the similarity of the consecutive cross-sections perpendicular to the vessel direction. It is in turn founded on the prior knowledge that a vessel segment is shaped like a cylinder in axial slices. To improve computational efficiency, an axial slice, which contains part of three main coronary arteries, is selected and regions of interest (ROIs) are extracted in the slice. Only for the voxels belonging to the ROIs, the proposed geometric feature is calculated. With the seed points, which are the centroids of the detected vessel regions, and their vessel directions, vessel tracking method can be used for artery extraction. Here a particle filtering-based tracking algorithm is tested. Using 19 clinical CCTA datasets, it is demonstrated that the proposed method detects seed points and can be used for full automatic coronary artery extraction. ROC (receiver operating characteristic) curve analysis shows the advantages of the proposed method.

Unraveling flow patterns through nonlinear manifold learning

From climatology to biofluidics, the characterization of complex flows relies on computationally expensive kinematic and kinetic measurements. In addition, such big data are difficult to handle in real time, thereby hampering advancements in the area of flow control and distributed sensing. Here, we propose a novel framework for unsupervised characterization of flow patterns through nonlinear manifold learning. Specifically, we apply the isometric feature mapping (Isomap) to experimental video data of the wake past a circular cylinder from steady to turbulent flows. Without direct velocity measurements, we show that manifold topology is intrinsically related to flow regime and that Isomap global coordinates can unravel salient flow features.

Vessel connectivity using Murray's hypothesis

We describe a new method for vascular image analysis that incorporates a generic physiological principle to estimate vessel connectivity, which is a key issue in reconstructing complete vascular trees from image data. We follow Murray's hypothesis of the minimum work principle to formulate the problem as an optimization problem. This principle reflects a global property of any vascular network, in contrast to various local geometric properties adopted as constraints previously. We demonstrate the effectiveness of our method using a set of microCT mouse coronary images. It is shown that the performance of our method has a statistically significant improvement over the widely adopted minimum spanning tree methods that rely on local geometric constraints.

Automatic centerline extraction of coronary arteries in coronary computed tomographic angiography

Coronary computed tomographic angiography (CCTA) is a non-invasive imaging modality for the visualization of the heart and coronary arteries. To fully exploit the potential of the CCTA datasets and apply it in clinical practice, an automated coronary artery extraction approach is needed. The purpose of this paper is to present and validate a fully automatic centerline extraction algorithm for coronary arteries in CCTA images. The algorithm is based on an improved version of Frangi's vesselness filter which removes unwanted step-edge responses at the boundaries of the cardiac chambers. Building upon this new vesselness filter, the coronary artery extraction pipeline extracts the centerlines of main branches as well as side-branches automatically. This algorithm was first evaluated with a standardized evaluation framework named Rotterdam Coronary Artery Algorithm Evaluation Framework used in the MICCAI Coronary Artery Tracking challenge 2008 (CAT08). It includes 128 reference centerlines which were manually delineated. The average overlap and accuracy measures of our method were 93.7% and 0.30 mm, respectively, which ranked at the 1st and 3rd place compared to five other automatic methods presented in the CAT08. Secondly, in 50 clinical datasets, a total of 100 reference centerlines were generated from lumen contours in the transversal planes which were manually corrected by an expert from the cardiology department. In this evaluation, the average overlap and accuracy were 96.1% and 0.33 mm, respectively. The entire processing time for one dataset is less than 2 min on a standard desktop computer. In conclusion, our newly developed automatic approach can extract coronary arteries in CCTA images with excellent performances in extraction ability and accuracy.

Segmentation of liver vasculature from contrast enhanced CT images using context-based voting

A novel vessel context-based voting is proposed for automatic liver vasculature segmentation in CT images. It is able to conduct full vessel segmentation and recognition of multiple vasculatures effectively. The vessel context describes context information of a voxel related to vessel properties, such as intensity, saliency, direction and connectivity. Voxels are grouped to liver vasculatures hierarchically based on vessel context. They are first grouped locally into vessel branches with the advantage of a vessel junction measurement, and then grouped globally into vasculatures, which is implemented using a multiple feature point voting mechanism. The proposed method has been evaluated on 10 clinical CT datasets. Segmentation of third-order vessel trees from CT images (0.76 × 0.76 × 2.0mm) of the portal venous phase takes less than 3 min on a PC with 2.0 GHz dual core processor and the average segmentation accuracy is up to 98%.

Vascular Tree Reconstruction by Minimizing A Physiological Functional Cost

The reconstruction of complete vascular trees from medical images has many important applications. Although vessel detection has been extensively investigated, little work has been done on how connect the results to reconstruct the full trees. In this paper, we propose a novel theoretical framework for automatic vessel connection, where the automation is achieved by leveraging constraints from the physiological properties of the vascular trees. In particular, a physiological functional cost for the whole vascular tree is derived and an efficient algorithm is developed to minimize it. The method is generic and can be applied to different vessel detection/segmentation results, e.g. the classic rigid detection method as adopted in this paper. We demonstrate the effectiveness of this method on both 2D and 3D data.

Coronary centerline extraction from CT coronary angiography images using a minimum cost path approach

PURPOSE

The application and large-scale evaluation of minimum cost path approaches for coronary centerline extraction from computed tomography coronary angiography (CTCA) data and the development and evaluation of a novel method to reduce the user-interaction time.

METHODS

A semiautomatic method based on a minimum cost path approach is evaluated for two different cost functions. The first cost function is based on a frequently used vesselness measure and intensity information, and the second is a recently proposed cost function based on region statistics. User interaction is minimized to one or two mouse clicks distally in the coronary artery. The starting point for the minimum cost path search is automatically determined using a newly developed method that finds a point in the center of the aorta in one of the axial slices. This step ensures that all computationally expensive parts of the algorithm can be precomputed.

RESULTS

The performance of the aorta localization procedure was demonstrated by a success rate of 100% in 75 images. The success rate and accuracy of centerline extraction was quantitatively evaluated on 48 coronary arteries in 12 images by comparing extracted centerlines with a manually annotated reference standard. The method was able to extract 88% and 47% of the vessel center-lines correctly using the vesselness/intensity and region statistics cost function, respectively. For only the proximal part of the vessels these values were 97% and 86%, respectively. Accuracy of centerline extraction, defined as the average distance from correctly automatically extracted parts of the centerline to the reference standard, was 0.64 mm for the vesselness/intensity and 0.51 mm for the region statistics cost function. The interobserver variability was 99% for the success rate measure and 0.42 mm for the accuracy measure. Qualitative evaluation using the best performing cost function resulted in successful centerline extraction for 233 out of the 252 coronaries (92%) in 63 additional CTCA images.

CONCLUSIONS

The presented results, in combination with minimal user interaction and low computation time, show that minimum cost path approaches can effectively be applied as a preprocessing step for subsequent analysis in clinical practice and biomedical research.

Standardized evaluation methodology and reference database for evaluating coronary artery centerline extraction algorithms

Efficiently obtaining a reliable coronary artery centerline from computed tomography angiography data is relevant in clinical practice. Whereas numerous methods have been presented for this purpose, up to now no standardized evaluation methodology has been published to reliably evaluate and compare the performance of the existing or newly developed coronary artery centerline extraction algorithms. This paper describes a standardized evaluation methodology and reference database for the quantitative evaluation of coronary artery centerline extraction algorithms. The contribution of this work is fourfold: (1) a method is described to create a consensus centerline with multiple observers, (2) well-defined measures are presented for the evaluation of coronary artery centerline extraction algorithms, (3) a database containing 32 cardiac CTA datasets with corresponding reference standard is described and made available, and (4) 13 coronary artery centerline extraction algorithms, implemented by different research groups, are quantitatively evaluated and compared. The presented evaluation framework is made available to the medical imaging community for benchmarking existing or newly developed coronary centerline extraction algorithms.