Hi Dan,<br>
<br>
Thank you for your suggestion. Unfortunately, this didn&#39;t help. I got the same binary output. Do you know any other possible reasons?<br><br><div class="gmail_quote">On Wed, Apr 22, 2009 at 2:01 AM, Dan Mueller <span dir="ltr">&lt;<a href="mailto:dan.muel@gmail.com">dan.muel@gmail.com</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;">Hi DING Feng,<br>
<br>
Try setting the output spacing and origin *before* updating the filter:<br>
<div class="im"><br>
&gt;    marcher-&gt;SetTrialPoints( points );<br>
&gt;    marcher-&gt;SetSpeedConstant( 1.0 );<br>
&gt;    marcher-&gt;SetStoppingValue( 100.0 );<br>
&gt;    marcher-&gt;SetOutputSize( size );<br>
<br>
</div>marcher-&gt;SetOutputSpacing( spacing );<br>
marcher-&gt;SetOutputOrigin( origin );<br>
<br>
This may be the source of your problem... FYI I use this filter<br>
regularly without issue.<br>
<br>
Hope this helps.<br>
<br>
Regards, Dan<br>
<br>
2009/4/21 Feng Ding &lt;<a href="mailto:carlding@gmail.com">carlding@gmail.com</a>&gt;:<br>
<div><div></div><div class="h5">&gt; Dear list,<br>
&gt;<br>
&gt; I&#39;m trying to use itkFastMarchingImageFilter by following the examples in<br>
&gt; the ITK source tree. Instead of getting a time arrival image which is<br>
&gt; similar to a distance map in this case, I only got a binary image with 0<br>
&gt; inside (a sphere) and 1.70141e+38 outside. I have tested my code on both<br>
&gt; 3.10 and 3.12 releases of ITK, and got the same results. Following is part<br>
&gt; of my code:<br>
&gt;<br>
&gt; // &gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; begin &gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;    index[0] = index[1] = index[2] = 100;<br>
&gt;    node.SetIndex( index );<br>
&gt;    node.SetValue( 0 );<br>
&gt;    points-&gt;InsertElement( 0, node );<br>
&gt;<br>
&gt;    marcher-&gt;SetTrialPoints( points );<br>
&gt;    marcher-&gt;SetSpeedConstant( 1.0 );<br>
&gt;    marcher-&gt;SetStoppingValue( 100.0 );<br>
&gt;    marcher-&gt;SetOutputSize( size );<br>
&gt;    marcher-&gt;DebugOn();<br>
&gt;    marcher-&gt;Update();<br>
&gt;    marcher-&gt;GetOutput()-&gt;SetOrigin( origin );<br>
&gt;    marcher-&gt;GetOutput()-&gt;SetSpacing( spacing );<br>
&gt;<br>
&gt;    volumeWriter-&gt;SetInput( marcher-&gt;GetOutput() );<br>
&gt;    volumeWriter-&gt;SetFileName( &quot;volume.mha&quot; );<br>
&gt;    volumeWriter-&gt;Write();<br>
&gt; // &lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;<br>
&gt;<br>
&gt; Anyone can give me some suggestion on that? Thanks!<br>
&gt;<br>
&gt; Best regards,<br>
&gt; DING Feng<br>
</div></div></blockquote></div><br>