This in an implementation of ITK's FastMarchingExtensionImageFilter - From http://www.itk.org/Doxygen/html/classitk_1_1FastMarchingExtensionImageFilter.html - Fast marching can be used to extend auxiliary variables smoothly from the zero level set. Starting from an initial position on the front, this class simultaneously calculate the signed distance and extend a set of auxiliary values. Implemenation of this class is based on Chapter 11 of "Level Set Methods and Fast Marching Methods", J.A. Sethian, Cambridge Press, Second edition, 1999.

fastMarchingExtension(speedImage, labelImage, valueImage)

Arguments

speedImage

defines the cost or distance function

labelImage

defines the known (value 1) and unknown (value 2) regions

valueImage

these values are extended into the unknown regions

Value

antsImage

Author

Duda, JT

Examples


img <- makeImage(c(7, 7), 0)
img[3:5, 3:5] <- 1:9
mask <- getMask(img)
mask[mask == 0] <- 2
speed <- smoothImage(mask, 1)
extendedImg <- fastMarchingExtension(speed, mask, img)