Hi All,<div><br></div><div>I am using itkFastMarchingImageFilter to create a distance map to a set of points, this seems like a very basic usage of ITK</div><div><br></div><div>//GetNodes</div><div><div>  seeds-&gt;Initialize();</div>
<div>  for (int i=0;i&lt;seedPoints.size();i++)</div><div>  {</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>    ImageType::IndexType seedIndex;</div><div>        input-&gt;TransformPhysicalPointToIndex( seedPoints[i],seedIndex);</div>
<div><span class="Apple-tab-span" style="white-space:pre">                </span>NodeType node;</div><div><span class="Apple-tab-span" style="white-space:pre">                </span>const double seedValue = 0.0;</div><div><span class="Apple-tab-span" style="white-space:pre">                </span>node.SetValue( seedValue );</div>
<div><span class="Apple-tab-span" style="white-space:pre">                </span>node.SetIndex( seedIndex );</div><div><span class="Apple-tab-span" style="white-space:pre">                </span>seeds-&gt;InsertElement( i, node );</div><div>  }</div>
</div><div>//execute fast marching</div><div><div>  FastMarchingFilterType::Pointer fastMarching = FastMarchingFilterType::New();</div><div>  fastMarching-&gt;SetInput( input);</div><div>  fastMarching-&gt;SetTrialPoints(seeds);</div>
<div>  fastMarching-&gt;SetStoppingValue(2*radius);</div><div>  fastMarching-&gt;SetSpeedConstant(1);</div><div>  fastMarching-&gt;SetNormalizationFactor(1);</div><div>  fastMarching-&gt;SetOutputSize( input-&gt;GetBufferedRegion().GetSize() );</div>
<div>  fastMarching-&gt;Update();</div></div><div><br></div><div>This works,but it is extremely slow......I think something went wrong but I do not know where.</div><div><br></div><div>May I ask did anyne of you also find out this problem?</div>
<div><br></div><div>Thanks a lot!</div><div><br></div><div>Best,</div><div><br></div><div>Hui</div>