<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear All: </DIV>
<DIV>&nbsp;</DIV>
<DIV>I tried to use ReinitializeLevelSetImageFilter.&nbsp;The output looks weird to me. </DIV>
<DIV>&nbsp;</DIV>
<DIV>The input to the filter is binary image binary image with float type (255&nbsp;as forground, 0 as background).</DIV>
<DIV>&nbsp;</DIV>
<DIV>I tried narroow band on&nbsp;and off. However, the output of this filter gives very high value (10^38) for voxles inside the binary mask.&nbsp;Almost&nbsp;all the outside mask voxels has value zero.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Below is my code, can somebody tell&nbsp;where I made mistake.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Baoyun</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>**************************code************************<BR>&nbsp;&nbsp; typedef typename itk::ReinitializeLevelSetImageFilter&lt; FeatureImageType &gt; DistanceFilterType; <BR>&nbsp;&nbsp; typename DistanceFilterType::Pointer DistanceFilter = DistanceFilterType::New(); <BR>&nbsp;&nbsp; DistanceFilter-&gt;SetInput( thresholder-&gt;GetOutput() ); <BR>&nbsp;&nbsp;// DistanceFilter-&gt;NarrowBandingOn(); <BR>&nbsp;&nbsp;// DistanceFilter-&gt;SetNarrowBandwidth( 10 ); <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; try<BR>&nbsp;&nbsp;&nbsp; {<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DistanceFilter-&gt;Update(); </DIV>
<DIV>&nbsp;&nbsp;&nbsp; }<BR>&nbsp;&nbsp;&nbsp; catch( itk::ExceptionObject &amp; excep )<BR>&nbsp;&nbsp;&nbsp; {<BR>&nbsp;&nbsp;&nbsp; std::cerr &lt;&lt; "Exception caught !" &lt;&lt; std::endl;<BR>&nbsp;&nbsp;&nbsp; std::cerr &lt;&lt; excep &lt;&lt; std::endl;<BR>&nbsp;&nbsp;&nbsp; }</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp; m_levelset2=DistanceFilter-&gt;GetOutput();</DIV>
<DIV>&nbsp;</DIV>
<DIV>****************code end*************************<BR></DIV></td></tr></table><br>