[Insight-users] MultiScaleHessianBasedMeasureImageFilter

[Juan Carlos Prieto]

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