Shape Analysis Using a Point-Based Statistical Shape Model Built on Correspondence Probabilities

A fundamental problem when computing statistical shape models is the determination of correspondences between the instances of the associated data set. Often, homologies between points that represent the surfaces are assumed which might lead to imprecise mean shape and variability results. We propose an approach where exact correspondences are replaced by evolving correspondence probabilities. These are the basis for a novel algorithm that computes a generative statistical shape model. We developed an unified MAP framework to compute the model parameters ('mean shape' and 'modes of variation') and the nuisance parameters which leads to an optimal adaption of the model to the set of observations. The registration of the model on the instances is solved using the Expectation Maximization--Iterative Closest Point algorithm which is based on probabilistic correspondences and proved to be robust and fast. The alternated optimization of the MAP explanation with respect to the observation and the generative model parameters leads to very efficient and closed-form solutions for (almost) all parameters. Experimental results on brain structure data sets demonstrate the efficiency and well-posedness of the approach. The algorithm is then extended to an automatic classification method using the k-means clustering and applied to synthetic data as well as brain structure classification problems.