Research Article  Open Access
An Improved Approach for Accurate and Efficient Measurement of Common Carotid Artery IntimaMedia Thickness in Ultrasound Images
Abstract
The intimamedia thickness (IMT) of common carotid artery (CCA) can serve as an important indicator for the assessment of cardiovascular diseases (CVDs). In this paper an improved approach for automatic IMT measurement with low complexity and high accuracy is presented. 100 ultrasound images from 100 patients were tested with the proposed approach. The ground truth (GT) of the IMT was manually measured for six times and averaged, while the automatic segmented (AS) IMT was computed by the algorithm proposed in this paper. The mean difference ± standard deviation between AS and GT IMT is 0.0231 ± 0.0348 mm, and the correlation coefficient between them is 0.9629. The computational time is 0.3223 s per image with MATLAB under Windows XP on an Intel Core 2 Duo CPU E7500 @2.93 GHz. The proposed algorithm has the potential to achieve realtime measurement under Visual Studio.
1. Introduction
In 2011, a report by the World Health Organization [1] revealed that cardiovascular diseases (CVDs) are the number one cause of death globally. An estimated 17.3 million people died from CVDs in 2008, representing 30% of all global deaths. The number of people who die from CVDs will increase to 23.3 million by 2030. Therefore, a growing body of studies is looking for an early diagnosis and treatment of cardiovascular diseases, which is crucial to prevent patients from suffering more serious pathologies.
Several large populationbased studies [2–5] have shown that the intimamedia thickness (IMT) of common carotid artery can serve as an important indicator for cardiovascular diseases at an early stage and can also be used to predict major cardiovascular events. In practice, ultrasound imaging, which has been widely used in medical diagnostic technique due to its noninvasive nature, low cost, and realtime examination, can be applied to measure IMT by visual assessment of the leading edges. The carotid arteries are most suitable for study because of their superficial location, size, and limited movement, and the CCA is easy to image for it is relatively close and parallel to the skin surface. As illustrated in Figure 1, IMT is the distance between the two approximately parallel lines, lumenintima interface (LII), and mediaadventitia interface (MAI), respectively. Conventionally, IMT was measured by manual tracing of the interfaces between tissue layers. This method requires substantial experience, and it is time consuming and varies according to the training, experience, and the subjective judgment of the experts. Therefore, the manual measurements suffer from considerable inter and intraobserver variability [6, 7].
(a)
(b)
During the past two decades, efforts have been made by numerous investigators worldwide to try to find an approach for interface detection that is less reliant on human operators; thus many computeraided methods have been developed to solve this problem. The earliest research on the feasibility study of using Bmode sonographic images to measure the CCA IMT was reported in 1986 [8]. Since then, a wide variety of techniques have been proposed for the IMT measurement, which can be classified into two groups: semiautomatic [9–14] and fully automatic measurement [15–20].
Almost all of the existing methods include two stages: estimating the initial boundaries (LII, MAI) and then adjusting them to their exact locations. For the purpose of the first stage, the semiautomatic techniques need human intervention to select a region of interest (ROI) [11–14] or some points of interest [9, 10]. By contrast, the fully automatic technique can complete the task without any user interaction. Petroudi et al. [18, 19] proposed an algorithm that used variational level set [21] method to segment the ultrasound image into two parts according to image intensity and then estimated the MAI and LII with the boundary of the segmented image. Delsanto et al. [15, 16] developed a completely userindependent algorithm to locate the CCA tunicae as a ROI considering the mean value and variance of image intensity. Molinari et al. [17] described a novel technique called completely automated multiresolution edge snapper (CAMES), which recognized CCA based on a combination of scalespace and statistical classification in a multiresolution framework. The main difference between semi and fully automatic groups is whether the technique can automatically complete the task of the first stage. In order to achieve the results in the second stage, different methods have been developed for the adjustment, such as active contours [9, 10, 18, 19], WilliamsShah snake [11], dual snakes [14], dynamic programming [12], and dual dynamic programming [13]. A stateoftheart review on IMT measurement and wall segmentation is presented in [22, 23], and the latest review article of this field is presented in [24].
The aim of this paper is to present an automatic technique for CCA IMT measurement with low computational cost and high accuracy. More specifically, a ROI is extracted at the beginning of the processing, and then the LII and MAI boundaries are estimated within the extracted ROI; finally the estimated boundaries are adjusted to the accurate locations by using improved dynamic programming. The rest of this paper is structured as follows. In Section 2, we will describe the approach in detail. Section 3 shows the results of our experiments on 100 test images. Finally, Section 4 discusses the accuracy, efficiency, limitations, and advantage of the proposed approach.
2. Materials and Methods
2.1. Image Dataset
Our dataset consisted of 100 Bmode longitudinal ultrasound images of the CCA from Cyprus Institute of Neurology of Nicosia (Cyprus). These images are acquired by the ATL HDI3000 ultrasound scanner (Advanced Technology Laboratories, Seattle, USA), which is equipped with 64element, finepitch, highresolution, 38 mm broadband array, a multielement ultrasound scan head with an operating frequency range of 4–7 MHZ, an acoustic aperture of , and a transmission focal range of 0.8–1.1 cm. Digital images were resized using the bicubic method to standard pixel density of 16.66 pixels/mm. The images were logarithmically compressed and were recorded digitally on a magnetooptical drive at size of pixels with 256 gray levels. The images were recorded at Cyprus Institute of Neurology and Genetics, in Nicosia, Cyprus, from 42 female and 58 male symptomatic patients aged between 26 and 95 years, with a mean age of 54 years. These subjects were at risk of atherosclerosis and have already developed clinical symptoms, such as a stroke or a transient ischemic attack [11].
In order to remove the textual markers of the images and facilitate the following processing, the original ultrasound images were automatically cropped into a size of pixels (from 100 to 500 lines and from 200 to 600 columns).
2.2. Overview of the Proposed Approach
The flowchart of the proposed approach is shown in Figure 2. Our approach can be divided into six cascaded steps: (1) the method of templatebased matching is utilized for ROI extraction; (2) bilateral filtering is applied to remove noise and artifacts in the ultrasound images; (3) initial LII boundary is estimated within the extracted ROI; (4) initial MAI is estimated based on the edge map and the LII boundary; (5) the estimated boundaries are adjusted to their exact locations with improved dynamic programming; (6) IMT measurement and validation are carried out.
2.3. ROI Extraction
The method of templatebased matching was applied to extract the ROI. Firstly, a template containing part of the IMC and some lumen and adventitia was constructed according to the pixel density and image contrast. Then, the template was used to find the positions of the IMCs by searching a series of ultrasound images. We call the search image , where represents the coordinates of each pixel in the ultrasound image and the template with denotes the coordinates.
In our algorithm, we simply moved the center of the template over each point on a vertical line in the image and calculated the sum of the absolute differences (SAD) between their pixel intensities. Thus the SAD can be defined as where and denote the number of rows and columns of the template, which are assumed as odd integers. In this method the lowest SAD score gives the best estimation for the IMC on that vertical line. Actually, we find seven such points on seven evenly spaced vertical lines. Figure 3 shows the seven estimated points.
If the image is damaged by speckle noises or artifacts to a large extent, some points may be located at wrong positions. As long as four points or more are positioned correctly, we can achieve a ROI by excluding the wrong points. In practice, the difference in vertical position of the points with each other is computed. If the vertical position of one point has a difference within 2 mm with other three points or more, this point is assumed to be correctly positioned since the CCA IMC almost locates in a horizontal position in the longitudinal CCA ultrasound image.
Among the correctly located points, a rectangular ROI can be extracted with top left and bottom right points, which can be denoted by and , respectively, where is a positive constant integer.
2.4. Bilateral Filtering
Speckle noises always affect the boundary detection on ultrasound images. Therefore we need to reduce the effect of noise over the extracted ROI at the first place. The bilateral filter [25] has been chosen because it is a simple, noniterative tool for edge preserving as well as noise reducing.
The basic idea underlying bilateral filtering is to combine domain and range filtering, thereby enforcing both geometric and photometric locality. Bilateral filtering can be described as follows: With the normalization where measures the geometric closeness between the neighborhood center and a near point and measures the intensity similarity. Both the closeness function and the similarity function are Gaussian functions of the Euclidean distance between their arguments. The filtered ROI can be seen in Figure 4(a).
(a)
(b)
(c)
(d)
(e)
2.5. LII Estimation
After filtering, the intensity of the ROI is normalized into the scale of 0 to 1. When examining the histogram of the ROI intensity, we could find that the histogram consists of three peaks and two valleys, as shown in Figure 5. The two valleys are located at and , respectively. A pixel with intensity between and may belong to the CCA intima, whereas pixels with intensity higher than probably belong to the CCA adventitia. The value can serve as the threshold to divide the pixels in the ROI into two clusters according to the intensity of each point. Suppose that the extracted ROI has a size of lines and columns; we need to find a point in each column in order to get a rough LII boundary. Inspired by Liu et al. [12], the rough LII was searched as follows.(1)Let ; increase from 1 to until ; is the starting point on the rough boundary.(2)Increase to ; search the first point in the next column.(3)If , set as the next point on the rough boundary, if , set the next point at , and if , select as the next point and then set vertical position of the new point as .(4)If , output the rough boundary; otherwise, return to (2).
After the above operations, we can get an eightneighborhood continuous line standing for the rough boundary of LII and a binary image segmented by the threshold . In some cases, intima layer may disappear at some places because of artifacts; thus, there would be a gap in the binary image, as shown in Figure 4(b). Morphological operations consisting of dilation followed by erosion were implemented to fill the gap. Finally we can get a flat line standing for the estimated LII boundary, as in Figure 4(c).
2.6. MAI Estimation
The boundary of MAI is estimated according to the LII, since the two lines are approximately paralleled to each other. Additionally, the MAI lies on the border of media and adventitia where a strong edge locates.
In order to estimate MAI boundary using the initial LII, we proposed a simple method. Firstly, the edge map of the ROI was computed; secondly, the initial LII boundary was displaced downward for several pixels taking into account the normal IMT range (0.4–1.4 mm) and the pixel density of the images; thus we got a series of parallel lines; and then the value of the edge strength on each line was summed up; at last, the line with the highest value was chosen as the estimated MAI boundary. The estimated LII and MAI can be seen in Figure 4(d).
2.7. Boundaries Adjustment with Dynamic Programming
Dynamic programming [26] is a technique for solving optimization problems when not all variables in the evaluation function are interrelated simultaneously. To formulate the boundarysnapping procedure as dynamic programming, we must define an evaluation function that embodies a notion of the “best boundary.” In our situation, both LII and MAI boundaries are nearly straight lines at the position of strong edges. Therefore, the “best boundary” can be denoted by a weighted sum of high cumulative edge strength and low cumulative curvature. For an segment curve boundary, the following formula can be established: where is a negative constant, is the point on the estimated boundary and can move upward or downward for a pixel, and the function is the edge strength at position and the function is the curvature of the polyline defined by the 3 points: . , the solution of the largest , can be found by applying dynamic programming, which consists of the points of the “best boundary.”
Canny [27] showed that Gaussian derivatives yield good compromise between localization and detection. Therefore, we convolved the ROI with a vertical directional firstorder derivative of Gaussian kernel to get the edge map, which was previously used in Section 2.6. The Gaussian kernel we applied has a size of with standard deviation of 1 in both horizontal and vertical directions. The vertical directional derivative of the Gaussian kernel has the advantage for detecting edges lying approximately in the horizontal direction.
Practically, curvature has become a basic tool to formulate a contour function. Williams and Shah’s work [28] presented several curvature approximation methods for discrete curves and evaluated their efficiency of computation, accuracy of the estimation, and presence of anomalies. In a situation where a curve consists of three not evenly spaced points (Figure 6), their work proved that taking the difference between the normalized vectors could give a good estimation for the curvature; that is, the formula of the curvature is given by After the dynamic programming procedure, we could get the accurate boundaries both for LII and MAI as shown in Figure 4(e).
2.8. IMT Measurement and Validation
The IMC lies almost horizontally in the ultrasound image, and the detected boundaries of LII and MAI have the same number of points. Therefore, the IMT can be measured by calculating the mean absolute distance between the two detected boundaries: where is the number of points constituting the two boundaries and is the index spanning the columns of the image. The maximum and minimum values of can be used to test whether the detected IMT is reasonable. If they are not within the range of 0.2 mm to 1.6 mm, the IMT measurement is failed.
When processing a large number of ultrasound images, the automatic segmented (AS) IMT values and ground truth (GT) IMT values should be compared. There are two important criteria for the comparisons: correlation and BlandAltman plot.
The two sets of IMT values are correlated and the Pearson coefficient is used to give an estimate of the measurement agreement: where is the number of subjects for IMT measurement, and are the automatic segmented IMT values and ground truth IMT values for subject, and and are the mean value of the two sets of IMT.
The mean and standard deviation (SD) between AS and GT IMT are plotted in the BlandAltman plot. This method is very effective in pointing out possible biases in IMT measurement.
3. Results
In our experiment, the ground truth (GT) of the IMT was measured by expert vascular clinician, which was manually computed by delineating the LII and MAI boundaries on the ultrasound images for six times and then averaged. The automatic segmented (AS) IMT is measured by applying the proposed algorithm. Table 1 shows the parameters’ list in our experiment, where is the positive constant integer in Section 2.3; is the threshold in Section 2.5; is the size of the structuring element applied to morphological operations in Section 2.5; is the weights of curvature in Section 2.7.

100 ultrasound images from 100 patients were tested by our algorithm. Seven images were failed for IMT measurement with our approach because of the absence of the intima layer. Among the 93 correctly detected images, the mean IMT ± standard deviation for GT and AS is mm and mm; the mean difference ± standard deviation between AS and GT is mm. The correlation coefficient is 0.9629 for AS and GT IMT. The average processing time for automated measurement is 0.3223 s with MATLAB under Windows XP on an Intel Core 2 Duo CPU E7500 @2.93 GHz. Figure 7 shows the scatter diagram of the AS IMT with respect to GT IMT. Figure 8 shows the BlandAltman plots for AS IMT versus GT IMT. Figure 9 shows another six pairs of the experiment results; the left sides of the pictures are the extracted ROI and the right halves are the final boundaries. The first five examples illustrate the correctly detected IMT and the last one shows the case of failed detection.
(a)
(b)
(c)
(d)
(e)
(f)
4. Discussions
In our experiment, the mean difference of AS and GT is 0.0231 mm among the 93 correctly tested images; that is, the IMT computed by our technique is more likely thicker than that manually measured by expert. This may be due to the ambiguous boundary understanding of MAI, while human did not delineate the concave contours according to gradient maxima. Figures 9(b) and 9(c) show the concave contours automatically detected by our approach.
The IMC template used in Section 2.3 should contain an IMC strip with a medium IMT value (about 0.7 mm). It can be cropped from one ultrasound image of the dataset and then smoothed by bilateral filtering. For a set of ultrasound images with the same pixel density and similar contrast, the template is quite robust for detecting IMC regions correctly. Even if some points (no more than four out of seven) are positioned at wrong places because of speckle noises or artifacts, the negative effects can be ignored.
The threshold in Section 2.5 can be adaptive according to the contrast of ultrasound images. The dataset used in our experiment comes from the same machine with the same parameters setting. Therefore the 93 images were correctly detected with . Among the seven failed images, the intima layer has completely disappeared. The ROI of IMC can be correctly detected, but the final LII and MAI boundaries are almost overlaid (Figure 9(f)). The intensity histogram of the failed ROI was shown in Figure 10. Compared to Figure 5, peaks and valleys are ambiguous in this histogram image.
Our approach is compared with some representative researches in Table 2. The results show low mean IMT bias and short processing time with other researches. However, there are three limitations related to the proposed approach. Firstly, the IMC template should be changed when the algorithm applied to ultrasound images from other machines because of different pixel density or image contrast. Secondly, the results of LII estimation are sensitive to plaques. Thirdly, the mean bias of LII and MAI were not investigated in our experiment for the lack of information on LII and MAI boundaries of manual tracing.

The proposed algorithm in this paper is especially suitable for medical uses, where some ultrasound equipment has been set to automatic IMT measurement for different patients. There are only two parameters that need to be set for a new machine, that is, the IMClike template and the intensity threshold for LII detection. After that the machine can be used for automatic IMT measurement from one patient to another. Since the processing time is less than 0.5 s per image in MATLAB, it has the potential to achieve realtime IMT measurement when using C++ under Visual Studio.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 “Global status report on noncommunicable diseases,” Tech. Rep., World Health Organization, Geneva, Switzerland, 2011. View at: Google Scholar
 M. L. Bots, A. W. Hoes, P. J. Koudstaal, A. Hofman, and D. E. Grobbee, “Common carotid intimamedia thickness and risk of stroke and myocardial infarction: the Rotterdam Study,” Circulation, vol. 96, no. 5, pp. 1432–1437, 1997. View at: Publisher Site  Google Scholar
 J. K. Virtanen, S. Voutilainen, P. Happonen et al., “Serum homocysteine, folate and risk of stroke: kuopio Ischaemic Heart Disease Risk Factor (KIHD) Study,” European Journal of Cardiovascular Prevention and Rehabilitation, vol. 12, no. 4, pp. 369–375, 2005. View at: Publisher Site  Google Scholar
 M. W. Lorenz, H. S. Markus, M. L. Bots, M. Rosvall, and M. Sitzer, “Prediction of clinical cardiovascular events with carotid intimamedia thickness: a systematic review and metaanalysis,” Circulation, vol. 115, no. 4, pp. 459–467, 2007. View at: Publisher Site  Google Scholar
 C. M. Robertson, F. Gerry, R. Fowkes, and J. F. Price, “Carotid intimamedia thickness and the prediction of vascular events,” Vascular Medicine, vol. 17, no. 4, pp. 239–248, 2012. View at: Publisher Site  Google Scholar
 R. Salonen, A. Haapanen, and J. T. Salonen, “Measurement of intimamedia thickness of common carotid arteries with highresolution bmode ultrasonography—interobserver and intraobserver variability,” Ultrasound in Medicine and Biology, vol. 17, no. 3, pp. 225–230, 1991. View at: Google Scholar
 D. Frost, A. Friedl, and W. Beischer, “Evaluation of the carotid artery intimamediathickness: influences of different methods, characteristics of the subjects, and investigators,” Ultraschall in der Medizin, vol. 19, no. 4, pp. 168–173, 1998. View at: Publisher Site  Google Scholar
 P. Pignoli, E. Tremoli, and A. Poli, “Intimal plus medial thickness of the arterial wall: a direct measurement with ultrasound imaging,” Circulation, vol. 74, no. 6, pp. 1399–1406, 1986. View at: Publisher Site  Google Scholar
 D. C. Cheng, A. SchmidtTrucksäss, K. Cheng, and H. Burkhardt, “Using snakes to detect the intimal and adventitial layers of the common carotid artery wall in sonographic images,” Computer Methods and Programs in Biomedicine, vol. 67, no. 1, pp. 27–37, 2002. View at: Publisher Site  Google Scholar
 J. L. IzquierdoZaragoza, L. C. BastidaJumilla, M. C. VerduMonedero et al., “Segmentation of the carotid artery in ultrasound images using frequencydesignen Bspline active contour,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '01), pp. 35–49, Prague Congress Ctr, Prague, Czech Republic, 2001. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 G. Liu, B. Wang, and D. C. Liu, “Detection of intimamedia layer of common carotid artery with dynamic programming based active contour model,” in Proceedings of the Chinese Conference on Pattern Recognition (CCPR '08), pp. 369–374, October 2008. View at: Publisher Site  Google Scholar
 D. 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. View at: Publisher Site  Google Scholar
 X. Xu, Y. Zhou, X. Cheng, E. Song, and G. Li, “Ultrasound intimamedia segmentation using Hough transform and dual snake model,” Computerized Medical Imaging and Graphics, vol. 36, no. 3, pp. 248–258, 2012. View at: Publisher Site  Google Scholar
 S. Delsanto, F. Molinari, P. Giustetto, W. Liboni, and S. Badalamenti, “CULEXcompletely userindependent layers EXtraction: ultrasonic carotid artery images segmentation,” in Proceeding of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEEEMBS '06), pp. 6468–6471, Shanghai, China, Jan 2006. View at: Publisher Site  Google Scholar
 S. Delsanto, F. Molinari, P. Giustetto, W. Liboni, S. Badalamenti, and J. S. Suri, “Characterization of a completely userindependent algorithm for carotid artery segmentation in 2D ultrasound images,” IEEE Transactions on Instrumentation and Measurement, vol. 56, no. 4, pp. 1265–1274, 2007. View at: Publisher Site  Google Scholar
 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 multiinstitutional database,” IEEE Transactions on Image Processing, vol. 21, no. 3, pp. 1211–1222, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Petroudi, C. Loizou, M. Pantziaris, M. Pattichis, and C. Pattichis, “A fully automated method using active contours for the evaluation of the intimamedia 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. View at: Publisher Site  Google Scholar
 S. Petroudi, C. Loizou, M. Pantziaris, and C. Pattichis, “Segmentation of the common carotid intimamedia complex in ultrasound images using active contours,” IEEE Transactions on Biomedical Engineering, vol. 59, no. 11, pp. 3060–3069, 2012. View at: Publisher Site  Google Scholar
 F. Molinari, U. Rajendra Acharya, G. Zeng, K. M. Meiburger, and J. S. Suri, “Completely automated robust edge snapper for carotid ultrasound IMT measurement on a multiinstitutional database of 300 images,” Medical and Biological Engineering and Computing, vol. 49, no. 8, pp. 935–945, 2011. View at: Publisher Site  Google Scholar
 R.M. MenchonLara, B.J. MaríaConsuelo, M.S. Juan, and S.G. JoséLuis, “Automatic detection of the intimamedia thickness in ultrasound images of the common carotid artery using neural networks,” Medical & Biological Engineering & Computing, vol. 52, no. 2, pp. 169–181, 2014. View at: Publisher Site  Google Scholar
 F. Molinari, G. Zeng, and J. S. Suri, “A state of the art review on intimamedia thickness (IMT) measurement and wall segmentation techniques for carotid ultrasound,” Computer Methods and Programs in Biomedicine, vol. 100, no. 3, pp. 201–221, 2010. View at: Publisher Site  Google Scholar
 F. Molinari and J. S. Suri, “Automated measurement of carotid artery intimamedia thickness,” in Ultrasound and Carotid Bifurcation Atherosclerosis, A. Nicolaides and K. W. Beach, Eds., pp. 177–192, Springer, London, UK, 2012. View at: Google Scholar
 V. Naik, R. S. Gamad, and P. P. Bansod, “Carotid artery segmentation in ultrasound images and measurement of intimamedia thickness,” BioMed Research International, vol. 2013, Article ID 801962, 10 pages, 2013. View at: Publisher Site  Google Scholar
 C. Tomasi and R. Manduchi, “Bilateral filtering for gray and color images,” in Proceedings of the 6th IEEE International Conference on Computer Vision, pp. 839–846, New Delhi, India, January 1998. View at: Google Scholar
 R. E. Bellman and S. E. Dreyfus, Applied Dynamic Programming, Princeton University Press, Princeton, NJ, USA, 1962. View at: MathSciNet
 J. Canny, “A computational approach to edge detection,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 8, no. 6, pp. 679–698, 1986. View at: Google Scholar
 D. J. Williams and M. Shah, “A fast algorithm for active contours and curvature estimation,” CVGIP: Image Understanding, vol. 55, no. 1, pp. 14–26, 1992. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Qiang Li 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.