Automated tracking of muscle fascicle orientation in B-mode ultrasound images

B-mode ultrasound can be used to non-invasively image muscle fascicles during both static and dynamic contractions. Digitizing these muscle fascicles can be a timely and subjective process, and usually studies have used the images to determine the linear fascicle lengths. However, fascicle orientations can vary along each fascicle (curvature) and between fascicles. The purpose of this study was to develop and test two methods for automatically tracking fascicle orientation. Images were initially filtered using a multiscale vessel enhancement (a technique used to enhance tube-like structures), and then fascicle orientations quantified using either the Radon transform or wavelet analysis. Tests on synthetic images showed that these methods could identify fascicular orientation with errors of less than 0.06°. Manual digitization of muscle fascicles during a dynamic contraction resulted in a standard deviation of angle estimates of 1.41° across ten researchers. The Radon transform predicted fascicle orientations that were not significantly different from the manually digitized values, whilst the wavelet analysis resulted in angles that were 1.35° less, and reasons for these differences are discussed. The Radon transform can be used to identify the dominant fascicular orientation within an image, and thus used to estimate muscle fascicle lengths. The wavelet analysis additionally provides information on the local fascicle orientations and can be used to quantify fascicle curvatures and regional differences with fascicle orientation across an image.     

Quantitative diffusion tensor MRI-based fiber tracking of human skeletal muscle.

Diffusion-tensor magnetic resonance imaging (DT-MRI) offers great potential for understanding structure-function relationships in human skeletal muscles. The purposes of this study were to demonstrate the feasibility of using in vivo human DT-MRI fiber tracking data for making pennation angle measurements and to test the hypothesis that heterogeneity in the orientation of the tibialis anterior (TA) muscle's aponeurosis would lead to heterogeneity in pennation angle. Eight healthy subjects (5 male) were studied. T(1)-weighted anatomical MRI and DT-MRI data were acquired of the TA muscle. Fibers were tracked from the TA's aponeurosis by following the principal eigenvector. The orientations of the aponeurosis and muscle fiber tracts in the laboratory frame of reference and the orientation of the fiber tracts with respect to the aponeurosis [i.e., the pennation angle (theta)] were determined. The muscle fiber orientations, when expressed relative to the laboratory frame of reference, did not change as functions of superior-to-inferior position. The sagittal and coronal orientations of the aponeurosis did not change in practically significant manners either, but the aponeurosis' axial orientation changed by approximately 40 degrees . As a result, the mean value for theta decreased from 16.3 (SD 6.9) to 11.4 degrees (SD 5.0) along the muscle's superior-to-inferior direction. The mean value of theta was greater in the deep than in the superficial compartment. We conclude that pennation angle measurements of human muscle made using DT-MRI muscle fiber tracking are feasible and reveal that in the foot-head direction, there is heterogeneity in the pennation properties of the human TA muscle.     

Passive and dynamic muscle architecture during transverse loading for gastrocnemius medialis in man

External forces from our environment impose transverse loads on our muscles. Studies in rats have shown that transverse loads result in a decrease in the longitudinal muscle force. Changes in muscle architecture during contraction may contribute to the observed force decrease. The aim of this study was to quantify changes in pennation angle, fascicle dimensions, and muscle thickness during contraction under external transverse load.Electrical stimuli were elicited to evoke maximal force twitches in the right calf muscles of humans. Trials were conducted with transverse loads of 2, 4.5, and 10 kg. An ultrasound probe was placed on the medial gastrocnemius in line with the transverse load to quantify muscle characteristics during muscle twitches.Maximum twitch force decreased with increased transverse muscle loading. The 2, 4.5, and 10 kg of transverse load showed a 9, 13, and 16% decrease in longitudinal force, respectively. Within the field of view of the ultrasound images, and thus directly beneath the external load, loading of the muscle resulted in a decrease in the muscle thickness and pennation angle, with higher loads causing greater decreases. During twitches the muscle transiently increased in thickness and pennation angle, as did fascicle thickness. Higher transverse loads showed a reduced increase in muscle thickness. Smaller increases in pennation angle and fascicle thickness strain also occurred with higher transverse loads.This study shows that increased transverse loading caused a decrease in ankle moment, muscle thickness, and pennation angle, as well as transverse deformation of the fascicles.     

Muscle fascicle and series elastic element length changes along the length of the human gastrocnemius during walking and running

Ultrasound imaging has recently been used to distinguish the length changes of muscle fascicles from those of the whole muscle tendon complex during real life movements. The complicated three-dimensional architecture of pennate muscles can however cause heterogeneity in the length changes along the length of a muscle. Here we use ultrasonography to examine muscle fascicle length and pennation angle changes at proximal, distal and midbelly sites of the human gastrocnemius medialis (GM) muscle during walking (4.5 km/h) and running (7.5 km/h) on a treadmill. The results of this study have shown that muscle fascicles perform the same actions along the length of the human GM muscle during locomotion. However the distal fascicles tend to shorten more and act at greater pennation angles than the more proximal fascicles. Muscle fascicles acted relatively isometrically during the stance phase during walking, however during running the fascicles shortened throughout the stance phase, which corresponded to an increase in the strain of the series elastic elements (SEEs) (consisting of the Achilles tendon and aponeurosis). Measurement of the fascicle length changes at the midbelly level provided a good approximation of the average fascicle length changes across the length of the muscle. The compliance of the SEE allows the muscle fascicles to shorten at a much slower speed, more concomitant with their optimal speed for maximal power output and efficiency, with high velocity shortening during take off in both walking and running achieved by recoil of the SEE.     

Reliability and robustness of muscle architecture measurements obtained using diffusion tensor imaging with anatomically constrained tractography

For detailed analyses of muscle adaptation mechanisms during growth, ageing or disease, reliable measurements of muscle architecture are required. Diffusion tensor imaging (DTI) and DTI tractography have been used to reconstruct the architecture of human muscles in vivo. However, muscle architecture measurements reconstructed with conventional DTI techniques are often anatomically implausible because the reconstructed fascicles do not terminate on aponeuroses, as real muscle fascicles are known to do. In this study, we tested the reliability of an anatomically constrained DTI-based method for measuring three-dimensional muscle architecture. Anatomical magnetic resonance images and diffusion tensor images were obtained from the left legs of eight healthy participants on two occasions one week apart. Muscle volumes, fascicle lengths, pennation angles and fascicle curvatures were measured in the medial and lateral gastrocnemius, soleus and the tibialis anterior muscles. Averaged across muscles, the intraclass correlation coefficient was 0.99 for muscle volume, 0.81 for fascicle length, 0.73 for pennation angle and 0.76 for fascicle curvature. Measurements of muscle architecture obtained using conventional DTI tractography were highly sensitive to variations in the stopping criteria for DTI tractography. The application of anatomical constraints reduced this sensitivity significantly. This study demonstrates that anatomically constrained DTI tractography can provide reliable and robust three-dimensional measurements of whole-muscle architecture. The algorithms used to constrain tractography have been made publicly available.     

Reproducibility of gastrocnemius medialis muscle architecture during treadmill running

The purpose of this study was to assess the reproducibility of fascicle length (FL) and pennation angle (PA) of gastrocnemius medialis (GM) muscle during running in vivo. Twelve male recreational long distance runners (mean ± SD; age: 24 ± 3 years, mass: 76 ± 7 kg) ran on a treadmill at a speed of 3.0 m/s, wearing their own running shoes, for two different 10 min sessions that were at least 2 days apart. For each test day 10 acceptable trials were recorded. Ankle and knee joint angle data were recorded by a Vicon 624 system with three cameras operating at 120 Hz. B-mode ultrasonography was used to examine fascicle length and pennation angle of gastrocnemius medialis muscle. The ultrasound probe was firmly secured on the muscle belly using a lightweight foam fixation. The results indicated that fascicle length and pennation angle demonstrated high reproducibility values during treadmill running both for within and between test days. The root mean square scores between the repeated waveforms of pennation angle and fascicle length were small (∼2° and ∼3.5 mm, respectively). However, ∼14 trials for pennation angle and ∼9 trials for fascicle length may be required in order to record accurate data from muscle architecture parameters. In conclusion, ultrasound measurements may be highly reproducible during dynamic movements such as treadmill running, provided that a proper fixation is used in order to assure the constant location and orientation of the ultrasound probe throughout the movement.     

Determining skeletal muscle architecture with Laplacian simulations: a comparison with diffusion tensor imaging

Determination of skeletal muscle architecture is important for accurately modeling muscle behavior. Current methods for 3D muscle architecture determination can be costly and time-consuming, making them prohibitive for clinical or modeling applications. Computational approaches such as Laplacian flow simulations can estimate muscle fascicle orientation based on muscle shape and aponeurosis location. The accuracy of this approach is unknown, however, since it has not been validated against other standards for muscle architecture determination. In this study, muscle architectures from the Laplacian approach were compared to those determined from diffusion tensor imaging in eight adult medial gastrocnemius muscles. The datasets were subdivided into training and validation sets, and computational fluid dynamics software was used to conduct Laplacian simulations. In training sets, inputs of muscle geometry, aponeurosis location, and geometric flow guides resulted in good agreement between methods. Application of the method to validation sets showed no significant differences in pennation angle (mean difference \(0.5{^{\circ }})\) or fascicle length (mean difference 0.9 mm). Laplacian simulation was thus effective at predicting gastrocnemius muscle architectures in healthy volunteers using imaging-derived muscle shape and aponeurosis locations. This method may serve as a tool for determining muscle architecture in silico and as a complement to other approaches.     

Fascicle length of gastrocnemius muscles in monozygous twins

A large inter-individual variation is seen in muscle fascicle length of the athletes but the reasons for this phenomenon are unclear. The purpose of this study was to determine whether genetic factors contribute to the variances in muscle architectural characteristics. Nine monozygous twin pairs (3 males and 6 females), mean age 23 years (range 17-40) were studied. Fascicle length, pennation angle, and muscle thickness of the medial (MG) and lateral (LG) gastrocnemius muscles were measured in vivo by B-mode ultrasound. In the LG muscle intrapair resemblance (P < 0.01) for fascicle length (r = 0.98), pennation angle (r = 0.94) and muscle thickness (r = 0.86) were observed. In MG muscle, however, there was no intrapair resemblance for fascicle length (r = 0.66, P > 0.05), although pennation angle (r = 0.73, P < 0.05) and muscle thickness (r = 0.86, P < 0.01) were significant. Mean percent intrapair difference in LG and MG muscles were 1.8% and 5.1% for fascicle length, 11.3% and 12.3% for pennation angle and 12.4% and 9.9% for muscle thickness, respectively. There is intrapair difference between muscle thickness and pennation angle in both MG (r = 0.69, P < 0.05) and LG (r = 0.70, P < 0.05) muscles. However, no significant correlation was observed for intrapair difference between muscle thickness and fascicle length in both muscles (MG, r = 0.46; LG, r = 0.40). It appears that genetic predisposition is the predominant factor for the determination of muscle fascicle length. However, a lack of intrapair resemblance in MG fascicle length raises the possibility that fascicle length may be further influenced by external environmental factors such as physical training.     

The effect of ultrasound probe orientation on muscle architecture measurement.

The purpose of this paper was to examine how muscle architecture parameter (MAP) measurements made using brightness-mode ultrasonography (BMU) differ based on probe orientation. The human tibialis anterior muscle was imaged from nine different probe orientations during concentric contractions at four joint angles to determine the effect of probe orientation on the measurement of muscle architecture parameters. Ankle dorsi-flexion torque, tibialis anterior electromyography (EMG), and measures of MAP including fascicle length (FL), pennation angle (PA) and muscle thickness (MT) were collected. Statistically significant differences were found between joint angles for measures of FL and PA and between probe orientations for measures of FL and MT. A comparison of actual MAP values to a geometric model used by researchers to determine an ideal probe orientation show that the actual values vary to a greater extent and do not follow the trend predicted by the model. The results suggest that ultrasound probe orientation affects measures of MAP but the effect either cannot be predicted from a geometric model and/or the error in the measurement technique does not allow a comparison.     

Determination of fascicle length and pennation in a contracting human muscle in vivo

Fukunaga, Tetsuo, Yoshiho Ichinose, Masamitsu Ito, YasuoKawakami, and Senshi Fukashiro. Determination of fascicle lengthand pennation in a contracting human muscle in vivo.J. Appl. Physiol. 82(1): 354-358, 1997.We have developed a technique to determine fascicle length inhuman vastus lateralis muscle in vivo by using ultrasonography. Whenthe subjects had the knee fully extended passively from a position of110° flexion (relaxed condition), the fascicle length decreasedfrom 133 to 97 mm on average. During static contractions at 10% ofmaximal voluntary contraction strength (tensed condition), fascicleshortening was more pronounced (from 126 to 67 mm), especially when theknee was closer to full extension. Similarly, as the knee was extended, the angle of pennation (fascicle angle, defined as the angle between fascicles and aponeurosis) increased (relaxed, from 14 to 18°; tensed, from 14 to 21°), and a greater increase in the pennation angle was observed in the tensed than in the relaxed condition when theknee was close to extension (<40°). We conclude that there aredifferences in fascicle lengths and pennation angles when the muscle isin a relaxed and isometrically tensed conditions and that thedifferences are affected by joint angles, at least at thesubmaximal contraction level.

     

In vivo motion transmission in the inactive gastrocnemius medialis muscle-tendon unit during ankle and knee joint rotation.

The purposes of this study were: (a) to quantify the influence of passive ankle and knee joint angular displacement on the estimated mechanical and architectural properties of the gastrocnemius medialis (GM) muscle-tendon unit, and (b) to determine the strain distribution of separate structures (tendon, aponeurosis and fascicle) during passive lengthening of the GM muscle-tendon unit at rest. Ten male subjects participated in the study. The passive ankle and knee joint movements were performed on an isokinetic dynamometer. The kinematics of the left leg were recorded using the Vicon 624 system with 8 cameras. Two ultrasound probes were used to examine the elongation of the tendon, the aponeurosis, the fascicles and the angle of pennation of the GM. To calculate the elongation of the GM muscle-tendon unit the Achilles tendon path was reconstructed using a series of small reflective markers. The results show that the passive ankle joint angular displacement has a considerable influence on the elongation of the tendinous and architectural structures of the GM muscle-tendon unit. In contrast, the influence of knee joint angular displacement on the GM fascicle length and pennation angle becomes relevant only at knee angles greater than 144 degrees . The contribution of the tendon to the elongation of the GM muscle-tendon unit at rest is relevant because of its greater resting length in comparison to the resting length of the GM fascicles. The results indicate the existence of slackness in the inactive GM muscle-tendon unit between 121 degrees and 107 degrees ankle angle and between 65 degrees and 144 degrees knee angle.     

Ultrasound-based subject-specific parameters improve fascicle behaviour estimation in Hill-type muscle model

The estimation of muscle fascicle behaviour is decisive in a Hill-type model as they are related to muscle force by the force–length–velocity relationship and the tendon force–strain relationship. This study was aimed at investigating the influence of subject-specific tendon force–strain relationship and initial fascicle geometry (IFG) on the estimation of muscle forces and fascicle behaviour during isometric contractions. Ultrasonography was used to estimate the in vivo muscle fascicle behaviour and compare the muscle fascicle length and pennation angle estimated from the Hill-type model. The calibration–prediction process of the electromyography-driven model was performed using generic or subject-specific tendon definition with or without IFG as constraint. The combination of subject-specific tendon definition and IFG led to muscle fascicle behaviour closer to ultrasound data and significant lower forces of the ankle dorsiflexor and plantarflexor muscles compared to the other conditions. Thus, subject-specific ultrasound measurements improve the accuracy of Hill-type models on muscle fascicle behaviour.     

Validity of architectural properties of the hamstring muscles: Correlation of ultrasound findings with cadaveric dissection

The purpose of this study was to compare the architectural parameters of the long head of biceps femoris (BFlh) and semitendinosus (ST) muscles by comparing measurements from ultrasound (US) with those obtained from direct dissection. The BFlh and ST architectures were examined bilaterally in 6 legs from 3 male cadavers. The fascicle length, pennation angle, muscle thickness and muscle and tendon length were obtained from direct measurement and US scans along each muscle. Intraclass correlation coefficients between the two methods ranged from 0.905 to 0.913 for the BFlh variables and from 0.774 to 0.974 for the ST parameters. Compared with the direct measurements, the US method showed a mean typical error of 0.09–0.14 cm for muscle thickness, 1.01–1.31° for the pennation angle, 0.92–1.71 cm for fascicle length and muscle–tendon length measurements. The US method is a valid alternative tool for assessing basic architectural parameters of ST and BFlh components of the hamstring muscles.     

Behavior of human muscle fascicles during shortening and lengthening contractions in vivo.

The aim of the present study was to investigate the behavior of human muscle fascicles during dynamic contractions. Eight subjects performed maximal isometric dorsiflexion contractions at six ankle joint angles and maximal isokinetic concentric and eccentric contractions at five angular velocities. Tibialis anterior muscle architecture was measured in vivo by use of B-mode ultrasonography. During maximal isometric contraction, fascicle length was shorter and pennation angle larger compared with values at rest (P < 0.01). During isokinetic concentric contractions from 0 to 4.36 rad/s, fascicle length measured at a constant ankle joint angle increased curvilinearly from 49.5 to 69.7 mm (41%; P < 0.01), whereas pennation angle decreased curvilinearly from 14.8 to 9.8 degrees (34%; P < 0.01). During eccentric muscle actions, fascicles contracted quasi-isometrically, independent of angular velocity. The behavior of muscle fascicles during shortening contractions was believed to reflect the degree of stretch applied to the series elastic component, which decreases with increasing contraction velocity. The quasi-isometric behavior of fascicles during eccentric muscle actions suggests that the series elastic component acts as a mechanical buffer during active lengthening.     

Computational methods for quantifying in vivo muscle fascicle curvature from ultrasound images

Muscle fascicles curve during contraction, and this has been seen using B-mode ultrasound. Curvature can vary along a fascicle, and amongst the fascicles within a muscle. The purpose of this study was to develop an automated method for quantifying curvature across the entirety of an imaged muscle, to test the accuracy of the method against synthetic images of known curvature and noise, and to test the sensitivity of the method to ultrasound probe placement. Both synthetic and ultrasound images were processed using multiscale vessel enhancement filtering to accentuate the muscle fascicles, wavelet-based methods were used to quantify fascicle orientations and curvature distribution grids were produced by quantifying local curvatures for each point within the image. Ultrasound images of ramped isometric contractions of the human medial gastrocnemius were acquired in a test–retest study.The methods enabled distinct curvatures to be determined in different regions of the muscle. The methods were sensitive to kernel sizes during image processing, noise within the image and the variability of probe placements during retesting. Across the physiological range of curvatures and noise, curvatures calculated from validation grids were quantified with a typical standard error of less than 0.026 m^?1, and this is about 1% of the maximum curvatures observed in fascicles of contracting muscle.     

Changes of calf muscle-tendon biomechanical properties induced by passive-stretching and active-movement training in children with cerebral palsy

Biomechanical properties of calf muscles and Achilles tendon may be altered considerably in children with cerebral palsy (CP), contributing to childhood disability. It is unclear how muscle fascicles and tendon respond to rehabilitation and contribute to improvement of ankle-joint properties. Biomechanical properties of the calf muscle fascicles of both gastrocnemius medialis (GM) and soleus (SOL), including the fascicle length and pennation angle in seven children with CP, were evaluated using ultrasonography combined with biomechanical measurements before and after a 6-wk treatment of passive-stretching and active-movement training. The passive force contributions from the GM and SOL muscles were separated using flexed and extended knee positions, and fascicular stiffness was calculated based on the fascicular force-length relation. Biomechanical properties of the Achilles tendon, including resting length, cross-sectional area, and stiffness, were also evaluated. The 6-wk training induced elongation of muscle fascicles (SOL: 8%, P = 0.018; GM: 3%, P = 0.018), reduced pennation angle (SOL: 10%, P = 0.028; GM: 5%, P = 0.028), reduced fascicular stiffness (SOL: 17%, P = 0.128; GM: 21%, P = 0.018), decreased tendon length (6%, P = 0.018), increased Achilles tendon stiffness (32%, P = 0.018), and increased Young's modulus (20%, P = 0.018). In vivo characterizations of calf muscles and Achilles tendon mechanical properties help us better understand treatment-induced changes of calf muscle-tendon and facilitate development of more effective treatments.     

Myofascial force transmission between the human soleus and gastrocnemius muscles during passive knee motion

The plantarflexors of the lower limb are often assumed to act as independent actuators, but the validity of this assumption is the subject of considerable debate. This study aims to determine the degree to which passive changes in gastrocnemius muscle length, induced by knee motion, affect the tension in the adjacent soleus muscle. A second aim is to quantify the magnitude of myofascial passive force transmission between gastrocnemius and adjacent soleus. Fifteen healthy volunteers participated. Simultaneous ultrasound images of the gastrocnemius and soleus muscles were obtained during passive knee flexion (0-90°), while keeping the ankle angle fixed at either 70° or 115°. Image correlation analysis was used to quantify muscle fascicle lengths in both muscles. The data show that the soleus muscle fascicles elongate significantly during gastrocnemius shortening. The approximate change in passive soleus force as a result of the observed change in fascicle length was estimated and appears to be <5 N, but this estimate is sensitive to the assumed slack length of soleus.     

In vivo determination of fascicle curvature in contracting human skeletal muscles.

Fascicle curvature of human medial gastrocnemius muscle (MG) was determined in vivo by ultrasonography during isometric contractions at three (distal, central, and proximal) locations (n = 7) and at three ankle angles (n = 7). The curvature significantly (P < 0.05) increased from rest to maximum voluntary contraction (MVC) (0.4-5.2 m(-1)). In addition, the curvature at MVC became larger in the order dorsiflexed, neutral, plantar flexed (P < 0.05). Thus both contraction levels and muscle length affected the curvature. Intramuscular differences in neither the curvature nor the fascicle length were found. The direction of curving was consistent along the muscle: fascicles were concave in the proximal side. Fascicle length estimated from the pennation angle and muscle thickness, under the assumption that the fascicle was straight, was underestimated by ~6%. In addition, the curvature was significantly correlated to pennation angle and muscle thickness. These findings are particularly important for understanding the mechanical functions of human skeletal muscle in vivo.     

Comparison of fascicle behaviors between superficial and deeper muscles of triceps brachii during isometric contractions

PurposeWe assessed fascicle behaviors of the upper extremities during isometric contractions at different joint angles in this study.MethodsThirteen healthy men and women performed isometric elbow extension tasks at 50% and 75% of maximal voluntary contraction (MVC) at 60°, 90°, and 120° of elbow extension (full extension = 180°). Extended field-of-view B-mode ultrasonography was used to obtain sagittal plane panoramic images of the long head (TB-Long) and medial head (TB-Med) of the triceps brachii at rest and during contraction; fascicle length and pennation angle were measured.ResultsIn the TB-Long, significant fascicle shortening from rest was found during 50% and 75%MVC at 60° and during 75%MVC at 90° of extension. There was no significant fascicle shortening in the TB-Med muscle under any conditions. There was no significant pennation angle change from rest in either muscle. The pennation angle of the TB-Long was significantly greater than that of the TB-Med under all conditions.ConclusionsThese results suggest that fascicle shortening in the TB-Long muscle occurs in flexion; however, no change was found in the TB-Med. In the upper extremity muscle–tendon complex, the superficial and deeper muscles may have different force-transmission efficiency at flexed joint angles.     

In vivo fascicle length measurements via B-mode ultrasound imaging with single vs dual transducer arrangements

Ultrasonography is a useful technique to study muscle contractions in vivo, however larger muscles like vastus lateralis may be difficult to visualise with smaller, commonly used transducers. Fascicle length is often estimated using linear trigonometry to extrapolate fascicle length to regions where the fascicle is not visible. However, this approach has not been compared to measurements made with a larger field of view for dynamic muscle contractions. Here we compared two different single-transducer extrapolation methods to measure VL muscle fascicle length to a direct measurement made using two synchronised, in-series transducers. The first method used pennation angle and muscle thickness to extrapolate fascicle length outside the image (extrapolate method). The second method determined fascicle length based on the extrapolated intercept between a fascicle and the aponeurosis (intercept method). Nine participants performed maximal effort, isometric, knee extension contractions on a dynamometer at 10° increments from 50 to 100° of knee flexion. Fascicle length and torque were simultaneously recorded for offline analysis. The dual transducer method showed similar patterns of fascicle length change (overall mean coefficient of multiple correlation was 0.76 and 0.71 compared to extrapolate and intercept methods respectively), but reached different absolute lengths during the contractions. This had the effect of producing force–length curves of the same shape, but each curve was shifted in terms of absolute length. We concluded that dual transducers are beneficial for studies that examine absolute fascicle lengths, whereas either of the single transducer methods may produce similar results for normalised length changes, and repeated measures experimental designs.