Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Assessing Model Calibration Workflows and Outcomes.

Model reproducibility is a point of emphasis for the National Institutes of Health (NIH) and in science, broadly. As the use of computational modeling in biomechanics and orthopedics grows, so does the need to assess the reproducibility of modeling workflows and simulation predictions. The long-term goal of the KneeHub project is to understand the influence of potentially subjective decisions, thus the modeler's "art", on the reproducibility and predictive uncertainty of computational knee joint models. In this paper, we report on the model calibration phase of this project, during which five teams calibrated computational knee joint models of the same specimens from the same specimen-specific joint mechanics dataset. We investigated model calibration approaches and decisions, and compared calibration workflows and model outcomes among the teams. The selection of the calibration targets used in the calibration workflow differed greatly between the teams and was influenced by modeling decisions related to the representation of structures, and considerations for computational cost and implementation of optimization. While calibration improved model performance, differences in the postcalibration ligament properties and predicted kinematics were quantified and discussed in the context of modeling decisions. Even for teams with demonstrated expertise, model calibration is difficult to foresee and plan in detail, and the results of this study underscore the importance of identification and standardization of best practices for data sharing and calibration.

3D trunk orientation measured using inertial measurement units during anatomical and dynamic sports motions.

Trunk motion is related to the performance and risk of injuries during dynamic sports motions. Optical motion capture is traditionally used to measure trunk motion during dynamic sports motions, but these systems are typically constrained to a laboratory environment. Inertial measurement units (IMUs) might provide a suitable alternative for measuring the trunk orientation during dynamic sports motions. The objective of the present study was to assess the accuracy of the three-dimensional trunk orientation measured using IMUs during dynamic sports motions and isolated anatomical trunk motions. The motions were recorded with two IMUs and an optical motion capture system (gold standard). Ten participants performed a total of 71 sports motions (19 golf swings, 15 one-handed ball throws, 19 tennis serves, and 18 baseball swings) and 125 anatomical trunk motions (42, 41, and 42 trials of lateral flexion, axial rotation, and flexion/extension, respectively). The root-mean-square differences between the IMU- and optical motion capture-based trunk angles were less than 5 degrees, and the similarity between the methods was on average across all trials "very good" to "excellent" (R ≥ 0.85; R2 ≥ 0.80). Across the dynamic sports motions, even higher measures of similarity were found (R ≥ 0.90; R2 ≥ 0.82). When aligned to the relevant segment, the current IMUs are a promising alternative to optical motion capture and previous presented IMU-based systems for the field-based measurement of the three-dimensional trunk orientation during dynamic sports motions and the anatomical trunk motions.

A Method to Compare Heterogeneous Types of Bone and Cartilage Meshes.

Accurately capturing the bone and cartilage morphology and generating a mesh remains a critical step in the workflow of computational knee joint modeling. Currently, there is no standardized method to compare meshes of different element types and nodal densities, making comparisons across research teams a significant challenge. The aim of this paper is to describe a method to quantify differences in knee joint bone and cartilages meshes, independent of bone and cartilage mesh topology. Bone mesh-to-mesh distances, subchondral bone boundaries, and cartilage thicknesses from meshes of any type of mesh are obtained using a series of steps involving registration, resampling, and radial basis function fitting after which the comparisons are performed. Subchondral bone boundaries and cartilage thicknesses are calculated and visualized in a common frame of reference for comparison. The established method is applied to models developed by five modeling teams. Our approach to obtain bone mesh-to-mesh distances decreased the divergence seen in selecting a reference mesh (i.e., comparing mesh A-to-B versus mesh B-to-A). In general, the bone morphology was similar across teams. The cartilage thicknesses for all models were calculated and the mean absolute cartilage thickness difference was presented, the articulating areas had the best agreement across teams. The teams showed disagreement on the subchondral bone boundaries. The method presented in this paper allows for objective comparisons of bone and cartilage geometry that is agnostic to mesh type and nodal density.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Variations in Model Development.

The use of computational modeling to investigate knee joint biomechanics has increased exponentially over the last few decades. Developing computational models is a creative process where decisions have to be made, subject to the modelers' knowledge and previous experiences, resulting in the "art" of modeling. The long-term goal of the KneeHub project is to understand the influence of subjective decisions on the final outcomes and the reproducibility of computational knee joint models. In this paper, we report on the model development phase of this project, investigating model development decisions and deviations from initial modeling plans. Five teams developed computational knee joint models from the same dataset, and we compared each teams' initial uncalibrated models and their model development workflows. Variations in the software tools and modeling approaches were found, resulting in differences such as the representation of the anatomical knee joint structures in the model. The teams consistently defined the boundary conditions and used the same anatomical coordinate system convention. However, deviations in the anatomical landmarks used to define the coordinate systems were present, resulting in a large spread in the kinematic outputs of the uncalibrated models. The reported differences and similarities in model development and simulation presented here illustrate the importance of the "art" of modeling and how subjective decision-making can lead to variation in model outputs. All teams deviated from their initial modeling plans, indicating that model development is a flexible process and difficult to plan in advance, even for experienced teams.

Remote Patient Monitoring with Wearable Sensors Following Knee Arthroplasty.

(Background) Inertial Measurement Units (IMUs) provide a low-cost, portable solution to obtain functional measures similar to those captured with three-dimensional gait analysis, including spatiotemporal gait characteristics. The primary aim of this study was to determine the feasibility of a remote patient monitoring (RPM) workflow using ankle-worn IMUs measuring impact load, limb impact load asymmetry and knee range of motion in combination with patient-reported outcome measures. (Methods) A pilot cohort of 14 patients undergoing primary knee arthroplasty for osteoarthritis was prospectively enrolled. RPM in the community was performed weekly from 2 up to 6 weeks post-operatively using wearable IMUs. The following data were collected using IMUs: mobility (Bone Stimulus and cumulative impact load), impact load asymmetry and maximum knee flexion angle. In addition, scores from the Oxford Knee Score (OKS), EuroQol Five-dimension (EQ-5D) with EuroQol visual analogue scale (EQ-VAS) and 6 Minute Walk Test were collected. (Results) On average, the Bone Stimulus and cumulative impact load improved 52% (p = 0.002) and 371% (p = 0.035), compared to Post-Op Week 2. The impact load asymmetry value trended (p = 0.372) towards equal impact loading between the operative and non-operative limb. The mean maximum flexion angle achieved was 99.25° at Post-Operative Week 6, but this was not significantly different from pre-operative measurements (p = 0.1563). There were significant improvements in the mean EQ-5D (0.20; p = 0.047) and OKS (10.86; p < 0.001) scores both by 6 weeks after surgery, compared to pre-operative scores. (Conclusions) This pilot study demonstrates the feasibility of a reliable and low-maintenance workflow system to remotely monitor post-operative progress in knee arthroplasty patients. Preliminary data indicate IMU outputs relating to mobility, impact load asymmetry and range of motion can be obtained using commercially available IMU sensors. Further studies are required to directly correlate the IMU sensor outputs with patient outcomes to establish clinical significance.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Overall Strategy.

Recent explorations of knee biomechanics have benefited from computational modeling, specifically leveraging advancements in finite element analysis and rigid body dynamics of joint and tissue mechanics. A large number of models have emerged with different levels of fidelity in anatomical and mechanical representation. Adapted modeling and simulation processes vary widely, based on justifiable choices in relation to anticipated use of the model. However, there are situations where modelers' decisions seem to be subjective, arbitrary, and difficult to rationalize. Regardless of the basis, these decisions form the "art" of modeling, which impact the conclusions of simulation-based studies on knee function. These decisions may also hinder the reproducibility of models and simulations, impeding their broader use in areas such as clinical decision making and personalized medicine. This document summarizes an ongoing project that aims to capture the modeling and simulation workflow in its entirety-operation procedures, deviations, models, by-products of modeling, simulation results, and comparative evaluations of case studies and applications. The ultimate goal of the project is to delineate the art of a cohort of knee modeling teams through a publicly accessible, transparent approach and begin to unravel the complex array of factors that may lead to a lack of reproducibility. This manuscript outlines our approach along with progress made so far. Potential implications on reproducibility, on science, engineering, and training of modeling and simulation, on modeling standards, and on regulatory affairs are also noted.

Application of a Simple Surveillance Method for Detecting the Prevalence and Impact of Overuse Injuries in Professional Men's Basketball.

The aim of this study was to use the Oslo Sports Trauma Research Center (OSTRC) Overuse Injury Questionnaire to record overuse injuries over a single season for a men's professional basketball team to (a) assess the prevalence and severity of overuse injuries and (b) determine the efficacy of this method in identifying overuse injuries in comparison with the team physiotherapist's detection of these injuries. Thirteen athletes from a men's professional basketball team participated in this study. The self-reported, OSTRC injury questionnaire was used to record overuse conditions of the ankle, knee, and lower back over an entire 24-week season. Standard time-loss injury registration methods were also used to record overuse conditions by the physiotherapist. Overuse injury rates per 1,000 hours of athlete exposure and average weekly prevalence of overuse injuries were calculated using the results of the questionnaire. A total of 183 overuse conditions were identified by the questionnaire, whereas only 28 overuse conditions were identified by the physiotherapist. The team's average weekly prevalence of all overuse conditions was 63% (95% confidence interval [CI]: 60-66), with the highest prevalence of injury affecting the lower back (25.9% [95% CI: 19.7-32.1]). The overuse injury rate per 1,000 hours of athlete exposure was 6.4. The OSTRC overuse injury questionnaire captures many more overuse injuries in basketball than standard time-loss methods. The prevalence of lower back injuries is higher than that previously reported in basketball. This additional method of overuse injury surveillance may more accurately quantify the overuse injury problem in basketball and aid earlier intervention and management of these conditions.

Muscle synergies may improve optimization prediction of knee contact forces during walking.

The ability to predict patient-specific joint contact and muscle forces accurately could improve the treatment of walking-related disorders. Muscle synergy analysis, which decomposes a large number of muscle electromyographic (EMG) signals into a small number of synergy control signals, could reduce the dimensionality and thus redundancy of the muscle and contact force prediction process. This study investigated whether use of subject-specific synergy controls can improve optimization prediction of knee contact forces during walking. To generate the predictions, we performed mixed dynamic muscle force optimizations (i.e., inverse skeletal dynamics with forward muscle activation and contraction dynamics) using data collected from a subject implanted with a force-measuring knee replacement. Twelve optimization problems (three cases with four subcases each) that minimized the sum of squares of muscle excitations were formulated to investigate how synergy controls affect knee contact force predictions. The three cases were: (1) Calibrate+Match where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously matched, (2) Precalibrate+Predict where experimental knee contact forces were predicted using precalibrated muscle model parameters values from the first case, and (3) Calibrate+Predict where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously predicted, all while matching inverse dynamic loads at the hip, knee, and ankle. The four subcases used either 44 independent controls or five synergy controls with and without EMG shape tracking. For the Calibrate+Match case, all four subcases closely reproduced the measured medial and lateral knee contact forces (R2 ≥ 0.94, root-mean-square (RMS) error < 66 N), indicating sufficient model fidelity for contact force prediction. For the Precalibrate+Predict and Calibrate+Predict cases, synergy controls yielded better contact force predictions (0.61 < R2 < 0.90, 83 N < RMS error < 161 N) than did independent controls (-0.15 < R2 < 0.79, 124 N < RMS error < 343 N) for corresponding subcases. For independent controls, contact force predictions improved when precalibrated model parameter values or EMG shape tracking was used. For synergy controls, contact force predictions were relatively insensitive to how model parameter values were calibrated, while EMG shape tracking made lateral (but not medial) contact force predictions worse. For the subject and optimization cost function analyzed in this study, use of subject-specific synergy controls improved the accuracy of knee contact force predictions, especially for lateral contact force when EMG shape tracking was omitted, and reduced prediction sensitivity to uncertainties in muscle model parameter values.

An articulated shape model to predict paediatric lower limb bone geometry using sparse landmarks.

Creating musculoskeletal models in a paediatric population currently involves either creating an image-based model from medical imaging data or a generic model using linear scaling. Image-based models provide a high level of accuracy but are time-consuming and costly to implement, on the other hand, linear scaling of an adult template musculoskeletal model is faster and common practice, but the output errors are significantly higher. An articulated shape model incorporates pose and shape to predict geometry for use in musculoskeletal models based on existing information from a population to provide both a fast and accurate method. From a population of 333 children aged 4-18 years old, we have developed an articulated shape model of paediatric lower limb bones to predict bone geometry from eight bone landmarks commonly used for motion capture. Bone surface root mean squared errors were found to be 2.63 ± 0.90 mm, 1.97 ± 0.61 mm, and 1.72 ± 0.51 mm for the pelvis, femur, and tibia/fibula, respectively. Linear scaling produced bone surface errors of 4.79 ± 1.39 mm, 4.38 ± 0.72 mm, and 4.39 ± 0.86 mm for the pelvis, femur, and tibia/fibula, respectively. Clinical bone measurement errors were low across all bones predicted using the articulated shape model, which outperformed linear scaling for all measurements. However, the model failed to accurately capture torsional measures (femoral anteversion and tibial torsion). Overall, the articulated shape model was shown to be a fast and accurate method to predict lower limb bone geometry in a paediatric population, superior to linear scaling.

A System for Reproducible 3D Ultrasound Measurements of Skeletal Muscles.

In 3D freehand ultrasound imaging, operator dependent variations in applied forces and movements can lead to errors in the reconstructed images. In this paper, we introduce an automated 3D ultrasound system, which enables acquisitions with controlled movement trajectories by using motors, which electrically move the probe. Due to integrated encoders there is no need of position sensors. An included force control mechanism ensures a constant contact force to the skin. We conducted 8 trials with the automated 3D ultrasound system on 2 different phantoms with 3 force settings and 10 trials on a human tibialis anterior muscle with 2 force settings. For comparison, we also conducted 8 freehand 3D ultrasound scans from 2 operators (4 force settings) on one phantom and 10 with one operator on the tibialis anterior muscle. Both freehand and automated trials showed small errors in volume and length computations of the reconstructions, however the freehand trials showed larger standard deviations. We also computed the thickness of the phantom and the tibialis anterior muscle. We found significant differences in force settings for the operators and higher coefficients of variation for the freehand trials. Overall, the automated 3D ultrasound system shows a high accuracy in reconstruction. Due to the smaller coefficients of variation, the automated 3D ultrasound system enables more reproducible ultrasound examinations than the freehand scanning. Therefore, the automated 3D ultrasound system is a reliable tool for 3D investigations of skeletal muscle.

Update on grand challenge competition to predict in vivo knee loads.

Validation is critical if clinicians are to use musculoskeletal models to optimize treatment of individual patients with a variety of musculoskeletal disorders. This paper provides an update on the annual Grand Challenge Competition to Predict in Vivo Knee Loads, a unique opportunity for direct validation of knee contact forces and indirect validation of knee muscle forces predicted by musculoskeletal models. Three competitions (2010, 2011, and 2012) have been held at the annual American Society of Mechanical Engineers Summer Bioengineering Conference, and two more competitions are planned for the 2013 and 2014 conferences. Each year of the competition, a comprehensive data set collected from a single subject implanted with a force-measuring knee replacement is released. Competitors predict medial and lateral knee contact forces for two gait trials without knowledge of the experimental knee contact force measurements. Predictions are evaluated by calculating root-mean-square (RMS) errors and R(2) values relative to the experimentally measured medial and lateral contact forces. For the first three years of the competition, competitors used a variety of methods to predict knee contact and muscle forces, including static and dynamic optimization, EMG-driven models, and parametric numerical models. Overall, errors in predicted contact forces were comparable across years, with average RMS errors for the four competition winners ranging from 229 N to 312 N for medial contact force and from 238 N to 326 N for lateral contact force. Competitors generally predicted variations in medial contact force (highest R(2 )= 0.91) better than variations in lateral contact force (highest R(2 )= 0.70). Thus, significant room for improvement exists in the remaining two competitions. The entire musculoskeletal modeling community is encouraged to use the competition data and models for their own model validation efforts.

A viscoelastic constitutive model can accurately represent entire creep indentation tests of human patella cartilage.

Cartilage material properties provide important insights into joint health, and cartilage material models are used in whole-joint finite element models. Although the biphasic model representing experimental creep indentation tests is commonly used to characterize cartilage, cartilage short-term response to loading is generally not characterized using the biphasic model. The purpose of this study was to determine the short-term and equilibrium material properties of human patella cartilage using a viscoelastic model representation of creep indentation tests. We performed 24 experimental creep indentation tests from 14 human patellar specimens ranging in age from 20 to 90 years (median age 61 years). We used a finite element model to reproduce the experimental tests and determined cartilage material properties from viscoelastic and biphasic representations of cartilage. The viscoelastic model consistently provided excellent representation of the short-term and equilibrium creep displacements. We determined initial elastic modulus, equilibrium elastic modulus, and equilibrium Poisson's ratio using the viscoelastic model. The viscoelastic model can represent the short-term and equilibrium response of cartilage and may easily be implemented in whole-joint finite element models.

Deep-Learning Markerless Tracking of Infant General Movements using Standard Video Recordings.

Monitoring spontaneous General Movements (GM) of infants 6-20 weeks post-term age is a reliable tool to assess the quality of neurodevelopment in early infancy. Abnormal or absent GMs are reliable prognostic indicators of whether an infant is at risk of developing neurological impairments and disorders such as cerebral palsy (CP). Therapeutic interventions are most effective at improving neuromuscular outcomes if administered in early infancy. Current clinical protocols require trained assessors to rate videos of infant movements, a time-intensive task. This work proposes a simple, inexpensive, and broadly applicable markerless pose-estimation approach for automatic infant movement tracking using conventional video recordings from handheld devices (e.g., tablets and mobile phones). We leverage the enhanced capabilities of deep-learning technology in image processing to identify 12 anatomical locations (3 per limb) in each video frame, tracking a baby's natural movement throughout the recordings. We validate the capability of resnet152 and a mobile-net-v2-1 to identify body-parts in unseen frames from a full-term male infant, using a novel automatic unsupervised approach that fuses likelihood outputs of a Kalman filter and the deep-nets. Both deep-net models were found to perform very well in the identification of anatomical locations in the unseen data with high average Percentage of Correct Keypoints (aPCK) performances of >99.65% across all locations.Clinical relevance-Results of this research confirm the feasibility of a low-cost and publicly accessible technology to automatically track infants' GMs and diagnose those at higher risk of developing neurological conditions early, when clinical interventions are most effective.

Predicting ankle and knee sagittal kinematics and kinetics using an ankle-mounted inertial sensor.

The purpose of this study was to develop a machine learning model to reconstruct time series kinematic and kinetic profiles of the ankle and knee joint across six different tasks using an ankle-mounted IMU. Four male collegiate basketball players performed repeated tasks, including walking, jogging, running, sidestep cutting, max-height jumping, and stop-jumping, resulting in a total of 102 movements. Ankle and knee flexion-extension angles and moments were estimated using motion capture and inverse dynamics and considered 'actual data' for the purpose of model fitting. Synchronous acceleration and angular velocity data were collected from right ankle-mounted IMUs. A time-series feature extraction model was used to determine a set of features used as input to a random forest regression model to predict the ankle and knee kinematics and kinetics. Five-fold cross-validation was performed to verify the model accuracy, and statistical parametric mapping was used to determine the difference between the predicted and experimental time series. The random forest regression model predicted the time-series profiles of the ankle and knee flexion-extension angles and moments with high accuracy (Kinematics: R2 ranged from 0.782 to 0.962, RMSE ranged from 2.19° to 11.58°; Kinetics: R2 ranged from 0.711 to 0.966, RMSE ranged from 0.10 Nm/kg to 0.41 Nm/kg). There were differences between predicted and actual time series for the knee flexion-extension moment during stop-jumping and walking. An appropriately trained feature-based regression model can predict time series knee and ankle joint angles and moments across a wide range of tasks using a single ankle-mounted IMU.

Assessment of reporting practices and reproducibility potential of a cohort of published studies in computational knee biomechanics.

Reproducible research serves as a pillar of the scientific method and is a foundation for scientific advancement. However, estimates for irreproducibility of preclinical science range from 75% to 90%. The importance of reproducible science has not been assessed in the context of mechanics-based modeling of human joints such as the knee, despite this being an area that has seen dramatic growth. Framed in the context of five experienced teams currently documenting knee modeling procedures, the aim of this study was to evaluate reporting and the perceived potential for reproducibility across studies the teams viewed as important contributions to the literature. A cohort of studies was selected by polling, which resulted in an assessment of nine studies as opposed to a broader analysis across the literature. Using a published checklist for reporting of modeling features, the cohort was evaluated for both "reporting" and their potential to be "reproduced," which was delineated into six major modeling categories and three subcategories. Logistic regression analysis revealed that for individual modeling categories, the proportion of "reported" occurrences ranged from 0.31, 95% confidence interval (CI) [0.23, 0.41] to 0.77, 95% CI: [0.68, 0.86]. The proportion of whether a category was perceived as "reproducible" ranged from 0.22, 95% CI: [0.15, 0.31] to 0.44, 95% CI: [0.35, 0.55]. The relatively low ratios highlight an opportunity to improve reporting and reproducibility of knee modeling studies. Ongoing efforts, including our findings, contribute to a dialogue that facilitates adoption of practices that provide both credibility and translation possibilities.

Reproducibility in modeling and simulation of the knee: Academic, industry, and regulatory perspectives.

Stakeholders in the modeling and simulation (M&S) community organized a workshop at the 2019 Annual Meeting of the Orthopaedic Research Society (ORS) entitled "Reproducibility in Modeling and Simulation of the Knee: Academic, Industry, and Regulatory Perspectives." The goal was to discuss efforts among these stakeholders to address irreproducibility in M&S focusing on the knee joint. An academic representative from a leading orthopedic hospital in the United States described a multi-institutional, open effort funded by the National Institutes of Health to assess model reproducibility in computational knee biomechanics. A regulatory representative from the United States Food and Drug Administration indicated the necessity of standards for reproducibility to increase utility of M&S in the regulatory setting. An industry representative from a major orthopedic implant company emphasized improving reproducibility by addressing indeterminacy in personalized modeling through sensitivity analyses, thereby enhancing preclinical evaluation of joint replacement technology. Thought leaders in the M&S community stressed the importance of data sharing to minimize duplication of efforts. A survey comprised 103 attendees revealed strong support for the workshop and for increasing emphasis on computational modeling at future ORS meetings. Nearly all survey respondents (97%) considered reproducibility to be an important issue. Almost half of respondents (45%) tried and failed to reproduce the work of others. Two-thirds of respondents (67%) declared that individual laboratories are most responsible for ensuring reproducible research whereas 44% thought that journals are most responsible. Thought leaders and survey respondents emphasized that computational models must be reproducible and credible to advance knee M&S.

Comparison of MRI and ¹8F-NaF PET/CT in patients with patellofemoral pain.

PURPOSE: To determine whether bone metabolic activity corresponds to bone and cartilage damage in patients with patellofemoral pain. MATERIALS AND METHODS: We acquired magnetic resonance imaging (MRI) and (18) F-NaF positron emission tomography (PET) / computed tomography (CT) scans of the knees of 22 subjects. We compared locations of increased tracer uptake on the (18) F-NaF PET images to bone marrow edema and cartilage damage visualized on MRI. RESULTS: We found that increased bone activity on (18) F-NaF PET does not always correspond to structural damage in the bone or cartilage as seen on MRI. CONCLUSION: Our results suggest that (18) F-NaF PET/CT may provide additional information in patellofemoral pain patients compared to MRI.

Predicting the hip joint centre in children: New regression equations, linear scaling, and statistical shape modelling.

Determination of the hip joint centre (HJC) is important to accurately estimate hip joint motion, moments and muscle forces. The most accurate method for HJC estimation without medical imaging is an area of interest in the biomechanics community, especially in a paediatric population, which has not been widely evaluated. HJC locations were calculated by sphere-fitting to the acetabulum of three-dimensional pelvises segmented from 333 CT scans of children aged 4 to 18 years old. Three methods for determining the HJC were compared: regression equations, linear scaling, and shape model prediction. The new regression equations developed in this study produced Euclidean distance errors of 6.23 mm ± 2.90 mm. Linear scaling of paediatric bone produced errors of 3.90 mm ± 2.52 mm and adult bone scaling of 5.45 mm ± 3.26 mm. Prediction of the HJC using a paediatric statistical shape model produced mean Euclidian distance errors of 2.95 mm ± 1.65 mm. Overall, shape model prediction of the HJC produced the lowest errors, with linear scaling of a mean paediatric pelvis providing better estimates than regression equations.

Morphological variation in paediatric lower limb bones.

Available methods for generating paediatric musculoskeletal geometry are to scale generic adult geometry, which is widely accessible but can be inaccurate, or to obtain geometry from medical imaging, which is accurate but time-consuming and costly. A population-based shape model is required to generate accurate and accessible musculoskeletal geometry in a paediatric population. The pelvis, femur, and tibia/fibula were segmented from 333 CT scans of children aged 4-18 years. Bone morphology variation was captured using principal component analysis (PCA). Subsequently, a shape model was developed to predict bone geometry from demographic and linear bone measurements and validated using a leave one out analysis. The shape model was compared to linear scaling of adult and paediatric bone geometry. The PCA captured growth-related changes in bone geometry. The shape model predicted bone geometry with root mean squared error (RMSE) of 2.91 ± 0.99 mm in the pelvis, 2.01 ± 0.62 mm in the femur, and 1.85 ± 0.54 mm in the tibia/fibula. Linear scaling of an adult mesh produced RMSE of 4.79 ± 1.39 mm in the pelvis, 4.38 ± 0.72 mm in the femur, and 4.39 ± 0.86 mm in the tibia/fibula. We have developed a method for capturing and predicting lower limb bone shape variation in a paediatric population more accurately than linear scaling without using medical imaging.