Hi Luca,<br><br>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.<br>
<br>_hessianmultiscale->Update();<br>_hessianmultiscale->GetOutput(); //this is working fine <br><br>TypeConstRegionIterator constit( _hessianmultiscale->GetScalesOutput(),<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
_hessianmultiscale->GetScalesOutput()->GetRequestedRegion() );</blockquote><br> for ( constit.GoToBegin(); !constit.IsAtEnd(); ++constit);<br>
{<br>
cout<<constit.Get()<< " "; //this output produces "0" only<br>
}<br><br>I'm going to try saving the image and then reading it to see what happens.<br>I'll let you know.<br><br>Juan Carlos Prieto<br><br><div class="gmail_quote">2010/6/2 Luca Antiga <span dir="ltr"><<a href="mailto:luca.antiga@gmail.com">luca.antiga@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Juan Carlos,<br>
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.<br>
Can you try it out on your image?<br><font color="#888888">
<br>
Luca<br>
<br>
<br>
</font><br><br>
<br>
<br>
On May 31, 2010, at 8:18 PM, Juan Carlos Prieto wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello Luis,<br>
<br>
I'm trying to use the MultiScaleHessianBasedMeasureImageFilter.<br>
My problem is that I don't know how to get the "GetScalesOutput()" to work, I have both methods to On "GenerateHessianOutputOn()" and<br>
"GenerateScalesOutputOn()", but the ScalesOutput produces nothing.<br>
<br>
This is my code<br>
<br>
_hessianvessel->SetAlpha1 (0.5);<br>
_hessianvessel->SetAlpha2 (2);<br>
<br>
_hessianmultiscale->SetInput(image);<br>
_hessianmultiscale->SetHessianToMeasureFilter(_hessianvessel);<br>
_hessianmultiscale->SetSigmaMinimum(0.5);<br>
_hessianmultiscale->SetSigmaMaximum(4.0);<br>
_hessianmultiscale->SetNumberOfSigmaSteps(5);<br>
_hessianmultiscale->GenerateHessianOutputOn();<br>
_hessianmultiscale->GenerateScalesOutputOn();<br>
_hessianmultiscale->Update();<br>
<br>
_outputimagefloat = ImageTypeFloat::New();<br>
_outputimagefloat->SetRegions(_hessianmultiscale->GetOutput()->GetRequestedRegion());<br>
_outputimagefloat->Allocate();<br>
<br>
TypeConstRegionIterator constit( _hessianmultiscale->GetScalesOutput(),<br>
_hessianmultiscale->GetScalesOutput()->GetRequestedRegion() );<br>
<br>
TypeRegionIterator outputit( _outputimagefloat, _outputimagefloat->GetRequestedRegion() );<br>
<br>
for ( constit.GoToBegin(); !constit.IsAtEnd(); ++constit);<br>
{<br>
cout<<constit.Get()<< " "; //this output produces "0" 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>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</blockquote>
<br>
<br></blockquote></div><br>