Hamilton-Jacobi Skeletons

  • Authors:
  • Kaleem Siddiqi;Sylvain Bouix;Allen Tannenbaum;Steven W. Zucker

  • Affiliations:
  • School of Computer Science & Center for Intelligent Machines, McGill University, Montreal, QC H3A 2A7, Canada. siddiqi@cim.mcgill.ca;School of Computer Science & Center for Intelligent Machines, McGill University, Montreal, QC H3A 2A7, Canada. sbouix@cim.mcgill.ca;Department of Electrical & Computer Engineering and Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA 30082-0250, USA. tannenba@ece.gatech.edu;Departments of Computer Science and Electrical Engineering and Center for Computational Vision & Control, Yale University, New Haven, CT 06520-8285, USA. zucker-steven@cs.yale.edu

  • Venue:
  • International Journal of Computer Vision
  • Year:
  • 2002

Quantified Score

Hi-index 0.00

Visualization

Abstract

The eikonal equation and variants of it are of significant interest for problems in computer vision and image processing. It is the basis for continuous versions of mathematical morphology, stereo, shape-from-shading and for recent dynamic theories of shape. Its numerical simulation can be delicate, owing to the formation of singularities in the evolving front and is typically based on level set methods. However, there are more classical approaches rooted in Hamiltonian physics which have yet to be widely used by the computer vision community. In this paper we review the Hamiltonian formulation, which offers specific advantages when it comes to the detection of singularities or shocks. We specialize to the case of Blum's grassfire flow and measure the average outward flux of the vector field that underlies the Hamiltonian system. This measure has very different limiting behaviors depending upon whether the region over which it is computed shrinks to a singular point or a non-singular one. Hence, it is an effective way to distinguish between these two cases. We combine the flux measurement with a homotopy preserving thinning process applied in a discrete lattice. This leads to a robust and accurate algorithm for computing skeletons in 2D as well as 3D, which has low computational complexity. We illustrate the approach with several computational examples.