<br>Hi Serena,<br><br>When you say that the code <br><br>         &quot;doesn&#39;t write anything&quot;<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 &quot;diff&quot; 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">&lt;<a href="mailto:fabbri@u.washington.edu">fabbri@u.washington.edu</a>&gt;</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-&gt;GetNumberOfParameters();<br>
<br>
        typedef MetricType::ScalesType ScalesType;<br>
        ScalesType scales( numberOfParameters2 );<br>
        scales.Fill( 0.0001 );<br>
        metric-&gt;SetDerivativeStepLengthScales(scales);<br>
<br>
        MetricType::HistogramType::SizeType histogramSize;<br>
        histogramSize[0] = 25;<br>
        histogramSize[1] = 25;<br>
        metric-&gt;SetHistogramSize( histogramSize );<br>
<br>
      optimizer-&gt;MaximizeOn();<br>
<br>
the problem is the program seems blocked  as the observer associated to registration pipeline doesn&#39;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>