Morphological variation in paediatric lower limb bones

OpenHands: An Open-Source Statistical Shape Model of the Finger Bones

This paper presents statistical shape models of the four fingers of the hand, with an emphasis on anatomic analysis of the proximal and distal interphalangeal joints. A multi-body statistical shape modelling pipeline was implemented on an exemplar training dataset of computed tomography (CT) scans of 10 right hands (5F:5M, 27-37 years, free from disease or injury) imaged at 0.3 mm resolution, segmented, meshed and aligned. Model generated included pose neutralisation to remove joint angle variation during imaging. Repositioning was successful; no joint flexion variation was observed in the resulting model. The first principal component (PC) of morphological variation represented phalanx size in all fingers. Subsequent PCs showed variation in position along the palmar-dorsal axis, and bone breadth: length ratio. Finally, the models were interrogated to provide gross measures of bone lengths and joint spaces. These models have been published for open use to support wider community efforts in hand biomechanical analysis, providing bony anatomy descriptions whilst preserving the security of the underlying imaging data and privacy of the participants. The model describes a small, homogeneous population, and assumptions cannot be made about how it represents individuals outside the training dataset. However, it supplements anthropometric datasets with additional shape information, and may be useful for investigating factors such as joint morphology and design of hand-interfacing devices and products. The model has been shared as an open-source repository ( https://github.com/abel-research/OpenHands ), and we encourage the community to use and contribute to it.

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.

3D gait analysis in children using wearable sensors: feasibility of predicting joint kinematics and kinetics with personalized machine learning models and inertial measurement units

Introduction: Children's walking patterns evolve with age, exhibiting less repetitiveness at a young age and more variability than adults. Three-dimensional gait analysis (3DGA) is crucial for understanding and treating lower limb movement disorders in children, traditionally performed using Optical Motion Capture (OMC). Inertial Measurement Units (IMUs) offer a cost-effective alternative to OMC, although challenges like drift errors persist. Machine learning (ML) models can mitigate these issues in adults, prompting an investigation into their applicability to a heterogeneous pediatric population. This study aimed at 1) quantifying personalized and generalized ML models' performance for predicting gait time series in typically developed (TD) children using IMUs data, 2) Comparing random forest (RF) and convolutional neural networks (CNN) models' performance, 3) Finding the optimal number of IMUs required for accurate predictions. Methodology: Seventeen TD children, aged 6 to 15, participated in data collection involving OMC, force plates, and IMU sensors. Joint kinematics and kinetics (targets) were computed from OMC and force plates' data using OpenSim. Tsfresh, a Python package, extracted features from raw IMU data. Each target's ten most important features were input in the development of personalized and generalized RF and CNN models. This procedure was initially conducted with 7 IMUs placed on all lower limb segments and then performed using only two IMUs on the feet. Results: Findings suggested that the RF and CNN models demonstrated comparable performance. RF predicted joint kinematics with a 9.5% and 19.9% NRMSE for personalized and generalized models, respectively, and joint kinetics with an NRMSE of 10.7% for personalized and 15.2% for generalized models in TD children. Personalized models provided accurate estimations from IMU data in children, while generalized models lacked accuracy due to the limited dataset. Furthermore, reducing the number of IMUs from 7 to 2 did not affect the results, and the performance remained consistent. Discussion: This study proposed a promising personalized approach for gait time series prediction in children, involving an RF model and two IMUs on the feet.

On the estimation of hip joint centre location with incomplete bone ossification for foetus-specific neuromusculoskeletal modeling

Childbirth is a complex physiological process in which a foetal neuromusculoskeletal model is of great importance to develop realistic delivery simulations and associated complication analyses. However, the estimation of hip joint centre (HJC) in foetuses remains a challenging issue. Thus, this paper aims to propose and evaluate a new approach to locate the HJC in foetuses. Hip CT-scans from 25 children (F = 11, age = 5.5 ± 2.6 years, height = 117 ± 21 cm, mass = 26 kg ± 9.5 kg) were used to propose and evaluate the novel acetabulum sphere fitting process to locate the HJC. This new approach using the acetabulum surface was applied to a population of 57 post-mortem foetal CT scans to locate the HJC as well as to determine associated regression equations using multiple linear regression. As results, the average distance between the HJC located using acetabulum sphere fitting and femoral head sphere fitting in children was 1.5 ± 0.7 mm. The average prediction error using our developed foetal HJC regression equations was 3.0 ± 1.5 mm, even though the equation for the x coordinate had a poor value of R2 (R2 for the x coordinate = 0.488). The present study suggests that the use of the acetabulum sphere fitting approach is a valid and accurate method to locate the HJC in children, and then can be extrapolated to get an estimation of the HJC in foetuses with incomplete bone ossification. Therefore, the present paper can be used as a guideline for foetus specific neuromusculoskeletal modelling.

A Narrative Review of Personalized Musculoskeletal Modeling Using the Physiome and Musculoskeletal Atlas Projects

In this narrative review, we explore developments in the field of computational musculoskeletal model personalization using the Physiome and Musculoskeletal Atlas Projects. Model geometry personalization; statistical shape modeling; and its impact on segmentation, classification, and model creation are explored. Examples include the trapeziometacarpal and tibiofemoral joints, Achilles tendon, gastrocnemius muscle, and pediatric lower limb bones. Finally, a more general approach to model personalization is discussed based on the idea of multiscale personalization called scaffolds.

Sex differences in linear bone measurements occur following puberty but do not influence femoral or tibial torsion

Torsional, angular, and linear measurements in a paediatric population are clinically important but not well defined and understood. Different methods of measurement and discrepancies between assessors leads to a lack of understanding of what should be defined as typical or atypical for the growing skeleton. From a large dataset of 333 paediatric CT scans, we extracted three-dimensional torsional, angular, and linear measurements from the pelvis, femur, and tibia/fibula. Sex differences in linear measurements were observed in bones of children aged 13+ (around puberty), but femoral and tibial torsion were similar between males and females. The rotational profile (femoral anteversion minus tibial torsion) tended to increase with growth. Epicondylar, condylar, and malleolar widths were smaller in females than males for the same bone length after the age of 13 years, which could explain why females may be more at risk for sport injuries during adolescence. This rich dataset can be used as an atlas for researchers and clinicians to understand typical development of critical rotational profiles and linear bone measurements in children.

Roadmap for an imaging and modelling paediatric study in rural NZ

Our study methodology is motivated from three disparate needs: one, imaging studies have existed in silo and study organs but not across organ systems; two, there are gaps in our understanding of paediatric structure and function; three, lack of representative data in New Zealand. Our research aims to address these issues in part, through the combination of magnetic resonance imaging, advanced image processing algorithms and computational modelling. Our study demonstrated the need to take an organ-system approach and scan multiple organs on the same child. We have pilot tested an imaging protocol to be minimally disruptive to the children and demonstrated state-of-the-art image processing and personalized computational models using the imaging data. Our imaging protocol spans brain, lungs, heart, muscle, bones, abdominal and vascular systems. Our initial set of results demonstrated child-specific measurements on one dataset. This work is novel and interesting as we have run multiple computational physiology workflows to generate personalized computational models. Our proposed work is the first step towards achieving the integration of imaging and modelling improving our understanding of the human body in paediatric health and disease.

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.