Modelling human musculoskeletal functional movements using ultrasound imaging
- Michael Peolsson†3,
- Tommy Löfstedt†1, 2,
- Susanna Vogt1, 2,
- Hans Stenlund1, 2,
- Anton Arndt4, 5 and
- Johan Trygg1, 2Email author
© Peolsson et al; licensee BioMed Central Ltd. 2010
Received: 20 July 2009
Accepted: 21 May 2010
Published: 21 May 2010
A widespread and fundamental assumption in the health sciences is that muscle functions are related to a wide variety of conditions, for example pain, ischemic and neurological disorder, exercise and injury. It is therefore highly desirable to study musculoskeletal contributions in clinical applications such as the treatment of muscle injuries, post-surgery evaluations, monitoring of progressive degeneration in neuromuscular disorders, and so on.
The spatial image resolution in ultrasound systems has improved tremendously in the last few years and nowadays provides detailed information about tissue characteristics. It is now possible to study skeletal muscles in real-time during activity.
The ultrasound images are transformed to be congruent and are effectively compressed and stacked in order to be analysed with multivariate techniques. The method is applied to a relevant clinical orthopaedic research field, namely to describe the dynamics in the Achilles tendon and the calf during real-time movements.
This study introduces a novel method to medical applications that can be used to examine ultrasound image sequences and to detect, visualise and quantify skeletal muscle dynamics and functions.
This new objective method is a powerful tool to use when visualising tissue activity and dynamics of musculoskeletal ultrasound registrations.
Image analysis has been used for decades as a tool to analyse different kinds of tissues. It provides a rich resource when presenting medical and clinical information and the plethora of techniques using images is increasing rapidly.
Various pathologic conditions benefit from medical imaging of soft tissues in terms of detection, progression, remission (effects of treatment), etc. Tumours, for example, share important similarities with skeletal muscles. They are both inherently smooth masses with dynamic characteristics. In the tumour case, a change in the tissue's dynamic quality is an important indicator of diagnostic evaluations and is therefore of great clinical value. For instance, decreased tissue dynamics in a muscle can, in the musculoskeletal case, be an indication of inflammation  or ischemia . Muscle function can also be related to a wide variety of conditions, for example pain, ischemic and neurological disorders, exercise and injury. There is thus an urging need for imaging tools that can be used to study the dynamics of skeletal muscles, and especially to do so in real-time.
The quality of ultrasound imaging has evolved quickly and is today capable of producing high spatial resolution image sequences with a high frame rate. This makes ultrasound imaging suitable for analysis of both structural and functional aspects of muscles. Examples include tendon injuries in sports medicine  and inflammation processes [4, 5].
Speckle tracking is a technique that can be applied to ultrasound sequences. Speckle tracking relies on the scattered echoes that arise from inhomogeneities in the tissues. Speckle tracking provides studying tissue motions in real-time by following a specific pattern frame-by-frame in a sequence of images . The speckle tracking can be described in terms of e.g. displacement maps [7–11]. Our studies here are focused mainly on transverse displacements, but are performing 2D measurement; a recent study has shown that measurements of 2D displacements can be accurately performed . Skeletal muscle dynamics and muscle fatigue have also been studied by using speckle tracking [13, 14].
The aim of our study is to develop a complementary method that can be used separately for analysis of skeletal muscles, or used together with current and established methods—such as speckle tracking—to improve model precision. The method that we present here is based on two parts: the discrete wavelet transform and multivariate data analysis. The wavelet transform is able to extract position, size and shape information that is present in greyscale B-scan ultrasound image sequences of muscle tissues movements. Not only are single pixel intensities considered, but regions of pixel intensities are considered and extracted as a whole. Secondly, multivariate analysis is applied to the wavelet transformed images in order to get a global model of the captured events. Multivariate analysis allows a multitude of different tools to be used on the data, e.g. comparisons, cluster analysis, discrimination analysis, etc; and helps to detect, visualise and quantify the skeletal muscle dynamics and functionality. Multivariate methods can also analyse speckle tracking data, as we will demonstrate in the given examples.
The intensity value changes within an image represent movements directly observable (the effect of noise is heavily reduced by using the MACI method) and which can for example, in the speckle tracking case, be directly and actively followed, "observed", around within an image in an image sequence. Using the discrete wavelet transform means instead that these intensity value changes, which represent actual muscle tissue movements, are indirectly and passively observed as they pass by the "observer" (the variable, the wavelet coefficient), meaning a particular muscle tissue movement is not represented by any single variable, but by a multitude of variables simultaneously and together describing a particular phenomenon within the sequence.
The MACI method is applied to describe the dynamics in the Achilles tendon and the calf during real-time heel raise movements. It will be argued that this technique is a valuable tool to use when studying and visualising muscle tissue dynamics.
The strategic significance of a method like MACI is very high since it could be applied to muscle rehabilitation programmes (including sports medicine), neurological disorders (e.g. whiplash and fibromyalgia) and pain related conditions (e.g. back pain).
Modelling tissue dynamics
The basis for image analysis is that an M × N two-dimensional digital intensity image can be seen as a function, I(x, y), with an intensity value for every point, (x, y), of the image .
When several such 2D images are stacked upon each other, like in an ultrasound image sequence, a multivariate image space occurs. In such a multivariate image space, each location is represented by a multitude of values and a massive amount of highly correlated data is the result. Different unfolding procedures can then be used in order to analyse these 3D data structures [16–18].
The feature extraction procedure used here is based on the discrete 2D wavelet transform (2D-DWT). The 2D-DWT is beneficial to use here since when using it to transform and compress large sets of images, the important position, size and shape information in the images is not lost. The images are compressed by choosing the wavelet coefficients that hold the most information about the ultrasound sequence. This is done by measuring the variance of each coefficient. The selected wavelet coefficients are then put in a regular two-way data table where each row represents an image and each column contains the wavelet coefficients for that image. Principal component analysis (PCA) is then used to analyse the two-way data table of wavelet coefficients.
In the examples below, movements in the Achilles tendon and the calf were analysed by multivariate analysis of congruent images (MACI). The results were then compared to an analysis of the same movements in the corresponding areas using state of the art speckle tracking. Specifically, the absolute transversal displacement was measured using speckle tracking and the results were compared to the results of the MACI method. Also, a grid of square speckle tracking ROIs was placed equidistantly in the whole image plane in the third example and a PCA of the resulting tracks was made. The purpose of using these speckle tracking measures was to show that the MACI technique indeed does capture the actual movements in the sequence.
The investigated types of movements are: standardised passive movement, standardised active movement and functional movement.
The MACI method was used in the example cases below to extract information from ultrasound image sequences. The resulting two-way data table was examined with PCA and the movements verified by speckle tracking. This section explains speckle tracking, PCA and the concept of MACI.
The acoustic patterns in an ultrasound signal change when the muscle being scanned moves. These acoustic markers, or speckle patterns, remain relatively stable over time and can therefore be followed frame-by-frame in a sequence of images . A commercial software package (EchoPac, GE Healthcare, Horten, Norway) was used together with the ultrasonic system Vivid 7 (GE Healthcare, Horten, Norway) in all case examples in this study. An in-house developed speckle tracking software was also used in the third example.
Where I is the image intensity, x and y are the pixel coordinates at time, or frame, t and w is a weighting function, which can be 1 in the simplest case.
The speckles are not always followed correctly when there are rapid alterations in the tissue, when there are large deformations in the tissue or when there are out of plane motions present. However, current research indicates that speckle tracking indeed is well suited for skeletal muscle investigations [20, 21].
The speckle tracking algorithm used in the in-house developed speckle tracking software is the Lucas-Kanade optical flow method  using pyramids  (hierarchical version) as implemented in the open source computer vision library (OpenCV) version 1.1 http://opencv.willowgarage.com/wiki/.
Some more information regarding greyscale B-scan ultrasound and speckle tracking is given in the supplemental material in Additional file 1.
Principal component analysis
Principal component analysis (PCA) is commonly used to identify patterns in multi- and megavariate data sets .
PCA is one of the cornerstones of multivariate data analysis. It is used to extract systematic variation in a dataset by projecting the data to a lower-dimensional space in which they are more easily analysed. PCA achieves this by finding the best least squares approximating hyperplane in this lower-dimensional space.
where k is the number of latent variables, T contains the scores, P contains the loadings and E is a matrix of residuals. X is the matrix of data being analysed.
The scores give information about the relation between the observations, for example the variation between the frames of an ultrasound sequence. Plotting score vectors against each other results in a score plot that gives information about trends, clusters and outliers in the data.
The loadings are seen as a measure of the importance of the variables. The loadings can be used to reconstruct the ultrasound sequence, resulting in a loading image. This image gives information about what parts of the frames are important for describing the movement.
PCA was performed using SIMCA-P+ 11.5 (Umetrics, Umeå, Sweden).
Multivariate analysis of congruent images
Multivariate analysis of congruent images (MACI) is used to find and express patterns over multivariate image spaces for the purpose of classification of, or finding relationships between, the images . The main goals of multivariate image analysis are: firstly, to compress the highly correlated data into terms of a few linear combinations of the intensity values; and secondly, to preserve the spatial information in the images.
By congruent is meant a set of images that are properly pre-processed, i.e. transformed, so that each image element corresponds to the same element across all images. When the images are not fully congruent initially, as they seldom are, they are made so by some means. In this case the discrete 2D wavelet transform (2D-DWT) of the images was used to make them congruent. The 2D-DWT was used to extract spatial (position) and frequency (shape and size) based features from the images. The wavelet basis Symlet 8 was used in the transforms.
The wavelet coefficients for each image are put into the rows of an ordinary two-way data table on which conventional multivariate methods, such as PCA, can be used to extract the information. In such methods, each variable must describe the same phenomenon for each observation; this is solved by using the congruent wavelet coefficients.
The development of the wavelet texture analysis-based methods have, in recent years, opened the door to this new field of multivariate data analysis .
Three case examples will be described in this manuscript, and they will be used to illustrate the usefulness of our method. The first example is of a standardised passive ankle movement; the second of a standardised active ankle movement; and the third of a functional ankle movement. The reason for these choices of case examples is that treatment and rehabilitation of patients with ruptured Achilles tendon surgery covers the described chosen movement types in the post-surgery rehabilitation phase.
The movement examined in the examples is dorsal-plantar flexion. The reason for this choice of movement is that it can be standardised by using an isokinetic dynamometer (Isomed 2000, D&R Ferstl GmbH, Henau, Germany). The isokinetic dynamometer controlled the muscular activity (EMG and force) and the range of motion (ROM) of the ankle. The EMG and dynamometer could, for instance, confirm that there was no active muscle contribution present in the passive ankle movement.
The passive dynamic and active isometric movements were performed using the isokinetic dynamometer, with the subject laying face down with the ankle fixated onto the dynamometer platform. The functional movement was performed with the subject in a standing position, holding a hand towards the wall in order to balance the body. Two heel raises were performed from this position using one leg at the time.
A commercial speckle tracking post-process software package from GE Medical (Echopac, GE Healthcare, Horten, Norway) was applied to the ultrasound loops in order to provide a reference analysis of the different phases of the ankle movements. An ROI was specified in a vertical position in the tendon, 2 cm proximal to the malleolus. The ROI was specified in the first frame and the speckle tracking method was applied in order to capture the tendon movements in the following frames in a frame-by-frame approach. The ROI was seen to move in accordance with the tissue during the movements. The absolute transversal displacement of the movements was measured, the distance (mm) that each segment of the ROI moved. Information regarding how to interpret the speckle tracking displacement plots is given in the supplemental material, in Additional file 1.
Rectangular ROIs were also specified equidistantly in a 100 × 100 grid over the entire image plane in the third example. This was done using an in-house developed speckle tracking software. The reason for this is to provide a more similar reference to the whole-image analysis that the MACI method performs on the image sequence in the third example.
EMG was used in the active foot ankle movement to confirm baseline rest (not shown).
Case study 1: Passive foot ankle movement
The objective here was to analyse the tissue response during a non-active dorsal-plantar flexion movement. This movement is the first type of movement in the clinical rehabilitation phase.
The dorsal-plantar flexions were provided by fixating the foot onto a platform that was connected to an isokinetic dynamometer. The ankle angle was varied at a constant speed of 30 degrees per second between 20 degrees of dorsal flexion and 20 degrees of plantar flexion, forming the range of motion (ROM).
A second objective was to study the repeatability of the method when repeating the dorsal-plantar flexions; therefore two subsequent ROMs were performed.
Case study 2: Active foot ankle movement
The objective in this case was to analyse a movement where the subject activates the calf muscles voluntarily. In order to do so, a standardised active muscle contraction was performed. The reason for this choice was that voluntary muscle activation is the second phase in the rehabilitation program after surgery. In this case, the subject was asked to perform a maximum plantar flexion towards the dynamometer platform (maximum voluntary contraction, MVC). The foot angle in this test was set to 0 degrees.
Case study 3: functional movement
The third type of ankle movement studied was a functional movement, in which the ankle was not fixated onto an isokinetic apparatus. The reason for this choice was that functional movements are central in clinical rehabilitation programs in order to regain functionality and smoothness in movements and to prevent secondary lesions. Stabilisation training, by way of balancing, strength training and relearning to use injured muscle structures, is important in rehabilitation.
In this example, the subject performed three subsequent 30 degrees heel raises (plantar flexions) in a standing position. The reason for the choice of three repetitions was to see if it was possible to get an indication of the methodological repeatability when analysing the ultrasound images using the two image analysis techniques.
A speckle tracking measurement of the pennation angle in the Soleus muscle during the movement was performed in order to provide a measurement of the muscle deformation. The pennation angle was measured as the angle of fascicular insertion into the deep fascia of the Achilles tendon. The pennation angles were measured using the in-house developed software and each measure was found as the angle between three points. The pennation angle measurements were confirmed using visual confirmation [27, 28].
All experiments were performed on one of the authors of this manuscript. A copy of the written consent of this author for the publication of this case report and any accompanying images is available for review by the Editor-in-Chief of this journal.
Case study 1: Passive foot ankle movement
The output from the commercial ultrasound software is shown in Figure 3(a). The original greyscale loop is presented in the upper left part of the image and an ROI is specified vertically within the Achilles tendon. The right part of this and subsequent such images presents the displacement as a function of time for each of the coloured segments of the ROI.
Figure 3(a) presents the speckle tracking result from the tissue movements during two subsequent dorsal flexions. The turning phases of the dorsal and plantar movements are clearly visible.
In Figure 3(b) and Figure 3(c) the MACI result is presented for a sub image covering the Achilles tendon of the same image sequence as in Figure 3(a). PCA was performed on the wavelet compressed images and the result is presented as a score plot and as a distance to model plot, DModX—which is equal to the residual standard deviation . The score plot for this MACI model is presented in Figure 3(b) and the DModX plot in Figure 3(c). The movement phases are highlighted in Figure 3(a) and the corresponding ones in Figure 3(b) and Figure 3(c).
The score plot in Figure 3(b) shows the relationship between the frames of the loop. The score plot shows that each movement starts and ends near to the origin. The first principal component separates the plantar flexions from the dorsal flexions and the second principal component indicates a difference between the two dorsal flexions. The difference between the dorsal flexions is probably due to a slightly unrelaxed muscle in these phases of the movement. The turning phases are easily identified since the observations, the frames, change direction in the score plot.
Critical points in the case example movements
Turn, dorsal flexion 1:
Turn, plantar flexion 1:
Turn, dorsal flexion 2:
Turn, plantar flexion 2:
Static contraction phase:
Turn, repetition 1:
Turn, repetition 2:
Turn, repetition 3:
Overall, the sequence was modelled well—the two score vectors, created from an initial set of 3000 variables, account for 13% of all variation in the compressed sequence and the observations are within the 5% level in DModX—and all movement phases can be identified and correlated. Similar movement pattern appears in the two subsequent repetitions, both when using speckle tracking and when using the MACI method.
The Euclidean distance between adjacent pairs of images in the PCA score plot in Figure 3(b) reveals important information. Images in the start and end zones and in the turning phase areas are tethered, implying that they are similar. Between the start and end zones and the turning phase areas there are larger distances between the images, implying greater differences. Thus, larger differences are seen where the muscle tissue is either compressed or stretched.
Case study 2: Active foot ankle movement
Figure 4(a) presents the speckle tracking reference analysis result of the movement, the displacement as a function of time for each of the coloured segments of the ROI. The first two phases—active concentric phase (isokinetic phase), active static contraction (isometric phase)—are clearly visible in the displacement plot. The last phase, the active eccentric phase (isokinetic phase), cannot be plotted together with the two first phases because of a time limitation in the commercial ultrasound software, but the displacement plot does return close to the baseline after the last phase of the movement. This is confirmed if the contraction phase is cropped and only the stationary and eccentric phases are examined, or if the in-house developed speckle tracking software is used.
The MACI result of this case study is presented as a PCA score plot in Figure 4(b). The MACI method was applied to a sub image of the same sequence, a sub image covering the Achilles tendon, just like in the previous case example. All three phases can be identified in the score plot and correlated to the displacement result from speckle tracking, as is also seen in Table 1. This PCA model managed to model the movement very well; the two score vectors, created from an initial set of 3000 variables, that was presented here account for 28% of all variation in the compressed sequence.
The differences in the Euclidian distances between the phases in Figure 4(b) indicate a difference between two kinds of contractions: dynamic isokinetic and static isometric, respectively. In the isometric MVC phase, in the second PCA component direction, it can be seen that clusters of images are separated by images farther apart. Hence, the tissue portrayed in these images both contains images showing hardly any variation (the clusters) and images that are separated by some distance, which suggests dynamic contributions in the tissue in the isometric phase as well. This is a likely result because of the difficulty of maintaining a constant force during the MVC, and MACI thus actually confirms research already made by for example . The distances between the images are larger in the dynamic isokinetic phases, implying more variation, which indicates more dynamics in the tissue during these phases.
Case study 3: Functional movement
The results of the functional movement analysis are presented in Figure 5. In this case, the speckle tracking ROI was placed in the vertical direction crossing all tissue layers present in the image, instead of just within the Achilles tendon as in the two previous case examples. The right part of Figure 5(a) shows clearly a rhythmic pattern when studying direction, magnitude and internal relationships between investigated segments in the calf.
A rhythmic pattern can also be seen in the PCA score plot in Figure 5(b) for the three sequential movements. All repetitions start in the right part of the score plot and move to the left, they turn and return to the right, close to the start area. The start zone, plantar flexion, turning zone, dorsal flexion and end zones can all be identified and correlated to the movements identified by the speckle tracking method. This can be seen in Figure 5(b) and in Table 1.
The first PCA component in Figure 5(b) separates the plantar flexions from the dorsal flexions. The second component describes inter-repetition variation and the isokinetic turning phase of the plantar flexions—we will return to this in a moment. Thus, the MACI method provides information about both progression characteristics and different phases of the movements as well as information about variations between the three repetitions. It is worth noting that each component carries unique information about the movements and could thus be investigated further.
A comparison of the three repetitions reveals a strong resemblance between them. The PCA scores corresponding to the first and the second repetition have a correlation of 0.64, the first and the third repetition a correlation of 0.72 and the second and third repetition have a correlation of 0.90. The first repetition has a slightly shorter duration than the other two, giving rise to the lower correlation with them.
It is clear that Figure 6(c) is very similar to Figure 5(b); the two methods have thus captured the same properties of the movements. The first score vector for the three repetitions have almost the same correlation between the repetitions as the MACI method had. However, for the second component, one can see that the start and the end zones of each repetition are much farther apart when using the speckle tracking method than when using the MACI method; something that it is reasonable to believe that they shouldn't be. The reason for this could be due to drift in the speckle tracking, or that the MACI method is better at capturing the isokinetic movement phase found in the second component in the left part of Figure 5(b), before and after the plantar flexion turn.
The frames at which the isometric phase begins and ends are visible in the pennation angle plot, Figure 7(b) and Figure 7(c), as frames at which the angle derivative is slightly smaller. The start and end frames of the isometric phase have been marked in Figure 7(b) and Figure 7(c) for the three repetitions. When the pennation angle is plotted together with the first and second principal component it is clear that the concentric and eccentric contraction phases are mainly caught by the first component and that the isokinetic phases are mainly caught by the second component. This is particularly apparent in the combined two-component score plot in Figure 5(b), where the isokinetic phase corresponds to the left-most, upwards-downwards moving phases of the score trajectory.
Ultrasound images of muscles and tendons captured in real-time during activity carry information about tissue movements and deformations in several different parts of the image. Since different muscles cooperate when performing different movements, changes within the images are usually not independent, but are coordinated together. The MACI method exploits these dependent coordination patterns naturally, and uses the whole array of stacked images when extracting its information. It was confirmed by speckle tracking that the MACI method can accurately identify the movement patterns that are present in ultrasound sequences, which is encouraging.
The MACI method can in fact also be used to find regions of interest, and not only using them. This is done by utilising the loading vectors, and reconstructing what we call loading images . Such images give information of what areas where important for building for example a PCA model, and can thus be used to see which areas are active in different phases of a movement. However, this methodology has not yet been properly developed to be used in muscle tissue application. This is some of the future work that will be performed in order to improve the MACI method.
The MACI approach also allows comparative studies to be made between patients. Different groups of patients, for example healthy and control, could be modelled individually or together and compared for similarities and differences. That has not been done in this study, however, but is one of the aims of our future studies.
The MACI method examines image sequences inductively, without being told what to look for. Instead, a data driven method extracts the information that is the most important for describing a particular sequence. It is thus not the analyst who decides what parts of the image contribute the most information, but the method that does so. This unbiased data exploration approach therefore has great potential in clinical situations.
The MACI method was able to correctly capture each passive reference movement and to describe the different phases of the ROM. It also captured the tissue dynamics during the different phases, partly in terms of the directions of the frames in the score plot, and partly in terms of distance relationships between images during the various movement phases. Also, both standardised and functional movements could be analysed with the MACI method.
To conclude, all sequences were modelled well with the MACI method and the different phases of the movements (turning phases, dorsal and plantar flexions) could all be identified and correlated in the three types of ankle movements. The same kind of pattern also arose when repeated flexions were compared; the repeatability provided by the MACI method is thus very encouraging.
Applying multivariate techniques based on whole image sequences is thus a powerful way to study variations and relationships between images in an ultrasound image sequence.
We argue that a multitude of applications are possible. The functionality of healthy muscles, injured muscles or training could be a focus. But also neuro-motor related diseases, such as multiple sclerosis, could be a possible area of application. Future studies will be made that will give more insight into what application areas are possible.
We would like to thank both reviewers (Antoine Nordez and Russell Witte) for their very constructive and detailed comments that significantly contributed to improve this paper.
The authors would like to acknowledge financial support from MKS Umetrics AB (TL), the Swedish Research Council (JT) and the Swedish Foundation for Strategic Research (JT).
- Kuo GP, Carrino JA: Skeletal muscle imaging and inflammatory myopathies. Current opinion in rheumatology. 2007, 19: 530-535. 10.1097/BOR.0b013e3282efdc66.View ArticlePubMedGoogle Scholar
- Elvin A, Siösteen AK, Nilsson A, Kosek E: Decreased muscle blood flow in fibromyalgia patients during standardised muscle exercise: A contrast media enhanced colour doppler study. European Journal of Pain. 2006, 10: 137-144. 10.1016/j.ejpain.2005.02.001.View ArticlePubMedGoogle Scholar
- Nørregaard J, Larsen CC, Bieler T, Langberg H: Eccentric exercise in treatment of Achilles tendinopathy. Scandinavian journal of medicine & science in sports. 2007, 17: 133-138.Google Scholar
- Grassi W, Filippucci E, Busilacchi P: Musculoskeletal ultrasound. Best Practice & Research Clinical Rheumatology. 2004, 18: 813-826.View ArticlePubMedGoogle Scholar
- Kane D, Balint PV, Sturrock R, Grassi W: Musculoskeletal ultrasound - a state of the art review in rheumatology. Part 1: Current controversies and issues in the development of musculoskeletal ultrasound in rheumatology. Rheumatology. 2004, 43: 823-828. 10.1093/rheumatology/keh214.View ArticlePubMedGoogle Scholar
- Bohs LN, Trahey GE: A novel method for angle independent ultrasonic imaging of blood flow and tissue motion. IEEE Transactions on Biomedical Engineering. 1991, 38: 280-286. 10.1109/10.133210.View ArticlePubMedGoogle Scholar
- Deffieux T, Gennisson JL, Tanter M, Fink M: Assessment of the mechanical properties of the musculoskeletal system using 2-D and 3-D very high frame rate ultrasound. IEEE Trans Ultrason Ferroelectr Freq Control. 2008, 55: 2177-2190. 10.1109/TUFFC.917.View ArticlePubMedGoogle Scholar
- Deffieux T, Gennisson JL, Tanter M, Fink M, Nordez A: Ultrafast imaging of in vivo muscle contraction using ultrasound. Applied Physics Letters. 2006, 89: 3-10.1063/1.2378616.View ArticleGoogle Scholar
- Nordez A, Gallot T, Catheline S, Guevel A, Cornu C, Hug F: Electromechanical delay revisited using very high frame rate ultrasound. J Appl Physiol. 2009, 106: 1970-1975. 10.1152/japplphysiol.00221.2009.View ArticlePubMedGoogle Scholar
- Revell J, Mirmehdi M, McNally D: Computer vision elastography: Speckle adaptive motion estimation for elastography using ultrasound sequences. IEEE Trans Med Imaging. 2005, 24: 755-766. 10.1109/TMI.2005.848331.View ArticlePubMedGoogle Scholar
- Farron J, Varghese T, Thelen DG: Measurement of tendon strain during muscle twitch contractions using ultrasound elastography. IEEE Trans Ultrason Ferroelectr Freq Control. 2009, 56: 27-35. 10.1109/TUFFC.2009.1002.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopata RGP, Nillesen MM, Hansen HHG, Gerrits IH, Thijssen JM, de Korte CL: Performance of two dimensional displacement and strain estimation techniques using a phased array transducer. Ultrasound in Medicine & Biology. 2009, 35: 2031-2041.View ArticleGoogle Scholar
- Witte RS, Dow DE, Olafsson R, Shi Y, O'Donnell M: High resolution ultrasound imaging of skeletal muscle dynamics and effects of fatigue. 2004 IEEE Ultrasonics Symposium. Edited by: Yuhas MP. 2004, New York: Ieee, 1-3: 764-767. full_text. Ultrasonics SymposiumView ArticleGoogle Scholar
- Witte RS, Kim K, Martin BJ, O'Donnell M: Effect of fatigue on muscle elasticity in the human forearm using ultrasound strain imaging. 2006, New York: IEEEView ArticleGoogle Scholar
- Gonzalez RC, Woods RE: Digital Image Processing. 2002, Prentice Hall, SecondGoogle Scholar
- Esbensen K, Geladi P: Strategy of multivariate image analysis (MIA). Chemometrics Intell Lab Syst. 1989, 7: 67-86. 10.1016/0169-7439(89)80112-1.View ArticleGoogle Scholar
- Geladi P: Some special topics in multivariate image analysis. Chemometrics Intell Lab Syst. 1992, 14: 375-390. 10.1016/0169-7439(92)80119-O.View ArticleGoogle Scholar
- Geladi P, Isaksson H, Lindqvist L, Wold S, Esbensen K: Principal component analysis of multivariate images. Chemometrics Intell Lab Syst. 1989, 5: 209-220. 10.1016/0169-7439(89)80049-8.View ArticleGoogle Scholar
- Korinek J, Wang JW, Sengupta PP, Miyazaki C, Kjaergaard J, McMahon E, Abraham TP, Belohlavek M: Two-dimensional strain - a Doppler-independent ultrasound method for quantitation of regional deformation: Validation in vitro and in vivo. J Am Soc Echocardiogr. 2005, 18: 1247-1253. 10.1016/j.echo.2005.03.024.View ArticlePubMedGoogle Scholar
- Amundsen BH, Helle-Valle T, Torp H, Crosby J, Støylen A, Ihlen H, Smiseth O, Slørdahl SA: Ultrasound speckle tracking reduces angle dependency of myocardial strain estimates - Validation by sonomicrometry. Euroecho 8. 2004, European Association of EchocardiographyGoogle Scholar
- Helle-Valle T, Crosby J, Edvardsen T, Lyseggen E, Amundsen BH, Smith HJ, Rosen BD, Lima JAC, Torp H, Ihlen H, Smiseth OA: New noninvasive method for assessment of left ventricular rotation - Speckle tracking echocardiography. Circulation. 2005, 112: 3149-3156. 10.1161/CIRCULATIONAHA.104.531558.View ArticlePubMedGoogle Scholar
- Lucas BD, Kanade T: An Iterative Image Registration Technique with an Application to Stereo Vision. International Joint Conference on Artificial Intelligence (IJCAI). 1981, 674-679.Google Scholar
- Bouguet JY: Pyramidal Implementation of the Lucas Kanade Feature Tracker: Description of the algorithm. Intel Corporation Microprocessor Research Labs. 2000Google Scholar
- Pearson K: On lines and planes of closest fit to systems of points in space. Philosophical Magazine. 1901, 2: 559-572.View ArticleGoogle Scholar
- Eriksson L, Wold S, Trygg J: Multivariate analysis of congruent images (MACI). J Chemometr. 2005, 19: 393-403. 10.1002/cem.944.View ArticleGoogle Scholar
- Bharati MH, Liu JJ, MacGregor JF: Image texture analysis: methods and comparisons. Chemometrics Intell Lab Syst. 2004, 72: 57-71. 10.1016/j.chemolab.2004.02.005.View ArticleGoogle Scholar
- Aagaard P, Andersen JL, Dyhre-Poulsen P, Leffers AM, Wagner A, Magnusson SP, Halkjær-Kristensen J, Simonsen EB: A mechanism for increased contractile strength of human pennate muscle in response to strength training: changes in muscle architecture. Journal of Physiology. 2001, 534: 613-623. 10.1111/j.1469-7793.2001.t01-1-00613.x.View ArticlePubMedPubMed CentralGoogle Scholar
- De Monte G, Arampatzis A: Influence of different shortening velocities preceding stretch on human triceps surae moment generation in vivo. Journal of Biomechanics. 2008, 41: 2272-2278. 10.1016/j.jbiomech.2008.04.008.View ArticlePubMedGoogle Scholar
- Wold S, Albano C, Dunn WJ, Edlund U, Esbensen K, Geladi P, Hellberg S, Johansson E, Lindberg W, Sjöström M: Multivariate data analysis in chemistry. Chemometrics - Mathematics and Statistics in Chemistry. Edited by: Kowalski BR. 1984, Dordrecht, The Netherlands: D. Reidel Publishing Company, 17-81.Google Scholar
- Ito L, 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.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/10/9/prepub