Landmark-free morphometric analysis of knee osteoarthritis using joint statistical models of bone shape and articular space variability

Purpose: Osteoarthritis (OA) is a common degenerative disease involving a variety of structural changes in the affected joint. In addition to narrowing of the articular space, recent studies involving statistical shape analysis methods have suggested that specific bone shapes might be associated with the disease. We aim to investigate the feasibility of using the recently introduced framework of functional shapes (Fshape) to extract morphological features of OA that combine shape variability of articular surfaces of the tibia (or femur) together with the changes of the joint space. Approach: Our study uses a dataset of three-dimensional cone-beam CT volumes of 17 knees without OA and 17 knees with OA. Each knee is then represented as an object (Fshape) consisting of a triangulated tibial (or femoral) articular surface and a map of joint space widths (JSWs) measured at the points of this surface (joint space map, JSM). We introduce a generative atlas model to estimate a template (mean) Fshape of the sample population together with template-centered variables that model the transformations from the template to each subject. This approach has two potential advantages compared with other statistical shape modeling methods that have been investigated in knee OA: (i) Fshapes simultaneously consider the variability in bone shape and JSW, and (ii) Fshape atlas estimation is based on a diffeomorphic transformation model of surfaces that does not require a priori landmark correspondences between the subjects. The estimated atlas-to-subject Fshape transformations were used as input to principal component analysis dimensionality reduction combined with a linear support vector machine (SVM) classifier to identify the morphological features of OA. Results: Using tibial articular surface as the shape component of the Fshape, we found leave-one-out cross validation scores of ≈91.18 % for the classification based on the bone surface transformations alone, ≈91.18 % for the classification based on the residual JSM, and ≈85.29 % for the classification using both Fshape components. Similar results were obtained using femoral articular surfaces. The discriminant directions identified in the statistical analysis were associated with medial narrowing of the joint space, steeper intercondylar eminence, and relative deepening of the medial tibial plateau. Conclusions: The proposed approach provides an integrated framework for combined statistical analysis of shape and JSPs. It can successfully extract features correlated to OA that appear consistent with previous studies in the field. Although future large-scale study is necessary to confirm the significance of these findings, our results suggest that the functional shape methodology is a promising new tool for morphological analysis of OA and orthopedics data in general.