Multi-rigid image segmentation and registration for the analysis of joint motion from three-dimensional magnetic resonance imaging

Subject-specific material properties of the heel pad: An inverse finite element analysis

Individuals with diabetes are at a higher risk of developing foot ulcers. To better understand internal soft tissue loading and potential treatment options, subject-specific finite element (FE) foot models have been used. However, existing models typically lack subject-specific soft tissue material properties and only utilize subject-specific anatomy. Therefore, this study determined subject-specific hindfoot soft tissue material properties from one non-diabetic and one diabetic subject using inverse FE analysis. Each subject underwent cyclic MRI experiments to simulate physiological gait and to obtain compressive force and three-dimensional soft tissue imaging data at 16 phases along the loading-unloading cycles. The FE models consisted of rigid bones and nearly-incompressible first-order Ogden hyperelastic skin, fat, and muscle (resulting in six independent material parameters). Then, calcaneus and loading platen kinematics were computed from imaging data and prescribed to the FE model. Two analyses were performed for each subject. First, the skin, fat, and muscle layers were lumped into a single generic soft tissue material and optimized to the platen force. Second, the skin, fat, and muscle material properties were individually determined by simultaneously optimizing for platen force, muscle vertical displacement, and skin mediolateral bulging. Our results indicated that compared to the individual without diabetes, the individual with diabetes had stiffer generic soft tissue behavior at high strain and that the only substantially stiffer multi-material layer was fat tissue. Thus, we suggest that this protocol serves as a guideline for exploring differences in non-diabetic and diabetic soft tissue material properties in a larger population.

A Cadaveric Comparison of the Kinematic and Anatomical Axes and Arthrokinematics of the Metatarsosesamoidal and First Metatarsophalangeal Joints

Presently, developments in weightbearing computed tomography and biplanar fluoroscopy technologies offer exciting avenues for investigating normative and pathologic foot function with increasing precision. Still, data quantifying sesamoid bone and proximal phalange motion are currently sparse. To express joint kinematics and compare various clinical cohorts, future studies of first ray motion will necessitate robust coordinate frames that respect the variations in underlying anatomy while also aligning closely with the functional, physiological axes of motion. These activity-dependent functional axes may be represented by a mean helical axis of the joint motion. Our cadaveric study quantified joint kinematics from weightbearing computed tomography scans during simulated toe lift and heel rise tasks. We compared the spatial orientations of the mean finite helical axes of the metatarsosesamoidal and metatarsophalangeal joints to the primary joint axis of two relevant methods for defining metatarsal coordinate frames: inertial axes and fitting of geometric primitives. The resultant kinematics exhibited less crosstalk when using a metatarsal coordinate system based on fitting cylindrical primitives to the bony anatomy compared to using principal component axes. Respective metatarsophalangeal and metatarsosesamoidal arthrokinematic contact paths and instantaneous centers of rotation were similar between activities and agree well with currently published data. This study outlines a methodology for quantitatively assessing the efficacy and utility of various anatomical joint coordinate system definitions. Improvements in our ability to characterize the shape and motion of foot bones in the context of functional tasks will elucidate their biomechanical roles and aid clinicians in refining treatment strategies.

Hind- and midfoot bone morphology varies with foot type and sex

Foot type has been associated with pain, injury, and altered gait mechanics. Morphological variations in foot bones due to foot type variation may impact surgical and therapeutic treatments. The purpose of this study was to utilize principal component analysis (PCA) to determine how morphology of the hind- and midfoot bones differs among foot types and sex. The calcaneus, talus, navicular, and cuboid were segmented using previously obtained computed tomography (CT) scans and converted to surface models. The CTs were sorted into four foot types-cavus, neutrally aligned, asymptomatic planus, and symptomatic planus. Morphometric shape analysis software (Geomorph) was used to perform a PCA to determine which components varied between foot types and between sexes. The calcaneus showed planus feet of both types to have calcanei that have decreased height and increased length compared to neutrally aligned feet. The talus demonstrated increased posterior mass for cavus feet compared to neutrally aligned feet. For the navicular, symptomatic planus had a more posteriorly positioned tuberosity and were wider than asymptomatic planus feet. The cuboid did not exhibit any differences between foot types. Sex differences, found only at the talus and navicular, were subtle. PCA is an objective technique that helped elucidate differences in bone morphology between foot types and sex without needing to determine the features of interest before comparing groups. Understanding these variations can help inform diagnosis of foot pathologies and surgical protocols as well as improve computer models of the foot. Published 2018. This article is a U.S. Government work and is in the public domain in the USA. J Orthop Res 9999:1-16, 2019.

Model-based tracking of the bones of the foot: A biplane fluoroscopy validation study

Measuring foot kinematics using optical motion capture is technically challenging due to the depth of the talus, small bone size, and soft tissue artifact. We present a validation of our biplane X-ray system, demonstrating its accuracy in tracking the foot bones directly. Using an experimental linear/rotary stage we imaged pairs of tali, calcanei, and first metatarsals, with embedded beads, through 30 poses. Model- and bead-based algorithms were employed for semi-automatic tracking. Translational and rotational poses were compared to the experimental stage (a reference standard) to determine registration performance. For each bone, 10 frames per pose were analyzed. Model-based: The resulting overall translational bias of the six bones was 0.058 mm with a precision of ± 0.049 mm. The overall rotational bias of the six bones was 0.42° with a precision of ± 0.41°. Bead-based: the overall translational bias was 0.037 mm with a precision of ± 0.032 mm and for rotation was 0.29° with a precision of ± 0.26°. We validated the accuracy of our system to determine the spatial position and orientation of isolated foot bones, including the talus, calcaneus, and first metatarsal over a range of quasi-static poses. Although the accuracy of dynamic motion was not assessed, use of an experimental stage establishes a reference standard.

Strategies towards rapid generation of forefoot model incorporating realistic geometry of metatarsals encapsulated into lumped soft tissues for personalized finite element analysis

Use of finite element (FE) foot model as a clinical diagnostics tool is likely to improve the specificity of foot injury predictions in the diabetic population. Here we proposed a novel workflow for rapid construction of foot FE model incorporating realistic geometry of metatarsals encapsulated into lumped forefoot's soft tissues. Custom algorithms were implemented to perform unsupervised segmentation and mesh generation to directly convert CT data into a usable FE model. The automatically generated model provided higher efficiency and comparable numerical accuracy when compared to the model constructed using a traditional solid-based mesh process. The entire procedure uses MATLAB as the main platform, and makes the present approach attractive for creating personalized foot models to be used in clinical studies.

Metatarsal Shape and Foot Type: A Geometric Morphometric Analysis

Planus and cavus foot types have been associated with an increased risk of pain and disability. Improving our understanding of the geometric differences between bones in different foot types may provide insights into injury risk profiles and have implications for the design of musculoskeletal and finite-element models. In this study, we performed a geometric morphometric analysis on the geometry of metatarsal bones from 65 feet, segmented from computed tomography (CT) scans. These were categorized into four foot types: pes cavus, neutrally aligned, asymptomatic pes planus, and symptomatic pes planus. Generalized procrustes analysis (GPA) followed by permutation tests was used to determine significant shape differences associated with foot type and sex, and principal component analysis was used to find the modes of variation for each metatarsal. Significant shape differences were found between foot types for all the metatarsals (p < 0.01), most notably in the case of the second metatarsal which showed significant pairwise differences across all the foot types. Analysis of the principal components of variation showed pes cavus bones to have reduced cross-sectional areas in the sagittal and frontal planes. The first (p = 0.02) and fourth metatarsals (p = 0.003) were found to have significant sex-based differences, with first metatarsals from females shown to have reduced width, and fourth metatarsals from females shown to have reduced frontal and sagittal plane cross-sectional areas. Overall, these findings suggest that metatarsal bones have distinct morphological characteristics that are associated with foot type and sex, with implications for our understanding of anatomy and numerical modeling of the foot.

Talonavicular joint coverage and bone morphology between different foot types

This study explored three dimensional (3D) talonavicular joint (TNJ) coverage/orientation and bone morphology to reveal parameters that could classify and identify predispositions to cavus and planus feet. 3D models of 65 feet from 40 subjects were generated from computed tomography images classified as pes cavus, neutrally aligned, or asymptomatic/symptomatic pes planus. We calculated the talar and navicular overlap (TNJ coverage). We also measured orientation of the navicular, morphological parameters of the talus and navicular, and angular position of the talar head to body. Pes cavus showed significantly less talonavicular coverage (58 ± 2% talus and 86 ± 2% navicular) compared to asymptomatic pes planus (63 ± 2% and 95 ± 2%) and neutrally aligned feet (98 ± 2% navicular), and significantly more navicular dorsiflexion and adduction relative to the talus (p < 0.0083). The talar head in cavus feet was inverted relative to the body compared to planus feet (p < 0.0083). For symptomatic pes planus, significant abduction was measured for the navicular relative to the talus and the talar head was plantar flexed relative to the body (p < 0.0083). The talar head in planus feet was everted relative to the body compared to neutrally aligned feet. Both intrinsic (bone morphology) and extrinsic (bone position) differences exist in groups of feet described as cavus and planus.

Finite element analysis of the foot: model validation and comparison between two common treatments of the clawed hallux deformity

BACKGROUND

Clawed hallux is defined by first metatarsophalangeal joint extension and first interphalangeal joint flexion; it can increase plantar pressures and ulceration risk. We investigated two corrective surgical techniques, the modified Jones and flexor hallucis longus tendon transfer.

METHODS

A finite element foot model was modified to generate muscle overpulls, including extensor hallucis longus, flexor hallucis longus and peroneus longus. Both corrective procedures were simulated, predicting joint angle and plantar pressure changes.

FINDINGS

The clawed hallux deformity was generated by overpulling: 1) extensor hallucis longus, 2) peroneus longus + extensor hallucis longus, 3) extensor hallucis longus + flexor hallucis longus and 4) all three together. The modified Jones reduced metatarsophalangeal joint angles, but acceptable hallux pressure was found only when there was no flexor hallucis longus overpull. The flexor hallucis longus tendon transfer reduced deformity at the metatarsophalangeal and interphalangeal joints but may extended the hallux due to the unopposed extensor hallucis longus. Additionally, metatarsal head pressure increased with overpulling of the extensor hallucis longus + flexor hallucis longus, and all three muscles together.

INTERPRETATION

The modified Jones was effective in correcting clawed hallux deformity involving extensor hallucis longus overpull without flexor hallucis longus overpull. The flexor hallucis longus tendon transfer was effective in correcting clawed hallux deformity resulting from the combined overpull of both extensor and flexor hallucis longus, but not with isolated extensor hallucis longus overpull. An additional procedure to reduce the metatarsal head pressure may be required concomitant to the flexor hallucis longus tendon transfer. However this procedure avoids interphalangeal joint fusion.

Evaluating foot kinematics using magnetic resonance imaging: from maximum plantar flexion, inversion, and internal rotation to maximum dorsiflexion, eversion, and external rotation

The foot consists of many small bones with complicated joints that guide and limit motion. A variety of invasive and noninvasive means [mechanical, X-ray stereophotogrammetry, electromagnetic sensors, retro-reflective motion analysis, computer tomography (CT), and magnetic resonance imaging (MRI)] have been used to quantify foot bone motion. In the current study we used a foot plate with an electromagnetic sensor to determine an individual subject's foot end range of motion (ROM) from maximum plantar flexion, internal rotation, and inversion to maximum plantar flexion, inversion, and internal rotation to maximum dorsiflexion, eversion, and external rotation. We then used a custom built MRI-compatible device to hold each subject's foot during scanning in eight unique positions determined from the end ROM data. The scan data were processed using software that allowed the bones to be segmented with the foot in the neutral position and the bones in the other seven positions to be registered to their base positions with minimal user intervention. Bone to bone motion was quantified using finite helical axes (FHA). FHA for the talocrural, talocalcaneal, and talonavicular joints compared well to published studies, which used a variety of technologies and input motions. This study describes a method for quantifying foot bone motion from maximum plantar flexion, inversion, and internal rotation to maximum dorsiflexion, eversion, and external rotation with relatively little user processing time.