Module: Thickness map

Description:

Assigns as greylevel to each pixel in the object the radius of the largest sphere contained in the object that also contains the pixel. Other pixels get a value of 0. This module implements the algorithm proposed by Coeurjolly and Montanvert ¹.
Example.

Connections:

Image

[required]
The input image comes in through this port. For now, the only input type that has been implemented is HxUniformScalarField3.

Ports:

Threshold

Threshold used for defining the pixels of the object.

Comparison

Comparison function used for defining the pixels of the object.

Compute

Push this button to compute.

Commands:

Additional options can be accessed when typing in the console Thickness_Map COMMAND_NAME. Typing the command again usually reverts back to original settings.

`useRDMA`

Uses the reduced discrete medial axis to compute the thickness map¹, the default. If not set, then the map is computed from the upper enveloppe of elliptic paraboloids (sub-optimal).

`verbose`

Displays timing information after the computation. Retype to hide info.

`create`

Runs the computation. Returns the name of the output, so it can be used in a script, such as set RESULT [Thickness_Map create].

References:

¹ Coeurjolly, D. & Montanvert, A. Optimal separable algorithms to compute the reverse euclidean distance transformation and discrete medial axis in arbitrary dimension, Ieee Transactions On Pattern Analysis and Machine Intelligence, 2007, 29, 437-448.