[Insight-users] Use Initial Segmentation as Input Level-Set

[Abayiz]

Hello, 


Thank you for your reply.

I just used the "GeodesicActiveContourImageFilter.cxx" in ITK to perform the level-set segmentation. And I used a reader to read the initially segmented data, then give that data as the input to the "GeodesicActiveContourFilter", that is, in the original code, it was like: 


*****************
  typedef   float           InternalPixelType;
  const     unsigned int    Dimension = 3;
  typedef itk::Image< InternalPixelType, Dimension >  InternalImageType;


typedef  itk::FastMarchingImageFilter< 
                              InternalImageType, 
                              InternalImageType >    FastMarchingFilterType;

FastMarchingFilterType::Pointer  fastMarching = FastMarchingFilterType::New();
...

geodesicActiveContour->SetInput(  fastMarching->GetOutput() );
*****************

But, I changed it as: 


*****************

typedef  itk::ImageFileReader< InternalImageType > ReaderType;
ReaderType::Pointer reader2 = ReaderType::New();
reader2->SetFileName( argv[10] );

geodesicActiveContour->SetInput(  reader2->GetOutput() );
*****************
The code compiles, but the result is nothing but exactly the same as the initially segmented data, which was read and given to the "GeodesicActiveContourFilter", even if I changed several kinds of parameters. 


Am I doing something wrong here??

Looking forward to receive your suggestions. Thanks in advance. 


Regards,
Abayiz









________________________________
From: Juan Cardelino <juan.cardelino at gmail.com>
To: Abayiz <abayiz at yahoo.com>
Cc: itk <insight-users at itk.org>
Sent: Monday, August 15, 2011 2:31 AM
Subject: Re: [Insight-users] Use Initial Segmentation as Input Level-Set

On Sun, Aug 14, 2011 at 9:29 AM, Abayiz <abayiz at yahoo.com> wrote:
> Hello,
> I want to perform level-set segmentation using a set of already achieved
> initially segmented 3D data. What I tried was giving the initial segmented
> data as the output of the level-set filter  (for instance instead of giving
> FastMarchingFilter->GetOutput(), I gave reader->GetOutput, after reading the
> data) . But it didn't work. How can I do this??

Well I do exactly that and works, it is just a matter of details. If
you provide more information about how exactly do you implement it
maybe we can help you.
In addition, it would help if you also tell us what exactly "didn't
work means". It doesn't compile? The program crashes? It gives
unexpected output?
Please try to be more precise, we are not seeing your screen so we
can't possibly guess what is going on.


> Sincerely thank you in advance!
> Regards,
> Abayiz
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20110815/5d3a6654/attachment.htm>