Automatic segmentation of the lumen of the carotid artery in ultrasound B-mode images

A new algorithm is proposed for the segmentation of the lumen and bifurcation boundaries of the carotid artery in B-mode ultrasound images. It uses the hipoechogenic characteristics of the lumen for the identification of the carotid boundaries and the echogenic characteristics for the identification of the bifurcation boundaries. The image to be segmented is processed with the application of an anisotropic diffusion filter for speckle removal and morphologic operators are employed in the detection of the artery. The obtained information is then used in the definition of two initial contours, one corresponding to the lumen and the other to the bifurcation boundaries, for the posterior application of the Chan-vese level set segmentation model. A set of longitudinal B-mode images of the common carotid artery (CCA) was acquired with a GE Healthcare Vivid-e ultrasound system (GE Healthcare, United Kingdom). All the acquired images include a part of the CCA and of the bifurcation that separates the CCA into the internal and external carotid arteries. In order to achieve the uppermost robustness in the imaging acquisition process, i.e., images with high contrast and low speckle noise, the scanner was adjusted differently for each acquisition and according to the medical exam. The obtained results prove that we were able to successfully apply a carotid segmentation technique based on cervical ultrasonography. The main advantage of the new segmentation method relies on the automatic identification of the carotid lumen, overcoming the limitations of the traditional methods

Vessel Axis Tracking Using Topology Constrained Surface Evolution

An approach to three-dimensional vessel axis tracking based on surface evolution is presented. The main idea is to guide the evolution of the surface by analyzing its skeleton topology during evolution, and imposing shape constraints on the topology. For example, the intermediate topology can be processed such that it represents a single vessel segment, a bifurcation, or a more complex vascular topology. The evolving surface is then re-initialized with the newly found topology. Re-initialization is a crucial step since it creates probing behavior of the evolving front, encourages the segmentation process to extract the vascular structure of interest and reduces the risk on leaking of the curve into the background. The method was evaluated in two computed tomography angiography applications: (i) extracting the internal carotid arteries including the region in which they traverse through the skull base, which is challenging due to the proximity of bone structures and overlap in intensity values, and (ii) extracting the carotid bifurcations including many cases in which they are severely stenosed and contain calcifications. The vessel axis was found in 90% (18/20 internal carotids in ten patients) and 70% (14/20 carotid bifurcations in a different set of ten patients) of the cases

Anatomical shape reconstruction and manufacturing: solving topological changes of lumen vessel trough geometric approach

Over the last years there has been an increasing growth of interest in Rapid Prototyping (RP) techniques applied to various fields of medicine. RP makes it possible, in vascular surgery, to produce accurate anatomic replicas of patient vessels. These replicas can help the customization of surgical invasive interventions such as in situ stent-graft insertion in carotid region. The main goal of this work is to obtain high quality in lumen reconstruction and manufacturing replicas by RP technique. This goal is achieved through the complete control of each phase of the generating process. We present a semi-automatic method for carotid lumen reconstruction based on Boundary Representation (BRep). All parameters influencing the quality of the shape reconstruction are presented and discussed: shape acquisition, shape reconstruction and shape manufacturing. The shape acquisition starts by extracting the points belonging to the boundary of the lumen vessel, from Computer Tomography (CT) images. These points, parameterised in a vector, are the input data of the shape reconstruction algorithm based on B-Spline interpolation. The B-Spline type for representing curves and surfaces were chosen because of their properties of continuity and local control. In the shape reconstruction stage we had to face problems due to the topological change on the vessel structure. For vessel regions where there are not changes of topology, we use the closed B-Spline curves (belonging to adjacent acquisition planes) as generating curves to build a B-Spline surface. For vessel regions with at least a change of topology (ex. bifurcation region) our algorithm split automatically the involved curves to obtain three rectangular B-Spline patches. Such patches are joined together to obtain the bifurcation vessel lumen. The set of lumen surfaces is then inserted in a Boundary Representation in order to get a valid solid. To analyse the quality of the reconstructed shapes, the final object is compared with the acquisition image. This solid is correctly tessellated in triangles to produce the data format used by the RP devices (STL)

Contrast-enhanced micro-CT imaging in murine carotid arteries : a new protocol for computing wall shear stress

Background: Wall shear stress (WSS) is involved in the pathophysiology of atherosclerosis. The correlation between WSS and atherosclerosis can be investigated over time using a WSS-manipulated atherosclerotic mouse model. To determine WSS in vivo, detailed 3D geometry of the vessel network is required. However, a protocol to reconstruct 3D murine vasculature using this animal model is lacking. In this project, we evaluated the adequacy of eXIA 160, a small animal contrast agent, for assessing murine vascular network on micro-CT. Also, a protocol was established for vessel geometry segmentation and WSS analysis. Methods: A tapering cast was placed around the right common carotid artery (RCCA) of ApoE(-/-) mice (n = 8). Contrast-enhanced micro-CT was performed using eXIA 160. An innovative local threshold-based segmentation procedure was implemented to reconstruct 3D geometry of the RCCA. The reconstructed RCCA was compared to the vessel geometry using a global threshold-based segmentation method. Computational fluid dynamics was applied to compute the velocity field and WSS distribution along the RCCA. Results: eXIA 160-enhanced micro-CT allowed clear visualization and assessment of the RCCA in all eight animals. No adverse biological effects were observed from the use of eXIA 160. Segmentation using local threshold values generated more accurate RCCA geometry than the global threshold-based approach. Mouse-specific velocity data and the RCCA geometry generated 3D WSS maps with high resolution, enabling quantitative analysis of WSS. In all animals, we observed low WSS upstream of the cast. Downstream of the cast, asymmetric WSS patterns were revealed with variation in size and location between animals. Conclusions: eXIA 160 provided good contrast to reconstruct 3D vessel geometry and determine WSS patterns in the RCCA of the atherosclerotic mouse model. We established a novel local threshold-based segmentation protocol for RCCA reconstruction and WSS computation. The observed differences between animals indicate the necessity to use mouse-specific data for WSS analysis. For our future work, our protocol makes it possible to study in vivo WSS longitudinally over a growing plaque