Dear list,<br><br>I'm trying to use itkFastMarchingImageFilter by following the examples in the ITK source tree. Instead of getting a time arrival image which is similar to a distance map in this case, I only got a binary image with 0 inside (a sphere) and 1.70141e+38 outside. I have tested my code on both 3.10 and 3.12 releases of ITK, and got the same results. Following is part of my code:<br>
<br>// >>>>>>>>>> begin >>>>>>>>>><br> index[0] = index[1] = index[2] = 100;<br> node.SetIndex( index );<br> node.SetValue( 0 );<br> points->InsertElement( 0, node );<br>
<br> marcher->SetTrialPoints( points );<br> marcher->SetSpeedConstant( 1.0 );<br> marcher->SetStoppingValue( 100.0 );<br> marcher->SetOutputSize( size );<br> marcher->DebugOn();<br> marcher->Update();<br>
marcher->GetOutput()->SetOrigin( origin );<br> marcher->GetOutput()->SetSpacing( spacing );<br><br> volumeWriter->SetInput( marcher->GetOutput() );<br> volumeWriter->SetFileName( "volume.mha" );<br>
volumeWriter->Write();<br>// <<<<<<<<<<<<<<<<<<<<<<<<<<br><br>Anyone can give me some suggestion on that? Thanks!<br><br>Best regards,<br>DING Feng<br>