- Open Access
Automatic thickness estimation for skeletal muscle in ultrasonography: evaluation of two enhancement methods
BioMedical Engineering OnLine volume 12, Article number: 6 (2013)
Ultrasonography is a convenient technique to investigate muscle properties and has been widely used to look into muscle functions since it is non-invasive and real-time. Muscle thickness, a quantification which can effectively reflect the muscle activities during muscle contraction, is an important measure for musculoskeletal studies using ultrasonography. The traditional manual operation to read muscle thickness is subjective and time-consuming, therefore a number of studies have focused on the automatic estimation of muscle fascicle orientation and muscle thickness, to which the speckle noises in ultrasound images could be the major obstacle. There have been two popular methods proposed to enhance the hyperechoic regions over the speckles in ultrasonography, namely Gabor Filtering and Multiscale Vessel Enhancement Filtering (MVEF).
A study on gastrocnemius muscle is conducted to quantitatively evaluate whether and how these two methods could help the automatic estimation of the muscle thickness based on Revoting Hough Transform (RVHT). The muscle thickness results obtained from each of the two methods are compared with the results from manual measurement, respectively. Data from an aged subject with cerebral infarction is also studied.
It’s shown in the experiments that, Gabor Filtering and MVEF can both enable RVHT to generate comparable results of muscle thickness to those by manual drawing (mean ± SD, 1.45 ± 0.48 and 1.38 ± 0.56 mm respectively). However, the MVEF method requires much less computation than Gabor Filtering.
Both methods, as preprocessing procedure can enable RVHT the automatic estimation of muscle thickness and MVEF is believed to be a better choice for real-time applications.
Ultrasonography, as an important medical image modality in the study of the musculoskeletal system, has been widely used to measure changes in muscle geometry, such as muscle thickness, muscle pennation angle, fascicle length and cross-sectional area[1–8], because it is versatile, inexpensive and radiation-free. More important, the change of muscle properties over time can form signals, representing architectural muscle behavior under contraction, namely sonomyography (SMG)[9–12]. Together with electromyography (EMG)[13–15], SMG can provide quantitative information of muscle contraction, help correlate the actions of several muscles, either agonistic or antagonistic, and subsequently advance medical and engineering efforts to understand, diagnose or even treat musculoskeletal disorders[5, 16–20].
Specifically speaking, muscle thickness (MT), which can effectively reflect the muscle activities during muscle contraction[21–24], is an important measure for musculoskeletal studies using ultrasonography. There have been reports on relationships between MT and many other measures, such as level of physical activity, age, muscle stiffness and response to exercise. Meanwhile, for major muscles such as quadriceps and gastrocnemius muscle (GM), another important parameter, fascicle length (FL) can hardly be measured directly since it’s beyond the dimension of most ultrasound images. Therefore FL is usually based on MT and pennation angle (for example,[29, 30]). However, MT measurement is traditionally based on manual selection of 2 reference points at the superficial and deep aponeuroses and usage of an on-screen caliper to obtain the MT frame by frame[31–34], which is time-consuming and subjective. There have been several studies targeting at automatic estimation of MT, such as the Revoting Hough Transform (RVHT) method[35, 36], Radon Transform, Region of Interest (ROI) and cross-correlation algorithm. Taking RVHT as an example, with assumption that the superficial and deep aponeuroses could be located using RVHT as the very first 2 lines detected, the MT could be computed readily as the distance between them. However, besides the fact that the speckle noise associated with ultrasound image, sarcopenia or loss of muscle mass due to ageing, can affect the quality of the images obtained by ultrasound from deconditioned muscles. This can make it difficult to accurately locate the borders of the muscles to obtain an accurate measure of the MT. There have been two popular methods proposed to enhance the hyperechoic regions over the speckles in ultrasonography, namely Gabor Filtering[36, 40, 41] and Multiscale Vessel Enhancement Filtering (MVEF)[42, 43]. In this paper, to investigate whether image enhancement can make sufficient preparation for automatic estimation of MT and which method is more appropriate for MT estimation, we will evaluate both Gabor Filtering and MVEF on gastrocnemius muscle images.
The automatic MT estimation procedure, based on RVHT method proposed previously, could include four steps: 1) image enhancement, 2) thresholding operation to generate black-white image and 3) locating of superficial and deep aponeuroses by RVHT and 4) computation of the distance between aponeuroses. The RVHT first computes accumulator matrix of Hough transform based on a black-white image called an “edge map”. This edge map represents the meaningful image contents. The RVHT method then locates the global maximum in the accumulator matrix of Hough transform that corresponds to the most dominant line-shaped feature points globally, using the standard Hough transform. Then the pixels close to the detected line are removed from the edge map and the Hough transform accumulator matrix is calculated again. The same procedure could be executed to search for another line. For ultrasound images of skeletal muscles, usually the very first two lines detected would be the superficial and deep aponeuroses according to a priori knowledge. Finally the mean distance between each two lines detected by RVHT is computed, and the two lines which have maximum distance between them are recognized as the superficial and deep aponeuroses, and finally the mean distance between them is computed as the MT. Diagram for the procedures is shown in Figure1.
Taking into the account that the fascicles and fibroadipose spetas in a sonogram are tubular and include coherent orientation tendencies, the Gabor Filter bank method can be used to implement the image enhancement. The method includes three steps: the orientation filed estimation, frequency map computation and Gabor filtering confined by the orientation reliability. More details can be consulted in[36, 44].
Multiscale vessel enhancement filtering
The Multiscale Vessel Enhancement Filtering method is based on the second order local structure, with excellent noise and background suppression performance. The method also includes three steps: the Hessian matrix estimation (including the choice of Gaussian kernels), computation of eigenvector for each scale and processing for the maximum vesselness response. More details can be found in[42, 43].
Three healthy male subjects (mean ± SD, age = 28.6 ± 0.6 years; body weight 67.0 ± 1.7 kg; height = 1.72 ± 0.01 m) volunteered to participate in this study. No participant had a history of neuromuscular disorders, and all were aware of experimental purposes and procedures. Human subject ethical approval was obtained from the relevant committee from the Hong Kong Polytechnic University, and informed consent was obtained from the subject prior to the experiment.
Experiment protocol and data acquisition
The testing position of the subject was in accordance with the User’s Guide of a Norm dynamometer (Humac/Norm Testing and Rehabilitation System, Computer Sports Medicine, Inc., Massachusetts, USA). The subject was instructed to generate dorsiflexion/plantar-flexion movements in prone position. The torque was measured by the aforementioned dynamometer and muscle contracting in a range from 0 to 90% maximal voluntary contraction (MVC) is imaged by ultrasonography. 90% MVC is set as the highest value to avoid muscle fatigue.
A real-time B-mode ultrasonic scanner (EUB-8500, Hitachi Medical Corporation, Tokyo, Japan) with an electronic linear array probe (EUP-L53L, 5.0-10.0 MHz Hitachi Medical Corporation, Tokyo, Japan) was used to obtain ultrasound images of muscles. The long axis of the ultrasound probe was arranged parallel to the long axis of the GM and on its muscle belly, decided accozrding to operator’s palpation. The ultrasound probe was fixed by a custom-designed foam container with fixing straps, and a very generous amount of ultrasound gel was applied to secure acoustic coupling between the probe and skin during muscle contractions, as shown in Figure2. The probe was adjusted to optimize the contrast of muscle fascicles in ultrasound images. Then the B-mode ultrasound images were digitized by a video card (NI PCI-1411, National Instruments, Austin, USA) at a rate of 25 frame/s for later analysis.
The experiments detailed in last section ended up in 300 ultrasound images and all images were cropped to keep the image content only. All data were processed off-line using programs written in Matlab 7.11 (MathWorks, Inc., Massachusetts, USA) and the correlation and agreement analysis were performed using the computer package SPSS/PC 4.0 (SPSS Inc, Chicago, USA).
Experiment setup for an aged subject with cerebral infarction
One male subject with unilateral limb dysfunction caused by cerebral infarction (age = 68 years; body weight = 71 kg; height = 1.72 m; right leg dysfunctional) volunteered to participate in this study. The human subject ethical approval was obtained from the relevant committee in Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences before carrying out the experiment. The subject was briefed about the procedure of experiment and written consents were collected prior to the experiment. During measurement, the subject performed plantarflextion both in left leg (Normal) and in right leg (Abnormal). A laptop ultrasound system (SS-10, Sonostar Technologies Co., Limited, Guangzhou, China), with a 7.5 MHz electronic linear array probe, was used to obtain the dynamic ultrasonic image sequences. The long axis of the probe was arranged parallel to the long axis of the gastrocnemius muscle. 64 consecutive frames from both legs are captured.
In this section, experimental results from the comparison between the performance of the Gabor Filtering method and the one of the MVEF method over a series of clinical gastrocnemius muscle ultrasound images were first presented. The originally acquired images (640 × 480 pixels) were first cropped to keep the image content only, and the cropped images (403 × 373 pixels) were then enhanced by Gabor Filtering and MVEF methods respectively and binarized, as shown in Figure3a-c. It should be noted that, in Gabor Filtering, empirically the default value of the reliability of the local orientation of image was set to 50%, which means that the corresponding location would be ignored unless when the reliability of the orientation is higher than 50%. The superficial and deep aponeuroses were then located by RVHT and MT was estimated automatically, as shown in Figure3b-c. For all 300 frames, after either Gabor Filtering or MVEF method, the superficial and deep aponeuroses had been detected by RVHT as the very first two lines, as expected.
In order to make a quantitative evaluation of the performance of the MT estimated by RVHT after the two methods, two operators were asked to manually draw lines and read MT of the 300 images digitally using a software (NIH Image, National Institutes of Health, USA). The superficial and deep aponeuroses were located independently and blindly in each image. And then the mean of the two operators’ results, named as “Op”, was considered to be a close estimation of the true value of MT. The P-P plots of the results obtained by manual drawing and the two automatic methods are shown in Figure4a-c, suggesting that roughly estimated, the results data were all normally distributed. Therefore, based on the current recommendations for reliability analysis[45, 46], inter-class correlation coefficient, intra-class correlation coefficient (ICC) and Bland-Altman plots were used to further explore the results. The correlation analysis of Op and the corresponding values estimated by RVHT after the two enhancement methods are all shown in Figure5a-c. The square of the Pearson product moment inter-class correlation coefficients of the above three comparison study were 91.3%, 91.3% and 87.8% respectively. Figure6 presents the details of the MT statistics using Bland-Altman plots. The differences between the manually-obtained results and the corresponding values obtained by using RVHT after Gabor Filtering and MVEF method were 1.45 ± 0.48 and 1.38 ± 0.56 mm, respectively. The difference between the results obtained by RVHT after Gabor Filtering and MVEF method was 0.07 ± 0.54 mm. The computation of ICCs was conducted using the SPSS, based on a single rating, absolute-agreement, 2-way mixed effect model. The standard error of the measurements (SEMs) was also calculated to quantify the measurement precision (i.e., where s x is the standard deviation). The minimal detectable changes (MDCs) were calculated as, which represents the minimal change in MT that must occur to be 95% confident that a true change occurred. The mean values ± SDs, ICC, SEMs, and MDCs in each subject are presented in Table1. The measurements based on RVHT after Gabor Filtering decreased the SEM by 12.9% as that of MVEF method. The MDCs based on RVHT after Gabor Filtering and MVEF method were 2.3% vs. 2.2%, 1.3% vs. 2.3%, 1.5% vs. 1.3% of the MT, respectively. The use of the mean among 3 subjects as an outcome measure resulted in MDCs of 1.8% (Gabor Filtering) and 2.1% (MVEF method) of the MT, respectively.
Two typical frames from left and right legs of the aged subject are shown in Figure7a-b. MTs estimated using the proposed method with MVEF are shown in Figure7c-d. For both group of data, in 5 out of total 128 (64 for each leg) frames, the MTs estimated are significantly different from the neighboring values. To have a better observation on the transitions of MT for both legs during the plantarflextion excercises, simple median filters were used to render the smooth versions. Frame number, corresponding to MTs with difference > 0.5 mm between the original and smooth versions are, are also marked out in Figure7. MT for the left and right legs are 15.06 ± 0.86 and 14.30 ± 0.76 mm respectively.
As both Gabor Filtering and MVEF are established methods for skeletal muscle enhancement[36, 43], in this work we focus on their performance as a pre-processing step of automatic estimation of MT. It’s hoped that this study would provide consultative guide for widespread application of the computerized MT estimation, thus to replace the time-consuming and subjective manual drawing.
First of all, to estimate the MT automatically, we use RVHT to detect lines in ultrasound images of skeletal muscle, and the superficial and deep aponeuroses are expected to be the very first two lines detected. Without enhancement procedure before RVHT, the performance of aponeuroses detection is quite poor, indicated by the fact that the superficial and deep aponeuroses could not be located as the very first 2 lines in 110 of 300 images. However, as mentioned before, it becomes 300 out of 300 images after either image enhancement method. In other words, for the automatic estimation of MT, we must drop the assumption that the superficial and deep aponeuroses are the strongest lines in ultrasound images, unless a proper image enhancement procedure is used.
As for the quantitative performance of the two methods, MT estimated by RVHT after the two methods have high correlation (R2 = 91.3% and 87.8% respectively) with small differences (1.45 ± 0.48 and 1.38 ± 0.56 mm respectively, when MT is about 38-46 mm) to manual drawing measurement. And Bland-Altman plots in Figure6 also show the good agreement between the MT estimated by RVHT after the two methods and manual results.
It should be noted that, as shown in Figure6, the differences of MT estimated by RVHT after either Gabor Filtering or MVEF method and that obtained by manual drawing is about 1.4 mm, and only in 3 out of 300 images the manual results were larger than RVHT method after MVEF, and none is larger after Gabor Filtering. It’s believed that this bias is due to the fact that in manual method, MT was calculated as the border/nearest distance between the superficial and deep aponeuroses in the middle of the ultrasound image at a 90-degree angle from the deep aponeurosis, while in automatic method using RVHT, the MT will be estimated as roughly the middle-line distances between superficial and deep aponeuroses. In other words, by nature there is this bias due to the different acting definitions of MT between the manual and automatic method using RVHT. We have calculated the thickness of superficial aponeuroses using active contour for the 300 images and the result is about 1.51 ± 0.05 mm, which could further reduce the difference/bias and add a piece of proof to the deduction that the bias is mainly caused by the mean thickness of aponeuroses themselves. On the other hand, we’d like to point out that although the mean results of two operator is taken to be the ‘true’ thickness, the two operators also, not surprisingly as shown in Figure5a, do not output identical results of MT and the difference is 0.44 ± 0.48 mm.
The time costs of the two methods to process one image were about 5.1 s and 0.3 s in Matlab, respectively, suggesting that MVEF is more suitable for real-time applications where MT transitions are often interested.
The results on the aged subject with cerebral infarction further demonstrated the feasibility of the proposed method as an automatic method for muscle MT estimation. Interestingly the dysfunctional right leg exhibited smaller MT compared to left leg, by 94.9% and 89.1%, in terms of mean and standard deviation respectively. We’d like to note that the current report is the first one to prove the feasibility of automatic measurement of MT in longitudinal direction, with the help of image enhancement techniques. The automatic measurement of MT then can make possible the continuous monitoring of the time course of MT along with muscle contraction. It’s objective and not labor-intensive, especially when the images to analyze are numerous.
In the further research, more musculoskeletal images should be collected including those under pathological conditions. Problems are expected to arise, for example, acute patients usually can not adopt positions such as standing, or prone lying to allow optimal access to measurement sites, which may increase the likelihood of measurement errors. And the corresponding further improvements on the automatic estimation of MT would certainly broaden the application area of ultrasound imaging in clinical musculoskeletal system, such as, study of the relationship between MT changes and tumbling for Tumble Monitoring and Identification (TMI).
In conclusion, this study aims at automatic estimation of skeletal MT using RVHT, and an evaluation of two image enhancement methods in ultrasound images, as a preprocessing step to make sure that superficial and deep aponeuroses can be detected by RVHT as the very first 2 line-shaped features. Experimental results from 300 images in total showed that, both methods can provide effective preparation for RVHT and to further estimate MT. Specifically speaking, the superficial and deep aponeuroses could not be located by RVHT as the very first 2 lines in 110 of 300 images without image enhancement. However, as mentioned before, it becomes 300 out of 300 images after either image enhancement method. Meanwhile, the results from both methods have high correlation and great agreement with manually-obtained results, while the MVEF runs much faster. This study demonstrated that with the help of proper image enhancement preprocessing procedures, RVHT can be used for automatic estimation of skeletal MT.
Revoting hough transform
Multiscale vessel enhancement filtering
Intra-class correlation coefficient.
Farris DJ, Sawicki GS: Human medial gastrocnemius force–velocity behavior shifts with locomotion speed and gait. Proc Natl Acad Sci 2012, 109: 977–982. 10.1073/pnas.1107972109
Fukunaga T, Kawakami Y, Kuno S, Funato K, Fukashiro S: Muscle architecture and function in humans. J Biomech 1997, 30: 457–463. 10.1016/S0021-9290(96)00171-6
Hodges P, Pengel L, Herbert R, Gandevia S: Measurement of muscle contraction with ultrasound imaging. Muscle Nerve 2003, 27: 682–692. 10.1002/mus.10375
Ito M, Kawakami Y, Ichinose Y, Fukashiro S, Fukunaga T: Nonisometric behavior of fascicles during isometric contractions of a human muscle. J Appl Physiol 1998, 85: 1230–1235.
Maganaris CN, Baltzopoulos V, Sargeant AJ: In vivo measurements of the triceps surae complex architecture in man: implications for muscle function. J Physiol 1998, 512: 603–614. 10.1111/j.1469-7793.1998.603be.x
Miyoshi T, Kihara T, Koyama H, Yamamoto SI, Komeda T: Automatic detection method of muscle fiber movement as revealed by ultrasound images. Medical engineering & physics 2009, 31: 558–564. 10.1016/j.medengphy.2008.11.004
Narici MV, Maganaris CN, Reeves ND, Capodaglio P: Effect of aging on human muscle architecture. J Appl Physiol 2003, 95: 2229–2234.
Reeves ND, Maganaris CN, Narici MV: Ultrasonographic assessment of human skeletal muscle size. Eur J Appl Physiol 2004, 91: 116–118. 10.1007/s00421-003-0961-9
Shi J, Zheng Y, Huang Q, Chen X: Continuous monitoring of sonomyography, electromyography and torque generated by normal upper arm muscles during isometric contraction: sonomyography assessment for arm muscles. Biomedical Engineering, IEEE Transactions on 2008, 55: 1191–1198.
Xie HB, Zheng YP, Guo JY, Chen X, Shi J: Estimation of wrist angle from sonomyography using support vector machine and artificial neural network models. Medical engineering & physics 2009, 31: 384–391. 10.1016/j.medengphy.2008.05.005
Zheng YP, Chan M, Shi J, Chen X, Huang QH: Sonomyography: Monitoring morphological changes of forearm muscles in actions with the feasibility for the control of powered prosthesis. Medical engineering & physics 2006, 28: 405–415. 10.1016/j.medengphy.2005.07.012
Zhou Y, Li J, Zhou G, Zheng Y: Dynamic measurement of pennation angle of gastrocnemius muscles during contractions based on ultrasound imaging. Biomed Eng Online 2012, 11: 63. 10.1186/1475-925X-11-63
Cram JR, Kasman GS, Holtz J: Introduction to surface electromyography. MD: Aspen Publishers Gaithersburg; 1998.
De Luca CJ: The use of surface electromyography in biomechanics. J Appl Biomech 1997, 13: 135–163.
Hodges PW, Bui BH: A comparison of computer-based methods for the determination of onset of muscle contraction using electromyography. Electroencephalography and Clinical Neurophysiology/Electromyography and Motor Control 1996, 101: 511–519. 10.1016/S0921-884X(96)95190-5
Chehrehrazi M, Arab AM, Karimi N, Zargham M: Assessment of pelvic floor muscle contraction in stress urinary incontinent women: comparison between transabdominal ultrasound and perineometry. International Urogynecology Journal 2009, 20: 1491–1496. 10.1007/s00192-009-0977-8
Kim YJ, Brox T, Feiden W, Weickert J: Fully automated segmentation and morphometrical analysis of muscle fiber images. Cytometry A 2007, 71: 8–15.
Manshadi FD, Parnianpour M, Sarrafzadeh J, Kazemnejad A: Abdominal hollowing and lateral abdominal wall muscles’ activity in both healthy men & women: An ultrasonic assessment in supine and standing positions. J Bodyw Mov Ther 2011, 15: 108–113. 10.1016/j.jbmt.2009.10.004
Shi J, Yan Z, Xue L, Zheng Y: Evaluation of the Muscle Fatigue Based on Ultrasound Images. The Six World Congress on Intelligent Control and Automation 2006, 2: 10332–10336.
Walker FO, Cartwright MS, Wiesler ER, Caress J: Ultrasound of nerve and muscle. Clin Neurophysiol 2004, 115: 495–507. 10.1016/j.clinph.2003.10.022
Blumhagen J, Noble H: Muscle thickness in hypertrophic pyloric stenosis: sonographic determination. Am J Roentgenol 1983, 140: 221–223. 10.2214/ajr.140.2.221
Campbell IT, Watt T, Withers D, England R, Sukumar S, Keegan MA, Faragher B, Martin DF: Muscle thickness, measured with ultrasound, may be an indicator of lean tissue wasting in multiple organ failure in the presence of edema. Am J Clin Nutr 1995, 62: 533–539.
Common MTD: Changes in deep abdominal muscle thickness during common trunk-strengthening exercises using ultrasound imaging. J Orthop Sports Phys Ther 2008, 38: 596–605.
Wong A, Gallagher K, Callaghan J: Computerised system for measurement of muscle thickness based on ultrasonography. Computer methods in biomechanics and biomedical engineering 2012. 10.1080/10255842.2011.617003
Ichinose Y, Kanehisa H, Ito M, Kawakami Y, Fukunaga T: Morphological and functional differences in the elbow extensor muscle between highly trained male and female athletes. Eur J Appl Physiol Occup Physiol 1998, 78: 109–114. 10.1007/s004210050394
Kubo K, Kanehisa H, Azuma K, Ishizu M, Kuno S, Okada M, Fukunaga T: Muscle architectural characteristics in young and elderly men and women. Int J Sports Med 2003, 24: 125–130. 10.1055/s-2003-38204
Ikezoe T, Asakawa Y, Fukumoto Y, Tsukagoshi R, Ichihashi N: Associations of muscle stiffness and thickness with muscle strength and muscle power in elderly women. Geriatr Gerontol Int 2012, 12(1):86–92. 10.1111/j.1447-0594.2011.00735.x
Starkey DB, Pollock ML, Ishida Y, Welsch MA, Brechue WF, Graves JE, Feigenbaum MS: Effect of resistance training volume on strength and muscle thickness. Medicine & Science in Sports & Exercise 1996, 28: 1311. 10.1097/00005768-199610000-00016
Blazevich AJ, Gill ND, Deans N, Zhou S: Lack of human muscle architectural adaptation after short‐term strength training. Muscle Nerve 2007, 35: 78–86. 10.1002/mus.20666
Kumagai K, Abe T, Brechue WF, Ryushi T, Takano S, Mizuno M: Sprint performance is related to muscle fascicle length in male 100-m sprinters. J Appl Physiol 2000, 88: 811–816.
Gill NW, Teyhen DS, Lee IE: Improved contraction of the transversus abdominis immediately following spinal manipulation: a case study using real-time ultrasound imaging. Man Ther 2007, 12: 280. 10.1016/j.math.2006.06.014
Koppenhaver SL, Hebert JJ, Fritz JM, Parent EC, Teyhen DS, Magel JS: Reliability of rehabilitative ultrasound imaging of the transversus abdominis and lumbar multifidus muscles. Arch Phys Med Rehabil 2009, 90: 87–94. 10.1016/j.apmr.2008.06.022
Moreau NG, Simpson KN, Teefey SA, Damiano DL: Muscle architecture predicts maximum strength and is related to activity levels in cerebral palsy. Phys Ther 2010, 90: 1619–1630. 10.2522/ptj.20090377
Raney NH, Teyhen DS, Childs JD: Observed changes in lateral abdominal muscle thickness after spinal manipulation: a case series using rehabilitative ultrasound imaging. J Orthop Sports Phys Ther 2007, 37: 472–479.
Zhou Y, Zheng Y: Estimation of Muscle Fiber Orientation in Ultrasound Images Using Revoting Hough Transform (RVHT). Ultrasound in Medicine & Biology 2008, 34: 1474–1481. 10.1016/j.ultrasmedbio.2008.02.009
Zhou Y, Zheng YP: Longitudinal Enhancement of the Hyperechoic Regions in Ultrasonography of Muscles Using a Gabor Filter Bank Approach: A Preparation for Semi-Automatic Muscle Fiber Orientation Estimation. Ultrasound Med Biol 2011, 37: 665–673. 10.1016/j.ultrasmedbio.2010.12.011
Zhao H, Zhang LQ: Automatic tracking of muscle fascicles in ultrasound images using localized radon transform. IEEE Trans Biomed Eng 2011, 58: 2094–2101.
Koo TKK, Wong C, Zheng Y: Reliability of sonomyography for pectoralis major thickness measurement. J Manipulative Physiol Ther 2010, 33: 386–394. 10.1016/j.jmpt.2010.05.009
English CK, Thoirs KA, Fisher L, McLennan H, Bernhardt J: Ultrasound Is a Reliable Measure of Muscle Thickness in Acute Stroke Patients, for Some, but Not All Anatomical Sites: A Study of the Intra-Rater Reliability of Muscle Thickness Measures in Acute Stroke Patients. Ultrasound Med Biol 2012, 38: 368–376. 10.1016/j.ultrasmedbio.2011.12.012
Jain AK, Farrokhnia F: Unsupervised texture segmentation using Gabor filters. Pattern recognition 1991, 24: 1167–1186. 10.1016/0031-3203(91)90143-S
Zhou Y, Zheng YP: Enhancement of muscle fibers in ultrasound images using Gabor filters. IEEE International Ultrasonics Symposium Proceedings 2009, 2296–2299. 101109/ULTSTYM.2009.5441467
Frangi A, Niessen W, Vincken K, Viergever M: Multiscale vessel enhancement filtering. Medical Image Computing and Computer-Assisted Interventation—MICCAI’98 1998, 1496: 130–7. 10.1007/BFb0056195
Rana M, Hamarneh G, Wakeling JM: Automated tracking of muscle fascicle orientation in B-mode ultrasound images. J Biomech 2009, 42: 2068–2073. 10.1016/j.jbiomech.2009.06.003
Hong L, Wan Y, Jain A: Fingerprint image enhancement: Algorithm and performance evaluation. Pattern Analysis and Machine Intelligence, IEEE Transactions on 1998, 20: 777–789. 10.1109/34.709565
Karanicolas PJ, Bhandari M, Kreder H, Moroni A, Richardson M, Walter SD, Norman GR, Guyatt GH: Evaluating agreement: conducting a reliability study. The Journal of Bone and Joint Surgery (American) 2009, 91: 99–106.
Rankin G, Stokes M: Reliability of assessment tools in rehabilitation: an illustration of appropriate statistical analyses. Clin Rehabil 1998, 12: 187–199. 10.1191/026921598672178340
Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986, 327: 307–310. 10.1016/S0140-6736(86)90837-8
McGraw KO, Wong S: Forming inferences about some intraclass correlation coefficients. Psychol Methods 1996, 1: 30.
Barger TS, Brown DE, Alwan M: Health-status monitoring through analysis of behavioral patterns. Systems, Man and Cybernetics, Part A: Systems and Humans, IEEE Transactions on 2005, 35: 22–27. 10.1109/TSMCA.2004.838474
This work is supported in part by grants from National Natural Science Foundation of China (NSFC: 81171402, 61103165), the next generation communication technology Major project of National S&T(2013ZX03005013), Guangdong Innovative Research Team Program (GIRTF-LCHT, No. 2011S013), Low-cost Healthcare Programs of Chinese Academy of Sciences and International Science and Technology Cooperation Program of Guangdong Province (2012B050200004) and Shenzhen Key Laboratory for Low-cost Healthcare (CXB201005260056A).
The authors declare that they have no competing interests.
PH: analysed the data and composed the manuscript together with YZ, YZ: proposed the idea. LO, YC, HL, GX and LW performed experiments, processed the data and made the discussions. 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.
About this article
Cite this article
Han, P., Chen, Y., Ao, L. et al. Automatic thickness estimation for skeletal muscle in ultrasonography: evaluation of two enhancement methods. BioMed Eng OnLine 12, 6 (2013). https://doi.org/10.1186/1475-925X-12-6
- Ultrasonography/Ultrasound imaging
- Image enhancement
- Gabor filtering
- Multiscale vessel enhancement filtering
- Muscle thickness