Employ a topologically constrained label propagation tool in order to thin an existing segmentation into a skeleton. The function will perform better if the initial segmentation is well-composed.

skeletonize(x, laplacianThreshold = 0.25, propagationOption = 1)

Arguments

x

input binary image

laplacianThreshold

between zero and one, controls the speed image

propagationOption

either 0 (none), 1 (well-composed) or 2 (topological) constraint

Value

output binary image

Author

Avants BB

Examples


img <- antsImageRead(getANTsRData("r16")) %>% thresholdImage("Otsu", 3)
skel <- skeletonize(thresholdImage(img, 3, 3))