A Computational Method for Segmenting Topological Point-Sets and Application to Image Analysis

• Authors: Stiliyan N. Kalitzin;Joes Staal;Bart M. ter Haar Romeny;Max A. Viergever
• Venue: IEEE Transactions on Pattern Analysis and Machine Intelligence
• Year: 2001

We propose a new computational method for segmenting topological subdimensional point-sets in scalar images of arbitrary spatial dimensions. The technique is based on calculating the homotopy class defined by the gradient vector in a subdimensional neighborhood around every image point. This neighborhood is defined as the linear envelope spawned over a given subdimensional vector frame. In the simplest case where the rank of this frame is maximal, we obtain a technique for localizing the critical points, i.e., extrema and saddle points. We consider, in particular, the important case of frames formed by an arbitrary number of the first largest by absolute value principal directions of the Hessian. The method then segments positive and and negative ridges as well as other types of critical surfaces of different dimensionalities. The signs of the eigenvalues associated to the principal directions provide a natural labeling of the critical subsets. The result, in general, is a constructive definition of a hierarchy of point-sets of different dimensionalities linked by inclusion relations. Because of its explicit computational nature, the method gives a fast way to segment height ridges or edges in different applications. The defined topological point-sets are connected manifolds and, therefore, our method provides a tool for geometrical grouping using only local measurements. We have demonstrated the grouping properties of our construction by presenting two different cases where an extra image coordinate is introduced. In one of the examples, we considered the image analysis in the framework of the linear scale-space concept, where the topological properties are gradually simplified through the scale parameter. This scale parameter can be taken as an additional coordinate. In the second example, a local orientation parameter was used for grouping and segmenting elongated structures.