- Technical note
- Open Access
Three-dimensional kinematic change of hindfoot during full weightbearing in standing: an analysis using upright computed tomography and 3D-3D surface registration
Journal of Orthopaedic Surgery and Research volume 14, Article number: 355 (2019)
Weightbearing of the hindfoot affects positional changes of the ankle joint and subtalar joint (ankle-joint complex [AJC]). However, it is difficult to assess the kinematic changes in the hindfoot in a natural full weightbearing condition using conventional CT or cone beam computed tomography (CT) due to limitations of acquiring foot images under a physiological weightbearing condition using those imaging modalities. Analysis of AJC kinematics using fluoroscopy and 2D-3D registration technique requires data on the number of steps and amount of time to build and match the bones. This study aimed to analyze the effect of full weightbearing on hindfoot motion when standing using upright CT and 3D-3D surface registration.
Forty-eight AJCs of 24 asymptomatic volunteers (13 women, 11 men) were examined under no weightbearing, 50% weightbearing, and single leg full weightbearing conditions while standing. The CT images were acquired from the distal femur to the whole foot using a 320-row upright CT scanner. The condition of each weightbearing stance was measured using a pressure mat. Bone-to-bone rotations of the talus relative to the tibia and calcaneus relative to the talus were evaluated using the surface registration technique. Image quality of the CT and intra- and interobserver reliabilities of the rotation angle were also evaluated.
All CT images were excellent or good quality and the intra- and interobserver correlation coefficients for the angle were 0.996 and 0.995, respectively. The motion of the ankle joint and subtalar joint under 50% and 100% weightbearing were as follows (in degrees); the talus plantarflexed (5.1 ± 4.5 and 6.8 ± 4.8), inverted (1.3 ± 1.4 and 2.0 ± 1.6), and internally rotated (2.4 ± 4.2 and 4.3 ± 4.6) relative to the tibia, and the calcaneus dorsiflexed (2.8 ± 1.4 and 3.8 ± 1.7), everted (5.3 ± 2.6 and 8.0 ± 3.6), and externally rotated (3.0 ± 2.0 and 4.1 ± 2.4) relative to the talus, respectively.
The effect of weightbearing was clearly identified using an upright CT and the 3D-3D registration technique. Three-dimensional kinematics under static full weightbearing were opposite between the ankle and subtalar joints on their respective axes.
The joint of the hindfoot, i.e., the ankle-joint complex (AJC), consists of the ankle joint and the subtalar joint. One of the major functions of the AJC is adjusting the lower limb alignment while weightbearing. It is difficult to assess the kinematic change of the AJC under the weightbearing condition using 2D radiography because of the complex bony shapes of the joint. Thus, the 2D-3D registration technique was developed to measure foot bone and AJC kinematics of patients [1, 2] or asymptomatic volunteers [3,4,5,6]. Although the accuracy of the 2D-3D registration technique is below 1.0 mm and 1.0°, the number of steps is required to achieve accuracy. First, computed tomography (CT) scans are crucial to build the bone model. Second, 2D fluoroscopic images require calibration to adjust enlarged images projected from radiographs. Finally, a matching algorithm involves a number of mathematical calculations and optimization. Therefore, a method to analyze the kinematics of AJC without such time and cost is needed.
Many studies have analyzed the effect of weightbearing on the hindfoot using magnetic resonance imaging  or conventional CT with loading devices [7,8,9,10,11,12,13,14,15] or upright cone beam CT [16,17,18,19,20,21,22,23]. However, the effect of natural full weightbearing in a standing position has not been evaluated due to limitations of acquiring foot images under a physiological weightbearing condition using those imaging modalities. Major limitations of these modalities include low resolution and longer scanning time (MRI), prone position and non-physiological weightbearing (conventional CT), and motion artifact and partial weightbearing (cone beam CT). Approximately 1 to 2° of bony motion were observed in the AJC under weightbearing [10, 11, 15], while no study has reported the kinematic change in the AJC between non weightbearing and full weightbearing positions. Recently, we developed an upright CT with an area detector with Canon Medical Systems , in which a CT scan under full weightbearing in a natural standing condition can be acquired.
The present study aimed to analyze the effect of full weightbearing on hindfoot motion of asymptomatic feet using an upright CT with a 320-row multidetector and 3D-3D surface registration technique. We hypothesized that the upright CT and 3D-3D registration technique clearly reveals kinematic change in the AJC due to natural full weightbearing.
Materials and methods
A total of 48 AJCs of 24 healthy volunteers (13 women, 11 men) with no history of a foot or ankle injury and no obvious foot deformations were enrolled in the present study. The mean (± standard deviation) age, body weight, and body mass index (BMI) of the participants were 28.3 ± 4.0 (range, 23–39) years, 59.0 ± 9.7 (range, 45.0–78.0) kg, and 21.4 ± 2.0 (range, 17.6–26.4) kg/m2, respectively. Each participant provided written informed consent, and the study protocol (ID#20150293) was approved by our ethical committee.
The CT images were acquired from the distal femur to the entire foot using the 320-row upright CT scanner (prototype TSX-401R; Canon Medical Systems, Otawara, Japan) (Fig. 1) . The CT examinations were performed using the following parameters: peak tube voltage, 100 kV; tube current, 10 to 350 mA (using a noise index of 15 for a slice thickness of 5 mm); rotation speed, 0.5 s; and slice thickness, 0.5 mm. No weightbearing, standing (50% weightbearing), and single leg full weightbearing (100% weightbearing) were assessed for each participant . The condition of each weightbearing stance was measured using a pressure mat (BIG-MAT; NITTA Corporation, Osaka, Japan) and pressure calculation system (FootMat; Tekscan, South Boston, MA, USA). In an upright CT scanner, all participants stood in a relaxed position with their bare feet shoulder width apart. Similar to the previous report, 2-kg weightbearing with the ankle in a neutral position was defined as the no weightbearing condition in the present study , while the opposite side was defined as full weightbearing. The CT data were accumulated using the Digital Imaging and Communication in Medicine (DICOM) data format. Two orthopedic surgeons with 8 and 19 years of experience diagnosing musculoskeletal CT images independently evaluated the image quality of all the images using the nine-level Likert scale, as described in the previous study : a score of 5 indicates “excellent” image quality without any artifacts; score of 4 indicates “good” image quality with minor artifacts; score of 3 indicates “fair” image quality with moderate artifacts; score of 2 indicates “poor” image quality with non-diagnostic quality and the visualization task can be identified; score of 1 indicates “very poor” image quality with non-diagnostic quality and the visualization task cannot be identified. A half score was added between each score and a total of nine grades were evaluated.
3D-3D surface registration
Three-dimensional surface data of the tibia, talus, and calcaneus were extracted from the CT DICOM data using 3D visualization software (AVIZO 6.4; Thermo Fisher Scientific, Tokyo, Japan). We matched the 3D surface of the talus in each weightbearing condition for each participant using the iterative closest point algorithm using Visualization Toolkit 8.1.0 (Kitware Inc., Clifton Park, NY, USA) for the 3D surface registration technique in which point data are superimposed onto another 3D surface by iterative steps to reach the closest points.
We modified and used the coordinate system of the tibia using the method defined by Sato et al.  and the International Society of Biomechanics . The coordinate systems of the talus and the calcaneus were defined using the method described by Gutekunst et al.  (Fig. 2).
The tibia coordinate system consisted of the following: the line connecting the center of the most medial point on the border of the medial tibial condyle (MC) and the most lateral point on the border of the lateral tibial condyle (LC) and the center of the tibia plafond pointing upward was defined as the Y axis. The line connecting the MC and LC was defined as the T (temporary) axis. The Z axis was the line perpendicular to the Y and T axes pointing laterally. The X axis was defined as the line perpendicular to the Y and Z axis pointing anteriorly.
The talus coordinate system consisted of the following: the line connecting the mediolateral midline of the posterior aspect of the talus and the center of the convex surface of the talar head centered both mediolaterally and vertically and pointing anteriorly was defined as the X axis. The line connecting the dorsal maximum of the medial edge of talar trochlea articular surface and the dorsal maximum of the lateral surface was defined as the T axis. The Y axis was the line perpendicular to the T and X axes pointing upward. The Z axis was defined as the line perpendicular to the X and Y axes pointing laterally.
The calcaneus coordinate system consisted of the following: the line connecting the midpoint of the posterior surface of the calcaneal tuberosity, centered both mediolaterally and vertically, and the center of the anterior surface of the calcaneus where it articulates with the cuboid pointing anteriorly was defined as the X axis. The line connecting the mediolateral midline of the posterior surface of the calcaneal tuberosity along the inferior border and the mediolateral midline of the posterior surface of the calcaneal tuberosity along the superior border was defined as the T axis. The Z axis was the line perpendicular to the T and X axes pointing upward. The Y axis was defined as the line perpendicular to the X and Y axes pointing laterally.
Analysis of joint motion
The ankle (talocrural) joint: the articulation formed between the talus and the tibia/fibula. The subtalar joint: the articulation formed between the talus and the calcaneus.
Rotation around the X axis was defined as inversion/eversion, rotation around the Y axis was defined as internal/external rotation, and rotation around the Z axis was defined as dorsiflexion/plantarflexion. Bone-to-bone rotations of the talus relative to the tibia and the calcaneus relative to the talus around each axis are described by the Euler/Cardan angles representing three sequential rotations about the anatomical axis of the proximal bone. The rotation sequence “Z-X-Y” was used.
The intra- and interobserver reliabilities were independently assessed by two orthopedic specialists and by reassessment of the data with an interval of longer than 2 months, respectively. The correlation coefficients were calculated to assess the intra- and interobserver reliabilities using SPSS ver. 24.0 (IBM, Armonk, NY, USA).
Image qualities of the 144 AJC scans of 24 subjects were good (diagnostic quality with minor artifacts) or excellent (diagnostic quality without any artifacts)  (Table 1). The intra- and interobserver correlation coefficients for the present study were 0.996 (95% confidence interval, 0.994–0.998) and 0.995 (95% confidence interval, 0.992–0.997). These data indicated that the present measurement was highly reliable.
In the ankle joint, the talus plantarflexed (50%/100% weightbearing, 5.07 ± 4.52/6.77 ± 4.84 degrees), inverted (50%/100% weightbearing, 1.28 ± 1.37/2.01 ± 1.58 degrees), and internally rotated (50%/100% weightbearing, 2.40 ± 4.18/4.30 ± 4.64°) relative to the tibia as the weight load increased. Conversely, at the subtalar joint, the calcaneus dorsiflexed (50%/100% weightbearing, 2.76 ± 1.42/3.82 ± 1.68°), everted (50%/100% weightbearing, 5.29 ± 2.56/7.99 ± 3.55°), and externally rotated (50%/100% weightbearing, 2.96 ± 1.95/4.13 ± 2.43°) relative to the talus as the weight load increased (Figs. 3 and 4). Three-dimensional kinematics were opposite between the ankle joint and the subtalar joint on their respective axes, and each angle increased as the weight load increased. Regarding the absolute value, sagittal and axial plane movements were larger in the ankle joint, while the coronal plane movement was larger in the subtalar joint.
Our approach using the upright CT and 3D-3D registration technique clearly described the effect of full weightbearing in AJC kinematics, and the results support our hypothesis.
Our method to analyze AJC kinematics has several advantages over the methods using fluoroscopy or other imaging modalities (Table 3). First, 3D-3D registration on CT images requires fewer steps to match the bone and evaluate AJC kinematics and it is easier to match 3D to 3D than 2D to 3D models. Analysis of foot bone and AJC kinematics using fluoroscopy and the 2D-3D registration technique has been reported [1,2,3,4,5,6]; however, its major limitation is the complex nature of the steps required to build and match the bones. The 2D images taken by fluoroscopic imaging are shadow pictures, and a 3D bone model based on CT images is required to accurately match the bones on the 2D images. Image calibration is also required to adjust enlarged images when using the X-ray system. Several matching algorithms have been developed, but the significant time and cost required to analyze the kinematics of the bones limit its use. The accuracy of the 3D-3D registration was below 0.2° in rotation . Second, only minor motion artifacts were found in AJC images with upright CT in the present study (Table 1). Changes in hindfoot alignment have been assessed using upright cone beam CT [16,17,18,19,20,21,22,23], but it takes as long as 20 to 48 s to acquire images, and it is necessary for participants to support the body to reduce artifacts. In fact, moderate to severe motion artifacts were observed in the cone beam CT images of the knee and ankle . In addition, participants must put their foot in a small tube of the cone beam CT, and thus the participants must set their contralateral foot somewhere aside from the tube or stabilize their body using supportive tools such as a pole. This position is not a natural standing position, and only partial weight is loaded on the foot. Third, physiological weightbearing while standing can be acquired in the upright CT, while simulated weight with loading devices was applied in the studies using conventional CT [7,8,9,10,11,12,13,14,15]. In those studies, the hip, shoulder, or knee must be fixed to reproduce the hypothetical loading conditions, and the lower limb muscles used to maintain the standing position was not active in the prone position. Those limit the representation of physiological loading and tarsal bone alignment while standing.
Several limitations of the present study should be noted. First, there were no patient data, and only asymptomatic subjects were included. However, our method using an upright CT and 3D-3D registration technique can be a powerful tool to investigate kinematic change in the AJC of the patients. The clinical relevance of the hindfoot motion during natural full weightbearing should be studied in the near future. Second, the imaging was divided into three categories, i.e., no weightbearing, 50% weightbearing, and full weightbearing, and static imaging was performed. Although continuous imaging in 4D was possible using an upright CT with 320-row multidetector, the image quality of 4D CT was insufficient to capture the tarsal bones; thus, we separately scanned the three loading conditions. To analyze the continuous dynamics of the hindfoot, we need to increase the observation points under different weightbearing conditions in a future study.
An upright CT and 3D-3D registration technique clearly described the kinematics of the AJC in a static full weightbearing condition. Our findings demonstrated that 3D motions were opposite between the ankle and subtalar joints on their respective axes.
Availability of data and materials
The datasets of the present study are available from the corresponding author on reasonable request.
Digital Imaging and Communication in Medicine
Lateral tibial condyle
Medial tibial condyle
Goto A, Moritomo H, Itohara T, Watanabe T, Sugamoto K. Three-dimensional in vivo kinematics of the subtalar joint during dorsi-plantarflexion and inversion-eversion. Foot Ankle Int. 2009;30(5):432–8.
Yamaguchi S, Sasho T, Kato H, Kuroyanagi Y, Banks SA. Ankle and subtalar kinematics during dorsiflexion-plantarflexion activities. Foot Ankle Int. 2009;30(4):361–6.
Ito K, Hosoda K, Shimizu M, Ikemoto S, Kume S, Nagura T, et al. Direct assessment of 3D foot bone kinematics using biplanar X-ray fluoroscopy and an automatic model registration method. J Foot Ankle Res. 2015;8:21.
Ito K, Hosoda K, Shimizu M, Ikemoto S, Nagura T, Seki H, et al. Three-dimensional innate mobility of the human foot bones under axial loading using biplane X-ray fluoroscopy. R Soc Open Sci. 2017;4(10):171086.
Kobayashi T, Saka M, Suzuki E, Yamazaki N, Suzukawa M, Akaike A, et al. In vivo kinematics of the talocrural and subtalar joints during weightbearing ankle rotation in chronic ankle instability. Foot Ankle Spec. 2014;7(1):13–9.
Roach KE, Foreman KB, Barg A, Saltzman CL, Anderson AE. Application of high-speed dual fluoroscopy to study in vivo tibiotalar and subtalar kinematics in patients with chronic ankle instability and asymptomatic control subjects during dynamic activities. Foot Ankle Int. 2017;38(11):1236–48.
Ananthakrisnan D, Ching R, Tencer A, Hansen ST Jr, Sangeorzan BJ. Subluxation of the talocalcaneal joint in adults who have symptomatic flatfoot. J Bone Joint Surg Am. 1999;81(8):1147–54.
Apostle KL, Coleman NW, Sangeorzan BJ. Subtalar joint axis in patients with symptomatic peritalar subluxation compared to normal controls. Foot Ankle Int. 2014;35(11):1153–8.
Ferri M, Scharfenberger AV, Goplen G, Daniels TR, Pearce D. Weightbearing CT scan of severe flexible pes planus deformities. Foot Ankle Int. 2008;29(2):199–204.
Kido M, Ikoma K, Imai K, Maki M, Takatori R, Tokunaga D, et al. Load response of the tarsal bones in patients with flatfoot deformity: in vivo 3D study. Foot Ankle Int. 2011;32(11):1017–22.
Kido M, Ikoma K, Imai K, Tokunaga D, Inoue N, Kubo T. Load response of the medial longitudinal arch in patients with flatfoot deformity: in vivo 3D study. Clin Biomech (Bristol, Avon). 2013;28(5):568–73.
Ledoux WR, Rohr ES, Ching RP, Sangeorzan BJ. Effect of foot shape on the three-dimensional position of foot bones. J Orthop Res. 2006;24(12):2176–86.
Van Bergeyk AB, Younger A, Carson B. CT analysis of hindfoot alignment in chronic lateral ankle instability. Foot Ankle Int. 2002;23(1):37–42.
Yoshioka N, Ikoma K, Kido M, Imai K, Maki M, Arai Y, et al. Weight-bearing three-dimensional computed tomography analysis of the forefoot in patients with flatfoot deformity. J Orthop Sci. 2016;21(2):154–8.
Zhang Y, Xu J, Wang X, Huang J, Zhang C, Chen L, et al. An in vivo study of hindfoot 3D kinetics in stage II posterior tibial tendon dysfunction (PTTD) flatfoot based on weight-bearing CT scan. Bone Joint Res. 2013;2(12):255–63.
Burssens A, Peeters J, Buedts K, Victor J, Vandeputte G. Measuring hindfoot alignment in weight bearing CT: a novel clinical relevant measurement method. Foot Ankle Surg. 2016;22(4):233–8.
Cody EA, Williamson ER, Burket JC, Deland JT, Ellis SJ. Correlation of talar anatomy and subtalar joint alignment on weightbearing computed tomography with radiographic flatfoot parameters. Foot Ankle Int. 2016;37(8):874–81.
Hirschmann A, Pfirrmann CW, Klammer G, Espinosa N, Buck FM. Upright cone CT of the hindfoot: comparison of the non-weight-bearing with the upright weight-bearing position. Eur Radiol. 2014;24(3):553–8.
Krahenbuhl N, Tschuck M, Bolliger L, Hintermann B, Knupp M. Orientation of the subtalar joint: measurement and reliability using weightbearing CT scans. Foot Ankle Int. 2016;37(1):109–14.
Lintz F, Welck M, Bernasconi A, Thornton J, Cullen NP, Singh D, et al. 3D biometrics for hindfoot alignment using weightbearing CT. Foot Ankle Int. 2017;38(6):684–9.
Richter M, Lintz F, Zech S, Meissner SA. Combination of PedCAT weightbearing CT with pedography assessment of the relationship between anatomy-based foot center and force/pressure-based center of gravity. Foot Ankle Int. 2018;39(3):361–8.
Richter M, Seidl B, Zech S, Hahn S. PedCAT for 3D-imaging in standing position allows for more accurate bone position (angle) measurement than radiographs or CT. Foot Ankle Surg. 2014;20(3):201–7.
Richter M, Zech S, Hahn S, Naef I, Merschin D. Combination of pedCAT(R) for 3D imaging in standing position with pedography shows no statistical correlation of bone position with force/pressure distribution. J Foot Ankle Surg. 2016;55(2):240–6.
Jinzaki M, Yamada Y, Nagura T, Nakahara T, Yokoyama Y, Narita K, et al. Development of upright CT with area detector for whole body scans: phantom study, efficacy on workflow, effect of gravity on human body, and potential clinical impact. Investig Radiol. in press
Ota T, Nagura T, Yamada Y, Yamada M, Yokoyama Y, Ogihara N, et al. Effect of natural full weight-bearing during standing on the rotation of the first metatarsal bone. Clin Anat. 2019;32(5):715–21.
Watanabe K, Ikeda Y, Suzuki D, Teramoto A, Kobayashi T, Suzuki T, et al. Three-dimensional analysis of tarsal bone response to axial loading in patients with hallux valgus and normal feet. Clin Biomech (Bristol, Avon). 2017;42:65–9.
Sisniega A, Thawait GK, Shakoor D, Siewerdsen JH, Demehri S, Zbijewski W. Motion compensation in extremity cone-beam computedtomography. Skeletal Radiol. 2019;48(12):1999–2007.
Sato T, Koga Y, Sobue T, Omori G, Tanabe Y, Sakamoto M. Quantitative 3-dimensional analysis of preoperative and postoperative joint lines in total knee arthroplasty: a new concept for evaluation of component alignment. J Arthroplast. 2007;22(4):560–8.
Wu G, Siegler S, Allard P, Kirtley C, Leardini A, Rosenbaum D, et al. ISB recommendation on definitions of joint coordinate system of various joints for the reporting of human joint motion--part I: ankle, hip, and spine. International Society of Biomechanics. J Biomech. 2002;35(4):543–8.
Gutekunst DJ, Liu L, Ju T, Prior FW, Sinacore DR. Reliability of clinically relevant 3D foot bone angles from quantitative computed tomography. J Foot Ankle Res. 2013;6:38.
Ochia RS, Inoue N, Renner SM, Lorenz EP, Lim TH, Andersson GB, et al. Three-dimensional in vivo measurement of lumbar spine segmental motion. Spine (Phila Pa 1976). 2006;31(18):2073–8.
The authors would like to thank Yu Iwama, Shu Kobayashi, and Yasuo Niki from the Department of Orthopedic Surgery, School of Medicine, for clinical advice. Further support was provided by Sumi Yamashita and Hiroko Arai of the Department of Orthopedic Surgery, School of Medicine, Keio University.
The following funding was received for this study: Uehara Memorial Foundation (N/A), recipient: Dr. Masahiro Jinzaki; Ministry of Health, Labour and Welfare (17H04266), recipient: Dr. Masahiro Jinzaki.
Ethics approval and consent to participate
The present study was approved by the Institutional Review Board of the School of Medicine, Keio University (ID#20150293), and informed consent was obtained from all participants.
Consent for publication
The approval of the Institutional Review Board included the consent for publication of any individual data.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kaneda, K., Harato, K., Oki, S. et al. Three-dimensional kinematic change of hindfoot during full weightbearing in standing: an analysis using upright computed tomography and 3D-3D surface registration. J Orthop Surg Res 14, 355 (2019) doi:10.1186/s13018-019-1443-z
- Upright computed tomography
- Surface registration