Automatic generation of personalised skeletal models of the lower limb from three-dimensional bone geometries.

The generation of personalised and patient-specific musculoskeletal models is currently a cumbersome and time-consuming task that normally requires several processing hours and trained operators. We believe that this aspect discourages the use of computational models even when appropriate data are available and personalised biomechanical analysis would be beneficial. In this paper we present a computational tool that enables the fully automatic generation of skeletal models of the lower limb from three-dimensional bone geometries, normally obtained by segmentation of medical images. This tool was evaluated against four manually created lower limb models finding remarkable agreement in the computed joint parameters, well within human operator repeatability. The coordinate systems origins were identified with maximum differences between 0.5 mm (hip joint) and 5.9 mm (subtalar joint), while the joint axes presented discrepancies between 1° (knee joint) to 11° (subtalar joint). To prove the robustness of the methodology, the models were built from four datasets including both genders, anatomies ranging from juvenile to elderly and bone geometries reconstructed from high-quality computed tomography as well as lower-quality magnetic resonance imaging scans. The entire workflow, implemented in MATLAB scripting language, executed in seconds and required no operator intervention, creating lower extremity models ready to use for kinematic and kinetic analysis or as baselines for more advanced musculoskeletal modelling approaches, of which we provide some practical examples. We auspicate that this technical advancement, together with upcoming progress in medical image segmentation techniques, will promote the use of personalised models in larger-scale studies than those hitherto undertaken.

Influence of material properties and boundary conditions on patient-specific models.

Patient-specific finite element models (PSFEM) are becoming more and more used. Different methods for assigning their material properties were studied on PSFEMs of 9 tibias along with the minimal required length of the CT acquisition window. Material properties are generally attributed to the PSFEM using relationships linking the grayscale of CT scans to the elasticity moduli. Using cortical-specific and trabecular-specific relationships or a generic one, did not result in significant differences. However, the use of homogeneous elastic moduli in the cortical and trabecular regions led to considerable differences. The result highlight that the PSFEM must comprise at least 40% of the tibia to ensure consistent results in the proximal 20%.

Tibial subchondral trabecular bone micromechanical and microarchitectural properties are affected by alignment and osteoarthritis stage.

At advanced knee osteoarthritis (OA) stages subchondral trabecular bone (STB) is altered. Lower limb alignment plays a role in OA progression and modify the macroscopic loading of the medial and lateral condyles of the tibial plateau. How the properties of the STB relate to alignment and OA stage is not well defined. OA stage (KL scores 2-4) and alignment (HKA from 17° Varus to 8° Valgus) of 30 patients were measured and their tibial plateau were collected after total knee arthroplasty. STB tissue elastic modulus, bone volume fraction (BV/TV) and trabecula thickness (Tb.Th) were evaluated with nanoindentation and µCT scans (8.1 µm voxel-size) of medial and lateral samples of each plateau. HKA and KL scores were statistically significantly associated with STB elastic modulus, BV/TV and Tb.Th. Medial to lateral BV/TV ratio correlated with HKA angle (R = -0.53, p = 0.016), revealing a higher ratio for varus than valgus subjects. STB properties showed lower values for KL stage 4 patients. Tissue elastic modulus ratios and BV.TV ratios were strongly correlated (R = 0.81, p < 0.001). Results showed that both micromechanical and microarchitectural properties of STB are affected by macroscopic loading at late stage knee OA. For the first time, a strong association between tissue stiffness and quantity of OA STB was demonstrated.

Articular-surface-based automatic anatomical coordinate systems for the knee bones.

Increasing use of patient-specific surgical procedures in orthopaedics means that patient-specific anatomical coordinate systems (ACSs) need to be determined. For knee bones, automatic algorithms constructing ACSs exist and are assumed to be more reliable than manual methods, although both approaches are based on non-unique numerical reconstructions of true bone geometries. Furthermore, determining the best algorithms is difficult, as algorithms are evaluated on different datasets. Thus, in this study, we developed 3 algorithms, each with 3 variants, and compared them with 5 from the literature on a dataset comprising 24 lower-limb CT-scans. To evaluate algorithms' sensitivity to the operator-dependent reconstruction procedure, the tibia, patella and femur of each CT-scan were each reconstructed once by three different operators. Our algorithms use principal inertia axis (PIA), cross-sectional area, surface normal orientations and curvature data to identify the bone region underneath articular surfaces (ASs). Then geometric primitives are fitted to ASs, and the ACSs are constructed from the geometric primitive points and/or axes. For each bone type, the algorithm displaying the least inter-operator variability is identified. The best femur algorithm fits a cylinder to posterior condyle ASs and a sphere to the femoral head, average axis deviations: 0.12°, position differences: 0.20 mm. The best patella algorithm identifies the AS PIAs, average axis deviations: 0.91°, position differences: 0.19 mm. The best tibia algorithm finds the ankle AS center and the 1st PIA of a layer around a plane fitted to condyle ASs, average axis deviations: 0.38°, position differences: 0.27 mm.

Three-dimensional deformation and transverse rotation of the human free Achilles tendon in vivo during isometric plantarflexion contraction.

Freehand three-dimensional ultrasound (3DUS) was used to investigate longitudinal and biaxial transverse deformation and rotation of the free Achilles tendon in vivo during a voluntary submaximal isometric muscle contraction. Participants (n = 8) were scanned at rest and during a 70% maximal voluntary isometric contraction (MVIC) of the plantarflexors. Ultrasound images were manually digitized to render a 3D reconstruction of the free Achilles tendon for the computation of tendon length, volume, cross-sectional area (CSA), mediolateral diameter (MLD), anteroposterior diameter (APD), and transverse rotation. Tendon longitudinal and transverse (CSA, APD, and MLD) deformation and strain at 70% MVIC were calculated relative to the resting condition. There was a significant main effect of contraction on tendon length and mean CSA, MLD, and APD (P < 0.05), but no effect on tendon volume (P = 0.70). Group mean transverse strains for CSA, MLD, and APD averaged over the length of the tendon were -5.5%, -8.7% and 8.7%, respectively. Peak CSA, MLD, and APD transverse strains all occurred between 40% and 60% of tendon length. Transverse rotation of the free tendon was negligible at rest but increased under load, becoming externally rotated relative to the calcaneal insertion. The relationship between longitudinal and transverse strains of the free Achilles tendon during muscle-induced elongation may be indicative of interfascicle reorganization. The finding that transverse rotation and strain peaked in midportion of the free Achilles tendon may have important implications for tendon injury mechanisms and estimation of tendon stress in vivo.