[Insight-users] MultiScaleHessianBasedMeasureImageFilter

[Luca Antiga]

Hi Juan Carlos,
  I ran the test for the filter in Testing/Code/Review on one of my  
images and it worked correctly, the scales image was filled with scale  
values. Find the code I ran attached.
Can you try it out on your image?

Luca


-------------- next part --------------
A non-text attachment was scrubbed...
Name: foo.cxx
Type: application/octet-stream
Size: 9177 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20100602/b9f2bf9f/attachment.obj>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: CMakeLists.txt
URL: <http://www.itk.org/pipermail/insight-users/attachments/20100602/b9f2bf9f/attachment.txt>
-------------- next part --------------



On May 31, 2010, at 8:18 PM, Juan Carlos Prieto wrote:

> Hello Luis,
>
> I'm trying to use the MultiScaleHessianBasedMeasureImageFilter.
> My problem is that I don't know how to get the "GetScalesOutput()"  
> to work, I have both methods to On "GenerateHessianOutputOn()" and
>         "GenerateScalesOutputOn()", but the ScalesOutput produces  
> nothing.
>
> This is my code
>
> _hessianvessel->SetAlpha1 (0.5);
>         _hessianvessel->SetAlpha2 (2);
>
>         _hessianmultiscale->SetInput(image);
>         _hessianmultiscale->SetHessianToMeasureFilter(_hessianvessel);
>         _hessianmultiscale->SetSigmaMinimum(0.5);
>         _hessianmultiscale->SetSigmaMaximum(4.0);
>         _hessianmultiscale->SetNumberOfSigmaSteps(5);
>         _hessianmultiscale->GenerateHessianOutputOn();
>         _hessianmultiscale->GenerateScalesOutputOn();
>         _hessianmultiscale->Update();
>
>         _outputimagefloat = ImageTypeFloat::New();
>         _outputimagefloat->SetRegions(_hessianmultiscale- 
> >GetOutput()->GetRequestedRegion());
>         _outputimagefloat->Allocate();
>
>         TypeConstRegionIterator constit( _hessianmultiscale- 
> >GetScalesOutput(),
>                          _hessianmultiscale->GetScalesOutput()- 
> >GetRequestedRegion() );
>
>         TypeRegionIterator outputit( _outputimagefloat,  
> _outputimagefloat->GetRequestedRegion() );
>
>         for ( constit.GoToBegin(); !constit.IsAtEnd(); ++constit);
>         {
>             cout<<constit.Get()<< " "; //this output produces "0" only
>         }
>         /*for ( constit.GoToBegin(), outputit.GoToBegin(); ! 
> constit.IsAtEnd(); ++constit, ++outputit)
>         {
>             outputit.Set(constit.Get()); //does not work
>         }*/
>
>
> I hope this is enough information
>
> Thank you for your time
>
> Juan Carlos Prieto
>
> _____________________________________
> 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