Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dynamic MRI of plantar flexion: A comprehensive repeatability study of electrical stimulation-gated muscle contraction standardized on evoked force

  • Xeni Deligianni ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Supervision, Validation, Visualization, Writing – original draft

    xeni.deligianni@unibas.ch

    Affiliations Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland, Department of Biomedical Engineering, University of Basel, Allschwil, Switzerland

  • Anna Hirschmann,

    Roles Investigation, Writing – review & editing

    Affiliation Department of Radiology, University Hospital Basel, Basel, Switzerland

  • Nicolas Place,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Institute of Sport Sciences, University of Lausanne, Lausanne, Switzerland

  • Oliver Bieri,

    Roles Resources, Writing – review & editing

    Affiliations Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland, Department of Biomedical Engineering, University of Basel, Allschwil, Switzerland

  • Francesco Santini

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Software, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliations Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland, Department of Biomedical Engineering, University of Basel, Allschwil, Switzerland

Abstract

Quantification of skeletal muscle contraction in Magnetic Resonance Imaging (MRI) is a non-invasive method for studying muscle motion and deformation. The aim of this study was to evaluate the repeatability of quantitative measures such as strain, based on single slice dynamic MRI synchronized with neuromuscular electrical stimulation (NMES) and standardized to a similar relative force level across various individuals. Unilateral electrical stimulation of the triceps surae muscles was applied in eight volunteers during single-slice, three-directional phase contrast MRI acquisition at a 3T MRI scanner. To assess repeatability, the same process was executed on two different days by standardizing the stimulation aiming at evoking a fixed percentage of their maximal voluntary force in the same position. Except from the force, the effect of using the current as reference was evaluated on day two as a secondary acquisition. Finally, the presence of fatigue induced by NMES was assessed (on day one) by examining the difference between consecutive measurements. Strain maps were derived from the acquired slice at every time point; distribution of strain in the muscle and peak strain over the muscle of interest were evaluated for repeatability. It was found that fatigue did not have an appreciable effect on the results. The stimulation settings based on evoked force produced more repeatable results with respect to using the current as the only reference, with an intraclass correlation coefficient between different days of 0.95 for the former versus 0.88 for the latter. In conclusion, for repeatable strain imaging it is advisable to record the force output of the evoked contraction and use that for the standardization of the NMES setup rather than the current.

Introduction

Quantification of skeletal muscle contraction in magnetic resonance imaging (MRI) can give valuable insight into muscle motion. The contraction can be either voluntary or evoked through neuromuscular electrical stimulation (NMES), which are different not only in terms of physiological response of the muscle (e.g. motor unit recruitment), but also in terms of practical details of their use in muscle diagnostics, rehabilitation and training. While voluntary contraction follows the Henneman size principle (i.e., small diameter motor units are recruited at lower force levels as compared to larger motor units) [1, 2], standard NMES induces a non-selective and mostly superficial random motor unit recruitment, allowing a greater contribution of type II muscle fiber at low force levels [1, 36]. Given the peculiarities of motor unit recruitment under NMES [4], NMES is a useful tool both for patients and for athletes in rehabilitation. Firstly, if a patient is unable to produce force levels high enough for a dynamic MRI muscle examination, NMES use can be a viable alternative [7, 8]. Secondly, since it follows a different activation pathway than voluntary contraction, it can yield complimentary information to examinations during voluntary exercise even in healthy individuals [9].

Muscle exercise in MRI can be visualized with various methods. Phase contrast (PC) MRI [1012], amongst other techniques such as spin-tagging [1315] and DENSE imaging [16], allows the dynamic acquisition of velocity and displacement images of various tissues, including the skeletal muscle. As imaging is usually not fast enough to be performed in real time (although there are reports of usage of this method, at the cost of lower spatial and temporal resolution [17]), a prospective or retrospective gating is usually implemented, under the assumption that the motion being quantified is occurring in a repetitive cycle through the same arc of motion. In this case, a portion of the k-space is acquired for every cycle and then the acquired datasets are reordered in a single period. For this reason, PC imaging requires the acquisition of several flow-sensitized images for each temporal frame of a periodic contraction.

Although there have been many promising studies [10, 11, 1822], the reproducibility of PC imaging over spin tagging has been questioned [23]. The error of phase contrast sequence was reduced with the proper settings such as the phase encoding direction [24, 25], but there has been no relevant study for motion guided by NMES.

PC velocity images can be used for the calculation of displacement maps and subsequently other biomechanical properties of the skeletal muscle, which describe the deformation. In this way, a three-dimensional quantification of muscle motion is obtained in a non-invasive way. Information about the deformation behavior of the muscle during contraction is valuable, because it gives input for the condition of the muscle such as e.g. existence of stiffness or atrophy [13]. Two parameters are mainly used for the final contractility evaluation: the strain [26] and the strain rate [13, 19]. The strain is a tensorial quantity, which is defined as change of length per unit length in each spatial direction of a material under stress with respect to its length at rest. The strain rate is the temporal derivative of the strain, and it is therefore an instantaneous measure that requires no reference state. It can either be calculated from the strain, or directly from the velocity field.

It has been shown that strain/ strain rate in leg muscles during voluntary or evoked isometric plantar flexion expresses age-related dependence [18, 27]. Lee et al. showed that unilateral lower limb suspension significantly influenced strain values [28]. Recently, a 3D PC scan protocol was presented for imaging of voluntary contraction [11]. Even though it has been shown that it is possible to visualize or yield meaningful quantitative contraction characteristics with MRI, the standardization of these measurements is needed before using these methods for diagnostic purposes.

The challenges of scanning standardization during voluntary contraction are the requirement for synchronization of the feedback loop with the MRI acquisition and the compliance of the subject with the exercise, which becomes even more important in case of patients [1012, 1921]. An alternative is to use NMES for external controlled triggering.

Recently, it has been shown that MRI of the thigh muscles during NMES is feasible and the resulting strain/strain rate significantly depend on the level of the applied current [29]. These results were shown on single slice acquisitions. Since it is well accepted that muscle motion is three dimensional, the general aim of the current study was to investigate if the results of a two-dimensional acquisition are sufficiently robust under rescan conditions and which NMES settings can improve reproducibility. Although the current used is a decisive factor, it might not be a sufficient reference to ensure similar repeatability conditions. Therefore, the induced force was used in this study as the main reference parameter [30] (in case of the isometric plantar flexion muscles, the force is expressed as resistance). The main goal of the present study was to investigate whether the force is a more robust reference parameter as the current for inducing a reliable periodic neuromuscular stimulation during MRI acquisition and finally if there is a muscle fatigue effect already during the acquisition. Finally, a preliminary control of how the results are affected by muscle size will be performed. This overall assessment is relevant for the translation of this method into clinical practice as it is crucial to quantify the inherent variability of the measurements induced by the acquisition method itself.

Materials and methods

The Ethics Commission of Northwest and Central Switzerland (Ethikkommission Nordwest- und Zentralschweiz-EKNZ) approved the study and all participants gave informed consent. Eight healthy volunteers (5 male/ 3 female, mean age: 31.9 ± 8.1 years, mean weight: 80.7 ± 15.1 kg, mean height: 1.79 ± 0.11 m) who usually performed low to moderate physical activity according to the short international physical activity questionnaire (IPAQ) were included in this study. The volunteers were asked not to perform any vigorous activity the day of the examination. The scans were performed on one side only and exclusion criteria, additional to the standard MRI contraindications, were recent injuries or operations in the examined leg and systemic pathologies with an involvement of the musculoskeletal system (e.g. arthritis, neuromuscular/neurodegenerative diseases, etc.). Experiments were performed on two different days within a time period of maximum four months.

An MR(Magnetic Resonance)-compatible force sensor was used to record the NMES-induced force during the dynamic scan at the height of the forefoot with a sampling rate of 10 Hz [30]. The sensor was attached to an MR-compatible supportive construction, which was fixed to the scanner bed and in addition kept the foot in a fixed position, as this setup is meant for isometric motion. After placing the volunteer on the scanner bed, in the same supine position as during the scan and just before moving the bed inside the scanner bore, the maximum voluntary plantar flexion force (MVF) was measured by the sensor. The leg was examined in a neutral position with the knee and hip in full extension, while the ankle angle was 90°. Measurements were repeated two to four times with at least 20 s recovery time until the MVF variation was less than 10% between the last two trials.

A summary of the various dynamic acquisitions is given in Table 1.

thumbnail
Table 1. Summary of the experimental protocol and the various scans of day 1 and day 2.

https://doi.org/10.1371/journal.pone.0241832.t001

The MRI acquisition was performed on a clinical 3T MRI scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany). The synchronization of the dynamic image acquisition and the electrical stimulation was accomplished as in [29]. In short, the NMES device was set to “synchronous mode,” where its two channels deliver stimulation at the same time. While the device remained outside the scanner room, one channel was guided to the room and attached through two electrodes to the subject’s skin for stimulation, while the other was fed as input to an Arduino custom made circuit [29] and used to generate the trigger signal for the MRI data acquisition.

For the triceps surae muscle belly stimulations, two rectangular electrodes were positioned over the gastrocnemii (below the popliteal fossa) and soleus (above the calcaneus) muscles. In order to identify the position of the electrodes (5.1 x 8.9 cm2 rectangular self-adhesive gel-based NMES electrodes [TensUnits.com, USA]) on the muscle belly of the calf, the subjects were asked to perform plantar flexion while standing prior to entering the scanner room. For localization purposes, glycerin markers were placed on the surface of both electrodes and a parasagittal slice through the markers was acquired.

The plane of data acquisition was decided based on the position of the two markers and in tendentially sagittal orientation. An out-of-phase gradient echo single slice image (voxel size: 1.125 x 1.125 x 3 mm3) was acquired to be used as a reference for precise segmentation.

The two-dimensional (2D) three-directional phase contrast images (i.e., single slice with velocity sensitivities in the three spatial directions, in-plane and through-plane) were acquired with velocity encoding of 30 cm/s on all directions, TR = 8.1 ms, TE = 5.63 ms and resolution 1.6 x 1.6 x 5 mm3, bandwidth/ pixel = 765 Hz/ Px, flip angle = 7°, FOV = 225 x 300 mm2, matrix: 192 x 144, 1 k-space line per segment, parallel imaging factor 2 (Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) [31], 24 integrated reference lines), acquisition time 2.1 min and 43 temporal phases. It should be noted that there is some additional dead time after release for triggering.

A standard digital NMES unit (EM49, Beurer GmbH, Germany) was used to induce periodic contractions. The NMES protocol consisted of 400 us, 80 Hz bipolar rectangular pulses. The stimulation current was raised until to the point that it induced a force equal to 15% of the participant’s maximum force, unless if the volunteer felt uncomfortable (i.e., in which case the stimulation current was reduced). Periodic stimulation pulses (750 ms contraction, 750 ms pause) were applied during the whole duration of dynamic image acquisition (appr. 2.1 min of repeated stimulation pulses for the acquisition of one slice). The stimulation current was kept constant during the whole acquisition.

For every PC dataset, the average force evolution over all contraction cycles was calculated. Subsequently, the percentage of difference with respect to the average force of the reference scan D0 was also calculated.

Post-processing

For strain calculation, a region of interest (ROI) was drawn on an out-of-phase gradient echo image including both regions of the gastrocnemius and the soleus muscle. The velocity images were elaborated offline with Matlab (The Mathworks, Inc., Natick, MA, USA, R2018b). An example of velocity images at the first velocity peak (beginning of contraction) and at the maximum point for strain (approximately at the middle of the contraction) is given on Fig 1. Velocity images were corrected for phase shading, the displacement was calculated with forward/backward integration and the Langrangian strain was computed [29]. To increase precision the image grid was interpolated (factor of 9.9) spatially with a cubic interpolation using not-a-knot end conditions. The analysis was performed with the assumption that the pixels of the acquired slice do not move out of the slice during the acquisition. The tensors were diagonalized and the positive eigenvalue was considered as pixel-wise principal strain. The principal strain maps were calculated from the displacement maps and then visualized as described before [29], but the method of displacement derivative calculation was replaced with the addition of MaxPol for smoother numerical differentiation [32, 33]. In this case, a Selesnick differentiator was used [34, 35] and a two-dimensional steerable derivative kernel (15th order tap polynomial, 3rd degree polynomial controling the cutoff threshold at x- or y- axes, 1st order of differentiation).

thumbnail
Fig 1. Example of raw velocity images with velocity sensitivity in the left-right (left) and top-bottom (right) directions (relative to the image) for two different frames at the first velocity peak (top) and at the time point where strain is maximum (bottom) are given.

https://doi.org/10.1371/journal.pone.0241832.g001

As a representative value of strain, the spatial average over the ROI (see Fig 2) was calculated for every time frame of the reconstructed single contraction period and the maximum of this time curve was considered as peak strain. The rates at which the strain reached the maximum (positive or buildup rate) and relaxed to zero (negative or release rate) were also calculated by fitting a sigmoid curve to the corresponding portions of the strain curve (see Fig 2). The differences of the various rates were normalized to the maximum calculated rate amongst all data points.

thumbnail
Fig 2. Example of a region-of-interest, a strain map with the MaxPlot algorithm and a strain map with the original postprocessing.

On the strain graph, strain is given as a function of time, with overlayed explanation graphs explaining positive and negative rate of the strain.

https://doi.org/10.1371/journal.pone.0241832.g002

Statistical tests—Evaluation of repeatability & reliability

Similarity index (SI) was calculated in Matlab to compare strain maps on the middle of the contraction plateau (middle frame: 19th frame). Affine registration was performed beforehand and then the calculation of the structural SI [36].

All statistical calculations and plots were produced using Matlab and R-studio [37]. The representative strain values of all eight volunteers were compared with scatter and boxplots. Three comparisons were performed following the scan protocol presented in Table 1; the basis scan D0 was compared to: a) Dfor, b) Dfat, c) Dcur.

To check for the repeatability by using the force as reference, D0 was compared to Dfor and respectively to check for the current as reference D0 was compared to Dcur.

The two-way intraclass correlation coefficient (ICC) [38] was calculated to compare the first scans of each day D0 vs Dfor and D0 to Dcur (observations of exchangeable order). The ICC was considered a measure of reliability since the measurements are exchangeable in order [39]. The repeatability coefficient (RC) was calculated as (39): [1] where SD is the within-subject standard deviation. According to the respective theory, if the differences between two measurements are normally distributed (e.g., the differences of strain between successive scans) the absolute difference between two measurements on 95% of occasions is expected to not differ more than RC.

Finally, to check for potential muscle fatigue induced by NMES, D0 and Dfat were compared by using the concordance correlation coefficient since the measurements are non-exchangeable (cccrm package of R for Concordance Correlation Coefficient (CC) for Repeated Measurements by U-statistics was used).

Results

To evoke a similar force, the current on day 2 was typically adapted by 1.5 mA (i.e., median difference), with the exception of two cases where the difference was 13 and 26 mA.

The strain maps of a time frame at the middle of contraction from the scan-rescan test (D0 vs Dfor) are presented in Fig 3. From a quantitative comparison, the average SI of all volunteers for scan-rescan strains with similar force output, denoted as SI(D0, Dfor), is in average 0.81 ± 0.02 comparable to 0.80 ± 0.03 of SI(D0,Dcur) (i.e., identical current setup in different days). This is only slightly worse than 0.90 ± 0.03 of the SI(D0, Dfat) (i.e., fatigue evaluation). With a few exceptions (e.g. subject 6 with SI(D0, Dfor) = 0.79 & SI(D0,Dcur) = 0.78) the strain maps are qualitatively very similar (see Fig 3).

thumbnail
Fig 3. Strain-maps of a central frame from the first scan of day 1 (D0) from all the subjects from subject 1 (S1) to subject 8 (S8).

https://doi.org/10.1371/journal.pone.0241832.g003

For the direct analysis of summarized strain values (i.e. maximum of the average value of the map, see Table 2) ICC/CC, RC and SD were compared. In order of increasing correlation and decreasing error, better repeatability and reliability is shown between consecutive measurements. A smaller ICC/CC and a larger error are found between scan-rescan measurements. The largest error and smallest ICC/CC are calculated between scan-rescan measurements setup based on current and not force. Finally, regarding the absolute strain values the SD was 0.01, which also gives the limit of precision.

thumbnail
Table 2. Reliability and repeatability coefficient of strain maps.

https://doi.org/10.1371/journal.pone.0241832.t002

The individual maximum strains were compared with boxplots (Fig 4). As expected, overall the minimum difference was observed between consecutive acquisitions (i.e. no fatigue from the stimulation). Between rescans based on force or current, the differences are qualitatively similar. However, if instead of the absolute values, we visualize the differences for each individual (see Fig 5, top) there appears to be slightly less variability when the stimulation intensity is regulated according to exerted force in comparison to current (see Fig 5, top, central vs right column). In addition, the positive and negative rates of strain are less variable in the case of using the force as standardization factor (see Fig 5, bottom, central vs right column). However, the negative strain rate shows higher variability than the positive in the comparison of consecutive scans (see Fig 5, bottom, left column).

thumbnail
Fig 4. Boxplots of the maximum strain.

(top) first scan of day 1 in comparison to first scan of day 2 which are setup at similar evoked force (D0-Dfor), (middle) the first scan of day 1 and a subsequent one to estimate fatigue (D0-Dfat), (bottom) the first scans of day 1 and the second of day 2 which are setup at the same current (D0-Dcur).

https://doi.org/10.1371/journal.pone.0241832.g004

thumbnail
Fig 5. Scatter plots of differences.

a. individual strain differences, b. normalized differences of the positive rate of the strain, c. normalized differences of negative rate of the strain between the repeated scan on day 1 (“fatigue effect”– D0-Dfat), the scans on day 1 and day 2 with same force output (“same force”– D0-Dfor), and the scans on day 1 and day 2 with same stimulation current (“same current”– D0-Dcur).

https://doi.org/10.1371/journal.pone.0241832.g005

Ultimately, for every acquisition (D0, Dfat, Dfor, Dcur) the average of all rates across volunteers was calculated, which were lower for the positive rates (137.0 s-1, 130.5 s-1, 135.3 s-1, 133.6 s-1) than for the negative (163.1 s-1, 164.5 s-1, 169.3 s-1, 163.8 s-1). Overall, it was observed that there is little difference between the various sessions.

Finally, examples of the curves of induced force are presented in Fig 6. The mean difference of the average scan force was maximum (43.8% (Q1/Q3 = 17.5/78.7%)) for D0-Dcur, minimum for D0-Dfat (10.7% (Q1/Q3 = 8.7/32.5%)). For the scan-rescan tests D0-Dfor the variation was in between (21.3% (Q1/Q3 = 14.3/44.3%)).

thumbnail
Fig 6. Evoked force during stimulation for three different cases.

(top) higher resulting force, (middle) average, (bottom) low. In every window, the force curves for every stimulation are overlayed and on the top the average force during the scan.

https://doi.org/10.1371/journal.pone.0241832.g006

Discussion

Established biomarkers for muscle characterization such as fat fraction and T2 relaxation are often indicators of an already progressed disease status and early detection of muscle alterations is very important [40]. The examined method of characterizing the condition of the muscle is simple, non-invasive, and of low additional cost on top of the cost of the MRI scan itself. However, since the settings of NMES are fundamental for the standardization of the evoked contraction [3, 7] and therefore of the MRI-based results, it is crucial to control the stimulation in the most repeatable way. The force output is employed for dynamic MRI of voluntary contraction [13], but to our knowledge there is no systematic study of a comparable setup for dynamic MRI during NMES.

In this study, we considered the effect of taking the evoked force or the applied current as a reference, as well as the influence of fatigue from continuously applying NMES. The results show that the evoked force during stimulation is the best reference for repeatable results. The effect of fatigue induced by this NMES protocol was minimal. It should be noted that in most cases only a minimal adaptation of the current was necessary to obtain the same force output, rendering the two control parameters de facto analogous. Nevertheless, in two cases the difference was large, which supports the fact the force is an important parameter to consider.

It should be noted that the sampling rate of the force is different (approximately two times longer) compared to the MR velocity acquisition. However, the two signals are fundamentally different, the force being recorded in real-time, whereas the MR velocity is acquired in a triggered fashion. The lower temporal resolution of the real-time force is in fact no limitation, as the actual dynamic of the force signal does not enter the analysis. It is not used for triggering the MR acquisition (in contrast to what was done in [41]), and it was in fact only recorded to ascertain the periodicity of the movement and the repeatability of the force output in the contracted state.

The remaining measurement variation of the quantitative strain measurements could be partly attributed to the three-dimensional nature of muscle motion. Undoubtedly muscle contraction is a three-dimensional effect and eventually a 3D analysis is important, but nevertheless this two-dimensional acquisition is certainly valuable since it is not time demanding. In the future, this limitation could be overcome by the application of advanced multi-planar or 3D imaging methods [11].

An alternative approach, not considered in this present study, could be the usage of spin tagging [42] methods instead of PC-MRI, which have shown good reliability for the detection of muscle motion [23]. However, Sinha et al [13] demonstrated good accordance between these two methods, and the fading of the tagging following T1 relaxation might render the choice of NMES stimulation protocols less flexible.

This study has some additional limitations. A factor that adds variability is the repositioning of the volunteer and of the electrodes. To properly evaluate repeatability, we would have had to repeat the scans, with the electrodes at the same position and with the use of longer breaks (i.e., to exclude fatigue). However, this was not possible due to time restrictions and we considered the implemented setup to be more representative of potential realistic conditions.

In this study including only young healthy volunteers the data showed quite homogeneous and repeatable results. In future investigations, various populations (e.g. in terms of age, sex and training) will be scanned as well, which may show a more stratified distribution of the assessed parameters. For example in a similar study [27], the response of senior volunteers was shown to be lower than in the younger cohort. In addition, strain needs to be evaluated in degenerative muscles to establish whether the proposed method can indeed differentiate between healthy and diseased muscle. Finally, dynamic muscle data from different age groups of healthy volunteers with normal daily activity is an essential first step in order to proceed on studying abnormal or pathophysiological muscle conditions.

In conclusion, this study shows that strain values produced with MRI of NMES-controlled evoked muscle contraction are highly repeatable when the induced electrically-evoked force during the scan is used as a reference. In this case the strain maps were not only qualitatively similar, but also the maximum strain values and the rates of strain increase, and release were very repeatable.

Acknowledgments

We would like thank Sebastian Manneck for his help with muscle segmentation testing.

References

  1. 1. Gregory CM, Bickel CS. Recruitment Patterns in Human Skeletal Muscle During Electrical Stimulation. Phys Ther. 2005 Apr 1;85(4):358–64. pmid:15794706
  2. 2. Henneman E, Somjen G, Carpenter DO. Functional significance of cell size in spinal motoneurons. J Neurophysiol. 1965 May 1;28(3):560–80. pmid:14328454
  3. 3. Maffiuletti NA. Physiological and methodological considerations for the use of neuromuscular electrical stimulation. Eur J Appl Physiol. 2010 Sep 1;110(2):223–34. pmid:20473619
  4. 4. Bickel CS, Gregory CM, Dean JC. Motor unit recruitment during neuromuscular electrical stimulation: a critical appraisal. Eur J Appl Physiol. 2011 Aug 26;111(10):2399. pmid:21870119
  5. 5. Babault N, Cometti G, Bernardin M, Pousson M, Chatard J-C. Effects of electromyostimulation training on muscle strength and power of elite rugby players. J Strength Cond Res. 2007 May;21(2):431–7. pmid:17530954
  6. 6. Kamen G, Knight CA. Training-related adaptations in motor unit discharge rate in young and older adults. J Gerontol Biol Sci Med Sci. 2004 Dec;59(12):1334–8. pmid:15699535
  7. 7. Maffiuletti NA, Gondin J, Place N, Stevens-Lapsley J, Vivodtzev I, Minetto MA. Clinical Use of Neuromuscular Electrical Stimulation for Neuromuscular Rehabilitation: What Are We Overlooking? Arch Phys Med Rehabil. 2018;99(4):806–12. pmid:29233625
  8. 8. Gondin J, Vilmen C, Cozzone PJ, Bendahan D, Duhamel G. High-field (11.75T) multimodal MR imaging of exercising hindlimb mouse muscles using a non-invasive combined stimulation and force measurement device. NMR Biomed. 2014 Aug;27(8):870–9. pmid:24890578
  9. 9. Gandevia SC. Spinal and Supraspinal Factors in Human Muscle Fatigue. Physiol Rev. 2001 Jan 10;81(4):1725–89. pmid:11581501
  10. 10. Sinha S, Sinha U. Dynamic MR Imaging of the Skeletal Musculature: From Static Measures to a Dynamic Assessment of the Muscular (Loco-) Motion. In: Weber M-A, editor. Magnetic Resonance Imaging of the Skeletal Musculature. Springer Berlin Heidelberg; 2013. p. 165–84.
  11. 11. Mazzoli V, Schoormans J, Froeling M, Sprengers AM, Coolen BF, Verdonschot N, et al. Accelerated 4D self-gated MRI of tibiofemoral kinematics. NMR Biomed. 30(11):e3791. pmid:28873255
  12. 12. Drace JE, Pelc NJ. Elastic deformation in tendons and myotendinous tissue: measurement by phase-contrast MR imaging. Radiology. 1994 Jun;191(3):835–9. pmid:8184075
  13. 13. Sinha S, Hodgson JA, Finni T, Lai AM, Grinstead J, Edgerton VR. Muscle kinematics during isometric contraction: development of phase contrast and spin tag techniques to study healthy and atrophied muscles. J Magn Reson Imaging. 2004 Dec;20(6):1008–19. pmid:15558560
  14. 14. Englund EK, Elder CP, Xu Q, Ding Z, Damon BM. Combined diffusion and strain tensor MRI reveals a heterogeneous, planar pattern of strain development during isometric muscle contraction. Am J Physiol Regul Integr Comp Physiol. 2011 May;300(5):R1079–1090. pmid:21270344
  15. 15. Moerman KM, Sprengers AMJ, Simms CK, Lamerichs RM, Stoker J, Nederveen AJ. Validation of continuously tagged MRI for the measurement of dynamic 3D skeletal muscle tissue deformation. Med Phys. 2012 Apr;39(4):1793–810. pmid:22482602
  16. 16. Zhong X, Epstein FH, Spottiswoode BS, Helm PA, Blemker SS. Imaging two-dimensional displacements and strains in skeletal muscle during joint motion by cine DENSE MR. J Biomech. 2008 Jan 1;41(3):532–40. pmid:18177655
  17. 17. Asakawa DS, Nayak KS, Blemker SS, Delp SL, Pauly JM, Nishimura DG, et al. Real-time imaging of skeletal muscle velocity. J Magn Reson Imaging. 2003 Dec;18(6):734–9. pmid:14635159
  18. 18. Sinha U, Malis V, Csapo R, Moghadasi A, Kinugasa R, Sinha S. Age-related differences in strain rate tensor of the medial gastrocnemius muscle during passive plantarflexion and active isometric contraction using velocity encoded MR imaging: Potential index of lateral force transmission. Magn Reson Med. 2015 May 1;73(5):1852–63. pmid:25046255
  19. 19. Malis V, Sinha U, Sinha S. Compressed sensing velocity encoded phase contrast imaging: Monitoring skeletal muscle kinematics. Magn Reson Med. 2020;84(1):142–56. pmid:31828833
  20. 20. Drace JE, Pelc NJ. Measurement of skeletal muscle motion in vivo with phase-contrast MR imaging. J Magn Reson Imaging. 1994 Apr;4(2):157–63. pmid:8180454
  21. 21. Drace JE, Pelc NJ. Skeletal muscle contraction: analysis with use of velocity distributions from phase-contrast MR imaging. Radiology. 1994 Nov;193(2):423–9. pmid:7972757
  22. 22. Zhou H, Novotny JE. Cine phase contrast MRI to measure continuum Lagrangian finite strain fields in contracting skeletal muscle. J Magn Reson Imaging. 2007 Jan;25(1):175–84. pmid:17152055
  23. 23. Borotikar B, Lempereur M, Lelievre M, Burdin V, Salem DB, Brochard S. Quantifying musculoskeletal motion using dynamic MRI techniques: A systematic review of validity and perspectives for musculoskeletal disorder evaluation. Ann Phys Rehabil Med. 2017 Sep 1;60:e28.
  24. 24. Behnam AJ, Herzka DA, Sheehan FT. Assessing the accuracy and precision of musculoskeletal motion tracking using cine-PC MRI on a 3.0T platform. J Biomech. 2011 Jan 4;44(1):193–7. pmid:20863502
  25. 25. Jensen ER, Morrow DA, Felmlee JP, Odegard GM, Kaufman KR. Error analysis of cine phase contrast MRI velocity measurements used for strain calculation. J Biomech. 2015 Jan 2;48(1):95–103. pmid:25433567
  26. 26. Finni T, Hodgson JA, Lai AM, Edgerton VR, Sinha S. Nonuniform strain of human soleus aponeurosis-tendon complex during submaximal voluntary contractions in vivo. J Appl Physiol. 2003 Aug;95(2):829–37. pmid:12716873
  27. 27. Deligianni X, Klenk C, Place N, Garcia M, Pansini M, Hirschmann A, et al. Dynamic MR imaging of the skeletal muscle in young and senior volunteers during synchronized minimal neuromuscular electrical stimulation. Magn Reson Mater Phys Biol Med. 2020;(33):393–400. pmid:31650419
  28. 28. Lee H-D, Finni T, Hodgson JA, Lai AM, Edgerton VR, Sinha S. Soleus aponeurosis strain distribution following chronic unloading in humans: an in vivo MR phase-contrast study. J Appl Physiol. 2006 Jun;100(6):2004–11. pmid:16424072
  29. 29. Deligianni X, Pansini M, Garcia M, Hirschmann A, Schmidt-Trucksäss A, Bieri O, et al. Synchronous MRI of muscle motion induced by electrical stimulation. Magn Reson Med. 2017 Feb 1;77(2):664–72. pmid:26898990
  30. 30. Santini F, Bieri O, Deligianni X. OpenForce MR: A low-cost open-source MR-compatible force sensor. Concepts Magn Reson Part B Magn Reson Eng. 2018;48B(4):e21404.
  31. 31. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med. 2002 Jun;47(6):1202–10. pmid:12111967
  32. 32. Hosseini MS, Plataniotis KN. Derivative Kernels: Numerics and Applications. IEEE Trans Image Process. 2017 Oct;26(10):4596–611. pmid:28613176
  33. 33. Hosseini MS, Plataniotis KN. Finite Differences in Forward and Inverse Imaging Problems—MaxPol Design. ArXiv170908321 Math [Internet]. 2017 Sep 25 [cited 2019 Dec 4]; http://arxiv.org/abs/1709.08321
  34. 34. Selesnick IW, Burrus CS. Maximally flat low-pass FIR filters with reduced delay. IEEE Trans Circuits Syst II Analog Digit Signal Process. 1998 Jan;45(1):53–68.
  35. 35. Selesnick IW, Burrus CS. Generalized digital Butterworth filter design. IEEE Trans Signal Process. 1998 Jun;46(6):1688–94.
  36. 36. Wang Z, Bovik AC, Sheikh HR, Simoncelli EP. Image quality assessment: from error visibility to structural similarity. IEEE Trans Image Process. 2004 Apr;13(4):600–12. pmid:15376593
  37. 37. RStudio Team (2018). RStudio: Integrated Development for R. RStudio, Inc. [Internet]. Boston, MA; 2018. http://www.rstudio.com/
  38. 38. Shrout PE, Fleiss JL. Intraclass correlations: uses in assessing rater reliability. Psychol Bull. 1979 Mar;86(2):420–8. pmid:18839484
  39. 39. Bartlett JW, Frost C. Reliability, repeatability and reproducibility: analysis of measurement errors in continuous variables. Ultrasound Obstet Gynecol. 2008 Apr 1;31(4):466–75. pmid:18306169
  40. 40. Strijkers GJ, Araujo ECA, Azzabou N, Bendahan D, Blamire A, Burakiewicz J, et al. Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease—A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI. J Neuromuscul Dis. 2019 Jan 1;6(1):1–30. pmid:30714967
  41. 41. Sinha S, Shin DD, Hodgson JA, Kinugasa R, Edgerton VR. Computer-controlled, MR-compatible foot-pedal device to study dynamics of the muscle tendon complex under isometric, concentric, and eccentric contractions. J Magn Reson Imaging. 2012 Aug;36(2):498–504. pmid:22392816
  42. 42. Axel L, Dougherty L. MR imaging of motion with spatial modulation of magnetization. Radiology. 1989 Jun;171(3):841–5. pmid:2717762