Skeletonization of volumetric angiograms for display

The display of three-dimensional angiograms can benefit from the knowledge of quantitative shape features such as tangent and curvature of the centerline of vessels. These can be obtained from a curve-like skeleton representation. If connectivity and topology are preserved, and if geometrical constraints such as smoothness and centeredness are satisfied, it is possible to estimate length, orientation, curvature, and torsion. It is also required that no part of the original object be left unrepresented. An efficient method for the identification of such shape components is developed. First, a suitable representation is obtained using a voxel coding approach to yield connected and labeled unit-thick paths. The desired features are estimated from a smoothed version of the skeleton produced by a moving average filter. The computational cost is linear, of the order of N(object), the total number of object voxels contained in the binary volumetric data. The method is also shown to be robust to boundary noise. Examples are discussed.