Characterizing the shape of anatomical structures with Poisson's equation

Poisson's equation, a fundamental partial differential equation in classical physics, has a number of properties that are interesting for shape analysis. In particular, the equipotential sets of the solution graph become smoother as the potential increases. We use the displacement map, the length of the streamlines formed by the gradient field of the solution, to measure the "complexity" (or smoothness) of the equipotential sets, and study its behavior as the potential increases. We believe that this function complexity=f(potential), which we call the shape characteristic, is a very natural way to express shape. Robust algorithms are presented to compute the solution to Poisson's equation, the displacement map, and the shape characteristic. We first illustrate our technique on two-dimensional synthetic examples and natural silhouettes. We then perform two shape analysis studies on three-dimensional neuroanatomical data extracted from magnetic resonance (MR) images of the brain. In the first study, we investigate changes in the caudate nucleus in Schizotypal Personality Disorder (SPD) and confirm previously published results on this structure . In the second study, we present a data set of caudate nuclei of premature infants with asymmetric white matter injury. Our method shows structural shape differences that volumetric measurements were unable to detect

[1]  Fred L. Bookstein,et al.  “Voxel-Based Morphometry” Should Not Be Used with Imperfectly Registered Images , 2001, NeuroImage.

[2]  Ali Shokoufandeh,et al.  Shock Graphs and Shape Matching , 1998, International Journal of Computer Vision.

[3]  G. Pinder,et al.  Numerical solution of partial differential equations in science and engineering , 1982 .

[4]  D. Louis Collins,et al.  Hippocampal shape analysis using medial surfaces , 2005, NeuroImage.

[5]  Martin Styner,et al.  Medial Models Incorporating Object Variability for 3D Shape Analysis , 2001, IPMI.

[6]  Thomas E. Nichols,et al.  Nonparametric permutation tests for functional neuroimaging: A primer with examples , 2002, Human brain mapping.

[7]  Claes Johnson Numerical solution of partial differential equations by the finite element method , 1988 .

[8]  Olivier D. Faugeras,et al.  Statistical shape influence in geodesic active contours , 2000, Proceedings IEEE Conference on Computer Vision and Pattern Recognition. CVPR 2000 (Cat. No.PR00662).

[9]  D. Louis Collins,et al.  Hippocampal shape analysis using medial surfaces , 2001, NeuroImage.

[10]  W. Eric L. Grimson,et al.  Detection and analysis of statistical differences in anatomical shape , 2005, Medical Image Anal..

[11]  William M. Wells,et al.  An EM algorithm for shape classification based on level sets , 2005, Medical Image Anal..

[12]  Gábor Székely,et al.  Multiscale Medial Loci and Their Properties , 2003, International Journal of Computer Vision.

[13]  T. Inder,et al.  Defining the nature of the cerebral abnormalities in the premature infant: a qualitative magnetic resonance imaging study. , 2003, Jornal de Pediatria.

[14]  Karl J. Friston,et al.  Why Voxel-Based Morphometry Should Be Used , 2001, NeuroImage.

[15]  R. Kikinis,et al.  Periventricular white matter injury in the premature infant is followed by reduced cerebral cortical gray matter volume at term , 1999, Annals of neurology.

[16]  D. Louis Collins,et al.  Retrospective evaluation of intersubject brain registration , 2003, IEEE Transactions on Medical Imaging.

[17]  N. Andreasen,et al.  Caudate size in first-episode neuroleptic-naive schizophrenic patients measured using an artificial neural network , 1999, Biological Psychiatry.

[18]  D. Louis Collins,et al.  3D Moment Invariant Based Morphometry , 2003, MICCAI.

[19]  L J Abernethy,et al.  Quantitative magnetic resonance imaging of the brain in survivors of very low birth weight , 2002, Archives of disease in childhood.

[20]  Martin Styner,et al.  Boundary and Medial Shape Analysis of the Hippocampus in Schizophrenia , 2003, MICCAI.

[21]  Ronen Basri,et al.  Shape Representation and Classification Using the Poisson Equation , 2006, IEEE Trans. Pattern Anal. Mach. Intell..

[22]  Jerry L. Prince,et al.  An Eulerian PDE approach for computing tissue thickness , 2003, IEEE Transactions on Medical Imaging.

[23]  R. Kikinis,et al.  MRI study of caudate nucleus volume and its cognitive correlates in neuroleptic-naive patients with schizotypal personality disorder. , 2002, The American journal of psychiatry.

[24]  Paul M. Thompson,et al.  Mapping hippocampal and ventricular change in Alzheimer disease , 2004, NeuroImage.

[25]  Martin Styner,et al.  Statistical shape analysis of neuroanatomical structures based on medial models , 2003, Medical Image Anal..

[26]  Guido Gerig,et al.  Parametrization of Closed Surfaces for 3-D Shape Description , 1995, Comput. Vis. Image Underst..

[27]  Haissam Haidar,et al.  New Numerical Solution of the Laplace Equation for Tissue Thickness Measurement in Three-Dimensional MRI , 2005, J. Math. Model. Algorithms.

[28]  H. Blum Biological shape and visual science (part I) , 1973 .

[29]  Timothy F. Cootes,et al.  Active Shape Models-Their Training and Application , 1995, Comput. Vis. Image Underst..

[30]  Raúl San José Estépar,et al.  Shape of caudate nucleus and its cognitive correlates in neuroleptic-naive schizotypal personality disorder , 2003, Biological Psychiatry.

[31]  John Suckling,et al.  Global, voxel, and cluster tests, by theory and permutation, for a difference between two groups of structural MR images of the brain , 1999, IEEE Transactions on Medical Imaging.

[32]  J. Volpe Neurology of the Newborn , 1959, Major problems in clinical pediatrics.

[33]  Christopher J. Cannistraci,et al.  Regional brain volume abnormalities and long-term cognitive outcome in preterm infants. , 2000, JAMA.

[34]  J. Sweeney,et al.  Decreased caudate volume in neuroleptic-naive psychotic patients. , 1998, The American journal of psychiatry.

[35]  I. Aharon,et al.  Three‐dimensional mapping of cortical thickness using Laplace's Equation , 2000, Human brain mapping.

[36]  U. Grenander,et al.  Hippocampal morphometry in schizophrenia by high dimensional brain mapping. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[37]  William H. Press,et al.  Numerical recipes , 1990 .

[38]  E. Ziegel Permutation, Parametric, and Bootstrap Tests of Hypotheses (3rd ed.) , 2005 .

[39]  David G. Stork,et al.  Pattern Classification (2nd ed.) , 1999 .