<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Hello, <br></span></div><div><br><span></span></div><div><span>Thank you for your reply.</span></div><div><br><span></span></div><div><span>I just used the "GeodesicActiveContourImageFilter.cxx" in ITK to perform the level-set segmentation. And I used a reader to read the initially segmented data, then give that data as the input to the "GeodesicActiveContourFilter", that is, in the original code, it was like: <br></span></div><div><br></div><div><span>*****************</span></div><div>&nbsp; typedef&nbsp;&nbsp; float&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; InternalPixelType;<br>&nbsp; const&nbsp;&nbsp;&nbsp;&nbsp; unsigned int&nbsp;&nbsp;&nbsp; Dimension = 3;<br>&nbsp; typedef itk::Image&lt; InternalPixelType, Dimension &gt;&nbsp; InternalImageType;<br></div><div><br></div><div>typedef&nbsp;
 itk::FastMarchingImageFilter&lt; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; InternalImageType, <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; InternalImageType &gt;&nbsp;&nbsp;&nbsp; FastMarchingFilterType;</div><div><br></div><div>FastMarchingFilterType::Pointer&nbsp; fastMarching = FastMarchingFilterType::New();</div><div>...<br><span></span></div><div><span>geodesicActiveContour-&gt;SetInput(&nbsp; fastMarching-&gt;GetOutput() );</span></div><div><span>*****************</span></div><div><br><span></span></div><div><span>But, I changed it as: <br></span></div><div><span><br></span></div><div><span>*****************<br></span></div><div><span>typedef&nbsp; itk::ImageFileReader&lt;
 InternalImageType &gt; ReaderType;</span></div><div><span>ReaderType::Pointer reader2 = ReaderType::New();</span></div><div><span>reader2-&gt;SetFileName( argv[10] );<br></span></div><div><span>geodesicActiveContour-&gt;SetInput(&nbsp; reader2-&gt;GetOutput() );</span></div><div><span>*****************</span></div><div>The code compiles, but the result is nothing but exactly the same as the initially segmented data, which was read and given to the "GeodesicActiveContourFilter", even if I changed several kinds of parameters. <br></div><div><br></div><div>Am I doing something wrong here??</div><div><br></div><div>Looking forward to receive your suggestions. Thanks in advance. <br></div><div><br></div><div>Regards,</div><div>Abayiz<br></div><div><br></div><div><br></div><div><br><span></span></div><div><span><br></span></div><div><span><br></span></div><div><span><br></span></div><div><br></div><div style="font-family: times new roman,new york,times,serif;
 font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Arial" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Juan Cardelino &lt;juan.cardelino@gmail.com&gt;<br><b><span style="font-weight: bold;">To:</span></b> Abayiz &lt;abayiz@yahoo.com&gt;<br><b><span style="font-weight: bold;">Cc:</span></b> itk &lt;insight-users@itk.org&gt;<br><b><span style="font-weight: bold;">Sent:</span></b> Monday, August 15, 2011 2:31 AM<br><b><span style="font-weight: bold;">Subject:</span></b> Re: [Insight-users] Use Initial Segmentation as Input Level-Set<br></font><br>On Sun, Aug 14, 2011 at 9:29 AM, Abayiz &lt;<a ymailto="mailto:abayiz@yahoo.com" href="mailto:abayiz@yahoo.com">abayiz@yahoo.com</a>&gt; wrote:<br>&gt; Hello,<br>&gt; I want to perform level-set segmentation using a set of already achieved<br>&gt; initially segmented 3D data. What I tried was giving the initial
 segmented<br>&gt; data as the output of the level-set filter&nbsp; (for instance instead of giving<br>&gt; FastMarchingFilter-&gt;GetOutput(), I gave reader-&gt;GetOutput, after reading the<br>&gt; data) . But it didn't work. How can I do this??<br><br>Well I do exactly that and works, it is just a matter of details. If<br>you provide more information about how exactly do you implement it<br>maybe we can help you.<br>In addition, it would help if you also tell us what exactly "didn't<br>work means". It doesn't compile? The program crashes? It gives<br>unexpected output?<br>Please try to be more precise, we are not seeing your screen so we<br>can't possibly guess what is going on.<br><br><br>&gt; Sincerely thank you in advance!<br>&gt; Regards,<br>&gt; Abayiz<br>&gt;<br>&gt; _____________________________________<br>&gt; Powered by www.kitware.com<br>&gt;<br>&gt; Visit other Kitware open-source projects at<br>&gt; <a
 href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>&gt;<br>&gt; Kitware offers ITK Training Courses, for more information visit:<br>&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>&gt;<br>&gt; Please keep messages on-topic and check the ITK FAQ at:<br>&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>&gt;<br>&gt; Follow this link to subscribe/unsubscribe:<br>&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>&gt;<br>&gt;<br><br><br></div></div></div></body></html>