Statistical shape model-based prediction of tibiofemoral cartilage

Population-based in silico modeling of anatomical shape variation of the knee and its impact on joint loading in knee osteoarthritis

Anatomical knee joint features and osteoarthritis (OA) severity are associated, however confirming causals link to altered knee loading is challenging. This study leverages statistical shape models (SSM) to investigate the relationship between joint shape/alignment and knee loading during gait in knee OA (KOA) patients to understand their contribution to elevated medial knee loading in OA. Musculoskeletal (MSK) models were created for the mean as well as the first eight SSM principal modes of variation (-3,-2,-1, +1, +2, +3 standard deviations for each mode) and used as input to a MSK modeling framework. Using an identical KOA gait pattern (i.e., joint kinematics and ground reaction forces), we ran simulations for each MSK model and evaluated medial compartment loading magnitude and contact distribution at the instant of first and second peak of knee joint loading. An increase in external rotation, posterior tibia translation and a decrease in medial joint space and medial femoral condylar size predisposed the medial compartment knee joint to overloading during gait. This was coupled with an anterior and medial shift in contact location with increasing external rotated tibial position and increasing posterior tibial translation with respect to the femur. Next, results also highlighted a posterior shift of the medial compartment loading location with decreasing medial joint space. This study provides important population-based insights on how knee shape and alignment predispose individuals with KOA to elevated medial compartmental knee loading. This information can be crucial in assessing the risk for medial KOA development and progression.

Personalized statistical modeling of soft tissue structures in the knee

Background and Objective: As in vivo measurements of knee joint contact forces remain challenging, computational musculoskeletal modeling has been popularized as an encouraging solution for non-invasive estimation of joint mechanical loading. Computational musculoskeletal modeling typically relies on laborious manual segmentation as it requires reliable osseous and soft tissue geometry. To improve on feasibility and accuracy of patient-specific geometry predictions, a generic computational approach that can easily be scaled, morphed and fitted to patient-specific knee joint anatomy is presented.

Methods: A personalized prediction algorithm was established to derive soft tissue geometry of the knee, originating solely from skeletal anatomy. Based on a MRI dataset (n = 53), manual identification of soft-tissue anatomy and landmarks served as input for our model by use of geometric morphometrics. Topographic distance maps were generated for cartilage thickness predictions. Meniscal modeling relied on wrapping a triangular geometry with varying height and width from the anterior to the posterior root. Elastic mesh wrapping was applied for ligamentous and patellar tendon path modeling. Leave-one-out validation experiments were conducted for accuracy assessment.

Results: The Root Mean Square Error (RMSE) for the cartilage layers of the medial tibial plateau, the lateral tibial plateau, the femur and the patella equaled respectively 0.32 mm (range 0.14–0.48), 0.35 mm (range 0.16–0.53), 0.39 mm (range 0.15–0.80) and 0.75 mm (range 0.16–1.11). Similarly, the RMSE equaled respectively 1.16 mm (range 0.99–1.59), 0.91 mm (0.75–1.33), 2.93 mm (range 1.85–4.66) and 2.04 mm (1.88–3.29), calculated over the course of the anterior cruciate ligament, posterior cruciate ligament, the medial and the lateral meniscus.

Conclusion: A methodological workflow is presented for patient-specific, morphological knee joint modeling that avoids laborious segmentation. By allowing to accurately predict personalized geometry this method has the potential for generating large (virtual) sample sizes applicable for biomechanical research and improving personalized, computer-assisted medicine.

Statistical modeling: Assessing the anatomic variability of knee joint space width

Population-based knee joint space width (JSW) assessments are promising for the prevention and early diagnosis of osteoarthritis. This study aimed to establish the statistical shape and alignment model (SSAM) of knee joints for assessing anatomic variation in knee JSW in the healthy Chinese male population. CT scans of asymptomatic knee joints of healthy male participants (n = 107) were collected for manual segmentation to create mesh samples. The as-scanned positional error was reduced by a standard processing flow of deformable mesh registration. Principal component analysis (PCA) was performed to create a tibiofemoral SSAM that was trained on all mesh samples. The anatomic variability of the JSW in the healthy Chinese male population was then assessed using the SSAM with regression analysis and 3D analysis by color-coded mapping. Almost all PCA modes had a linear influence on the anatomic variation of the medial and lateral JSW. The JSW variability within the SSAM was mainly explained by mode 1 (45.1 % of variation), demonstrating that this mode had the greatest influence on JSW variation. 3D assessment of the JSW showed that the minimum medial JSW varied from 2.76 to 3.23 mm, and its site shifted a short distance on the medial tibial plateau. The root-mean-square fitting and generalization errors of the SSAM were below 1 mm. This study will benefit the design and optimization of prosthetic devices, and may be applicable to the prevention and early diagnosis of osteoarthritis.

Imaging in Osteoarthritis

Osteoarthritis (OA) is the most frequent form of arthritis with major implications on both individual and public health care levels. The field of joint imaging, and particularly magnetic resonance imaging (MRI), has evolved rapidly due to the application of technical advances to the field of clinical research. This narrative review will provide an introduction to the different aspects of OA imaging aimed at an audience of scientists, clinicians, students, industry employees, and others who are interested in OA but who do not necessarily focus on OA. The current role of radiography and recent advances in measuring joint space width will be discussed. The status of cartilage morphology assessment and evaluation of cartilage biochemical composition will be presented. Advances in quantitative three-dimensional morphologic cartilage assessment and semi-quantitative whole-organ assessment of OA will be reviewed. Although MRI has evolved as the most important imaging method used in OA research, other modalities such as ultrasound, computed tomography, and metabolic imaging play a complementary role and will also be discussed.

Artificial Intelligence Based Patient-Specific Preoperative Planning Algorithm for Total Knee Arthroplasty

Previous studies have shown that the manufacturer's default preoperative plans for total knee arthroplasty with patient-specific guides require frequent, time-consuming changes by the surgeon. Currently, no research has been done on predicting preoperative plans for orthopedic surgery using machine learning. Therefore, this study aims to evaluate whether artificial intelligence (AI) driven planning tools can create surgeon and patient-specific preoperative plans that require fewer changes by the surgeon. A dataset of 5409 preoperative plans, including the manufacturer's default and the plans corrected by 39 surgeons, was collected. Features were extracted from the preoperative plans that describe the implant sizes, position, and orientation in a surgeon- and patient-specific manner. Based on these features, non-linear regression models were employed to predict the surgeon's corrected preoperative plan. The average number of corrections a surgeon has to make to the preoperative plan generated using AI was reduced by 39.7% compared to the manufacturer's default plan. The femoral and tibial implant size in the manufacturer's plan was correct in 68.4% and 73.1% of the cases, respectively, while the AI-based plan was correct in 82.2% and 85.0% of the cases, respectively, compared to the surgeon approved plan. Our method successfully demonstrated the use of machine learning to create preoperative plans in a surgeon- and patient-specific manner for total knee arthroplasty.

Using Statistical Shape Models to Optimize TKA Implant Design

(1) TKA implants should well fit on each patient’s anatomy. Statistical Shape Models (SSM) statistically represent the anatomy of a given population. The aims of this study were to assess how to generate a valid SSM for implant design and provide guidelines and examples on how to use the SSMs to evaluate the anatomic fit of TKA components. (2) Methods: A Caucasian SSM was built from 120 anatomies (65 female, 55 male) and an Asian SSM was based on 112 patients (75 female, 37 male). These SSMs were used to generate a database of 20 bone models. The AP/ML dimensions of the bone models were compared to those of the input population. Design input parameters, such as the tibial contour, trochlea, and femur curvature were extracted from the SSMs. Femur and patella components were virtually implanted on the bone models. (3) Results: the dimensions of the generated bone models well represented the population. The overhang of the femoral component as well as the coverage and peak restoration of the patella component were visualized. (4) Conclusions: SSMs can be used to efficiently gain input into TKA design and evaluate the implant fit on the studied population.

Intelligent Segmentation Algorithm for Diagnosis of Meniere's Disease in the Inner Auditory Canal Using MRI Images with Three-Dimensional Level Set

This paper aimed to explore segmentation effects of the magnetic resonance imaging (MRI) images of the inner auditory canal of patients with Meniere's disease under the intelligent segmentation method of the inner ear based on three-dimensional (3D) level set (IS3DLS). The statistical shape model and the level set segmentation algorithm were combined to propose the IS3DLS. First, the shape training samples of the inner ear model were determined, and the results were manually segmented to further obtain region of interest (ROI) of the inner ear. The IS3DLS was employed to accurately segment MRI images of the inner auditory canal of patients with Meniere's disease. The segmentation performance of IS3DLS was compared with the expert manual segmentation method and the region growth level set-based segmentation algorithm. Results showed that Matthews correlation coefficient (MCC), Dice similarity coefficient (DSC), false positive rate (FPR), and false negative rate (FNR) of this algorithm were 0.9599, 0.9594, 0.0325, and 0.03655, respectively. Therefore, the IS3DLS could achieve good segmentation effect in MRI images of the inner auditory canal of patients with Meniere's disease, which was helpful for diagnosis and subsequent treatment of Meniere's disease.

Predicting Knee Joint Instability Using a Tibio-Femoral Statistical Shape Model

Statistical shape models (SSMs) are a well established computational technique to represent the morphological variability spread in a set of matching surfaces by means of compact descriptive quantities, traditionally called "modes of variation" (MoVs). SSMs of bony surfaces have been proposed in biomechanics and orthopedic clinics to investigate the relation between bone shape and joint biomechanics. In this work, an SSM of the tibio-femoral joint has been developed to elucidate the relation between MoVs and bone angular deformities causing knee instability. The SSM was built using 99 bony shapes (distal femur and proximal tibia surfaces obtained from segmented CT scans) of osteoarthritic patients. Hip-knee-ankle (HKA) angle, femoral varus-valgus (FVV) angle, internal-external femoral rotation (IER), tibial varus-valgus (TVV) angles, and tibial slope (TS) were available across the patient set. Discriminant analysis (DA) and logistic regression (LR) classifiers were adopted to underline specific MoVs accounting for knee instability. First, it was found that thirty-four MoVs were enough to describe 95% of the shape variability in the dataset. The most relevant MoVs were the one encoding the height of the femoral and tibial shafts (MoV #2) and the one representing variations of the axial section of the femoral shaft and its bending in the frontal plane (MoV #5). Second, using quadratic DA, the sensitivity results of the classification were very accurate, being all >0.85 (HKA: 0.96, FVV: 0.99, IER: 0.88, TVV: 1, TS: 0.87). The results of the LR classifier were mostly in agreement with DA, confirming statistical significance for MoV #2 (p = 0.02) in correspondence to IER and MoV #5 in correspondence to HKA (p = 0.0001), FVV (p = 0.001), and TS (p = 0.02). We can argue that the SSM successfully identified specific MoVs encoding ranges of alignment variability between distal femur and proximal tibia. This discloses the opportunity to use the SSM to predict potential misalignment in the knee for a new patient by processing the bone shapes, removing the need for measuring clinical landmarks as the rotation centers and mechanical axes.