- Open Access
Curvelet based automatic segmentation of supraspinatus tendon from ultrasound image: a focused assistive diagnostic method
BioMedical Engineering OnLine volume 13, Article number: 157 (2014)
Disorders of rotator cuff tendons results in acute pain limiting the normal range of motion for shoulder. Of all the tendons in rotator cuff, supraspinatus (SSP) tendon is affected first of any pathological changes. Diagnosis of SSP tendon using ultrasound is considered to be operator dependent with its accuracy being related to operator’s level of experience.
The automatic segmentation of SSP tendon ultrasound image was performed to provide focused and more accurate diagnosis. The image processing techniques were employed for automatic segmentation of SSP tendon. The image processing techniques combines curvelet transform and mathematical concepts of logical and morphological operators along with area filtering. The segmentation assessment was performed using true positives rate, false positives rate and also accuracy of segmentation. The specificity and sensitivity of the algorithm was tested for diagnosis of partial thickness tears (PTTs) and full thickness tears (FTTs). The ultrasound images of SSP tendon were taken from medical center with the help of experienced radiologists. The algorithm was tested on 116 images taken from 51 different patients.
The accuracy of segmentation of SSP tendon was calculated to be 95.61% in accordance with the segmentation performed by radiologists, with true positives rate of 91.37% and false positives rate of 8.62%. The specificity and sensitivity was found to be 93.6%, 94% and 95%, 95.6% for partial thickness tears and full thickness tears respectively. The proposed methodology was successfully tested over a database of more than 116 US images, for which radiologist assessment and validation was performed.
The segmentation of SSP tendon from ultrasound images helps in focused, accurate and more reliable diagnosis which has been verified with the help of two experienced radiologists. The specificity and sensitivity for accurate detection of partial and full thickness tears has been considerably increased after segmentation when compared with existing literature.
Ultrasound is a medical imaging modality preferred by radiologist and physicians due to its capability to produce real time imaging data, reduced time in diagnosis, and high patient acceptability. Musculoskeletal (MSK) ultrasound is rapidly gaining importance for the assessment of joints and soft tissue disorders [1, 2]. Occurrence of MSK disorders are mostly found in sports, adulthood (>30 years of age) or in accidents due to sudden impact which leads to reduced functions of daily life.
Supraspinatus (SSP) tendon is among the four muscles found in rotator cuff in shoulder; it runs from supraspinatus fossa superior of scapula to the greater tuberosity in humerus. Supraspinatus (SSP) tendon disorders are third most prevalent in MSK [3, 4]. The disorders in SSP tendon occur in the form of tear (loss of connective tissues binding collagen fibre), tendinosis (inflammation) and results in acute pain, insomnia and, reduced mobility of shoulder. SSP tendon as imaged using ultrasound equipment is shown in Figure 1. In ultrasound images, soft tissues such as muscle, fat, and other connective tissues reflect different echogenic pattern allowing the radiologist to differentiate between healthy and ailing tissue. Healthy muscle tissues have uniform and organized patterns which tend to absorb the ultrasound beam and appear hypoechoic compared to fat [5, 6]. However, tendons with pathological conditions are disorganized, diffused and have hypoechoic appearance compared to healthy tendons. Conventionally, health of tendons was visually examined based on image texture provided by ultrasound images which makes modality highly operator dependent. Recently, ultrasound applications are being trained to provide quantitative information about diagnostics of patient’s health and fitness [7–10]. To date, for supraspinatus tendon quantitative analysis is limited to cross sectional area and thickness calculation . In one study , quantitative ultrasound techniques were employed wherein it computes structural measurements and mean echogenicity to discriminate between muscle pathologies and healthy conditions. In another recently published research, musculoskeletal ultrasound was used to quantify muscle kinematics during dynamic activities such as drop landing in healthy subjects . Reliability study was conducted for the use of ultrasound as a quantitative information tool by a research group for inter-rater and intra-rater reliability, the study reveals that appropriately designed protocol will allow radiologist to identify structural changes within tendons [13, 14].
Despite growing interest and established techniques, few attempts are made to automate quantitative analysis of pathologies existing in tendon. In [15, 16], Horng et al. used feature extraction method based on gray level co-occurrence matrix (GLCM) and texture feature coding method (TFCM) to classify lesions in tendon. In , Nielsen et al. proposed method for tissue characterization using first order and higher order statistics for SSP tendon and thigh muscles. In , Michael et al. suggested a measurement method to estimate muscle dystrophy and sarcopenia via echogenicity measurement using grayscale analysis for musculoskeletal muscles. This work is focused towards automatic segmentation of SSP tendon from ultrasound image for focused analysis and serves as first step towards automated diagnosis of tendon.
Segmentation in ultrasound images can be difficult for numerous reasons such as contrast and resolution of image, speckle noise which is inherent property of ultrasound imaging modality, operator dependency of the modality. The paper addresses the issue of contrast enhancement, despeckling and issues occurring due to operator dependency for accurate segmentation of SSP tendon.
Delineating boundaries and region of interest segmentation is a challenging task due to low contrast, resolution and inheritance of high intrinsic noise (speckle) in ultrasound images. Researchers [19–40] have suggested application specific segmentation technique for segmentation of ultrasound images. In , a brief review about ultrasound image segmentation techniques has been covered. The feature based segmentation techniques are often used in segmentation of ultrasound images. Image features such as gray level distribution , phase information , gradient information , shape of anatomical region  and temporal information have been exploited for accurate segmentation. In [24, 25], Rayleigh model has been extensively used in numerous occasions with anisotropic diffusion for removal of speckle noise and segmentation of ultrasound images. Rayleigh mixture model  was used for segmentation of skin lesions from ultrasound image. Rodtook  and Wang  have incorporated gradient information using active contour method and level set method for extraction of region of interest from ultrasound images. Intensity gradients are used along with anisotropic diffusion , morphological filtering , total variation filtering  and other filtering techniques [31, 32] for ultrasound image segmentation. Using shape priors is also very popular technique for segmentation of ultrasound images and has been used for segmentation of kidney , intravascular ultrasound (IVUS) , cardiac , and prostate ultrasound images . Curvelet decomposition [36, 37] has recently been introduced in medical imaging for segmentation and enhancement of images. In , curvelet features are used for analysis of retinal images along with connected component analysis. Semi-supervised ultrasound image segmentation  technique is proposed using curvelet features and texture analysis. In , curvelet features are used in ultrasound images for segmentation of kidney images for better and enhanced diagnosis.
This is the pioneer work towards focused and automated segmentation of SSP tendon for accurate, focused and more reliable diagnosis. In this work, automatic segmentation of SSP tendon from ultrasound image is proposed. The method uses curvelet transform for feature extraction based on energy analysis of features followed by connected component analysis and morphological operations to accomplish the task.
The database of 116 images of SSP tendon ultrasound image with and without pathological condition was collected with the help of trained radiologist from Universiti Malaya Medical Centre. The data was collected over a period of three years with patients varying in age group from 19–54 years of age. The database was collected using Philips iU22 ultrasound and GE Logiqe BT2011 ultrasound machines. The database consists of images of patients with abnormalities such as single or multiples tears and tendinosis. 80 patients were observed to be suffering from only tendinosis, 100 patients were suffering from single or multiple tears in tendinosis or healthy tendon. The implementation of proposed methodology was done on MATLAB version 8.0.
Ultrasound image of the SSP tendon when observed using database of collected images acquired by two different radiologists on two different machines with around 51 patients found to contain variations such as:
The width of the tendon,
Topography or location of tendon in images,
Intensity level of SSP tendon in image,
Radius of curvature of convex SSP tendon.
The observed differences are due to inter and intra-operator variability, acquisition system and patients anatomy. Therefore, for a method to be able to segment automatically the SSP tendon, it should be robust towards above variations. Despite above variations, invariants or unique features in ultrasound image of SSP tendon irrespective of machines, patients and radiologist are:
Convex nature of SSP tendon.
Compressed between bursae at the top and cortical bone at the bottom.
High intensity of bursae and cortical bone compared to nearby structures.Convex nature of SSP tendon which is evident from Figure 1, is found in all the images of the database. The location of SSP tendon runs from supraspinatus fossa superior of scapula to insert in greater tuberosity in humerus which makes it suppressed between bursae and cortex in all the images. Bursae and Cortex due to its high echogenic property exhibits hyperechoic structure compared to neighboring areas. These three unique features for SSP tendon were exploited to propose a robust method invariant to above variations for segmentation of SSP tendon. The process flow for proposed methodology for automatic segmentation of SSP tendon is shown in Figure 2.
Consequently, the proposed algorithm for automatic segmentation of SSP tendon follows the steps:
Step 1: Select ultrasound image of SSP tendon from database.
Step 2: Apply image enhancement using Rayleigh adaptive contrast enhancement followed by despeckling using anisotropic diffusion method.
Step 3: Decomposition of enhanced image by real valued curvelet coefficients using wrapping function.
Step 4: Morphological, Logical and Area filtering operation to remove the regions that do not belong to SSP tendon or are outliers.
Step 5: Polynomial curve fitting to smooth tendon and recover lost boundary points and construct the mask as per radiologist requirement.
The detailed methodology for proposed work is discussed in subsequent section.
Image enhancement and feature extraction
The ultrasound image suffers from poor resolution and contrast which causes traditional segmentation algorithm to poorly segment region of interest from the given image. In proposed method, ultrasound image obtained is first contrast enhanced. The contrast enhancement is performed using Rayleigh distribution following the concept that speckle in ultrasound image follows Rayleigh pattern . The probability and cumulative density function of Rayleigh distribution is exploited to map the contrast enhanced pixel to image grid. The image pixels are mapped to the contrast enhanced image using the following equations.
Mapping for pixels in contrast enhanced image is performed using cumulative density function (CDF) of Rayleigh distribution which is described below
The probability density function (PDF) for Rayleigh distribution is given by,
The pixel wise mapping was done using,
nj gives the number of pixels in image with jth gray level and x is the gray value. The value of σ for enhancement was empirically selected to be 0.4.
Contrast enhancement makes the boundaries as well as speckle in ultrasound image more prominent. Therefore removal of speckle noise from the homogeneous region so as to keep the boundaries prominent is performed using anisotropic diffusion method . Number of iterations of anisotropic diffusion performed for smoothing image is 10. The resultant image from contrast enhancement and despeckling is shown in Figure 3(b) and Figure 3(c) respectively. The preprocessing performed enhances the features needed for automatic extraction of tendon. The feature extraction from preprocessed image is performed using curvelet transform and is discussed in subsequent section.
Directional edge feature extraction
The capability of curvelet transform to extract directional edge features at different orientations is exploited for extracting edge features from ultrasound image. Curvelet transform provides details such as spectral information successfully at different orientations with reduced complexity. Curvelet decomposition using wrapping function method and real coefficient is used because of reduced complexity and faster computation. The feature extraction from ultrasound image is performed using method illustrated in Figure 4. Detailed description about curvelet transform can be found in literature .
Curvelet transform decomposes the image in different direction and frequencies in the curvelet domain. Curvelet coefficients after decomposition are well adapted to contain edges from the image. SSP tendon ultrasound image consist of bright features, which comprise of Bursae and Cortex along with some outliers in the form of muscle fat mixture or abnormalities such as calcification. The brightest structures along with some outliers are likely to present in biggest curvelet coefficients after decomposition. To find best scale and percentage of coefficients, energy analysis is conducted at scales 2, 3, 4 and 5 with orientation 16 and 32. The percentage of curvelet coefficients were varied from 5 to 10. The method used for calculating the optimum number of features was based on the energy of the region delineated by radiologist. The energy of the region delineated by radiologist is correlated with the energy of the reconstructed image using curvelet coefficients at different scales and orientation.The manual delineation of SSP tendon was performed by radiologist and is shown in color outline in Figure 5. The segmentation region includes bursae and cortex. The biggest curvelet coefficients contains region with highest intensity values which in this case happens to be: bursae, cortex and outliers. The idea is to filter out coefficients of high intensity values and remove outliers generating a mask to segment SSP tendon. The energy of manually segmented bursae and cortex by radiologist was calculated using the formula
where, p(i, j) is the pixel intensity at (i, j) and m* n are the total number of pixels in an image. The fraction of the energy contained in the area segmented by radiologist is computed. The ratio is calculated for the amount of energy contained in manual segmented region to that of total energy of the image.
The amount of energy contained in region extracted by radiologist was found to be consistent at 40 ± 5%, with almost all the images of the database. The best scale and orientation from curvelet decomposition was selected based on the energy analysis performed over scales: 2, 3, 4 and 5 with different orientation. The results of the analysis are shown in Table 1.
Table 1, shows the ratio of energy when 5%-10% of the curvelet coefficients are used to reconstruct ultrasound image. The ratio is calculated between the image reconstructed using varying percentages from 5%-10% to that of total energy of the image. It suggests that scales 3, 4, 5 at different orientation provides energy ranging from 95%-99.9% with curvelet coefficients percentage from 5%-10% respectively.
From the analysis performed above, on segmented image by radiologist, it has been found that bursae and cortex together contains 40 ± 5% of energy. Therefore, the decomposition of curvelet at higher scales was avoided. The decomposition at scale 2 with orientation 16 and 32 was performed, as evident from Table 1, and it is found that the amount of energy contained, when two different orientations are chosen is almost same. Therefore, in order to save computational time scale 2 with 16 orientations were chosen. The energy level increases with increasing percentage of coefficient from 5-10%. The image reconstructed with percentage coefficient 6.5%-7% was found to contain the region manually segmented by radiologist along with outliers in the form of muscle fat-mixture.The reconstructed images with 6.5% curvelet coefficient are shown in Figure 6 and Figure 3(d). The image contains several outliers that are unavoidable along with bursae and cortex. If lower percentage of coefficients is chosen, then, the structure for bursae and cortex starts to deteriorate. The reconstructed image with 6.5% curvelet coefficients contains 64.06% of the energy, wherein the bursae and cortex along with approximately 15%-20% of outliers is found. The energy test was conducted on several images from dataset and the results were found to be consistent. The second level with 16 orientations was found to be giving the best approximation for energy and was closest to radiologist’s segmentation. The increase in the energy of the reconstructed image is consistent with the outliers that are generated in the process to automate the segmentation. The next step of the algorithm focuses on removing the outliers. In the subsequent steps, the removal of outliers is performed using mathematical concepts of morphological operations and area filtering.
Area filtering and connected component analysis
The outliers in the reconstructed image remain in the form of bright structures which could be either muscle fat mixture or some form of an abnormality, i.e. calcification. The reconstructed image is rescaled to gray levels (0 ~ 255) to perform connected component analysis and thresholding. The removal of intensity components with lower intensity levels is performed by applying threshold based on formula introduced in  that is,
where, μ is the mean and σ is the standard deviation respectively. The value of α (α < 1) is chosen small enough. The thresholding is performed based on
where, p(i, j) is the value of pixel intensity at location (i, j) in an image. To remove undesired pixels or blobs from image area filtering is applied. To perform this, connected component analysis (CCA) with eight neighborhood is performed where the connected pixel component that are identified below certain threshold are removed. The threshold value is specified empirically, which involves analyzing the number of pixels in bursae or tendon for reconstructed image. From the analysis of more than 100 images with the size (421x580), it has been found that blobs with total number of pixels values less than (~ < 5000) tend to be the outlier. The thresholding value is relative to the size of the image. The image formed after the thresholding is shown in Figure 3(e).
Although most of the undesired objects are removed in this step, but some of the components connected with bursae and tendon are still inevitable. Therefore, desired mask for segmentation of SSP tendon is obtained by removal of these remaining outliers.
The dilation, erosion, opening, closing are standard mathematical morphological operations that helps in the removal of distorted texture and noise from image using structure element. For this application, disk was used as structure element. Morphological operations leave the structures bigger than SE unchanged. The drawback of conventional dilation and erosion is that they do not preserve edge information perfectly. The new operator proposed by Bangham et al. , takes care for above drawback and emphasizes on size of structure but forgets the shape completely. Morphological operators that take care of above problem and consider both shape and size were proposed . If the image is I(x, y) and g is the mask image, then general equation for geodesic dilation is defined by equation
And, the equation for geodesic erosion is written as
Step by step execution and pointwise application of maximum (minimum) operator after each iteration of elementary geodesic erosion (dilation) is necessary to control the marker image. In general, after execution of some steps both the above equations becomes stable and no further changes occur that means
based on above definitions, opening and closing operation by reconstruction is performed as below:
Application of conventional morphological operations sometimes introduces new edges and contours which are not needed. Reconstruction by geodesic morphological operations yields results wherein these drawbacks are removed. Sample test of geodesic morphological operations is shown in Figure 7. The morphological operation shown uses disk as the structure element.The morphological operation is performed to remove the remaining objects. The method involves two steps: one is morphological dilation followed by erosion operation. The structure element used for morphological operations is a disc with the size 17x6. The result obtained from morphological operations is shown in Figure 3(f). In the final step, the area filtering is applied wherein the high threshold value is chosen so as to obtain the mask for SSP tendon. The sequence of steps for image segmentation is shown in Figure 3.
Therefore, the algorithm is scripted as follows:
A thresholding operation is performed to remove high intensity values. The threshold value is chosen based on Eq. (27).
Connected component analysis (CCA) with 8-neighborhood is performed for removal of blobs not belonging to SSP tendon which are considered as outliers.
Area filtering is performed to remove the blobs with small number of pixels which do not belong to tendon.
Morphological operation of opening is performed using structure element of 17×6.
In last step, again area fileting is used to obtain the clear structure of tendon.The results at this stage, shown in Figure 3(g), were presented to radiologist for assessment. After discussion, recommendation was put forward to extract smooth boundary for SSP tendon which enclose bursae and cortex.
Polynomial curve fitting and morphology
The boundary for extracted mask was found and polynomial curve fitting is applied to make the boundary of the mask smooth for better visualization of the tendon. Two different quadratic functions were estimated based on the boundary point of the extracted region. The curve fitting was performed by using quadratic polynomial function. The result for this step is shown in Figure 3(h).
The curve fitting also recovered some of the missing boundary points from images thereby increasing the accuracy of the results. The experimental results of proposed algorithm and radiologist assessment of segmentation is discussed in next section.
The proposed algorithm was assessed based on the quantitative and qualitative analysis. The quantitative analysis for the assessment of proposed methodology for segmentation was assessed using three metrics: 1) false positive rate (FPR); 2) true positive rate (TPR) ; 3) Accuracy (ACC). The radiologist were requested to delineate SSP tendon manually ultrasound image. The comparison of the area delineated by radiologist and automatic segmentation performed by proposed algorithm was done based on the area of true positive, false positive, true negative, false negative. The Pictorial representation for corresponding area is shown in Figure 8.
The formula used to calculate FPR, TPR and Accuracy is,
where , Am is the manual segmentation for area of SSP tendon performed by an experienced radiologist. Aa is the automatic segmentation of SSP tendon performed by proposed methodology.
In the second round of assessment, the segmented images were provided to radiologist and blind evaluation was performed for partial thickness tear (PTTs) and full thickness tear (FTTs). And the results obtained by radiologist assessment were modelled so as to give sensitivity and specificity value for the method. The true positives, false positives, false negatives and true negatives in this were defined as
True Positives- pathology existent and diagnosed correctly
False positives- pathology not present and diagnosed
False negatives- pathology present but not diagnosed
True Negatives- pathology not present and not diagnosed
The sensitivity and specificity were calculated using the formula below
The results for assessment of pathological condition were computed and compared with the diagnosis of PTTs and FTTs using ultrasound images.
The result of tendon segmentation as per radiologist assessment was performed with the set operations described in proposed methodology section. The results of each step are shown in Figure 9.Figure 10, shows qualitative comparison of the results using manual segmentation by radiologist and result obtained using proposed algorithm. Figure 10(a) shows the true positive results wherein accurate segmentation as per radiologist requirement was attained. Figure 10(b) it can be seen that the region segmented by radiologist is present in the result with an outlier region shown. The false positives arise because of the inaccurate detection of muscle fat outlier which plays important role in segmentation. In Figure 10(c), result shows inaccuracy of the algorithm in the form of detection of false negatives, wherein the reason for inaccuracy is poor visibility of bursae. No cases of true negatives were found (true negatives are cases when the segmentation area completely lies outside the region manually segmented by radiologists). The quantitative assessment of results obtained using proposed algorithm was performed using the above mentioned parameters.
Table 2 shows the accuracy of the results in the form of true positive rate, false positive rate and accuracy of the result.
The performance of proposed method was also tested for focused assessment with the help of two radiologist and result were compared with the existing literature [47–49]. The results were compared for the diagnosis of partial thickness and full thickness tear using sensitivity and specificity values.
The results of comparison are shown in Table 3. In , authors used 40 images to compute the sensitivity and specificity of the diagnosis of PTTs and FTTs in SSP tendon, whereas Singh et al. in  uses 36 images to quantify the diagnosis of PTTs and FTTs in SSP tendon and Rutten et al. in  uses a database of 68 images to quantify the diagnosis of PTTs and FTTs in SSP tendon using ultrasound images. In this work, authors have used a set of 116 ultrasound images to qualitatively evaluate the performance of proposed method for assessment of PTTs and FTTs in ultrasound images. The sensitivity and specificity values for PTTS computed using [47, 48], and  are 92.3% and 92.6%, 66.7% and 93.5%, 92% and 33.5% whereas using the proposed method the values are 94% and 93.6%. Similarly, for FTTs sensitivity and specificity for [47, 48], and  are 92.6% and 94%, 92.3% and 94.4%, 94% and 94% resp. whereas sensitivity and specificity for diagnosis of FTTs using segmented SSP tendon are 95.6% and 95% respectively.
Algorithm run-time for automatic segmentation of SSP tendon is 1.4 s. The algorithm was implemented using i7, 3.4 GHz processor and 8GB RAM. The enhancement of image for segmentation takes 0.55-0.65 s and segmentation part takes about 0.75-0.85 s. Therefore, the proposed approach is fast, reliable with scope of possible future application in real time diagnosis. The assessment parameter and computation time for the algorithm suggests high performance for proposed methodology.
Transportability of method
Ultrasound is a deterministic imaging modality which means images acquired under identical circumstances will yield similar result. Despite being deterministic in nature ultrasound image is highly operator dependent. In the proposed method, the attempt has been made to decrease the operator dependency of ultrasound machine for diagnosis of pathological conditions in SSP tendon. The image features are studied from images taken from different operator and machines and analyzed for unique features. The method uses image processing methods to extract features from the tendon and reconstruct SSP tendon. In all 116 images taken from different source and two different radiologists it was found that SSP tendon is convex in nature and is always compressed between bursae at the top and humeral cortex at the bottom. Both bursae and cortex have hyperechoic texture as compared with tendon. Since proposed method is based on the redundant nature of occurrence of SSP tendon in ultrasound image and the algorithm is tested for variability in equipment and operators visualization of SSP tendon. Therefore, proposed method can be used in clinical settings for post processing of SSP tendon ultrasound image for pathologies in SSP tendon with reduced inter and intra observer variability. The method can also be effectively used for training of medical officers for focused and effective diagnosis of SSP tendon pathology. The future work is intended to reduce the complexity of algorithm so that it can be effectively used as real time imaging tool embedded in ultrasound equipment for clinical analysis of SSP tendon.
Discussion and future work
Before starting the research, literature survey on the existing methods for diagnosis of SSP tendon was studied. Also, several rounds of discussions with radiologist were performed regarding the state of art methods for diagnosis of SSP tendon. The problems faced by patients and radiologist during examination were consulted. It was found that MSK ultrasound techniques provide subjective evaluation regarding existing pathologies in tendon and are also painful because of duration of examination. The initiative was taken to automate the pathologies in SSP tendon. To further the research and automate pathologies first challenge was to locate region of interest i.e. SSP tendon automatically so that focused and accurate diagnosis for the ailment can be performed.
In this paper, a novel method for automatic segmentation of SSP tendon from ultrasound image is proposed. The method involves image enhancement and feature extraction from ultrasound image. The image was contrast enhanced using statistically adaptive method followed by speckle removal using anisotropic diffusion method. The image was then decomposed using curvelet transform. The energy analysis of decomposition was performed to select the amount of curvelet features needed for mask generation. It was found, that 6.5% of curvelet features, at scale 2 and 16 orientations, provides best mask for segmentation. Images were reconstructed using extracted curvelet features and geodesic morphological operations were used to extract edges and remove outliers. Connected component analysis and area filtering were applied to remove the remaining false areas and perform accurate detection. There is a trade-off between selecting curvelet features and removal of false areas. High percentage of curvelet features results in increase of false positives. The polynomial curve fitting is used to smooth the area of SSP tendon as per radiologist’s recommendations. The segmented SSP tendon will assist the radiologist for focused and accurate diagnosis of abnormalities in the tendon. The quantitative assessment performed for segmentation and results of diagnosis for pathological conditions suggests the effectiveness of proposed algorithm. Also the computation time for algorithm shows the capability of the algorithm to be made available for real time diagnosis of pathologies in SSP tendon.
In future, the work will be focused to provide an automated system for pathology in SSP tendon. The computation time of the algorithm will be reduced by refining extracted coefficients so that possible implementation in real time diagnosis is possible.
The automatic segmentation of SSP tendon was successfully achieved and radiologist assessment for segmentation was performed. As per radiologist comments, the results help in enhancing the accuracy of diagnosed pathology because of focused assessment of tendon. The accuracy for the assessment of segmentation of SSP tendon is 95.61%. When diagnosing tendinosis, tear or calcification in SSP tendon segmented tendon provide focused and more reliable result with increased sensitivity and specificity. The proposed algorithm well suited for real time applications for musculoskeletal ultrasound in SSP tendon.
Eranki A, Cortes N, Ferencek ZG, Sikdar S: A novel application of musculoskeletal ultrasound imaging. J Vis Exp 2013, 79: e50595.
Zanetti M, Hodler J: Imaging of degenerative and posttraumatic disease in the shoulder joint with ultrasound. Eur J Radiol 2000,35(2):119–125. 10.1016/S0720-048X(00)00227-8
National Center for Health Statistics: National Health Interview Survey. Hyattsville, MD: US Department of Health and Human Services; 1995.
Thain LM, Adler RS: Sonography of the rotator cuff and biceps tendon: technique, normal anatomy, and pathology. J Clin Ultrasound 1999,27(8):446–458. 10.1002/(SICI)1097-0096(199910)27:8<446::AID-JCU6>3.0.CO;2-M
Heckmatt JZ, Pier N, Dubowitz V: Real-time ultrasound imaging of muscles. Muscle Nerve 1988,11(1):56–65. 10.1002/mus.880110110
Nielsen PK, Jensen BR, Darvann T, Jørgensen K, Bakke M: Quantitative ultrasound image analysis of the supraspinatus muscle. Clin Biomech 2000,15(Supplement 1(0)):S13-S16.
Shoma A, Moutamed A, Ameen M, Abdelwahab A: Ultrasound for accurate measurement of invasive breast cancer tumor size. Breast J 2006,12(3):252–256. 10.1111/j.1075-122X.2006.00249.x
Baroncelli GI: Quantitative ultrasound methods to assess bone mineral status in children: technical characteristics, performance, and clinical application. Pediatr Res 2008,63(3):220–228. 10.1203/PDR.0b013e318163a286
Binkley TL, Berry R, Specker BL: Methods for measurement of pediatric bone. Rev Endocr Metab Disord 2008,9(2):95–106. 10.1007/s11154-008-9073-5
ᅟ ACUSON SC2000 Volume Imaging Ultrasound System. Siemens: Automated Volumetric Cardiac Ultrasound Analysis (White Paper) 2008. [https://library.peirce.edu/html/citations/whitepaper]
Leung JL, Griffith JF: Sonography of chronic Achilles tendinopathy: a case–control study. J Clin Ultrasound 2008,36(1):27–32. 10.1002/jcu.20388
Scholten RR, Pillen S, Verrips A, Zwarts MJ: Quantitative ultrasonography of skeletal muscles in children: normal values. Muscle Nerve 2003,27(6):693–698. 10.1002/mus.10384
Collinger JL, Gagnon D, Jacobson J, Impink BG, Boninger ML: Reliability of quantitative ultrasound measures of the biceps and supraspinatus tendons. Acad Radiol 2009,16(11):1424–1432. 10.1016/j.acra.2009.05.001
Collinger JL, Fullerton B, Impink BG, Koontz AM, Boninger ML: Validation of grayscale-based quantitative ultrasound in manual wheelchair users: relationship to established clinical measures of shoulder pathology. Am J Phys Med Rehabil 2010,89(5):390–400. 10.1097/PHM.0b013e3181d8a238
Horng M-H: The glowworm swarm optimization for training the radial basis function network in ultrasonic supraspinatus image classification. Adv Sci Lett 2013,19(9):2724–2727. 10.1166/asl.2013.4977
Horng M-H, Chen S-M: Multi-class classification of ultrasonic supraspinatus images based on radial basis function neural network. J Med Biol Eng 2009,29(5):242–250.
Nielsen PK, Jensen BR, Darvann T, Jorgensen K, Bakke M: Quantitative ultrasound tissue characterization in shoulder and thigh muscles–a new approach. BMC Musculoskelet Disord 2006, 7: 2. 10.1186/1471-2474-7-2
Harris-Love MO, Monfaredi R, Ismail C, Blackman MR, Cleary K: Quantitative ultrasound: measurement considerations for the assessment of muscular dystrophy and sarcopenia. Front Aging Neurosci 2014, 6: 172.
Noble JA, Boukerroui D: Ultrasound image segmentation: a survey. Med Imag IEEE Trans On 2006,25(8):987–1010.
Seabra J, Ciompi F, Radeva P, Sanches J: A Rayleigh Mixture Model for IVUS Imaging. In Ultrasound Imaging. Edited by: Sanches JM, Laine AF, Suri JS. USA: Springer; 2012:25–47.
Belaid A, Boukerroui D, Maingourd Y, Lerallut JF: Phase-based level set segmentation of ultrasound images. IEEE Trans Inf Technol Biomed 2011,15(1):138–147.
Rodtook A, Makhanov SS: Multi-feature gradient vector flow snakes for adaptive segmentation of the ultrasound images of breast cancer. J Vis Comm Image Represent 2013,24(8):1414–1430. 10.1016/j.jvcir.2013.09.009
Huang J, Yang X, Chen Y, Tang L: Ultrasound kidney segmentation with a global prior shape. J Vis Comm Image Represent 2013,24(7):937–943. 10.1016/j.jvcir.2013.05.013
Sarti A, Corsi C, Mazzini E, Lamberti C: Maximum likelihood segmentation of ultrasound images with Rayleigh distribution. IEEE Trans Ultrason Ferroelectr Freq Control 2005,52(6):947–960.
Deng Y, Wang Y, Shen Y: Speckle reduction of ultrasound images based on Rayleigh-trimmed anisotropic diffusion filter. Pattern Recogn Lett 2011,32(13):1516–1525. 10.1016/j.patrec.2011.06.007
Pereyra M, Dobigeon N, Batatia H, Tourneret JY: Segmentation of skin lesions in 2-D and 3-D ultrasound images using a spatially coherent generalized Rayleigh mixture model. IEEE Trans Med Imaging 2012,31(8):1509–1520.
Wang G, Dong Q, Pan Z, Zhao X, Yang J, Liu C: Active contour model for ultrasound images with Rayleigh distribution. Math Probl Eng 2014, 2014: 12.
Zhang F, Yoo YM, Koh LM, Kim Y: Nonlinear diffusion in Laplacian pyramid domain for ultrasonic speckle reduction. IEEE Trans Med Imaging 2007,26(2):200–211.
Marsousi M, Eftekhari A, Kocharian A, Alirezaie J: Endocardial boundary extraction in left ventricular echocardiographic images using fast and adaptive B-spline snake algorithm. Int J Comput Assist Radiol Surg 2010,5(5):501–513. 10.1007/s11548-010-0404-0
Chen B, Yuen PC, Lai J-H, Chen W-S: Image segmentation and selective smoothing based on variational framework. J Signal Process Syst 2009,54(1–3):145–158.
Khare A, Khare M, Jeong Y, Kim H, Jeon M: Despeckling of medical ultrasound images using Daubechies complex wavelet transform. Signal Process 2010,90(2):428–439. 10.1016/j.sigpro.2009.07.008
Gupta D, Anand RS, Tyagi B: Despeckling of ultrasound medical images using ripplet domain nonlinear filtering. Signal Image Video Process 2013. doi:10.1007/s11760–013–0540–4:1–19
Unal G, Bucher S, Carlier S, Slabaugh G, Fang T, Tanaka K: Shape-driven segmentation of the arterial wall in intravascular ultrasound images. IEEE Trans Inf Technol Biomed 2008,12(3):335–347.
Yang X, Jin J, Xu M, Wu H, He W, Yuchi M, Ding M: Ultrasound common carotid artery segmentation based on active shape model. Comput Math Methods Med 2013, 2013: 11.
Stefan D: Prostate ultrasound image processing. Crossroads 2007,13(3):8. 10.1145/1247238.1247246
Candès E, Demanet L, Donoho D, Ying L: Fast discrete curvelet transforms. Multiscale Modeling Simulation 2006,5(3):861–899. 10.1137/05064182X
Candes EJ, Donoho DL: Curvelets—A Surprisingly Effective Non-Adaptive Representation for Objects With Edges. Curves and Surfaces Nashville, TN: Vanderbilt Univ Press; 1999:123–143.
Miri MS, Mahloojifar A: Retinal image analysis using curvelet transform and multistructure elements morphology by reconstruction. IEEE Trans Biomed Eng 2011,58(5):1183–1192.
Yun T, Xu YQ, Cao L: Semi-supervised ultrasound image segmentation based on curvelet features. Appl Mech Mater 2012, 239–240: 104–114.
Jokar E, Pourghassem H: Kidney Segmentation in Ultrasound Images Using Curvelet Transform and Shape Prior. 2495424: IEEE Computer Society: Proceedings of the 2013 International Conference on Communication Systems and Network Technologies; 2013:180–185.
Wagner RF, Smith SW, Sandrik JM, Lopez H: Statistics of speckle in Ultrasound B-Scans. Sonics Ultrasonics IEEE Trans On 1983,30(3):156–163.
Yongjian Y, Acton ST: Speckle reducing anisotropic diffusion. Image Process IEEE Trans On 2002,11(11):1260–1270. 10.1109/TIP.2002.804276
Hamadani N Master’s thesis, Air Force Inst Technol, WPAFB, OH. Automatic Target Cueing in IR Imagery 1981.
Bangham JA, Ling PD, Aldridge RV, Harvey RW: Morphological scale-space preserving transforms in many dimensions. J Electron Imag 1996,5(3):283–299. 10.1117/12.243349
Mukhopadhyay S, Chanda B: Multiscale morphological segmentation of gray-scale images. Image Process IEEE Trans On 2003,12(5):533–549. 10.1109/TIP.2003.810757
Madabhushi A, Metaxas DN: Combining low-, high-level and empirical domain knowledge for automated segmentation of ultrasonic breast lesions. IEEE Trans Med Imaging 2003,22(2):155–169. 10.1109/TMI.2002.808364
Abd-ElGawad EA, Ibraheem MA, Fouly EH: Evaluation of supraspinatus muscle tears by ultrasonography and magnetic resonance imaging in comparison with surgical findings. Egyptian J Radiol Nucl Med 2013,44(4):829–834. 10.1016/j.ejrnm.2013.08.001
Singh JP: Shoulder ultrasound: what you need to know. Indian J Radiol Imag 2012,22(4):284–292. 10.4103/0971-3026.111481
Rutten MJ, Jager GJ, Kiemeney LA: Ultrasound detection of rotator cuff tears: observer agreement related to increasing experience. AJR Am J Roentgenol 2010,195(6):W440-W446. 10.2214/AJR.10.4526
The authors would like to thank UTP for their assistance, and Ministry of Science, Technology and Innovation (MOSTI) and Ministry of Education of Malaysia for sponsoring the project under the E-Science Grant Scheme, E-Science-03-02-02-SF0109 and FRGS Grant Scheme, FRGS-0153AB-I55 respectively.
The authors declare that they have no competing interests.
RG, EI, SCD, JG contributed to the study concept and design. JG, FZ contributed in the acquisition of data and assessment of the algorithm. RG, EI, SCD were involved in the interpretation and analysis of the acquired assessment results. RG, EI, SCD contributed to the methodological development of algorithmic concepts and discussions. PV, IF participated in the statistical analyses of the results and critical revision of the manuscript for intellectual content. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gupta, R., Elamvazuthi, I., Dass, S.C. et al. Curvelet based automatic segmentation of supraspinatus tendon from ultrasound image: a focused assistive diagnostic method. BioMed Eng OnLine 13, 157 (2014). https://doi.org/10.1186/1475-925X-13-157
- Ultrasound image segmentation
- Supraspinatus (SSP) tendon
- Curvelet transform
- Morphological operations