Dual-energy X-ray absorptiometry derived knee shape may provide a useful imaging biomarker for predicting total knee replacement: Findings from a study of 37,843 people in UK Biobank.

Objective: We aimed to create an imaging biomarker for knee shape using knee dual-energy x-ray absorptiometry (DXA) scans and investigate its potential association with subsequent total knee replacement (TKR), independently of radiographic features of knee osteoarthritis and established risk factors. Methods: Using a 129-point statistical shape model, knee shape (expressed as a B-score) and minimum joint space width (mJSW) of the medial joint compartment (binarized as above or below the first quartile) were derived. Osteophytes were manually graded in a subset of images and an overall score was assigned. Cox proportional hazards models were used to examine the associations of B-score, mJSW and osteophyte score with TKR risk, adjusting for age, sex, height and weight. Results: The analysis included 37,843 individuals (mean age 63.7 years). In adjusted models, B-score was associated with TKR: each unit increase in B-score, reflecting one standard deviation from the mean healthy shape, corresponded to a hazard ratio (HR) of 2.25 (2.08, 2.43), while a lower mJSW had a HR of 2.28 (1.88, 2.77). Among the 6719 images scored for osteophytes, mJSW was replaced by osteophyte score in the most strongly predictive model for TKR. In ROC analyses, a model combining B-score, osteophyte score, and demographics outperformed a model including demographics alone (AUC â€‹= â€‹0.87 vs 0.73). Conclusions: Using statistical shape modelling, we derived a DXA-based imaging biomarker for knee shape that was associated with kOA progression. When combined with osteophytes and demographic data, this biomarker may help identify individuals at high risk of TKR, facilitating targeted interventions.

Prediction of the Proximal Humerus Morphology Based on a Statistical Shape Model with Two Parameters: Comparison to Contralateral Registration Method.

(1) Background: Complex proximal humerus fractures often result in complications following surgical treatment. A better understanding of the full 3D displacement would provide insight into the fracture morphology. Repositioning of fracture elements is often conducted by using the contralateral side as a reconstruction template. However, this requires healthy contralateral anatomy. The purpose of this study was to create a Statistical Shape Model (SSM) and compare its effectiveness to the contralateral registration method for the prediction of the humeral proximal segment; (2) Methods: An SSM was created from 137 healthy humeri. A prediction for the proximal segment of the left humeri from eight healthy patients was made by combining the SSM with parameters. The predicted proximal segment was compared to the left proximal segment of the patients. Their left humerus was also compared to the contralateral (right) humerus; (3) Results: Eight modes explained 95% of the variation. Most deviations of the SSM prediction and the contralateral registration method were below the clinically relevant 2 mm distance threshold.; (4) Conclusions: An SSM combined with parameters is a suitable method to predict the proximal humeral segment when the contralateral CT scan is unavailable or the contralateral humerus is unhealthy, provided that the fracture pattern allows measurements of these parameters.

Regional plantar forces and surface geometry variations of a chronic ankle instability population described by statistical shape modelling.

BACKGROUND: Understanding detailed foot morphology as well as regional plantar forces could provide insight into foot function and provide recommendation for footwear design for chronic ankle instability (CAI) people. RESEARCH QUESTION: This study presented 3-dimensional statistical shape models of feet from three different populations including CAI, copers and healthy individuals, with regional plantar forces also acquired. METHODS: Sixty-six males (22 participants per group) were included in this study to capture 3-dimensional foot shapes under a standing condition and regional plantar forces during a cutting maneuver. Principal component analysis was performed to generate a mean foot shape of each group as well as modes of variations. A generalized procrustes analysis was used to achieve rapid registration of mean shapes. Besides, regional plantar forces and contact duration among these three populations were compared. RESULTS: For 3-dimensional foot shapes, although no significant differences of the average distance between each mode and mean shape were found among three populations, there were subtle variations in mean shapes. The CAI population presented a more bulging of the lateral malleolus; copers were characterized by the flexion of the lesser toes, a more bulging of the medial foot in the sagittal plane; and healthy individuals showed a greater heel width and a more bulging of the heel in the sagittal plane. In terms of plantar forces, healthy individuals had significantly greater summated plantar forces and greater plantar forces in the lateral heel area during the early contact phase compared to copers and CAI participants. SIGNIFICANCE: Overall, this study suggested that repetitive ankle sprains may lead to the bulging of the lateral malleolus. Further, CAI and copers seem to stabilize the ankle joint by medially shifting the center of pressure compared to healthy individuals under the static and less challenging dynamic conditions.

Geometric variation of the human tibia-fibula: a public dataset of tibia-fibula surface meshes and statistical shape model.

Background: Variation in tibia geometry is a risk factor for tibial stress fractures. Geometric variability in bones is often quantified using statistical shape modelling. Statistical shape models (SSM) offer a method to assess three-dimensional variation of structures and identify the source of variation. Although SSM have been used widely to assess long bones, there is limited open-source datasets of this kind. Overall, the creation of SSM can be an expensive process, that requires advanced skills. A publicly available tibia shape model would be beneficial as it enables researchers to improve skills. Further, it could benefit health, sport and medicine with the potential to assess geometries suitable for medical equipment, and aid in clinical diagnosis. This study aimed to: (i) quantify tibial geometry using a SSM; and (ii) provide the SSM and associated code as an open-source dataset. Methods: Lower limb computed tomography (CT) scans from the right tibia-fibula of 30 cadavers (male n = 20, female n = 10) were obtained from the New Mexico Decedent Image Database. Tibias were segmented and reconstructed into both cortical and trabecular sections. Fibulas were segmented as a singular surface. The segmented bones were used to develop three SSM of the: (i) tibia; (ii) tibia-fibula; and (iii) cortical-trabecular. Principal component analysis was applied to obtain the three SSM, with the principal components that explained 95% of geometric variation retained. Results: Overall size was the main source of variation in all three models accounting for 90.31%, 84.24% and 85.06%. Other sources of geometric variation in the tibia surface models included overall and midshaft thickness; prominence and size of the condyle plateau, tibial tuberosity, and anterior crest; and axial torsion of the tibial shaft. Further variations in the tibia-fibula model included midshaft thickness of the fibula; fibula head position relative to the tibia; tibia and fibula anterior-posterior curvature; fibula posterior curvature; tibia plateau rotation; and interosseous width. The main sources of variation in the cortical-trabecular model other than general size included variation in the medulla cavity diameter; cortical thickness; anterior-posterior shaft curvature; and the volume of trabecular bone in the proximal and distal ends of the bone. Conclusion: Variations that could increase the risk of tibial stress injury were observed, these included general tibial thickness, midshaft thickness, tibial length and medulla cavity diameter (indicative of cortical thickness). Further research is needed to better understand the effect of these tibial-fibula shape characteristics on tibial stress and injury risk. This SSM, the associated code, and three use examples for the SSM have been provided in an open-source dataset. The developed tibial surface models and statistical shape model will be made available for use at: https://simtk.org/projects/ssm_tibia.

Anatomical variation in humeri: gender and side comparison using statistical shape modelling.

PURPOSE: The surgical management of proximal humeral fractures remains challenging. Anatomical reduction of the fracture has been reported as the keystone for a sufficient surgical fixation and successful outcome. However, mostly there is no example of its premorbid state. Literature suggests that the mirrored contralateral side can be used as a reconstruction template. But is this a correct technique to use? The purpose of this study is to define anatomical variation between humeri based on gender and side comparison. METHODS: Two different statistical shape models of the humerus were created and their modes of variation were described. One model contained 110 unpaired humeri. The other model consisted of 65 left and corresponding right humeri. RESULTS: The compactness of the statistical shape model containing 110 humeri showed that two principal components explain more than 95% of the variation and the generalization showed that a random humerus can be described with an accuracy of 0.39 mm. For only three parameters, statistically significant differences were observed between left and right. However, comparing the mean of the different metrics on the humeri of men and women, almost all were significant. CONCLUSION: Since there were only small differences between left and right humeri, using the mirrored contralateral side as a reconstruction template for fracture reduction can be defended. The variable anatomy between men and women could explain why locking plates not always fit to the bone.

The accuracy of statistical shape models in predicting bone shape: A systematic review.

BACKGROUND: This systematic review aims to ascertain how accurately 3D models can be predicted from two-dimensional (2D) imaging utilising statistical shape modelling. METHODS: A systematic search of published literature was conducted in September 2022. All papers which assessed the accuracy of 3D models predicted from 2D imaging utilising statistical shape models and which validated the models against the ground truth were eligible. RESULTS: 2127 papers were screened and a total of 34 studies were included for final data extraction. The best overall achievable accuracy was 0.45 mm (root mean square error) and 0.16 mm (average error). CONCLUSION: Statistical shape modelling can predict detailed 3D anatomical models from minimal 2D imaging. Future studies should report the intended application domain of the model, the level of accuracy required, the underlying demographics of subjects, and the method in which accuracy was calculated, with root mean square error recommended if appropriate.

Foot shape is related to load-induced shape deformations, but neither are good predictors of plantar soft tissue stiffness.

Modern human feet are considered unique among primates in their capacity to transmit propulsive forces and re-use elastic energy. Considered central to both these capabilities are their arched configuration and the plantar aponeurosis (PA). However, recent evidence has shown that their interactions are not as simple as proposed by the theoretical and mechanical models that established their significance. Using three-dimensional foot scans and statistical shape and deformation modelling, we show that the shape of the longitudinal and transverse arches varies widely among the healthy adult population, and that the former is subject to load-induced arch flattening, whereas the latter is not. However, longitudinal arch shape and flattening are only one of the various foot shape-deformation relationships. PA stiffness was also found to vary widely. Yet only a small amount of this variability (approx. 10-18%) was explained by variations in foot shape, deformation and their combination. These findings add to the mounting evidence showing that foot mechanics are complex and cannot be accurately represented by simple models. Especially the interactions between longitudinal arch and PA appear to be far less constrained than originally proposed, most likely due to the many degrees of freedom provided by the structural complexity of our feet.

Shape modelling of the oropharynx distinguishes associations with body morphology but not whiplash-associated disorder.

Characterization of the oropharynx, a subdivision of the pharynx between the soft palate and the epiglottis, is limited to simple measurements. Structural changes in the oropharynx in whiplash-associated disorder (WAD) cohorts have been quantified using two-dimensional (2D) and three-dimensional (3D) measures but the results are inconsistent. Statistical shape modelling (SSM) may be a more useful tool for systematically comparing morphometric features between cohorts. This technique has been used to quantify the variability in boney and soft tissue structures, but has not been used to examine a hollow cavity such as the oropharynx. The primary aim of this project was to examine the utility of SSM for comparing the oropharynx between WAD cohorts and control; and WAD severity cohorts. The secondary aim was to determine whether shape is associated with sex, height, weight and neck length. Magnetic resonance (MR) T1-weighted images were obtained from healthy control (n = 20), acute WAD (n = 14) and chronic WAD (n = 14) participants aged 18-39 years. Demographic, WAD severity (neck disability index) and body morphometry data were collected from each participant. Manual segmentation of the oropharynx was undertaken by blinded researchers between the top of the soft palate and tip of the epiglottis. Digital 3D oropharynx models were constructed from the segmented images and principal component (PC) analysis was performed with the PC weights normalized to z-scores for consistency. Statistical analyses were undertaken using multivariate linear models. In the first statistical model the independent variable was group (acute WAD, chronic WAD, control); and in the second model the independent variable was WAD severity (recovered/mild, moderate/severe). The covariates for both models included height, weight, average neck length and sex. Shape models were constructed to visualize the effect of perturbing these covariates for each relevant mode. The shape model revealed five modes which explained 90% of the variance: mode 1 explained 59% of the variance and primarily described differences in isometric size of the oropharynx, including elongation; mode 2 (13%) primarily described lateral (width) and AP (depth) dimensions; mode 3 (8%) described retroglossal AP dimension; mode 4 (6%) described lateral dimensions at the retropalatal-retroglossal junction and mode 5 (4%) described the lateral dimension at the inferior retroglossal region. There was no difference in shape (mode 1 p = 0.52; mode 2 p = 0.96; mode 3 p = 0.07; mode 4 p = 0.54; mode 5 p = 0.74) between control, acute WAD and chronic WAD groups. There were no statistical differences for any mode (mode 1 p = 0.12; mode 2 p = 0.29; mode 3 p = 0.56; mode 4 p = 0.99; mode 5 p = 0.96) between recovered/mild and moderate/severe WAD. Sex was not significant in any of the models but for mode 1 there was a significant association with height (p = 0.007), mode 2 neck length (p = 0.044) and in mode 3 weight (p = 0.027). Although SSM did not detect differences between WAD cohorts, it did detect associations with body morphology indicating that it may be a useful tool for examining differences in the oropharynx.

Statistical shape modelling of the thoracic spine for the development of pedicle screw insertion guides.

Spinal fixation and fusion are surgical procedures undertaken to restore stability in the spine and restrict painful or degenerative motion. Malpositioning of pedicle screws during these procedures can result in major neurological and vascular damage. Patient-specific surgical guides offer clear benefits, reducing malposition rates by up to 25%. However, they suffer from long lead times and the manufacturing process is dependent on third-party specialists. The development of a standard set of surgical guides may eliminate the issues with the manufacturing process. To evaluate the feasibility of this option, a statistical shape model (SSM) was created and used to analyse the morphological variations of the T4-T6 vertebrae in a population of 90 specimens from the Visible Korean Human dataset (50 females and 40 males). The first three principal components, representing 39.7% of the variance within the population, were analysed. The model showed high variability in the transverse process (~ 4 mm) and spinous process (~ 4 mm) and relatively low variation (< 1 mm) in the vertebral lamina. For a Korean population, a standardised set of surgical guides would likely need to align with the lamina where the variance in the population is lower. It is recommended that this standard set of surgical guides should accommodate pedicle screw diameters of 3.5-6 mm and transverse pedicle screw angles of 3.5°-12.4°.

Automated volumetric and statistical shape assessment of cam-type morphology of the femoral head-neck region from clinical 3D magnetic resonance images.

Background: Femoroacetabular impingement (FAI) cam morphology is routinely assessed using manual measurements of two-dimensional (2D) alpha angles which are prone to high rater variability and do not provide direct three-dimensional (3D) data on these osseous formations. We present CamMorph, a fully automated 3D pipeline for segmentation, statistical shape assessment and measurement of cam volume, surface area and height from clinical magnetic resonance (MR) images of the hip in FAI patients. Methods: The novel CamMorph pipeline involves two components: (I) accurate proximal femur segmentation generated by combining the 3D U-net to identify both global (region) and local (edge) features in clinical MR images and focused shape modelling to generate a 3D anatomical model for creating patient-specific proximal femur models; (II) patient-specific anatomical information from 3D focused shape modelling to simulate 'healthy' femoral bone models with cam-affected region constraints applied to the anterosuperior femoral head-neck region to quantify cam morphology in FAI patients. The CamMorph pipeline, which generates patient-specific data within 5 min, was used to analyse multi-site clinical MR images of the hip to measure and assess cam morphology in male (n=56) and female (n=41) FAI patients. Results: There was excellent agreement between manual and CamMorph segmentations of the proximal femur as demonstrated by the mean Dice similarity index (DSI; 0.964±0.006), 95% Hausdorff distance (HD; 2.123±0.876 mm) and average surface distance (ASD; 0.539±0.189 mm) values. Compared to female FAI patients, male patients had a significantly larger median cam volume (969.22 vs. 272.97 mm3, U=240.0, P<0.001), mean surface area [657.36 vs. 306.93 mm2, t(95)=8.79, P<0.001], median maximum-height (3.66 vs. 2.15 mm, U=407.0, P<0.001) and median average-height (1.70 vs. 0.86 mm, U=380.0, P<0.001). Conclusions: The fully automated 3D CamMorph pipeline developed in the present study successfully segmented and measured cam morphology from clinical MR images of the hip in male and female patients with differing FAI severity and pathoanatomical characteristics.

Generation of hemipelvis surface geometry based on statistical shape modelling and contralateral mirroring.

Personalised fracture plates manufactured using 3D printing offer an improved treatment option for unstable pelvic ring fractures that may not be adequately secured using off-the-shelf components. To design fracture plates that secure the bone fragments in their pre-fracture positions, the fractures must be reduced virtually using medical imaging-based reconstructions, a time-consuming process involving segmentation and repositioning of fragments until surface congruency is achieved. This study compared statistical shape models (SSMs) and contralateral mirroring as automated methods to reconstruct the hemipelvis using varying amounts of bone surface geometry. The training set for the geometries was obtained from pelvis CT scans of 33 females. The root-mean-squared error (RMSE) was quantified across the entire surface of the hemipelvis and within specific regions, and deviations of pelvic landmarks were computed from their positions in the intact hemipelvis. The reconstruction of the entire hemipelvis surfaced based on contralateral mirroring had an RMSE of 1.21 ± 0.29 mm, whereas for SSMs based on the entire hemipelvis surface, the RMSE was 1.11 ± 0.29 mm, a difference that was not significant (p = 0.32). Moreover, all hemipelvis reconstructions based on the full or partial bone geometries had RMSEs and landmark deviations from contralateral mirroring that were significantly lower (p < 0.05) or statistically equivalent to the SSMs. These results indicate that contralateral mirroring tends to be more accurate than SSMs for reconstructing unilateral pelvic fractures. SSMs may still be a viable method for hemipelvis fracture reconstruction in situations where contralateral geometries are not available, such as bilateral pelvic factures, or for highly asymmetric pelvic anatomies.

Personalised statistical modelling of soft tissue structures in the ankle.

BACKGROUND AND OBJECTIVE: Revealing the complexity behind subject-specific ankle joint mechanics requires simultaneous analysis of three-dimensional bony and soft-tissue structures. 3D musculoskeletal models have become pivotal in orthopedic treatment planning and biomechanical research. Since manual segmentation of these models is time-consuming and subject to manual errors, (semi-) automatic methods could improve the accuracy and enlarge the sample size of personalised 'in silico' biomechanical experiments and computer-assisted treatment planning. Therefore, our aim was to automatically predict ligament paths, cartilage topography and thickness in the ankle joint based on statistical shape modelling. METHODS: A personalised cartilage and ligamentous prediction algorithm was established using geometric morphometrics, based on an 'in-house' generated lower limb skeletal model (N = 542), tibiotalar cartilage (N = 60) and ankle ligament segmentations (N = 10). For cartilage, a population-averaged thickness map was determined by use of partial least-squares regression. Ligaments were wrapped around bony contours based on iterative shortest path calculation. Accuracy of ligament path and cartilage thickness prediction was quantified using leave-one-out experiments. The novel personalised thickness prediction was compared with a constant cartilage thickness of 1.50 mm by use of a paired sample T-test. RESULTS: Mean distance error of cartilage and ligament prediction was 0.12 mm (SD 0.04 mm) and 0.54 mm (SD 0.05 mm), respectively. No significant differences were found between the personalised thickness cartilage and segmented cartilage of the tibia (p = 0.73, CI [-1.60 .10-17, 1.13 .10-17]) and talus (p = 0.95, CI[ -1.35 .10-17, 1.28 .10-17]). For the constant thickness cartilage, a statistically significant difference was found in 89% and 92% of the tibial (p < 0.001, CI [0.51, 0.58]) and talar (p < 0.001, CI [0.33, 0.40]) cartilage area. CONCLUSIONS: In this study, we described a personalised prediction algorithm of cartilage and ligaments in the ankle joint. We were able to predict cartilage and main ankle ligaments with submillimeter accuracy. The proposed method has a high potential for generating large (virtual) sample sizes in biomechanical research and mitigates technological advances in computer-assisted orthopaedic surgery.

Correlation between head shape and volumetric changes following spring-assisted posterior vault expansion.

The aim of the study was to investigate whether different head shapes show different volumetric changes following spring-assisted posterior vault expansion (SA-PVE) and to investigate the influence of surgical and morphological parameters on SA-PVE. Preoperative three-dimensional skull models from patients who underwent SA-PVE were extracted from computed tomography scans. Patient head shape was described using statistical shape modelling (SSM) and principal component analysis (PCA). Preoperative and postoperative intracranial volume (ICV) and cranial index (CI) were calculated. Surgical and morphological parameters included skull bone thickness, number of springs, duration of spring insertion and type of osteotomy. In the analysis, 31 patients were included. SA-PVE resulted in a significant ICV increase (284.1 ± 171.6 cm3, p < 0.001) and a significant CI decrease (-2.9 ± 4.3%, p < 0.001). The first principal component was significantly correlated with change in ICV (Spearman ρ = 0.68, p < 0.001). Change in ICV was significantly correlated with skull bone thickness (ρ = -0.60, p < 0.001) and age at time of surgery (ρ = -0.60, p < 0.001). No correlations were found between the change in ICV and number of springs, duration of spring insertion and type of osteotomy. SA-PVE is effective for increasing the ICV and resolving raised intracranial pressure. Younger, brachycephalic patients benefit more from surgery in terms of ICV increase. Skull bone thickness seems to be a crucial factor and should be assessed to achieve optimal ICV increase. In contrast, insertion of more than two springs, duration of spring insertion or performing a fully cut through osteotomy do not seem to impact the ICV increase. When interpreting ICV increases, normal calvarial growth should be taken into account.

The reliability and validity of triceps surae muscle volume assessment using freehand three-dimensional ultrasound in typically developing infants.

This study assessed the intra-acquirer, intra- and inter-processor reliability, and validity of the in vivo assessment of the medial gastrocnemius (MG), lateral gastrocnemius (LG) and soleus (SOL) muscle volumes using freehand 3D ultrasound (3DUS) in typically developing infants. Reliability assessments of freehand 3DUS were undertaken in infants across three ages groups: three, six and twelve months of age, with validity testing completed against magnetic resonance imaging (MRI) in infants at 3 months of age. Freehand 3DUS scanning was carried out by a single acquirer, with two independent processors manually segmenting images to render volumes. MRI images were segmented independently by a separate processor, with the volumes compared to those obtained via freehand 3DUS. Reliability was assessed using intraclass correlation (ICC), coefficient of variance (CV) and minimal detectable change (MDC) across each assessment time point. Validity was assessed using the limits of agreement. ICCs for intra-acquirer reliability of the acquisition process for freehand 3DUS ranged from 0.91 to 0.99 across all muscles. ICCs for intra-processor and inter-processor reliability for the segmentation process of freehand 3DUS ranged from 0.80 to 0.98 across all muscles. Acceptable levels of agreement between muscle volume obtained by freehand 3DUS and MRI were found for all muscles; however, freehand 3DUS overestimated muscle volume of MG and LG and underestimate the SOL compared with MRI, with average absolute differences of MG = 0.3 ml, LG = 0.3 ml and Sol = 1.2 ml. Freehand 3DUS is a reliable method for measuring in vivo triceps surae muscle volume in typically developing infants. We conclude that freehand 3DUS is a useful tool to assess changes in muscle volume in response to growth and interventions in infants.

Statistical Shape Modelling the In Vivo Location of Acetabular Wear in Retrieved Hip Implants.

Edge-wear in acetabular cups is known to be correlated with greater volumes of material loss; the location of this wear pattern in vivo is less understood. Statistical shape modelling (SSM) may provide further insight into this. This study aimed to identify the most common locations of wear in vivo, by combining CT imaging, retrieval analysis and SMM. Shape variance was described in 20 retrieved metal-on-metal acetabular surfaces. These were revised after a mean of 90 months, from 13 female and seven male patients. They were positioned with a mean inclination and anteversion of 53° and 30°, respectively. Their orientation, in vivo, was established using their stabilising fins, visible in pre-revision CT imaging. The impact of wear volume, positioning, time, gender and size on the in vivo location of wear was investigated. These surfaces had a mean wear volume of 49.63 mm3. The mean acetabular surface displayed superior edge-wear centred 7° within the posterosuperior quadrant, while more of the volumetric wear occurred in the anterosuperior quadrant. Components with higher inclination had greater superior edge-wear scars, while a relationship was observed between greater anteversion angles and more posterosuperior edge-wear. This SSM method can further our understanding of hip implant function, informing future design and may help to refine the safe zone for implant positioning.

Statistical shape modelling to analyse the talus in paediatric clubfoot.

One fifth of idiopathic clubfoot deformities cannot be fully corrected by Serial Ponseti casting and deformity recurs in 20%-30% of cases. To avoid x-ray exposure, the joints with largely unossified bones are diagnosed with magnetic resonance images (MRI). Typically, geometric measurements are made in the MRI planes; however, this method is inaccurate compared to measurements on three-dimensional (3D) models of the joint. More accurate measurements using the 3D bone shapes may be better at identifying differences between groups; and therefore, improve diagnosis. The entire set of shape features from MRI can be analysed simultaneously through statistical shape modelling (SSM) which assesses bone morphology of clubfoot in a more sensitive way. A method for SSM of the talus is developed in this study and the shape of the normal talus is compared with the one in clubfeet with residual deformity through both geometric measurements and SSM. Significant differences between two groups were found by both methods; and therefore, might contribute to improve diagnosis of clubfoot.

Statistical shape modelling for the analysis of head shape variations.

The aim of this study is, firstly, to create a population-based 3D head shape model for the 0 to 2-year-old subjects to describe head shape variability within a normal population and, secondly, to test a combined normal and sagittal craniosynostosis (SAG) population model, able to provide surgical outcome assessment. 3D head shapes of patients affected by non-cranial related pathologies and of SAG patients (pre- and post-op) were extracted either from head CTs or 3D stereophotography scans, and processed. Statistical shape modelling (SSM) was used to describe shape variability using two models - a normal population model (MODEL1) and a combined normal and SAG population model (MODEL2). Head shape variability was described via principal components analysis (PCA) which calculates shape modes describing specific shape features. MODEL1 (n = 65) mode 1 showed statistical correlation (p < 0.001) with width (125.8 ± 13.6 mm), length (151.3 ± 17.4 mm) and height (112.5 ± 11.1 mm) whilst mode 2 showed correlation with cranial index (83.5 mm ± 6.3 mm, p < 0.001). The remaining 9 modes showed more subtle head shape variability. MODEL2 (n = 159) revealed that post-operative head shape still did not achieve full shape normalization with either spring cranioplasty or total calvarial remodelling. This study proves that SSM has the potential to describe detailed anatomical variations in a paediatric population.

Reliability and quality of statistical shape and deformation models constructed from optical foot scans.

The unique shape of modern human feet, and how they change shape when loaded are thought to be integral to effective upright gait. This unique shape, and the natural variations therein, have previously been analysed using a range of methods; from visual assessments, anthropometric measurements, and footprints, to x-ray, ultrasound and magnetic resonance images. However, these methods are often limited by their use of linear two-dimensional measures. Only recently have advances in three-dimensional (3D) scanning technology and statistical shape analysis been applied to studying 3D foot shape variations. Given their novelty, information regarding the reliability and repeatability of 3D foot scanning and shape modelling is lacking. To investigate whether repeated foot scans captured by two examiners give the same 3D shape and produce consistent statistical shape models, 17 healthy adults' left feet were scanned while bearing half and full bodyweight, as well as minimal weight. Surface to surface distances between corresponding foot meshes and differences between shape model quality criteria were both found to be small and insignificant. The only exception being the specificity criterion for minimally loaded foot scans. Furthermore, Euclidean vectors were used to model the magnitude and direction of deformation that feet undergo as a consequence of increased loading. The deformation models showed that loading a minimally loaded foot results in greater, but less consistent, shape changes than when increasing the load on an already loaded foot. These results show that the utilized methods offer a valuable, reliable and repeatable approach to analysing foot shape and deformation.

Development and application of the average pelvic shape in virtual pelvic fracture reconstruction.

BACKGROUND: With unilateral pelvic fractures, the contralateral hemipelvis can be used as a template in virtual reconstruction; however, this cannot be applied for bilateral fractures. Therefore, statistical shape modelling was used to build average pelvic shapes that can serve as templates when reconstructing bilaterally fractured pelvises. METHODS: Four average shape models were created for male and female, left and right hemipelves from 20 male and 20 female subjects. They were used as templates to reconstruct eight unilaterally fractured pelvises. RESULTS: The average root-mean-square of deviations between the reconstructed and intact hemipelves was 1.46 ± 0.32 mm, which is less than the 2 mm threshold for causing hip joint complications. CONCLUSION: This indicates that the reconstructions are reliable and the average shape models can be used to reconstruct bilaterally fractured pelvises. The proposed technique can potentially provide quick and accurate treatment plans for pelvic fracture patients, which is necessary for recovery.

Automated creation and tuning of personalised muscle paths for OpenSim musculoskeletal models of the knee joint.

Computational modelling is an invaluable tool for investigating features of human locomotion and motor control which cannot be measured except through invasive techniques. Recent research has focussed on creating personalised musculoskeletal models using population-based morphing or directly from medical imaging. Although progress has been made, robust definition of two critical model parameters remains challenging: (1) complete tibiofemoral (TF) and patellofemoral (PF) joint motions, and (2) muscle tendon unit (MTU) pathways and kinematics (i.e. lengths and moment arms). The aim of this study was to develop an automated framework, using population-based morphing approaches to create personalised musculoskeletal models, consisting of personalised bone geometries, TF and PF joint mechanisms, and MTU pathways and kinematics. Informed from medical imaging, personalised rigid body TF and PF joint mechanisms were created. Using atlas- and optimisation-based methods, personalised MTU pathways and kinematics were created with the aim of preventing MTU penetration into bones and achieving smooth MTU kinematics that follow patterns from existing literature. This framework was integrated into the Musculoskeletal Atlas Project Client software package to create and optimise models for 6 participants with incrementally increasing levels of personalisation with the aim of improving MTU kinematics and pathways. Three comparisons were made: (1) non-optimised (Model 1) and optimised models (Model 3) with generic joint mechanisms; (2) non-optimised (Model 2) and optimised models (Model 4) with personalised joint mechanisms; and (3) both optimised models (Model 3 and 4). Following optimisation, improvements were consistently shown in pattern similarity to cadaveric data in comparison (1) and (2). For comparison (3), a number of comparisons showed no significant difference between the two compared models. Importantly, optimisation did not produce statistically significantly worse results in any case.