Hello
I hope I have explained what I mean clearly.
I do not understand why the Extract Skeleton module cannot extract the skeleton.
Is this node (LabelMapVolume node) too complicated?
As can be seen in the following image, some parts of the vessels are disconnected and not connected to each other. Does this prevent this module from extracting the skeleton of vascular tree?
To do this, I did the following steps.
The following steps were performed to extract the skeleton of pulmonary vascular system, respectively.
1- Segment the lung region from the chest CT image
2- Segment the pulmonary vessels by applying the Frangi filter (argument) in the command of vmtkimagevesselenhancement
3- Apply Threshold Effect in the module of Segment Editor
4- Removal of the outer edge of the lung region and the common border between the lung and the heart by the Sobel filter in the module of Simple Fiters
5- Split segmented vessel tree by Island Effect (Split islands to segments) in the module of Segment Editor
6- Convert the result of the previous step to Model and save them in vtp format
7- Merge vessel component models by the MergeModels (CLI)
The following image is shown in the result of above steps.
No difference, they both refer to the same thing: medial line or surface.
1D enforces extraction of a line (medial line), 2D extracts a surface (medial surface).
Computation time may be very long for 2D mode.
Number of points specifies the maximum. If your image has lower resolution and the centerline is short then the line may be represented with less points than what you have requested. If you need more points then you can get arbitrarily resampling of the curve using markups module.
Result of “Extract skeleton” module is a labelmap volume. You can use volume rendering or you can import it into Segmentation to visualize in 3D views.
If the difference between spacing along different axis is more than a few ten percent then I would recommend to use Crop volume to crop and resample the volume using isotropic spacing.
This is unnecessarily complicated and slow way of removing small islands. You can achieve the same by using “Remove small islands” method in Islands effect.
Probably those branches are too small. If you share the scene that contains the segmentation that you created using Thresholding effect (but before Island effect applied or other further processing is done) then I may have a look.
Thank you so much for your detailed explanation. Excuse me, but I do not understand how to share with you the scene that contains the segmentation. Do you mean a scene with .mrml extension (MRML scene) or .mrb extension (Medical Reality Bundle)? Where should I upload this file?
