New removal operators for surface skeletonization

New 3×3×3 operators are introduced to compute the surface skeleton of a 3D object by either sequential or parallel voxel removal. We show that the operators can be employed without creating disconnections, cavities, tunnels and vanishing of object components. A final thinning process, aimed at obtaining a unit-thick surface skeleton, is also described.