Inter-subject variability is inherently present in patient anatomy and is apparent in differences in shape, size and relative alignment of the bony structures. Understanding the variability in patient anatomy is useful for distinguishing between pathologies and to assist in surgical planning. With the aim of supporting the development of stratified orthopaedic interventions, this work introduces an Articulated Statistical Shape Model (ASSM) of the lower limb. The model captures inter-subject variability and allows reconstructing ‘virtual’ knee joints of the lower limb shape while considering pose.

A training dataset consisting of 173 lower limbs from CT scans of 110 subjects (77 male, 33 female) was used to construct the ASSM of the lower limb. Each bone of the lower limb was segmented using ScanIP (Simpleware Ltd., UK), reconstructed into 3D surface meshes, and a SSM of each bone was created. A series of sizing and positioning procedures were carried out to ensure all the lower limbs were in full extension, had the same femoral length and that the femora were aligned with a coincident centre. All articulated lower limbs were represented as: (femur scale factor) × (full extension articulated lower limb + relative transformation of tibia, fibula and patella to femur). Articulated lower limbs were in full extension were used to construct a statistical shape model, representing the variance of lower limb morphology. Relative transformations of the tibia, fibula and patella versus the femur were used to form a statistical pose model. Principal component analysis (PCA) was used to extract the modes of changes in the model.

The first 30 modes of the shape model covered 90% of the variance in shape and the first 10 modes of the pose model covered 90% of the pose variance. The first mode captures changes of the femoral CCD angle and the varus/valgus alignment of the knee. The second mode represents the changes in the ratio of femur to tibia length. The third mode reflects change of femoral shaft diameter and patella size. The first mode characterising pose captures the medial/lateral translation between femur and tibia. The second mode represents variation in knee flexion. The third mode reflects variation in tibio-femoral joint space.

An articulated statistical modelling approach was developed to characterize inter-subject variability in lower limb morphology for a set of training specimens. This model can generate large sets of lower limbs to systematically study the effect of anatomical variability on joint replacement performance. Moreover, if a series of images of the lower limb during a dynamic activity are used as training data, this method can be applied to analyse variance of lower limb motion across a population.