[Insight-users] ConnectedThresholdImageFilter, intialize output

[Luis Ibanez]

Hi Bryn,

Currently this filter doesn't take an initial segmentation as input.

The simplest way to extend it, is to follow your first suggestion:


       To pass through all the pixels of you initial
       segmentation and to add them as seed points.


Grafting the output will not do the trick because that will not
update the internal structures of this filter.


I would strongly encourage you to write this as a separate filter,
that takes a segmentation as input and produces as output a list
of indices (seed points) that can then be passed to any of the
region growing filters.

Such initialization filter will be relatively easy to write,
and will be a great contribution to the Insight Journal   :-)
http://www.insight-journal.org/


    Regards,


        Luis


-----------------
Bryn Lloyd wrote:
> Hi,
> 
> I am using the ConnectedThresholdImageFilter to segment vessels from MRA 
> images (reader -> hessian -> hessian2vesselness -> region growing).
> http://www.itk.org/Doxygen/html/classitk_1_1ConnectedThresholdImageFilter.html 
> 
> 
> Is it possible to setup the filter to use an initial given segmentation.
> I guess I could iterate over all pixels and add each object pixel as a 
> seed point...
> 
> But could I also use the function 
> ConnectedThresholdImageFilter::GraftOutput to initialize the filter with 
> an intial solution?
> 
> Any ideas?
> 
> 
> Thanks
> Bryn
> 
> 
> 
> 
> 
> _____________________________________
> Powered by www.kitware.com
> 
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
> 
> Please keep messages on-topic and check the ITK FAQ at: 
> http://www.itk.org/Wiki/ITK_FAQ
> 
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>