Posterior shape models

We present a method to compute the conditional distribution of a statistical shape model given partial data. The result is a "posterior shape model", which is again a statistical shape model of the same form as the original model. This allows its direct use in the variety of algorithms that include prior knowledge about the variability of a class of shapes with a statistical shape model. Posterior shape models then provide a statistically sound yet easy method to integrate partial data into these algorithms. Usually, shape models represent a complete organ, for instance in our experiments the femur bone, modeled by a multivariate normal distribution. But because in many application certain parts of the shape are known a priori, it is of great interest to model the posterior distribution of the whole shape given the known parts. These could be isolated landmark points or larger portions of the shape, like the healthy part of a pathological or damaged organ. However, because for most shape models the dimensionality of the data is much higher than the number of examples, the normal distribution is singular, and the conditional distribution not readily available. In this paper, we present two main contributions: First, we show how the posterior model can be efficiently computed as a statistical shape model in standard form and used in any shape model algorithm. We complement this paper with a freely available implementation of our algorithms. Second, we show that most common approaches put forth in the literature to overcome this are equivalent to probabilistic principal component analysis (PPCA), and Gaussian Process regression. To illustrate the use of posterior shape models, we apply them on two problems from medical image analysis: model-based image segmentation incorporating prior knowledge from landmarks, and the prediction of anatomically correct knee shapes for trochlear dysplasia patients, which constitutes a novel medical application. Our experiments confirm that the use of conditional shape models for image segmentation improves the overall segmentation accuracy and robustness.

[1]  B. Stuyts,et al.  Trochleoplasty in dysplastic knee trochlea , 2005, Knee Surgery, Sports Traumatology, Arthroscopy.

[2]  Michael E. Tipping,et al.  Probabilistic Principal Component Analysis , 1999 .

[3]  Thomas Vetter,et al.  Reconstructing the Complete 3D Shape of Faces from Partial Information (Rekonstruktion der dreidimensionalen Form von Gesichtern aus partieller Information) , 2002, Informationstechnik Tech. Inform..

[4]  Ron Kikinis,et al.  3D Slicer , 2012, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[5]  Sam T. Roweis,et al.  EM Algorithms for PCA and SPCA , 1997, NIPS.

[6]  C. Taylor,et al.  Active shape models - 'Smart Snakes'. , 1992 .

[7]  Nicholas Ayache,et al.  Rigid, affine and locally affine registration of free-form surfaces , 1996, International Journal of Computer Vision.

[8]  Nicholas Ayache,et al.  Non-parametric Diffeomorphic Image Registration with the Demons Algorithm , 2007, MICCAI.

[9]  Thomas Vetter,et al.  Modeling the Remaining Flexibility of Partially Fixed Statistical Shape Models , 2008 .

[10]  Gábor Székely,et al.  Conditional Variability of Statistical Shape Models Based on Surrogate Variables , 2009, MICCAI.

[11]  Gábor Székely,et al.  Statistical model based shape prediction from a combination of direct observations and various surrogates: Application to orthopaedic research , 2012, Medical Image Anal..

[12]  Dinggang Shen,et al.  Predictive modeling of anatomic structures using canonical correlation analysis , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[13]  U. Grenander,et al.  Computational anatomy: an emerging discipline , 1998 .

[14]  Stefan Klein,et al.  Conditional Shape Models for Cardiac Motion Estimation , 2010, MICCAI.

[15]  Alejandro F Frangi,et al.  Automatic construction of 3-D statistical deformation models of the brain using nonrigid registration , 2003, IEEE Transactions on Medical Imaging.

[16]  Akinobu Shimizu,et al.  Relaxed Conditional Statistical Shape Models and Their Application to Non-contrast Liver Segmentation , 2012, Abdominal Imaging.

[17]  Nicholas Ayache,et al.  The Correlation Ratio as a New Similarity Measure for Multimodal Image Registration , 1998, MICCAI.

[18]  Dinggang Shen,et al.  Spatial normalization of spine MR images for statistical correlation of lesions with clinical symptoms. , 2002, Radiology.

[19]  Philippe C. Cattin,et al.  Statismo - A framework for PCA based statistical models , 2012, The Insight Journal.

[20]  Marleen de Bruijne,et al.  Confidence of model based shape reconstruction from sparse data , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[21]  Thomas Vetter,et al.  Probabilistic Modeling and Visualization of the Flexibility in Morphable Models , 2009, IMA Conference on the Mathematics of Surfaces.

[22]  Marco Zanetti,et al.  Femoral Trochlear Dysplasia : MR Findings 1 , 2000 .

[23]  Weiguo Lu,et al.  A simple fixed-point approach to invert a deformation field. , 2007, Medical physics.

[24]  J Romero,et al.  Femoral trochlear dysplasia: MR findings. , 2000, Radiology.

[25]  Nasser M. Nasrabadi,et al.  Pattern Recognition and Machine Learning , 2006, Technometrics.

[26]  Nikos Paragios,et al.  Matching Distance Functions: A Shape-to-Area Variational Approach for Global-to-Local Registration , 2002, ECCV.

[27]  Berthold K. P. Horn,et al.  Closed-form solution of absolute orientation using unit quaternions , 1987 .

[28]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[29]  Michael I. Miller,et al.  On The Geometry and Shape of Brain Sub-Manifolds , 1997, Int. J. Pattern Recognit. Artif. Intell..

[30]  Thomas Vetter,et al.  A morphable model for the synthesis of 3D faces , 1999, SIGGRAPH.

[31]  Mads Nielsen,et al.  Conditional Point Distribution Models , 2010, MCV.

[32]  Max A. Viergever,et al.  elastix: A Toolbox for Intensity-Based Medical Image Registration , 2010, IEEE Transactions on Medical Imaging.

[33]  K. Strimmer,et al.  Statistical Applications in Genetics and Molecular Biology A Shrinkage Approach to Large-Scale Covariance Matrix Estimation and Implications for Functional Genomics , 2011 .

[34]  Marleen de Bruijne,et al.  Quantitative vertebral morphometry using neighbor-conditional shape models , 2007, Medical Image Anal..

[35]  Timothy F. Cootes,et al.  Constrained active appearance models , 2001, Proceedings Eighth IEEE International Conference on Computer Vision. ICCV 2001.

[36]  G. Strang.  Introduction to Linear Algebra , 1993 .

[37]  Hans-Peter Meinzer,et al.  Statistical shape models for 3D medical image segmentation: A review , 2009, Medical Image Anal..

[38]  Gábor Székely,et al.  Confidence Regions for Statistical Model Based Shape Prediction From Sparse Observations , 2012, IEEE Transactions on Medical Imaging.