<div>&nbsp;Thanks, Luca,<br><br>&nbsp;&nbsp;&nbsp; I have found the problem. The LevelSetFunction should be initialized first,otherwise the radius will set 0, then such an error happens. Also, it seems that the SpeedImage should be calculated before sending to the parallelfilter. <br></div><div>&nbsp;<br>&nbsp;&nbsp;&nbsp; The code as follows:<br>&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; &nbsp; m_LevelSetFunction &nbsp; &nbsp;= &nbsp; &nbsp; ThresholdFunctionType::New();<br>
 &nbsp; &nbsp;&nbsp;&nbsp; m_LevelSetFunction &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetUpperThreshold(150);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; m_LevelSetFunction &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetLowerThreshold(100);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; m_LevelSetFunction &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetFeatureImage(SourceImage);<br>
&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;  <br>&nbsp;&nbsp;&nbsp; &nbsp; <font color="#800000">ThresholdFunctionType::RadiusType &nbsp;&nbsp; r;<br>&nbsp;&nbsp;&nbsp; &nbsp; r[0] = r[1] = r[2] =1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; m_LevelSetFunction-&gt;Initialize(r);<br>&nbsp;&nbsp;&nbsp; &nbsp; <br>&nbsp;&nbsp;&nbsp; &nbsp;</font><font color="#800000"> m_LevelSetFunction-&gt;AllocateSpeedImage();<br>&nbsp;&nbsp;&nbsp; &nbsp;</font><font color="#800000"><font color="#800000"> m_LevelSetFunction-&gt;CalculateSpeedImage();</font> </font>&nbsp;&nbsp; &nbsp;&nbsp;  <br><br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;= &nbsp; &nbsp; ParallelSparseFieldFilterType::New();<br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetInput(InitImage);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetNumberOfLayers(3);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetIsoSurfaceValue(0.0);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;-&gt; &nbsp; &nbsp;SetDifferenceFunction(m_LevelSetFunction);<br>
 &nbsp; &nbsp;&nbsp;&nbsp; parallelfilter &nbsp; &nbsp;-&gt; &nbsp; &nbsp;Update();<br>&nbsp;&nbsp;&nbsp; <br></div><div>&nbsp;<br>&nbsp;&nbsp;&nbsp; Appreciate and Thanks all.<br><br>Kun <br></div><blockquote id="isReplyContent" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;"><div dir="ltr"><div class="gmail_quote"><br>From: <b class="gmail_sendername">Luca Antiga</b> <span dir="ltr">&lt;<a href="mailto:luca.antiga@gmail.com">luca.antiga@gmail.com</a>&gt;</span><br>
Date: Tue, Jul 29, 2008 at 10:59 AM<br>Subject: Re: [Insight-users] problem with ParallelSparseFieldLevelSetImageFilter<br>To: ITK Users &lt;<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>&gt;<br><br><br>
Dear Kun,<br>
&nbsp;I'll try to reproduce the problem and get back to you.<br>
Keep in touch<br><font color="#888888">
<br>
Luca</font><div><div></div><div class="Wj3C7c"><br>
<br>
<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
 &nbsp; <br>
<br></blockquote><br>
</div></div></div>
</div>
</blockquote>