Quantitative ultrasound tissue characterization in shoulder and thigh muscles – a new approach
© Nielsen et al; licensee BioMed Central Ltd. 2006
Received: 14 January 2005
Accepted: 18 January 2006
Published: 18 January 2006
The echogenicity patterns of ultrasound scans contain information of tissue composition in muscles. The aim was: (1) to develop a quantitative ultrasound image analysis to characterize tissue composition in terms of intensity and structure of the ultrasound images, and (2) to use the method for characterization of ultrasound images of the supraspinatus muscle, and the vastus lateralis muscle.
Computerized texture analyses employing first-order and higher-order grey-scale statistics were developed to objectively characterize ultrasound images of m. supraspinatus and m. vastus lateralis from 9 healthy participants.
The mean grey-scale intensity was higher in the vastus lateralis muscle (p < 0.05) than in the supraspinatus muscle (average value of middle measuring site 51.4 compared to 35.0). Furthermore, the number of spatially connected and homogeneous regions (blobs) was higher in the vastus lateralis (p < 0.05) than in the supraspinatus (average for m. vastus lateralis: 0.092 mm-2 and for m. supraspinatus: 0.016 mm-2).
The higher intensity and the higher number of blobs in the vastus lateralis muscle indicates that the thigh muscle contained more non-contractile components than the supraspinatus muscle, and that the muscle was coarser. The image analyses supplemented each other and gave a more complete description of the tissue composition in the muscle than the mean grey-scale value alone.
Conventional real-time B-mode ultrasonography is a valuable tool for quantitative assessment of skeletal muscle thickness, and high positive correlations have been found between muscle thickness and muscle volume and strength , respectively. Further, different tissue types reflect ultrasound differently. Muscle tissue appears dark whereas fat and connective tissue are bright in ultrasound images [2–4]. Morphological changes and diversity regarding the distribution of contractile and non-contractile tissue has been identified in the human muscle due to age, exposure, training condition, miscellaneous diseases, oedema, and pathological alterations. Healthy muscles look dark with sharp, bright lines whereas muscles with different neuromuscular diseases like muscle dystrophy, myopathy or motor neuron disorders are brighter and more diffuse in the structure when evaluated by ultrasonography [2, 3, 5–8]. In elderly well-trained muscles, the components are depicted distinctly with dark muscle tissue and sharp bright lines of connective tissue, whereas young diseased muscles and elderly untrained muscles are more brightly speckled compared to young healthy muscles and elderly trained muscles [9, 10]. Because invasive biopsy technique is not suitable as a diagnostic procedure for larger groups, there is an obvious need for non-invasive screening methods applicable for evaluation of whole muscles.
Ultrasonographic evaluation of the human muscles has been performed in studies on the quadriceps muscle, the arm, low-back, and masticatory muscles [6, 11, 12]. Most studies used a subjective evaluation in terms of a visual inspection of the images or measurements of the muscle thickness, and only a few studies have used quantitative ultrasound to characterize the various components of the muscle tissue [4, 11, 13–15]. The quantitative evaluation included the distribution of the grey-scale levels of the pixels in the images to describe tissue characteristics. However, muscles with very different muscle structure in terms of the spatial distributions determined by quantitative ultrasound image analysis could result in identical grey-scale intensity. Even if diseased muscles typically look diffused or unstructured, a well-trained muscle could be characterized by almost the same mean grey-scale value because of a distinct intramuscular tendon, e.g. in the supraspinatus muscle. Thus, supplementary analyses are needed to characterize muscles effectively. Therefore, we developed a method, a so-called "blob analysis", related to the higher order grey-scale statistics of the image, i.e. the structures. The blob analysis aims at detecting spatially connected and homogeneous (in terms of grey level) regions ( = blobs) in the image. The frequency distribution (histogram) of the size (area, mm2) and number (mm-2) of blobs detected characterize the image in terms of the amount and extension of spatial structures present. It has been hypothesized that various forms of musculoskeletal diseases and overload may facilitate formation of non-contractile tissue in muscles [16, 17].
The aim of the present study was to develop a non-invasive method based on quantitative ultrasound image analysis to characterize tissue composition in terms of intensity and structure.
The right and left supraspinatus muscles, and the right vastus lateralis muscle were scanned in nine participants (7 females and 2 males) with the mean age of 35.8 (23–47) years. All participants except one were right-handed and all were without musculoskeletal disorders. None of them had performed physical exercise prior to the measurements. In addition, the same muscles were scanned in one right-handed female laboratory technician (56 years old) with severe work-related musculoskeletal symptoms in the right shoulder region related to the supraspinatus muscle. The local ethical committee approved the study, and all participants were fully informed about the protocol and procedures, before giving their consent.
Ultrasonography was performed with an ultrasound scanner (Diasonics VST Masters Series System, GE-Medical, Denmark: "Small parts") fitted with a 5 (3.5–7.0) MHz transducer (5.0 MHz 60R CLA Curved Linear Array HI-Density Probe, GE-Medical, Denmark). The digital recordings were stored on discs. The settings of the ultrasound system was standardized for all the participants and kept constant during all measurements to avoid potential changes in the images due to different setttings. We used a depth setting of 5 cm. Depth-gain compensation was built-in in the used scanner, meaning that the gain automatically was increased, while the depth increased to compensate for the loss through the tissue scanned. The acoustic signal received by the ultrasound transducer was converted into digital 8 bit intensity values making the ultrasound image. Each image consisted of pixels with a size of 0.13 mm × 0.13 mm = 0.017 mm2. As the approximate speed of sound through an average tissue is 1540 m/s and the transducer used was 5 MHz, the wavelength was approximately 0.3 mm [2, 3, 18]. During the ultrasound scanning, the participants were sitting upright in a chair with the upper arms relaxed and kept in a vertical position and with the forearms resting on the thighs. Contact gel was used for acoustic coupling and care was taken not to exert undue pressure on the imaged tissue, since hard pressure from the transducer may influence the grey scale intensity of the pixels in the images. The transducer was placed perpendicular to the surface of the muscles under examination and adjusted to get the brightest image.
M. supraspinatus: to determine the sites for measuring, the distance between the lateral margin of the acromion and the medial margin of the scapulae at the level of the spina scapulae was measured with a ruler, and a line was drawn parallel to and 2 cm cranial to the spina scapulae. Three measuring sites were marked. The middle site was defined as the midpoint of the drawn line. The medial and lateral sites were located on the line 2 cm medially and laterally to the middle measuring site, respectively, both on the right and left supraspinatus muscles. For each of the supraspinatus muscles, the scanning session consisted of measurements at each of the three measuring sites, longitudinally as well as transversely .
M. vastus lateralis: five measuring sites were marked on the right m. vastus lateralis. The middle measuring site was marked in the longitudinal direction by selecting the midpoint between trochanter major and epicondylus lateralis, and in the transverse direction by selecting the midpoint based on muscle palpation. The measuring sites proximal, distal, medial, and lateral, were located 2 cm proximally, distally, medially, and laterally to the middle measuring site.
Grey scale analysis
To evaluate the ultrasound images, two types of computerized image texture analyses programmed in Interactive Data Language (IDL) were used (Interactive Data Language, Research Systems, Inc., Boulder, Colorado, USA). The two analyses derive properties of the image content that are related to intensity, and spatial structure (blob analysis). The blob analysis was developed for the present study. The two types of analyses reflect our hypothesis that the muscle tissue composition can be described by a combination of intensity and structure measures.
The intensity analysis consisted of a computation of the frequency distribution of the grey-scale values of the pixels inside the ROI. All images were filtered by the use of a median filter (kernel size = 3 × 3 pixels) to remove single pixel noise. The reproducibility of this analysis has been described earlier .
Mean blob size, s, in the blob image at threshold T (1); Mean number of blobs, N, for a set of blob images computed at a range of thresholds (2); and Mean blob size over a range of thresholds, S, (3) was calculated as follows:
p j is the size of blob number j in the blob image,
n is the number of blobs,
n i is the number of blobs in the blob image at threshold i,
m is the number of thresholds in the range, and
s i is the mean blob size at threshold i.
Data are presented as group mean of 9 participants. In the tables average values, (range), and [quartiles] are used for the average values of the 9 participants, whereas the mean, the variance, the skewness, the kurtosis, and the median are used to describe the distribution of the pixels in the individual ROI-images.
Data was tested for normality, and non-parametric statistics in terms of the Friedman and Wilcoxon-tests were used. The level of significance was selected at p < 0.05 . The Friedman test was used when comparing different measuring sites, and the Wilcoxon test was used when comparing muscles. The Minitab Statistical Software – Release 14 (Minitab Inc, State College, PA, USA).
Grey scale intensity
First-order grey-scale statistics of m. supraspinatus. M. supraspinatus, dominant arm (D) and non-dominant arm (ND) for n = 9 healthy participants. Average values, (range), and [quartiles] are used for the average value of the 9 participants, whereas the mean, the variance, the skewness, the kurtosis, and the median are used to describe the distribution of the pixels in the individual ROI-images.
First-order grey-scale statistics of m. vastus lateralis. M. vastus lateralis for n = 9 healthy participants. Average values, (range), and [quartiles] are used for the average value of the 9 participants, whereas the mean, the variance, the skewness, the kurtosis, and the median are used to describe the distribution of the pixels in the individual ROI-images.
The intensities for the supraspinatus muscle from the shoulder patient appeared higher than from the healthy participants, i.e. exceeding the 25–75% quartiles. Similar deviations were also found for the vastus lateralis muscle (Tables 1 and 2).
Grey-scale structure (blob analysis)
Corresponding to the anatomical structure and composition each muscle seems to have its own specific image defined by the grey-scale intensity and structure of the muscle based on by the computerized ultrasonic B-scan texture analysis. Such quantitative ultrasound is also a potential examination method to detect possible changes associated with use, age and pathology in whole muscles  where invasive biopsy technique is not suitable, and for screening purposes or supplementary diagnostic procedures.
In this study the gain settings etc were standardized for all the participants and were kept constant during all measurements to avoid differences in the images due to different setttings of the ultrasound scanner. It is evident that different gain settings and transducers as well as different postprocessing of the images will all influence the absolute grey-scale values. Therefore, so far, each laboratory needs to compile their own reference values. In the future it would be preferable to develop image analyses with a higher level of independency of the equipment.
Furthermore, it is important to avoid undue hard pressure on the imaged tissue since pressure from the transducer may influence the grey scale intensity of the pixels in the images due to changes in the muscle shape and thereby minor changes in the insonation angles. In the present study care was taken that only little pressure was applied on the muscles. As the difference in grey-scale intensity between the shoulder and the thigh muscle was far larger than the potential small difference in grey-scale intensity associated with the pressure exerted on the muscle during scanning we do not think this had blurred the conclusion.
Thus, being aware of these limitations in the grey-scale analyses, an increased proportion of non-contractile tissue might be expected in musculoskeletal disorders [21–24]. Accordingly, the affected supraspinatus muscle differed markedly from the non-affected side and the muscles from the healthy participants, with reservations for the interpretation since only one patient participated. However, some differences could also be an effect of ageing where muscle fibres are often replaced by fat and connective tissue  as differences were also present between the patient's non-affected supraspinatus muscle and the vastus lateralis and the muscles in the healthy participants.
The grey-scale intensity of the vastus lateralis muscle was higher than the intensity of the supraspinatus muscle, indicating that vastus lateralis contained more non-contractile components than supraspinatus. This result is in accordance with the general anatomy of the two muscles . The vastus lateralis is a long muscle forming the major part of the largest muscle in the body, the quadriceps femoris, and has two large aponeuroses and a tendon. In contrast the supraspinatus muscle is small and triangular with a short tendinous structure. However, in other studies [14, 15, 27] using a similar method an increased echo intensity of the vastus lateralis muscle was assumed to reflect an increased relative proportion of non-contractile tissue. In the liver Layer et al.  demonstrated that connective tissue only lead to a weak increase in the grey-level intensity, whereas the addition of connective tissue to a given tissue lipid might reduce the image brightness. Thus the intensity is so far not able to distinguish between different textures of fat and connective tissue and whether an increased brightness of the images is associated with fat or connective tissue .
We found differences between the measuring sites of the supraspinatus muscle, indicating inhomogeneity of the muscle. This is in accordance with previous anatomical findings . Therefore more than one measuring site is needed to describe the muscle as a whole. No differences were found between measuring sites in vastus lateralis, but this fact does not necessarily indicate that the muscle is homogeneous as a whole. The thigh muscle is far larger than the shoulder muscles, and the five measuring sites were concentrated on the middle of the muscle.
Since the first-order statistics expresses the distribution of grey-scale values without taking their spatial arrangement into account, muscles with identical echogenicity may have very different tissue structure. Therefore, to improve the fidelity of the image analysis we developed the higher-order grey-scale statistics blob analysis, reflecting the structures in the images, to obtain a more complete description of the scanned muscles. Theoretically, a muscle with a significant tendon structure will result in few and large blobs, whereas a muscle with a diffuse structure will result in many small blobs. This means that the blob analysis may quantitatively separate scans of different muscle composition, even though the mean grey-scale intensity is the same.
The number of the blobs was higher and the size of the blobs tended to be larger in the vastus lateralis muscle compared to the supraspinatus muscle at the level corresponding to non-contractile tissue. These findings indicate as expected from anatomical descriptions that the structures in the vastus lateralis muscle are coarser and larger (Figure 1).
Not only the content of the non-contractile components but also the structure of the muscle are indicated with the newly developed computerized, blob texture analysis, based on the spatial distribution of the grey-scale intensities of the ultrasonic B-scan. Thus, the supraspinatus and the vastus lateralis muscle in the present study seemed to have both different muscle compositions and structures corresponding to their known anatomy. In combination the two quantitative analyses presented, i.e. of the grey-scale intensity and the grey-scale structure (blobs), supplemented each other and gave a more complete description of the tissue composition in the muscle than the mean grey-scale value alone.
The authors wish to thank Senior researcher MSc PhD Bjarne Laursen from the National Institute of Public Health, Copenhagen, Denmark, and Associate Professor MScEE PhD Carsten E. Thomsen from the Department of Oral Medicine, Clinical Oral Physiology, Oral Pathology & Anatomy, School of Dentistry, University of Copenhagen, Denmark, for their advice and technical assistance regarding data acquisition.
- Bakke M, Tuxen A, Vilmann P, Jensen BR, Vilmann A, Toft M: Ultrasound image of human masseter muscle related to bite force, electromyography, facial morphology, and occlusal factors. Scand J Dent Res. 1992, 100: 164-171.PubMedGoogle Scholar
- Heckmatt JZ, Leeman S, Dubowitz V: Ultrasound imaging in the diagnosis of muscle disease. J Pediatr. 1982, 101: 656-660.View ArticlePubMedGoogle Scholar
- Heckmatt JZ, Pier N, Dubowitz V: Real-time ultrasound imaging of muscles. Muscle & Nerve. 1988, 11: 56-65. 10.1002/mus.880110110.View ArticleGoogle Scholar
- Nielsen PK, Jensen BR, Darvann T, Jørgensen K, Bakke M: Quantitative ultrasound image analyses of the supraspinatus muscle. Clin Biomech. 2000, 15 (Suppl 1): S13-S16. 10.1016/S0268-0033(00)00053-X.View ArticleGoogle Scholar
- Lamminen A, Jääskeläinen J, Rapola J, Suramo I: High-frequency ultrasonography of skeletal muscle in children with neuromuscular disease. J Ultrasound Med. 1988, 7: 505-509.PubMedGoogle Scholar
- Kiliaridis S, Engvall M, Tzakis MG: Ultrasound imaging of the masseter muscle in myotonic dystrophy patients. J Oral Rehabil. 1995, 22: 619-625.View ArticlePubMedGoogle Scholar
- Kaschka WP, Druschky KF, Rott HD: Myotonic muscular dystrophy: structural changes visualized by ultrasound. J Neurol. 1987, 234: 122-123. 10.1007/BF00314117.View ArticlePubMedGoogle Scholar
- Heckmatt JZ, Dubowitz V, Leeman S: Detection of pathological change in dystrophic muscle with B-scan ultrasound imaging. The Lancet. 1980, June: 1389-1391. 10.1016/S0140-6736(80)92656-2.View ArticleGoogle Scholar
- Sipilä S, Suominen H: Ultrasound imaging of the quadriceps muscle in elderly athletes and untrained men. Muscle & Nerve. 1991, 14: 527-533. 10.1002/mus.880140607.View ArticleGoogle Scholar
- Sipilä S, Suominen H: Knee extension strength and walking speed in relation to quadriceps muscle composition and training in elderly women. Clin Physiol. 1994, 14: 433-442.View ArticlePubMedGoogle Scholar
- Eisele R, Schmid R, Kinzl L, Kramer P, Katzmaier P, Hartwig E: Soft tissue texture analysis by B-mode ultrasound in the evaluation of impairment in chronic low back pain. European Journal of Ultrasound. 1998, 8: 167-175. 10.1016/S0929-8266(98)00073-1.View ArticlePubMedGoogle Scholar
- Scholten RR, Pillen S, Verrips A, Zwarts MJ: Quantitative ultrasonography of skeletal muscles in children: Normal values. Muscle & Nerve. 2004, 2003: 693-698.Google Scholar
- Schedel H, Reimers CD, Nägele M, Witt TN, Pongratz DE, Vogl T: Imaging techniques in myotonic dystrophy. A comparative study of ultrasound, computed tomography and magnetic resonance imaging of skeletal muscles. Eur J Radiol. 1992, 15: 230-238. 10.1016/0720-048X(92)90113-N.View ArticlePubMedGoogle Scholar
- Sipilä S, Suominen H: Quantitative ultrasonography of muscle: detection of adaptions to training in elderly women. Arch Phys Med Rehabil. 1996, 77: 1173-1178. 10.1016/S0003-9993(96)90143-4.View ArticlePubMedGoogle Scholar
- Hicks JE, Shawker TH, Jones BL, Linzer M, Gerber LH: Diagnostic Ultrasound: Its use in the evaluation of muscle. Arch Phys Med Rehabil. 1984, 65: 129-131.PubMedGoogle Scholar
- Jensen BR, Bakke M: Prolonged work with shoulder muscles and other small muscle groups: Use, function, and pain. In Muscle atrophy: Disuse and disease Edited by: Capodaglio P and Narici MV. Pavia, Italy, Le Collane della Fondazione Salvatore Maugeri; 1998:149-161.Google Scholar
- Hargens AR, Schmidt DA, Evans KL, Gonsalves MR, Cologne JB, Garfin SR, Mubarak SJ, Hagan PL, Akeson WH: Quantitation of skeletal-muscle necrosis in a model compartment syndrome. Bone Joint Surg (Am). 1981, 63-A: 631-636.Google Scholar
- Shirley IM, Blackwell R, Cusick G, Farman DJ, Vicary FR: Ultrasound. A User's Guide to Diagnostic Ultrasound. Edited by: Shirley IM, Blackwell R, Cusick G, Farman DJ and Vicary FR. 1988, , Pitman Medical Publishing Co LtD, 32-39.Google Scholar
- Jensen BR, Jørgensen K, Sjøgaard G: The effect of prolonged isometric contractions on muscle fluid balance. Eur J Appl Physiol. 1994, 69: 439-444. 10.1007/BF00865409.View ArticleGoogle Scholar
- Conover WJ: Practical Nonparametric Statistics. 1980, New York, Chinchester, Brisbane, Toronto, John Wiley & Sons, 1-493. 2Google Scholar
- Kadi F, Hägg G, Håkansson R, Holmner S, Butler-Brown GS, Thornell LE: Structural changes in male trapezius muscle with work-related myalgia. Acta Neurolpathol (Berl). 1998, 95: 352-360. 10.1007/s004010050810.View ArticleGoogle Scholar
- Kadi F, Waling K, Ahlgren C, Sundelin G, Holmner S, Butler-Brown GS, Thornell LE: Pathological mechanisms implicated in localized female trapezius myalgia. Pain. 1998, 78: 191-196. 10.1016/S0304-3959(98)00126-2.View ArticlePubMedGoogle Scholar
- Ljung BO, Lieber RL, Fridén J: Wrist extensor muscle pathology in lateral epicondylitis. J Hand Surg. 1999, 24: 177-183.View ArticleGoogle Scholar
- Hägg GM: Human muscle fibre abnormalities related to occupational load. Eur J Appl Physiol. 2000, 83: 159-165. 10.1007/s004210000274.View ArticlePubMedGoogle Scholar
- Larsson L: Motor units: remodeling in aged animals. J Gerontol A Biol Sci Med Sci. 1995, 50 Spec No: 91-95.PubMedGoogle Scholar
- Rohen JW, Yokochi C: Color Atlas of Anatomy. 1988, New York, Igaku-Shoin Ldt., 1-469. 2Google Scholar
- Pillen S, Scholten RR, Zwarts MJ, Verrips A: Quantitative skeletal muscle ultrasonography in children with suspected neuromuscular disease. Muscle & Nerve. 2004, 2003: 699-705.Google Scholar
- Layer G, Zuna I, Lorenz A, Zerban H, Haberkorn U, Bannasch P, van Kaick G, Rath U: Computerized ultrasound B-scan texture analysis of experimental diffuse parenchymal liver disease: correlation with histopathology and tissue composition. J Clin Ultrasound. 1991, 19: 193-201.View ArticlePubMedGoogle Scholar
- Jensen BR, Jørgensen K, Huijing PA, Sjøgaard G: Soft tissue architecture and intramuscular pressure in the shoulder region. Eur J Morphol. 1995, 33: 205-220.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2474/7/2/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.