Cardan angle rotation sequence effects on first-metatarsophalangeal joint kinematics: implications for measuring hallux valgus deformity
© Glasoe et al.; licensee BioMed Central Ltd. 2014
Received: 3 January 2014
Accepted: 7 May 2014
Published: 14 May 2014
There currently are no recommended standards for reporting kinematics of the first-metatarsophalangeal joint. This study compared 2 different rotation sequences of Cardan angles, with implications for understanding the measurement of hallux valgus deformity.
Thirty-one women (19 hallux valgus; 12 controls) participated. All were scanned in an open-upright magnetic resonance scanner, their foot posed to simulate the gait conditions of midstance, heel-off, and terminal stance. Using computer processes, selected tarsals were reconstructed into virtual bone models and embedded with principal-axes coordinate systems, from which the rotation matrix between the hallux and first metatarsal was decomposed into Cardan angles. Joint angles were then compared using a within factors (rotation sequence and gait condition) repeated-measures analysis of variance (ANOVA).
Only the transverse plane-first sequence consistently output incremental increases of dorsiflexion and abduction across gait events in both groups. There was an interaction (F ≥ 25.1; p < 0.001). Follow-up comparisons revealed angles were different (p < 0.05) at terminal stance.
Different rotation sequences yield different results. Extracting the first rotation in the transverse plane allows for the resting alignment of the hallux to deviate from the sagittal plane. Therefore, representing first-metatarsophalangeal joint kinematics with the transverse plane-first rotation sequence may be preferred, especially in cases of hallux valgus deformity.
KeywordsBiomechanics measurement Bunion deformity Gait simulation
Three-dimensional kinematics are often calculated with Cardan (Euler) angle equations. Cardan angles are quantified as ordered rotations about 3 axes, and represent joint motion displacements in the anatomical body planes. Foot dorsiflexion/plantar flexion joint rotations occur about the mediolateral-axis in the sagittal plane; inversion/eversion about the longitudinal-axis in the frontal plane; adduction/abduction about the vertical axis in the transverse plane. To reduce the sequence dependent nature of the matrix calculation, the first rotation is typically extracted about the axis of largest joint motion and the last rotation about the long-axis of the distal segment . Following these basic guidelines, the International Society of Biomechanics (ISB) recommended researchers adopt the sagittal, frontal, transverse rotation sequence as the standard for reporting kinematics of the ankle joint . No standards have been recommended for the first-metatarsophalangeal (1-MTP) joint.
Hallux valgus is chiefly characterized by an offset in transverse plane alignment of the 1-MTP joint. The hallux is abducted in relation to the long-axis orientation of the first metatarsal and when deformity is severe, abduction of the hallux exceeds 40 degrees . Even absent of deformity, the resting posture of the hallux tends to align in some degree of abduction with respect to the first metatarsal and in relation to the midsagittal line of the body. To account for this offset in joint alignment, Glasoe et al. [10, 11] modeled the 1-MTP joint using a transverse plane-first sequence. Other researchers [6, 8, 9] have followed ISB recommendations for the ankle, and modeled the 1-MTP joint using a sagittal plane-first sequence. Both solutions are correct, so preference becomes a matter of clinical interpretation.
The purpose of this study was to compare Cardan angle rotation sequence effects on 1-MTP joint kinematics during gait with implications for the measurement of hallux valgus deformity. Citing data recorded in gait trials [6, 7], and premised on how the foot is placed during stepping (Figure 1) and the trajectory of body’s center of pressure which moves in response to ground reaction forces , dorsiflexion and abduction of the hallux should occur during heel rise. Due to the significance of the vertical force and because the collateral ligaments check motion of the hallux in the transverse plane, dorsiflexion should be largest. About 50 degrees of 1-MTP joint motion is needed for normal gait [6–8]. Deformity increases loading on the medial border of the hallux , so the direction of the respective external and internal foot-to-ground contact forces may reduce dorsiflexion while further encouraging abduction of the hallux.
In summary, this study compared “sagittal” and “transverse” plane-first rotation sequence effects of 1-MTP joint angles in women grouped based on the presence or absence of hallux valgus deformity. The data were analyzed under 2 hypotheses. First, that different Cardan angle rotation sequences would output different angles. Second, that both rotation sequences would output increasing amounts of hallux dorsiflexion and abduction across gait events regardless of group assignment.
The experiment was retrospective in design. The data were collected on 31 (of 31 total) adult women enrolled in a two-part investigation [10, 11] that used weightbearing imaging to study tarsal kinematics. A comprehensive description of the participants and research methods has been reported elsewhere [10, 11]. Nineteen of the women had hallux valgus. Measurement (on images) of a hallux angle larger than 15° served as the threshold for classifying deformity. The hallux angle quantifies the offset in transverse plane alignment of the 1-MTP joint. The hallux was abducted 33° (range, 18 to 64°) with deformity, compared to 8° (range, -2 to 14°) in controls. Participants were group matched for demographics: age 48 ± 17 years; BMI 26 ± 6 kg/m2. All gave informed consent in accordance with University of Minnesota Human Research Institutional Review Board (#0709 M16823) guidelines.
Following imaging, selected bones (Figure 3) including the proximal phalanx of the hallux and the first metatarsal were reconstructed and registered across the gait conditions within the MR laboratory frame-of-reference system using MIMICS (Materialize, Leuven, Belgium). Principal component inertial-axes systems were embedded to define the spatial orientation of each virtual bone model . Coordinate system axes were directed Z-lateral (dorsi/plantar flexion), the Y-up (add/abduction), and X-anterior (in/eversion). The coordinate systems defining the hallux and first metatarsal embedded consistently and aligned nearly coincident with the MR frame of reference. The error associated with data processing was < 3 degrees . The measurement of tarsal motion was shown reliable and valid across the simulated gait conditions [10, 11].
First metatarsophalangeal joint angles were extracted from two different transformation matrices calculated from one Matlab (Mathworks Inc, Natick, Massachusetts) code. The sagittal plane-first (Z-Y-X) sequence followed ISB guidelines for the ankle joint [6, 8, 9], while the transverse plane-first (Y-Z-X) sequence was thought to better represent kinematics of the 1-MTP joint . Joint angles were described as mean values across gait conditions, and examined using a 2-way within factors (rotation sequence and gait condition) repeated measure analysis of variance (ANOVA). In the presence of an interaction, follow-up pairwise comparisons were completed at each level of condition. The ANOVA model was run separately on the two groups (hallux valgus and controls). Significance was set at p < 0.05.
There was an interaction (F ≥ 25.1; p < 0.001) between variables. The follow-up comparisons revealed all mean angles were different (p < 0.05) at TS. Differences were ≥ 5.1° in the sagittal plane about the Z-axis, ≥ 12.5° in the transverse plane about the Y-axis, and ≥ 8.2° in the frontal plane about the X-axis. The frontal plane angles were also divergent in direction.
De Mits et al.  and Leardini et al.  used the sagittal plane-first sequence to model 1-MTP joint with Cardan angles on adults without foot deformity. Dorsiflexion and abduction started before midstance and continued unabated until the hallux lifted off the ground at the last instant of push-off. Dorsiflexion was dominant in both studies. De Mitts  recorded 55 degrees of dorsiflexion and 50 degrees of abduction. Leardini  recorded 40 degrees of dorsiflexion and 10 degrees of abduction. Although the difference in joint motion between trials is notable, these results indicate that dorsiflexion and abduction should occur at the 1-MTP joint into late stance phase of gait (Figures 1 and 2). Consistent with data reported in these past gait trials [6, 7], increasing amounts of dorsiflexion and abduction were expected in this study - where gait was simulated by posing the foot on a system of wedges (Figure 3). When interpreting these results, it is important to remember the participant’s foot was positioned the same for both sequences, so any difference in 1-MTP joint angles shows the effects of the rotation sequence and nothing else.
Regardless of group assignment, both sequences captured increasing dorsiflexion (Figure 4) as should happen during gait (Figure 1). Dorsiflexion surpassed 50 degrees, with 15 degrees the increment of increase of from MS to HO, and 35 degrees from HO to TS. This total amount as well as the incremental change in dorsiflexion between gait conditions closely approximate the angle of wedging used to advance the foot to simulate gait (Figure 3).
Transverse plane abduction was the next largest component of motion (Figure 4). Of the 2 tested sequences, only the transverse plane-first sequence (Y-Z-X) output increasing abduction (~20°) across gait conditions in both groups. The sagittal plane-first sequence, in comparison, output as abduction (<5°) from MS to HO, and adduction (<5°) from HO to TS in the group with hallux valgus. Based on the placement of the foot, it is doubtful the hallux could adduct while pressed into the ground as the heel lifts and weight advances forward and displacement in the direction of adduction becomes even more unlikely in the group having hallux valgus deformity. Although the kinetics were not simulated, terminal stance gait propulsion has an active component powered by the muscles that flex the toes and ankle, and a passive component attributed to the tightening of the plantar fascia . As a consequence of hallux valgus deformity, pull of the flexor muscles and plantar fascia acting on the hallux shift to the lateral side of the 1-MTP joint. Loading must therefore abduct the hallux, but meaningful abduction (~20°) was represented in those with hallux valgus only when the 1-MTP joint angles were modeled with the transverse plane-first rotation sequence.
This implausible Cardan angle group result merits further explanation. The greater disparity between sequences in the hallux valgus group (Figure 4) is due to greater deviation in initial resting alignment of the 1-MTP from the straight (sagittal plane) dorsi/plantar flexion axis. Because resting posture alignment is not even close to straight (hallux angle = 33°) with hallux valgus, use of the sagittal plane-first sequence results in the first rotation occurring about a non-physiologic joint axis, whereas use of the transverse plane-first sequence adjusts for the initial presence of the abduction deformity, and then records dorsi/plantar flexion about that position.
The frontal plane angles were calculated as the last order-of-rotation. Divergent in direction (Figure 4), this difference in 1-MTP joint angles best highlight the sequence dependency effects of the Cardan angle calculation. Compounding the challenge of modeling the loading behaviors of the hallux in the frontal plane, surface palpation as required by most motion analysis techniques concentrates error around the long-axis of the hallux and first metatarsal . Researchers choosing to calculate the frontal plane position of the 1-MTP joint as the last rotation should exercise caution when interpreting this result.
This retrospectively designed study had limitations. Gait was simulated with static conditions. Posing a person in a MR scanner does not replicate the loading demands of gait. Importantly, joint angles were not validated by other methods of direct measurement, though the computation of 1-MTP joint dorsiflexion did approximate the angle of wedging used to advance the foot in staging gait. The study was also narrow in scope, as it singularly compared the computation of 1-MTP joint motion using the 2 most common rotation sequences. This line of research could be broadened to compare possible matrix multiplication differences associated 3D helical axis computations, or projection errors associated with direction cosine. In light of these limitations, this study does not settle the issue but demonstrates why basic guidelines that govern order-of-rotation decisions at other joints may not apply to the 1-MTP joint, as the hallux typically rests in a position of relative abduction.
The rotation matrices describing the Cardan angles between the first metatarsal and hallux were decomposed for the purpose of comparing sagittal plane-first and transverse plane-first rotation sequence effects. Both sequences output increasing amounts of dorsiflexion, only the transverse plane-first rotation sequence consistently output increasing amounts of abduction, and angles of in/eversion were divergent when representing 1-MTP joint displacements across gait conditions. Representing 1-MTP joint kinematics with a transverse, sagittal, and frontal plane sequence (Y-Z-X) may be preferred, especially in cases of hallux valgus deformity.
The authors acknowledge the contribution of Dr. Paula Ludewig. She collaborated to collect the data and in editing the final copy of this manuscript. This project was funded by the Minnesota Medical Foundation (#3970-9238-09).
- Cole GK, Nigg BM, Ronsky JL, Yeadon MR: Application of the joint co-ordinate system to three-dimensional joint attitude and movement representation: a standardization proposal. J Biomech Eng. 1993, 115: 344-349. 10.1115/1.2895496.View ArticlePubMedGoogle Scholar
- Wu G, Siegler S, Allard P, Kirtley C, Leardini A, Rosenbaum D, Whittle M, D'Lima DD, Cristofolini L, Witte H, Schmid O, Stokes I: ISB recommendation on definitions of joint coordinate system of various joints for the reporting of human joint motion-part I: ankle, hip, and spine. J Biomech. 2002, 35: 543-548. 10.1016/S0021-9290(01)00222-6.View ArticlePubMedGoogle Scholar
- Coughlin MJ, Jones CP: Hallux valgus: demographics, etiology, and radiographic assessment. Foot Ankle Int. 2007, 28: 759-777. 10.3113/FAI.2007.0759.View ArticlePubMedGoogle Scholar
- Michaud TC: Human locomotion: the conservative management of gait-related disorders. 2011, Newton, MA: Newton Biomechanics, 29-30.Google Scholar
- Root ML, Orien WP, Weed JH: Normal and abnormal function of the foot. 1977, Los Angeles, CA: Clinical Biomechanics Corp, 55-60.Google Scholar
- De Mits S, Segers V, Woodburn J, Elewaut D, De Clercq D, Roosen P: A clinically applicable six-segmented foot model. J Orthop Res. 2012, 30: 655-661. 10.1002/jor.21570.View ArticlePubMedGoogle Scholar
- Leardini A, Benedetti MG, Cantani F, Simoncini L, Giannini S: An anatomically based protocol for the description of foot segment kinematics during gait. Clin Biomech. 1999, 14: 528-536. 10.1016/S0268-0033(99)00008-X.View ArticleGoogle Scholar
- Nawoczenski DA, Baumhauer JF, Umberger BR: Relationship between clinical measurements and motion of the first metatarsophalangeal joint during gait. J Bone Joint Surg Am. 1999, A-81: 370-376.Google Scholar
- Umberger B, Nawoczenski D, Baumhauer J: Reliability and validity of first metatarsophalangeal joint orientation measured with an electromagnetic tracking device. Clin Biomech. 1999, 14: 74-76. 10.1016/S0268-0033(98)00052-7.View ArticleGoogle Scholar
- Glasoe WM, Phadke V, Pena FA, Nuckley DJ, Ludewig PM: An image-based gait simulation study of tarsal kinematics in women with hallux valgus. Phys Ther. 2013, 93: 1551-1562. 10.2522/ptj.20130025.View ArticlePubMedGoogle Scholar
- Glasoe WM, Jensen DD, Kampa BB, Karg LK, Krych AR, Pena FA, Ludewig PM: First ray kinematics in women with rheumatoid arthritis and bunion deformity: a gait simulation imaging study. Arthritis Care Res. in pressGoogle Scholar
- Glasoe WM, Pena FA, Phadke V, Ludewig PM: Cardan sequences on 1-MTP joint kinematics: implications for foot bunion [abstract]. Proceedings of the Annual American Society of Biomechanics (ASB) Conference. 2013, Omaha, 4-7. http://www.openconf.org/asb2013/modules/request.php?module=oc_program%26action=summary.php%26id=2061,Google Scholar
- Chiu MC, Wu HC, Chang LY, Wu MH: Center of pressure progression characteristics under the plantar region for elderly adults. Gait Posture. 2013, 37: 408-412. 10.1016/j.gaitpost.2012.08.010.View ArticlePubMedGoogle Scholar
- Hirsch BE, Udupa JK, Samarasekera S: New method of studying joint kinematics from three-dimensional reconstructions of MRI data. J Am Podiatr Med Assoc. 1996, 86: 4-15.View ArticlePubMedGoogle Scholar
- Glasoe WM: PhD thesis. Kinematic study of hallux valgus foot deformity. 2011, University of Minnesota, Department of KinesiologyGoogle Scholar
- Hamel AJ, Donahue SW, Sharkey NA: Contributions of active and passive toe flexion to forefoot loading. Clin Orthop. 2001, 393: 326-334.View ArticlePubMedGoogle Scholar
- Brown KM, Bursey DE, Arneson LJ, Andrews CA, Ludewig PM, Glasoe WM: Consideration of digitization precision when building local coordinate axes for a foot model. J Biomech. 2009, 42: 1263-1239. 10.1016/j.jbiomech.2009.03.013.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.