BioMed Research International

Volume 2013 (2013), Article ID 801962, 10 pages

http://dx.doi.org/10.1155/2013/801962

Review Article

Carotid Artery Segmentation in Ultrasound Images and Measurement of Intima-Media Thickness

Vaishali Naik, R. S. Gamad, and P. P. Bansod

Department of Electronics and Instrumentation Engineering, Shri Govindram Seksaria Institute of Technology and Science, Indore 23, Park Road, Indore 452003, India

Received 17 April 2013; Accepted 28 May 2013

Academic Editor: Manuel F. Casanova

Copyright © 2013 Vaishali Naik et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. The segmentation of the common carotid artery (CCA) wall is imperative for the determination of the intima-media thickness (IMT) on B-mode ultrasound (US) images. The IMT is considered an important indicator in the evaluation of the risk for the development of atherosclerosis. In this paper, authors have discussed the relevance of measurements in clinical practices and the challenges that one has to face while approaching the segmentation of carotid artery on ultrasound images. The paper presents an overall review of commonly used methods for the CCA segmentation and IMT measurement along with the different performance metrics that have been proposed and used for performance validation. Summary and future directions are given in the conclusion.

1. Introduction

Cardiovascular disease (CVD) is one of the leading causes of deaths in the metropolitan cities. A recent survey by the World Health Organization revealed that up to 17.3 million people died from CVDs in 2008. Scientifically presaged the upcoming 2030, almost 23.3 million human deaths are resulting from CVDs [1]. CVD ailments are relative to atherosclerosis (arterial disease). Atherosclerosis is responsible for the thickening of the artery walls, and the IMT is used as a validated measure for the evaluation of atherosclerosis. Prominently, the increase in IMT jeopardizes the brain infarction or cardiac attack [2, 3].

Usually, the B-mode ultrasound scan of CCA in its common tract is used for the evaluation of the artery status and for the measurement of the IMT. Ultrasound methodology has manifest benefits of being real-time, noninvasive, low-cost, reliable, and absolutely safe for the patients. The essential drawbacks of this methodology are the B-mode image having low signal-to-noise ratio and ultrasounds are operator dependent [4]. Conventionally, the IMT is manually measured by the trained operator from the US scan images. This methodology is highly user dependent, time consuming, tedious, and infeasible in presence of large image databases [5, 6]. During the past 20 years, several computerized techniques have been developed for segmentation of CCA. These methods can be broadly classified into two categories: first category includes techniques that are completely automated, whereas second category includes those that require user interaction (semiautomated). This review will focus on the techniques that have been developed to perform CCA wall segmentation and IMT measurement in B-mode ultrasound images in both automated and semiautomated manners.

2. Clinical Significance of Vessel Wall Segmentation

2.1. Common Carotid Artery Intima-Media Thickness

The CCA longitudinal section is shown in Figure 1. It is characterized by a longitudinal tract (common carotid) that, after an enlargement (carotid sinus) containing a flow divider, bifurcates into two arteries, one internal and one external, on the basis of their position in relation to neck skin [7].

Figure 1: The carotid artery view with the interfaces: (1) periadventitia-adventitia (NW), (2) adventitia-media (NW), (3) intima-lumen (NW), (4) lumen-intima (FW), (5) media-adventitia (FW), and (6) adventitia-periadventitia (FW) [7].

The bifurcation and the internal carotid artery (ICA) are more threatening to atherosclerosis, due to stronger hemodynamic stresses in the bifurcation and branching zones, but it is difficult to visualize the “double-line” pattern in these locations. So an IMT measurement in the CCA is preferred in the development of segmentation algorithms and in clinical practice [8].

Classically, IMT is defined as a double-line pattern visualized by echotomography on both walls of the CCA in a longitudinal image. It consists of two parallel anatomical boundaries referred to as the lumen-intima and media-adventitia interfaces [9]. US waves are reflected differently by blood (vessel lumen) and wall layers because of their differences in density and elasticity. US waves are not reflected by vessel lumen and tunica media, thus allow detection of the lumen-intima (LI) and media-adventitia (MA) interfaces, as depicted in Figure 2 B-mode ultrasound CCA image [10].

Figure 2: B-mode CCA image. The media layer thickness (MLT) is defined as the distance between the intima-media and the media-adventitia interface [11].

2.2. Importance of Carotid Arterial Intima-Media Thickness

The increased IMT reflects early stages of atherosclerosis and cardiovascular risk. Higher blood pressure and changes in shear stress are the potential causes of intimal thickening. Changes in shear stress and blood pressure may cause a local delay in lumen transportation of potentially atherogenic particles, which favors the accumulation of particles in the arterial wall and consequent plaque formation [12]. It is found that type 1 diabetes is a significant risk factor for increased carotid IMT in children [13]. It is confirmed that the increase in IMT is directly associated with an increase risk of myocardial infarction and stroke in older adults without a history of cardiovascular disease [14]. The assessment of IMT in prediction of the degrees of atherosclerosis and the risk of stroke and CVDs has been demonstrated by a lot of studies [15–17]. It is found that obesity especially abdominal obesity in childhood and adolescence is closely related to IMT [18]. The studies confirmed that the increase of the IMT value above 0.9-1.0 mm is indicative of a significant increase of CVD risk when taking into account the population of healthy elderly [4, 7]. In the following sections, we will discuss the difficulties confronted in the intima and media detection and review the recent advancements in CCA segmentation.

3. Difficulties of Intima and Adventitia Detection

The difficulties encountered in detecting intima and adventitia layers are as follows [19]:

(i)presence of speckles in US image;(ii)the structure of IM complex or the intimal layer changes due to diseases such as atherosclerotic plaques; (iii)variation in echo characteristic on intima and adventitia on images with the variation in sonographic instrumentations.

4. Segmentation Techniques

Plethora of ultrasound-segmentation techniques have been reviewed in recent surveys by Noble and Boukerroui [29] and Molinari et al. [30]. Thereby authors are steadfast in reporting on recent techniques, used for the segmentation of the CCA on US images. These studies on the segmentation of the carotid artery boundaries include the application of dynamic programming, deformable snakes, hough transforms, and classification approaches to detect the carotid boundaries on longitudinally oriented images [27, 31–33]. For each methodology, the authors have described principles, performance, advantages, and limitations. Table 1 gives and summarizes the technique that is mentioned below.

Table 1: Overview of recent CCA IMC segmentation techniques for ultrasound imaging.

Name,

year, and

ref. no.
Common carotid artery IMT segmentation technique
Advantages & limitations of the methods
Selection method of ROI

(SA/FA)
Performance metric
Processing

time/frame 
IMT error

(mm)

Ilea et al., 2013

[20]
Model-based approach, video tracking procedure: spatially coherent algorithm
Advantages: Method can deal with inconsistencies in the appearance of the IMC over the cardiac cycle. Robustness with respect to data captured under different imaging conditions.
FA
MAD
80 sec

(8 sec in 1st frame & 72 sec in tracking 28 frames)
(−0.007) ± 0.176
40*

Xu et al.,

2012,

[21]
Hough transform and dual snake model
Advantages: It is less likely to be affected by noise, compensates the holes or missing boundaries, and can estimate the missing LI interface boundary.

Limitations: The method would work fine for early thickening of IMC but fails for irregular boundaries in the presence of plaques and eliminates minor details.
SA
MAD
0.465 sec
0.02 ± 0.03
50

Molinari et al.,

2012,

[22]
Multiresolution edge snapper
Advantages: Complete automation, robustness to noise, and real-time computation.

Limitations: Robustness with respect to noise, but the LI/MA representation is less accurate.
FA
MAD
Less than 15 sec
0.078 ± 0.112
365

Destrempes et al.,

2011,

[23]
Nakagami distributions, Bayesian model
Advantages: Robust to a reasonable variability in the initialization, lowest tracing error for LI & MA, method is not sensitive to the degree of stenosis or calcification.

Limitations: Depends on initial segmentation
SA
MAD

HD
38 sec
—
8988

Petroudi et al., 2011, [24]
Active contours & active contours without edges
Advantages: Fully automated, fast, does not require any user interaction, and works well for noisy images.
FA
MAD
—
0.09 ± 0.10
30

Destrempes et al.,

2009,

[25]
Nakagami distributions, stochastic optimization
Advantages: Reasonable average computation time, robust to the estimation procedures.

Limitations: Method suitable for healthy arteries, extensive tuning & training, so computational cost is high, for different scanner requires retraining & retuning.
SA
MAD

HD
24 sec
—
7283

Faita et al.,

2008, [26]
First-order absolute moment edge operator
Advantages: Suited for fast real-time implementation, operator can have immediate feedback on the quality of the images.

Limitations: Depends on ROI selection.
SA
MAD
—
0.001 ± 0.035
150

Liang et al.,

2000,

[27]
Multiscale dynamic programming
Advantages: No initial human setting, capable of processing images of different quality, ambiguous cases user can intervene, and reduced interobserver variability.

Limitations: Training required, for different scanner retraining needed, searched LI & MA interfaces may cross each other.
FA
MAD
0.7 min
0.042 ± 0.02
50

Gustavsson et al., 1994,

[28]
Dynamic programming
Advantages: Fully automated, low computational complexity; suitable for clinical purposes, human correction allowed.

Limitations: Initial human setting & training required, fails for slanting IMC with weak boundary.
SA
MAD
—
0.03 ± 0.032
22

: number of images/cases, SA: semiautomated, FA: fully automated, SD: standard deviation, HD: Hausdorff distance, MAD: mean absolute distance, *video sequences, : average processing time/frame or image.

4.1. Dynamic Programming Techniques (DP)

DP technique is concisely used to solve optimization problems, where desired segmentation minimizes the cost function defined by the particular application. Local measurements of echo intensity, edge strength, and boundary continuity are included as weighted terms in a cost function. All possible set of spatially consecutive points forming a polyline is being considered, and favor is given to that which minimizes the cost function [34].

The polylines are represented as a vector:

where is the horizontal pixel position, and are neighbor points, and is the horizontal length of the search region. A vertical search window of size pixels is used to scan the boundary from left to right at horizontal positions as shown in Figure 3. At scan position , the boundary point, , in (1) can be any pixel in this window. The optimized connection is searched for each point in this window and the cost accumulated. At the end of the scanning, the optimal polyline is the one that minimizes the cost function:

The local cost is a weighted sum of cost terms:

where , , and are weighting factors, , , and are the echo intensity, intensity gradient, and boundary continuity cost terms, respectively. Based on DP techniques in 1994 Gustavsson et al. [28] introduced a procedure for automatic ultrasonic measurements of the carotid artery, and lumen diameter (LD) and IMT were computed. Intermethod (auto versus manual) variability as well as inter- and intraobserver variability was studied by computing the conventional coefficient of variation (CV). A major advantage of this methodology was complete automation and low computational complexity, thus suitable for clinical purposes. This method requires interactive tools for manual tracing in order to correct the remaining detection errors. The major limitation of this technique is the need for training of the system. In 1997, Gustavsson et al. [34] have compared four algorithms: the dynamic programming, the maximum gradient, the model-based, and the matched filter algorithm and confirmed that the DP algorithm provides superior performance in terms of accuracy and robustness. In 2008 Liu et al. [35] proposed a segmentation method in which the energy definition of active contour model was used and DP was employed to search the shortest path. To reduce the effects of speckle noise, anisotropic diffusion method was adopted. It is advantageous as it requires less manual input. Holdfeldt et al. [36] proposed a method based on DP for boundary detection in ultrasound image sequences. According to the author, this method gives favorable results on both synthetic and real ultrasound data. Cheng and Jiang [37] proposed a novel dual dynamic programming (DDP) technique that detected intimal and adventitial layers of the CCA of the B-mode US images. In this, the robustness against the speckles was increased by embedding the anatomical knowledge into its structure. Therefore, the researcher reported that the DDP technique achieved a detection performance comparable to manual segmentation.

Figure 3: Detecting interfaces I2 and I7 in an artery image.

4.2. Hough Transform (HT)

HT technique used to detect straight lines. A straight line at a distance and orientation θ can be represented by

The HT of this line is just a point in the () plane; that is all the points on this line map into a single point. This fact is utilized to detect straight lines in a given set of boundary points. If given boundary points are (), for some selected quantized values of parameter and , map each () into the () space and count , the number of edge points that map into the location (); that is, set

Then the local maxima of give the different straight line segments through the edge points. Generalized HT can be used to detect curves other than the straight lines [38]. Segmentation algorithm based on the HT was demonstrated by Golemati et al. [39] in 2007 to segment both longitudinal and transverse images. The HT is effective in detecting lines (longitudinal images) or circles (transverse images), but it may fail in detecting curved vessels. In 2008, Stoitsis et al. [40] proposed the HT-initialized active contour methodology. In 2010, Petroudi et al. [41] proposed a fully automated method that was proposed for the delineation of the intima-media complex (IMC). In this technique after speckle removal and HT used for boundary detection followed by image normalization, the corresponding results were used to provide the initial statistical information needed for a Markov random field (MRF). In 2011, Matsakou et al. [42] proposed a method in which an HT-based methodology was used for the definition of the initial snake followed by a gradient vector flow (GVF) snake deformation for the final contour detection. The author reported that the sensitivity, specificity, and accuracy were 0.97, 0.99 and 0.98, respectively, for both diastolic and systolic cases. Recently Xu et al. [21] proposed a segmentation method using HT and dual snake model; two contours are initialized from line segments generated by HT. Author admits that the technique is not suitable for irregular boundaries and decimate minor details.

4.3. Nakagami Mixture Modelling

In 2009, Destrempes et al. [25] introduced a segmentation technique based on Nakagami mixture modeling and stochastic optimization. The echogenicity of the region of interest (ROI) comprising the intima-media layers, the lumen, and the adventitia in an ultrasonic B-mode image is modeled by a mixture of three Nakagami distributions. In a first step the expectation maximization (EM) algorithm was used to compute the maximum A posterior (MAP) estimator of the proposed model, then computes the optimal segmentation based on the estimated distributions as well as a statistical prior for disease-free IMT using a variant of the exploration/selection (ES) algorithm. This method requires minimal manual initialization. Destrempes et al. [43] proposed a method for segmentation of plaques in sequences of ultrasound B-mode images of carotid arteries based on motion estimation and Nakagami distributions. In it, a local geometrical smoothness constraint and an original spatiotemporal cohesion constraint were incorporated, envisaging the segmented plaque based on motion field estimation. In 2011, Destrempes et al. [23] proposed a method for the segmentation of plaques in the sequence of ultrasound B-mode images of carotid arteries based on motion estimation and a Bayesian model. Authors have reported that the algorithm was not sensitive to the degree of stenosis or calcification.

4.4. Active Contour

The basic concept of active contour model is to fit a contour to local image information, for example, gradient. There exist several implementations of this basic idea such as snakes [44], discrete dynamic contour model [45], and level sets [46]. Based on the involved feature image, they can be categorized as edge based [47], region based [48, 49], and higher level knowledge based [50, 51]. Several studies have adopted the traditional snake model as proposed by Williams and Shah [52]. Snakes are also called active parametric contours, which have been widely used in medical image segmentation. The major limitations of this method are sensitive to noise, depends on the initial contour that is provided by the user, need for optimization of the parameters. In 2008, Moursi and El-Sakka [53] proposed an active contour-based segmentation technique, in which user only requires to place seed points in the ROI with the aim of reducing user interaction. Author admitted that the computational time depends on the size of the carotid artery and the location of the seed point. In 2010, Bastida-Jumilla et al. [54] used geodesic active contours for IMC detection. In 2011, Petroudi et al. [24] proposed the fully automated segmentation algorithm based on active contours, and active contours without edges were proposed in which anatomical information was incorporated to achieve accurate segmentation. The segmented regions were used to automatically achieve image normalization, which is followed by speckle removal. The resulting smoothed LI boundary combined with anatomical information provides an excellent initialization for parametric active contours that provide the final IMC segmentation. No information about an inter- and intraobserver variability and its effect on segmentations was given by author.

4.5. Edge Detection and Gradient-Based Techniques

The edge detection methods could detect the variation of gray levels, but it is sensitive to noise and may suffer from the focusing artifact. In 2001, Liguori et al. [7] proposed the segmentation technique based on an edge detection, in which image gradient was used. In this method, for each column of the image the gradient of the intensity profile has been computed. It was assumed that pixels belonging to lumen were black and that the carotid wall layers originate with gradient transitions. It is a semiautomatic method; ROI is selected by the user. The pattern recognition, edge detection (PRED) algorithm, and the measurement algorithm were used for carotid IMT measurement. Its main task is to find out all the pixels belonging to the two required interfaces (LI and MA) for each wall. The measured intensity gradient was different from the theoretical one, due to noise. In order to reduce the effect of noise, a statistical thresholding was adopted before computing the image gradient. Robustness of the edge detection algorithm had been evaluated with respect to the ROI. The gradient-based segmentation mainly suffers from the problem of superimposed noise, which precludes a proper individuation of the LI and MA transitions. In 2008, Faita et al. [26] proposed a method in which the gradient performance was improved by the use of a first-order absolute moment edge operator (FOAM) and a pattern recognition approach. The overall performance of this methodology was very high: IMT measurement error was equal to  μm. Moreover, FOAM operator and intelligent procedure determines maxima, ensuring a good robustness to noise. As the technique is real-time, it suits well to clinical application. It is a semiautomatic technique. Recently, Mahmoud et al. [55] introduced a method, which employs a multistep gradient-based algorithm. This method principally uses intensity, intensity gradient, and interface continuity of pixels to determine the ultrasound interface. Author reported that this technique eliminates subjectivity associated with conventional manual tracing and semiautomated gradient methods that employ seed point selection.

4.6. Combined Approaches

Delsanto et al. [6] proposed a combined approach for classification and a snake-based segmentation to perform IMT measurement. Completely user-independent layer extraction based on signal analysis (CULEXsa) is a completely user independent algorithm for IMT measurement. Firstly ROI is identified automatically followed by gradient-based initial segmentation, and then active contour technique is used for segmentation refinement. For improving segmentation performance, Molinari et al. [56] have combined the three IMT segmentation methods: (i) signal processing approach, combined with snakes and fuzzy clustering, (ii) integrated approach based on seed and line detection, followed by probability-based connectivity and classification, and (iii) morphological approach and fused the resulting boundaries using a greedy method described by the “ball and basket” to minimize the system error. In 2010, Molinari et al. [57] proposed a completely automated layer extraction technique (named CALEXia). The IMT measurement error was equal to pixels ( mm). Author admitted that CALEXia showed limited performance in segmenting the LI interface. Meiburger et al. [58] introduced the Carotid Automated Double-Line Extraction System based on the Edge Flow (CADLES-EF). It is characterized and validated by comparing the output of the algorithm with CALEXia and CULEXsa. Validation was performed on a multi-institutional database of 300 longitudinal B-mode carotid images with normal and pathologic arteries. CADLES-EF showed an IMT bias of  mm in comparison with CALEXia and CULEXsa that showed  mm and  mm, respectively. The system’s Figure of Merit (FoM) showed an improvement when compared with CALEXia and CULEXsa, leading to values of 84.7%, 91.5%, while CADLES-EF performed the best with 94.8%. In 2011, Molinari et al. [59] proposed a method called CARES 3.0 (a patented technology). CARES 3.0 is completely automated and adopts an integrated approach for segmentation of carotid artery in the image frame. The FoM of CARES 3.0 was 97.4%. In 2012, Molinari et al. [22] proposed completely automated multiresolution edge snapper (named CAMES). In it carotid artery is recognized automatically using a scale-space and statistical classification in a multiresolution framework. Recently, Ilea et al. [20] proposed a fully automated segmentation and tracking of the intima-media thickness in ultrasound video sequences of the CCA. For the video tracking procedure, a spatially coherent algorithm is introduced, which prevents the tracking process from converging to wrong arterial interfaces. Author reported that method can deal with inconsistencies in the appearance of the IMC over the cardiac cycle.

5. Validation/Quantitative Performance Assessment

Validation experiments are necessary in order to quantify the performance of a segmentation method. Validation is usually performed using truth models such as phantom studies, animal model studies, simulation, comparing the automated segmentation method with manually obtained segmentations. Following are the most used performance metrics to validate IMT measurements and computer traced boundaries [30]:

(i)Mean absolute distance (MAD),(ii)Hausdorff distance (HD),(iii)Polyline distance metric (PDM),(iv)Percent statistic test,(v)Reproducibility of manual procedures. It is assessed by calculating intraoperator and interoperator variability using either of CV, MAD, HD, regression analysis, and Bland-Altman statistics,(vi)Manual and computer-measured IMT (intermethod). Comparison between the two sets was done using correlation or Bland-Altman plot.

6. Conclusion

This paper reports an extensive review of ultrasound carotid artery IMT segmentation techniques. Active contours, dynamic programming, and integrated approaches have been presented to segment the carotid wall and trace the boundaries of the LI and MA interfaces. None of the existing techniques were overwhelmingly good in all aspects. Characterization and validation studies will be required in order to carefully assess the effect of such variability on segmentation performance. Finally, we recognize that in the future, more work is likely to be done in segmentation based on adaptive segmentation for determination of IMT in ultrasound images of CCA with high IMT measurement accuracy, robustness, automation and reducing processing time.

As in the case of fully automatic techniques, detection is not reliable, since it may detect the jugular veins edges. In addition performance of semiautomatic segmentation techniques is better than fully automatic segmentation techniques. Therefore, in future we will develop techniques in which human operator will select an ROI manually, and methodology will be based on adaptive segmentation with the aim of high accuracy, great robustness, and with reducing processing time. The proposed accuracy of detection of IMT algorithm falls within the inter- and intraobserver variability for the manual determination.

References

World Health Organization, “Cardiovascular Disease,” http://www.who.int/cardiovascular_diseases/en/

L.-N. Pu, Z. Zhao, and Y.-T. Zhang, “Investigation on cardiovascular risk prediction using genetic information,” IEEE Information Technology in Biomedicine, vol. 16, pp. 795–808, 2012.

P. Touboul, A. Elbaz, C. Koller et al., “Common carotid artery intima-media thickness and brain infarction: the etude du profil genetique de l'infarctus cerebral (GENIC) case-control study,” Circulation, vol. 102, no. 3, pp. 313–318, 2000.

S. Delsanto, F. Molinari, P. Giustetto, W. Liboni, S. Badalamenti, and J. S. Suri, “Characterization of a completely user-independent algorithm for carotid artery segmentation in 2-D ultrasound images,” IEEE Transactions on Instrumentation and Measurement, vol. 56, no. 4, pp. 1265–1274, 2007.

N. Santhiyakumari and M. Madheswaran, “Extraction of intima- media layer of arteria-carotis and evaluation of its thickness using active contour approach,” in Proceedings of the International Conference on Intelligent and Advanced Systems (ICIAS '07), pp. 582–586, November 2007.

S. Delsanto, F. Molinari, P. Giustetto, W. Liboni, and S. Badalamenti, “CULEX-completely user-independent layers EXtraction: ultrasonic carotid artery images segmentation,” in Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEE-EMBS '05), pp. 6468–6471, September 2005.

C. Liguori, A. Paolillo, and A. Pietrosanto, “An automatic measurement system for the evaluation of carotid intima-media thickness,” IEEE Transactions on Instrumentation and Measurement, vol. 50, no. 6, pp. 1684–1691, 2001.

M. Halenka, “Noninvasive measurement of early atherosclerosis by high-resolution B-mode ultrasonography,” Acta Universitatis Palackianae Olomucensis Facultatis Medicae, vol. 142, pp. 7–11, 1999.

P.-J. Touboul, M. G. Hennerici, S. Meairs et al., “Mannheim carotid intima-media thickness consensus (2004–2006),” Cerebrovascular Diseases, vol. 23, no. 1, pp. 75–80, 2007.

F. Faita, V. Gemignani, E. Bianehini, C. Giannarelli, and M. Demi, “Real-time measurement system for the evaluation of the Intima Media Thickness with a new edge detector,” in Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '06), pp. 715–718, September 2006.

C. P. Loizou, C. S. Pattichis, A. N. Nicolaides, and M. Pantziaris, “Manual and automated media and intima thickness measurements of the common carotid artery,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 56, no. 5, pp. 983–994, 2009.

I. M. Van Der Meer, M. L. Bots, A. Hofman, A. I. Del Sol, D. A. M. Van Der Kuip, and J. C. M. Witteman, “Predictive value of noninvasive measures of atherosclerosis for incident myocardial infarction: The Rotterdam Study,” Circulation, vol. 109, no. 9, pp. 1089–1094, 2004.

M. J. Järvisalo, A. Putto-Laurila, L. Jartti et al., “Carotid artery intima-media thickness in children with type 1 diabetes,” Diabetes, vol. 51, no. 2, pp. 493–498, 2002.

D. H. O'Leary, J. F. Polak, R. A. Kronmal, T. A. Manolio, G. L. Burke, and S. K. Wolfson Jr., “Carotid-artery intima and media thickness as a risk factor for myocardial infarction and stroke in older adults,” The New England Journal of Medicine, vol. 340, no. 1, pp. 14–22, 1999.

C.-C. Huang, T. Chen, T.-H. Chao, and Y.-F. Juan, “The investigation of the relationship between carotid intima-media thickness and vascular compliance in patients with coronary artery disease,” Biomedical Engineering, vol. 16, no. 1, pp. 37–42, 2004.

M. W. Lorenz, H. S. Markus, M. L. Bots, M. Rosvall, and M. Sitzer, “Prediction of clinical cardiovascular events with carotid intima-media thickness: a systematic review and meta-analysis,” Circulation, vol. 115, no. 4, pp. 459–467, 2007.

H. N. Hodis, W. J. Mack, L. LaBree et al., “The role of carotid arterial intima—media thickness in predicting clinical coronary events,” Annals of Internal Medicine, vol. 128, no. 4, pp. 262–269, 1998.

J. Fang, J. P. Zhang, C. X. Luo, X. M. Yu, and L. Q. Lv, “Carotid intima-media thickness in childhood and adolescent obesity relations to abdominal obesity, high triglyceride level and insulin resistance,” International Journal of Medical Sciences, vol. 7, no. 5, pp. 278–283, 2010.

D.-C. Cheng, X. Jiang, A. Schmidt-Trucksäss, and K. Cheng, “Automatic intima-media thickness measurement of carotid artery wall in B-mode sonographic images,” in Proceedings of the 3rd IEEE International Symposium on Biomedical Imaging: from Nano to Macro, pp. 912–915, April 2006.

D. E. Ilea, C. Duffy, L. Kavanagh, A. Stanton, and P. F. Whelan, “Fully automated segmentation and tracking of the intima media thickness in ultrasound video sequences of the common carotid artery,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 60, pp. 158–177, 2013.

X. Xu, Y. Zhou, X. Cheng, E. Song, and G. Li, “Ultrasound intima-media segmentation using Hough transform and dual snake model,” Computerized Medical Imaging and Graphics, vol. 36, no. 3, pp. 248–258, 2012.

F. Molinari, C. S. Pattichis, G. Zeng et al., “Completely automated multiresolution edge snapper—a new technique for an accurate carotid ultrasound IMT measurement: clinical validation and benchmarking on a multi-institutional database,” IEEE Transactions on Image Processing, vol. 21, no. 3, pp. 1211–1222, 2012.

F. Destrempes, J. Meunier, M. Giroux, G. Soulez, and G. Cloutier, “Segmentation of plaques in sequences of ultrasonic B-mode images of carotid arteries based on motion estimation and a bayesian model,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 8, pp. 2202–2211, 2011.

S. Petroudi, C. Loizou, M. Pantziaris, M. Pattichis, and C. Pattichis, “A fully automated method using active contours for the evaluation of the intima-media thickness in carotid US images,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 8053–8057, September 2011.

F. Destrempes, J. Meunier, M. Giroux, G. Soulez, and G. Cloutier, “Segmentation in ultrasonic B-mode images of healthy carotid arteries using mixtures of Nakagami distributions and stochastic optimization,” IEEE Transactions on Medical Imaging, vol. 28, no. 2, pp. 215–229, 2009.

F. Faita, V. Gemignani, E. Bianchini, C. Giannarelli, L. Ghiadoni, and M. Demi, “Real-time measurement system for evaluation of the carotid intima-media thickness with a robust edge operator,” Journal of Ultrasound in Medicine, vol. 27, no. 9, pp. 1353–1361, 2008.

Q. Liang, I. Wendelhag, J. Wikstrand, and T. Gustavsson, “A multiscale dynamic programming procedure for boundary detection in ultrasonic artery images,” IEEE Transactions on Medical Imaging, vol. 19, no. 2, pp. 127–142, 2000.

T. Gustavsson, Q. Liang, I. Wendelhag, and J. Wikstrand, “A dynamic programming procedure for automated ultrasonic measurement of the carotid artery,” in Proceedings of the IEEE Computers in Cardiology Conference, pp. 297–300, 1994.

J. A. Noble and D. Boukerroui, “Ultrasound image segmentation: a survey,” IEEE Transactions on Medical Imaging, vol. 25, no. 8, pp. 987–1010, 2006.

F. Molinari, G. Zeng, and J. S. Suri, “A state of the art review on intima-media thickness (IMT) measurement and wall segmentation techniques for carotid ultrasound,” Computer Methods and Programs in Biomedicine, vol. 100, no. 3, pp. 201–221, 2010.

C. P. Loizou, C. S. Pattichis, M. Pantziaris, T. Tyllis, and A. Nicolaides, “Snakes based segmentation of the common carotid artery intima media,” Medical and Biological Engineering and Computing, vol. 45, no. 1, pp. 35–49, 2007.

S. Golemati, J. Stoitsis, T. Balkizas, and K. S. Nikita, “Comparison of B-mode, M-mode and Hough transform methods for measurement of arterial diastolic and systolic diameters,” in Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEE-EMBS '05), pp. 1758–1761, September 2005.

F. Molinari, G. Zeng, and J. S. Suri, “An integrated approach to computer-based automated tracing and its validation for 200 common carotid arterial wall ultrasound images: a new technique,” Journal of Ultrasound in Medicine, vol. 29, no. 3, pp. 399–418, 2010.

T. Gustavsson, R. Abu-Gharbieh, G. Hamarneh, and Q. Liang, “Implementation and comparison of four different boundary detection algorithms for quantitative ultrasonic measurements of the human carotid artery,” in Proceedings of the 24th Annual Meeting on Computers in Cardiology, pp. 69–72, September 1997.

G. Liu, B. Wang, and D. C. Liu, “Detection of intima-media layer of common carotid artery with dynamic programming based active contour model,” in Proceedings of the Chinese Conference on Pattern Recognition (CCPR '08), pp. 1–6, October 2008.

P. Holdfeldt, M. Viberg, and T. Gustavsson, “A new method based on dynamic programming for boundary detection in ultrasound image sequences,” Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society Conference, vol. 2008, pp. 3072–3074, 2008.

D.-C. Cheng and X. Jiang, “Detections of arterial wall in sonographic artery images using dual dynamic programming,” IEEE Transactions on Information Technology in Biomedicine, vol. 12, no. 6, pp. 792–799, 2008.

A. K. Jain, Fundamentals of Digital Image Processing, PHI.

S. Golemati, J. Stoitsis, E. G. Sifakis, T. Balkizas, and K. S. Nikita, “Using the hough transform to segment ultrasound images of longitudinal and transverse sections of the carotid artery,” Ultrasound in Medicine and Biology, vol. 33, no. 12, pp. 1918–1932, 2007.

J. Stoitsis, S. Golemati, S. Kendros, and K. S. Nikita, “Automated detection of the carotid artery wall in B-mode ultrasound images using active contours initialized by the Hough Transform,” Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 2008, pp. 3146–3149, 2008.

S. Petroudi, C. P. Loizou, and C. S. Pattichis, “Atherosclerotic carotid wall segmentation in ultrasound images using Markov random fields,” in Proceedings of the 10th International Conference on Information Technology and Applications in Biomedicine (ITAB '10), pp. 1–5, November 2010.

A. I. Matsakou, S. Golemati, J. S. Stoitsis, and K. S. Nikita, “Automated detection of the carotid artery wall in longitudinal B-mode images using active contours initialized by the Hough transform,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 571–574, September 2011.

F. Destrempes, G. Soulez, M. Giroux, J. Meunier, and G. Cloutier, “Segmentation of plaques in sequences of ultrasonic B-mode images of carotid arteries based on motion estimation and Nakagami distributions,” in Proceedings of the IEEE International Ultrasonics Symposium (IUS '09), pp. 2480–2483, September 2009.

M. Kass, A. Witkin, and D. Terzopoulos, “Snakes: active contour models,” International Journal of Computer Vision, vol. 1, no. 4, pp. 321–331, 1988.

S. Lobregt and M. A. Viergever, “Discrete dynamic contour model,” IEEE Transactions on Medical Imaging, vol. 14, no. 1, pp. 12–24, 1995.

J. A. Sethian, Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision and Materials Science, Cambridge University, 1999.

V. Caselles, R. Kimmel, and G. Sapiro, “Geodesic active contours,” International Journal of Computer Vision, vol. 22, no. 1, pp. 61–79, 1997.

C. Li, C. Kao, J. C. Gore, and Z. Ding, “Minimization of region-scalable fitting energy for image segmentation,” IEEE Transactions on Image Processing, vol. 17, no. 10, pp. 1940–1949, 2008.

S. C. Zhu, “Region competition: unifying snakes, region growing, and bayes/mdl for multiband image segmentation,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 18, no. 9, pp. 884–900, 1996.

D. Cremers, F. Tischhäuser, J. Weickert, and C. Schnörr, “Diffusion snakes: Introducing statistical shape knowledge into the Mumford-Shah functional,” International Journal of Computer Vision, vol. 50, no. 3, pp. 295–313, 2002.

E. L. Michael, W. E. L. Grimson, and F. Olivier, “Statistical shape influence in geodesic active contours,” Computer Vision and Image Understanding, vol. 1, pp. 316–323, 2000.

D. J. Williams and M. Shah, “A Fast algorithm for active contours and curvature estimation,” CVGIP, vol. 55, no. 1, pp. 14–26, 1992.

S. G. Moursi and M. R. El-Sakka, “Active contours initialization for ultrasound carotid artery images,” in Proceedings of the 6th IEEE/ACS International Conference on Computer Systems and Applications (AICCSA '08), pp. 629–636, April 2008.

M. C. Bastida-Jumilla, J. Morales-Sánchez, R. Verdú-Monedero, J. Larrey-Ruiz, and J. L. Sancho-Gómez, “Detection of the intima and media walls of the carotid artery with geodesic active contours,” in Proceedings of the 17th IEEE International Conference on Image Processing (ICIP '10), pp. 2213–2216, September 2010.

A. Mahmoud, A. Morsy, and E. De Groot, “A new gradient-based algorithm for edge detection in ultrasonic carotid artery images,” in Proceedings of the 32nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '10), pp. 5165–5168, September 2010.

F. Molinari, G. Zeng, and J. S. Suri, “Intima-media thickness: setting a standard for a completely automated method of ultrasound measurement,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 57, no. 5, pp. 1112–1124, 2010.

F. Molinari, G. Zeng, and J. S. Suri, “Inter-Greedy technique for fusion of different carotid segmentation boundaries leading to high-performance IMT measurement,” in Proceedings of the 32nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '10), pp. 4769–4772, September 2010.

K. M. Meiburger, F. Molinari, G. Zeng, L. Saba, and J. S. Suri, “Carotid automated ultrasound double line extraction system (CADLES) via Edge-Flow,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 575–578, September 2011.

F. Molinari, K. M. Meiburger, U. R. Acharya et al., “CARES 3.0: a two stage system combining feature-based recognition and edge-based segmentation for CIMT measurement on a multi-institutional ultrasound database of 300 images,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 5149–5152, September 2011.