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