On the computation of the 〈3, 4, 5〉 curve skeleton of 3D objects

An algorithm to compute the curve skeleton of 3D objects in voxel images is presented. The skeleton is stable under isometric transformations of the object, since the algorithm is based on the use of the 〈3, 4, 5〉 weighted distance transform, which is a good approximation of the Euclidean distance transform. The 〈3, 4, 5〉 weighted distance transform is used both to identify suitable anchor points, and to efficiently check object voxels according to their distance to the background. The curve skeleton is symmetrically placed within the object, is topologically equivalent to the object, is unit-wide and provides a satisfactory representation of the object. Though the size of the object reconstructed from the curve skeleton via the reverse distance transformation is not as thick as that of the input, shape information is mostly retained by the skeleton, since all regions perceived as significant in the input can still be found in the reconstructed object.

Research highlights
► Selection as anchor points of a suitable subset of centers of maximal balls.
► Voxel classification in a set consisting of surfaces and curves.
► Final thinning preserving the structure of surfaces and curves.