Demo

Femoroacetabular impingement (FAI), is a syndrome when the femoral head ball runs abnormally or prohibits a normal range of motion in the acetabular socket. The cam-type FAI is illustrated in the left subfigure above. Right two subfigures above show a group of femur shapes and the constructed Statistical Shape Model (SSM) in the form of modes.

Figure above shows the first three modes of the SSM for the HEALTHY femur group, where the major varation patterns are captured as: Mode 1: femur head ball enlarging and the trochanter major shrinking (4.3mm); Mode 2: additional thickening pattern of intertrochanteric valley (1.2 mm); Mode 3: bone material depositing at femur neck near trochanter major(0.7mm)

Figure above shows the first three modes of the SSM for the UNHEALTHY femur group, where the major varation patterns are captured as: Mode 1: femur head ball enlarging and the trochanter major shrinking (3.9mm); Mode 2: head ball enlarging (1.7 mm); Mode 3: bone material depositing at femur neck near head (1mm)

Figure above shows the mean shape of the UNHEALTHY femur group on top of that of the HEALTHY group, in both superimposed view and shell distance field view. It is observed that the bone material 0.7mm in thickness of the unhealthy mean on top of the healthy mean shape is a probable cause of cam impingement among the unhealthy group.

Abstract

This dissertation proposes an efficient optimization approach for obtaining shape correspondence across a group of objects for statistical shape modeling. With each shape represented in a B-spline based parametric form, the correspondence across the shape population is cast as an issue of seeking a reparametrization for each shape so that a quality measure of the resulting shape correspondence across the group is optimized. The quality measure is the description length of covariance matrix of the shape population, with landmarks sampled on each shape. The movement of landmarks on each B-spline shape is controlled by the reparameterization of the B-spline shape. The reparameterization itself is also represented with B-splines and B-spline coefficients are used as optimization parameters. We have developed formulations for ensuring the bijectivity of the reparameterization. A gradient-based optimization approach is developed, including techniques such as constraint aggregation and adjoint senstivity for efficient, direct di↵eomorphic reparameterization of landmarks to improve the group-wise shape correspondence. Numerical experiments on both synthetic and real 2D and 3D data sets demonstrate the efficiency and e↵ectiveness of the proposed approach.