<br>Hi Serena,<br><br>When you say that the code <br><br> "doesn't write anything"<br><br>do you mean that it:<br><br>A) hangs ?<br>B) seg faults ?<br>C) throws an exception ?<br>D) terminate without errors, but without running any iteration ?<br>
<br><br>Also,<br>the value of<br><br> scales.Fill( 0.0001 );<br><br>is a bit extreme....<br><br>Why are you using such a small value ?<br><br><br> Please let us know,<br><br><br> Thanks<br><br><br> Luis<br>
<br><br><br>PS. It will be useful is you post a "diff" between<br>the original file and yours.<br><br>-----------------------------------------------------<br><div class="gmail_quote">On Mon, Jul 13, 2009 at 1:01 PM, Serena Fabbri <span dir="ltr"><<a href="mailto:fabbri@u.washington.edu">fabbri@u.washington.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
Hi All,<br>
<br>
I am registering MRI and CT image and I use Normalized Mutual Information.<br>
I am using DeformableRegistration8 (BSpline Tranformtion) and I have replaced in the code only the line about the metric.<br>
So I have written:<br>
<br>
<br>
const unsigned int numberOfParameters2 = transform->GetNumberOfParameters();<br>
<br>
typedef MetricType::ScalesType ScalesType;<br>
ScalesType scales( numberOfParameters2 );<br>
scales.Fill( 0.0001 );<br>
metric->SetDerivativeStepLengthScales(scales);<br>
<br>
MetricType::HistogramType::SizeType histogramSize;<br>
histogramSize[0] = 25;<br>
histogramSize[1] = 25;<br>
metric->SetHistogramSize( histogramSize );<br>
<br>
optimizer->MaximizeOn();<br>
<br>
the problem is the program seems blocked as the observer associated to registration pipeline doesn't write anything.<br>
Did i miss to write some instructions? any idea?<br>
I appreciate any suggestion.<br>
<br>
Thank you.<br>
<br>
Serena.<br>
<br>
<br>
<br>
<br>
<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>
Please keep messages on-topic and check the ITK FAQ at: <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></div><br>