## Hamilton-Jacobi Skeletons (1999)

### Cached

### Download Links

- [www.cs.mcgill.ca]
- [www.cim.mcgill.ca]
- [www.cim.mcgill.ca]
- [pages.cs.wisc.edu]
- [www.cim.mcgill.ca]
- [www.cim.mcgill.ca]
- DBLP

### Other Repositories/Bibliography

Citations: | 119 - 12 self |

### BibTeX

@MISC{Siddiqi99hamilton-jacobiskeletons,

author = {Kaleem Siddiqi and Sylvain Bouix and Allen Tannenbaum and Steven W. Zucker},

title = {Hamilton-Jacobi Skeletons},

year = {1999}

}

### Years of Citing Articles

### OpenURL

### 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 grass fire flow and measure the average outward ux 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 ux 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.

### Citations

873 | Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations
- Osher, Sethian
- 1988
(Show Context)
Citation Context ...omputer vision since, e.g., it is they which denote the skeleton (see Figure 3). To continue the evolution while preserving shocks, the technology of level set methods introduced by Osher and Sethian =-=[15]-=-, has proved to be extremely powerful. The approach relies on the notion of a weak solution, developed in viscosity theory [6], and the introduction of an appropriate entropy condition to select it. T... |

628 | User’s Guide to Viscosity Solutions of Second Order Partial Differential Equations, volume 27
- Crandall, Ishii, et al.
- 1992
(Show Context)
Citation Context ...cks, the technology of level set methods introduced by Osher and Sethian [15], has proved to be extremely powerful. The approach relies on the notion of a weak solution, developed in viscosity theory =-=[6]-=-, and the introduction of an appropriate entropy condition to select it. The representation of the evolving front as a level set of a hypersurface allows topological changes to be handled in a natural... |

427 | A fast marching level set method for monotonically advancing fronts
- Sethian
- 1996
(Show Context)
Citation Context ...The representation of the evolving front as a level set of a hypersurface allows topological changes to be handled in a natural way, and robust, efficient implementations have recently been developed =-=[20]-=-. Level set methods are Eulerian in nature because computations are restricted to grid points whose locations are fixed. For such methods, the question of computing the locus of shocks for dynamically... |

358 | A level set approach for computing solutions to incompressible two-phase flow - Sussman, Smereka, et al. - 1994 |

318 | Shape from shading - Brooks - 1989 |

204 | Digital topology: introduction and survey - Kong, Rosenfeld - 1989 |

203 | Differential equations and dynamical systems - Perko - 2001 |

203 | Shock graphs and shape matching - Siddiqi, Shokoufandeh, et al. - 1999 |

197 |
A viscosity solutions approach to shape-from-shading
- Rouy, Tourin
- 1992
(Show Context)
Citation Context ...'s fast marching method [20], which relies on an interpretation of Huygens's principle to efficiently propagate the solution from the initial curve, and Rouy and Tourin's viscosity solutions approach =-=[17]-=-. However, neither of these methods address the issue of shock detection explicitly, and more work has to be done to track shocks. A different approach, which is related to the solution surfacesT (x; ... |

196 | Distance transformations in arbitrary dimensions - Borgefors - 1984 |

190 | Set Methods: Evolving Interfaces - Level - 1996 |

181 |
The Variational Principles Of Mechanics
- Lanczos
- 1957
(Show Context)
Citation Context ...roach with several numerical examples including skeletons of complex 2D and 3D shapes. 1. Introduction Variational principles emerged naturally from considerations of energy minimization in mechanics =-=[11]-=-. We consider these in the context of the eikonal equation, which arises in geometrical optics and, recently, which has become of great interest for problems in computer vision [4]. It is the basis fo... |

146 |
High-order essentially nonoscillatory schemes for Hamilton–Jacobi equations
- Osher, Shu
- 1991
(Show Context)
Citation Context ...ction methods which rely on interpolation of the underlying hypersurface are computationally very expensive. Numerical thresholds are introduced and high order accurate numerical schemes must be used =-=[14, 22]-=-. On the other hand, there are more classical methods rooted in Hamiltonian physics, which can also be used to study shock theory. To the best of our knowledge, these have not been considered in the c... |

144 |
Biological shape and visual science
- Blum
(Show Context)
Citation Context ...ntly, which has become of great interest for problems in computer vision [4]. It is the basis for continuous versions of mathematical morphology [3, 18, 25], as well as for Blum's grassfire transform =-=[2]-=- and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, 5], shape-fromshading [10] and stereo [8]. The ... |

122 | Hierarchic voronoi skeletons
- Ogniewicz, Kubler
- 1995
(Show Context)
Citation Context ...ing surface evolution based segmentation techniques on volumetric MR data. Whereas alternate approaches based on Voronoi techniques provide a topologically organized set of skeletal branches or faces =-=[13, 12]-=-, their complexity increases with the number of points on the bounding curve or surface. More significantly, heuristic measures of significance are used for pruning, in order to obtain reasonable resu... |

104 | Complete dense stereovision using level set methods
- Faugeras, Keriven
- 1998
(Show Context)
Citation Context ...sform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, 5], shape-fromshading [10] and stereo =-=[8-=-]. The numerical simulation of this equation is non-trivial, School of Computer Science & Center for Intelligent Machines McGill University, Montreal, PQ, Canada H3A 2A7 y School of Computer Science ... |

76 |
Level set methods: evolving interfaces in geometry, fluid mechanics, computer vision and material science
- Sethian
- 1996
(Show Context)
Citation Context ..., 25], as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis =-=[19, 5-=-], shape-fromshading [10] and stereo [8]. The numerical simulation of this equation is non-trivial, School of Computer Science & Center for Intelligent Machines McGill University, Montreal, PQ, Canad... |

76 | Symmetry-Based Indexing of Image Databases - Sharvit, Chan, et al. - 1998 |

72 | Simulating the grassfire transform using an active contour model - Leymarie, Levine - 1992 |

69 | A common framework for curve evolution, segmentation, and anisotropic diffusion
- Shah
- 1996
(Show Context)
Citation Context ...of shock detection explicitly, and more work has to be done to track shocks. A different approach, which is related to the solution surfacesT (x; y) viewed as a graph, has been proposed by Shah et al =-=[21, 23]-=-. Here the key idea is to use an edge strength q o q q . o q S (q) = t D S p = Figure 2: Direction of a ray _ q and the direction of motion of the wave front p. From [1]. functional v in place of the ... |

63 | Shape, shocks, and deformations {I}: The components of two-dimensional shape and the reaction-diffusion space - Kimia, Tannenbaum, et al. - 1995 |

62 | Building skeleton models via 3-D medial surface/axis thinning algorithms - Lee, Kashyap, et al. - 1994 |

60 | Principles of quantum mechanics - Shankar - 1994 |

59 | Mathematical Methods of Classical Mechanics, 2nd edn - Arnold - 1989 |

57 |
Evolution equations for continuous-scale morphology
- Brockett, Maragos
- 1992
(Show Context)
Citation Context ... equation, which arises in geometrical optics and, recently, which has become of great interest for problems in computer vision [4]. It is the basis for continuous versions of mathematical morphology =-=[3, 18, 25]-=-, as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, ... |

55 | Distance-ordered homotopic thinning: A skeletonization algorithm for 3d digital images - Pudney - 1998 |

52 | Approximate tree matching and shape similarity - Liu, Geiger - 1999 |

49 | Topological segmentation of discrete surfaces - Malandain, Bertrand, et al. - 1993 |

48 |
The eikonal equation: Some results applicable to computer vision
- Bruss
- 1982
(Show Context)
Citation Context ...ation in mechanics [11]. We consider these in the context of the eikonal equation, which arises in geometrical optics and, recently, which has become of great interest for problems in computer vision =-=[4]-=-. It is the basis for continuous versions of mathematical morphology [3, 18, 25], as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It h... |

46 | An optimal control formulation and related numerical methods for a problem in shape reconstruction - Dupuis, Oliensis - 1994 |

42 |
Implementing continuousscale morphology via curve evolution
- Sapiro, Kimmel, et al.
- 1993
(Show Context)
Citation Context ... equation, which arises in geometrical optics and, recently, which has become of great interest for problems in computer vision [4]. It is the basis for continuous versions of mathematical morphology =-=[3, 18, 25]-=-, as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, ... |

39 | S.: Euclidean Skeletons - Malandain, Fernández-Vidal - 1998 |

38 | An algorithm for the medial axis transform of 3d polyhedral solids - Sherbrooke, Patrikalakis, et al. - 1996 |

34 | Representation and self-similarity of shapes - Geiger, Liu, et al. |

33 | A parallel thinning algorithm for medial surfaces - Bertrand - 1995 |

32 | Robust and efficient skeletal graphs - Dimitrov, Phillips, et al. - 2000 |

31 | C.-W.: Geometric shock-capturing ENO schemes for subpixel interpolation, computation and curve evolution. Graphical Models and Image Processing 59
- SIDDIQI, KIMIA, et al.
- 1997
(Show Context)
Citation Context ...ction methods which rely on interpolation of the underlying hypersurface are computationally very expensive. Numerical thresholds are introduced and high order accurate numerical schemes must be used =-=[14, 22]-=-. On the other hand, there are more classical methods rooted in Hamiltonian physics, which can also be used to study shock theory. To the best of our knowledge, these have not been considered in the c... |

30 |
The morphological structure of images: the differential equations of morphological scale–space
- Boomgaard, Smeulders
- 1994
(Show Context)
Citation Context ... equation, which arises in geometrical optics and, recently, which has become of great interest for problems in computer vision [4]. It is the basis for continuous versions of mathematical morphology =-=[3, 18, 25]-=-, as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, ... |

29 | Assisted articulation of closed polygonal models - Teichmann, Teller - 1998 |

29 | 3d skeleton and centerline generation based on an approximate minimum distance field,” The Visual - Zhou, Kaufman, et al. - 1998 |

26 |
Extraction of shape skeletons from grayscale images
- Tari, Shah, et al.
- 1997
(Show Context)
Citation Context ...er vision [4]. It is the basis for continuous versions of mathematical morphology [3, 18, 25], as well as for Blum's grassfire transform [2] and new dynamic theories of shape representation including =-=[9, 23-=-]. It has also been widely used for applications in image processing and analysis [19, 5], shape-fromshading [10] and stereo [8]. The numerical simulation of this equation is non-trivial, School of C... |

24 | di Baja, A width independent fast thinning algorithm - Arcelli, S - 1985 |

24 | Discrete Voronoi Skeletons, Hartung-Gorre Verlag - Ogniewicz - 1993 |

24 | Shape description by medial surface construction - Sheehy, Armstrong, et al. - 1996 |

23 | G.S.: A width-independent fast thinning algorithm - Arcelli, Baja - 1985 |

22 | ªTopological Numbers and Singularities in Scalar Images Scale-Space Evolution Properties,º Gaussian Scale-Space Theory - Kalitzin - 1997 |

22 | B.: Medial faces from a concise 3D thinning algorithm - Manzanera, Bernard, et al. - 1999 |

21 | M-reps: A new object representation for graphics - Pizer, Thall, et al. |

20 | Shape from shading: Level set propagation and viscosity solutions
- Kimmel, Siddiqqi, et al.
- 1995
(Show Context)
Citation Context ...s grassfire transform [2] and new dynamic theories of shape representation including [9, 23]. It has also been widely used for applications in image processing and analysis [19, 5], shape-fromshading =-=[10-=-] and stereo [8]. The numerical simulation of this equation is non-trivial, School of Computer Science & Center for Intelligent Machines McGill University, Montreal, PQ, Canada H3A 2A7 y School of Co... |

20 | Symmetry maps of freeform curve segments via wave propagation - Tek, Kimia - 2003 |