[Insight-users] MultiScaleHessianBasedMeasureImageFilter

[Juan Carlos Prieto]

Hi Luca,

Thank you for the response, I ran some other tests on the filter like
printing values directly on the template
(itkMultiScaleHessianBasedMeasureImageFilter.txx) and I think is working
fine because the filter managed to show them. But my problem is that I can't
get an iterator to go through the values after the update of the filter,
like the example I sent you.

_hessianmultiscale->Update();
_hessianmultiscale->GetOutput(); //this is working fine

TypeConstRegionIterator constit( _hessianmultiscale->GetScalesOutput(),

>
> _hessianmultiscale->GetScalesOutput()->GetRequestedRegion() );


 for ( constit.GoToBegin(); !constit.IsAtEnd(); ++constit);
       {
           cout<<constit.Get()<< " "; //this output produces "0" only
       }

I'm going to try saving the image and then reading it to see what happens.
I'll let you know.

Juan Carlos Prieto

2010/6/2 Luca Antiga <luca.antiga at gmail.com>

> 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
>
>
>
>
>
>
> 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
>>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20100602/470b2acd/attachment.htm>