<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7654.12">
<TITLE>RE: [Insight-users] itkScalarChanAndVeseSparseLevelSetImageFilter</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->

<P><FONT SIZE=2>Hello,<BR>
thank you very much for the help. Unfortunately I work with medical data that I am not supposed to give.<BR>
Can you explain a little more point 2. You mean that I have to define my level set as inside positive and outside negative to have the opposite propagation direction?<BR>
I tried it but it did not work. I both cases, if I initialize the LS inside the boundary of my data, it shrinks. I use mu=0, nu=0, l1=1 and l2=1. Do I need negative l1?When I initialize almost the whole LS outside my data it seems to work in some cases.Generally, I get different results from the insight journal that I was getting correct.<BR>
<BR>
Thank you,<BR>
Eleni<BR>
<BR>
<BR>
-----Original Message-----<BR>
From: Kishore Mosaliganti [<A HREF="mailto:kishoreraom@gmail.com">mailto:kishoreraom@gmail.com</A>]<BR>
Sent: Fri 6/5/2009 1:02 PM<BR>
To: Sgouritsa, Eleni<BR>
Cc: Luis Ibanez; insight-users@itk.org<BR>
Subject: Re: [Insight-users] itkScalarChanAndVeseSparseLevelSetImageFilter<BR>
<BR>
Hi Eleni,<BR>
<BR>
Thank you for pointing out these problems. You are right that the Review<BR>
code is not very stable as yet. We are working on fixing them. During the<BR>
transfer from the IJ paper to the Review directory, a lot of modifications<BR>
were made. We might have inadvertently introduced new errors.<BR>
<BR>
Regarding your problems:<BR>
<BR>
2. The level-set function is defined by convention as inside is negative and<BR>
outside is positive. So the negative values become your foreground and<BR>
positive values become the background. This controls the propagation<BR>
direction.<BR>
<BR>
3. I am looking at the Halt function. Looks like the boundary case of 0<BR>
iterations is not coded right.<BR>
<BR>
4.&nbsp; Regarding the high RMS values, the level-set function is re-initialized<BR>
to a signed distance function after every few iterations. It so happens that<BR>
the default setting is 1 iteration. Hence, the RMS change is being<BR>
calculated after the reinitialization while it needs to be done before. We<BR>
are also working on correcting this. Thank you for pointing this out.<BR>
<BR>
<BR>
Is it possible for you to post a small sample of your data for us to work on<BR>
the errors.<BR>
<BR>
<BR>
Kishore<BR>
<BR>
On Tue, Jun 2, 2009 at 11:24 AM, Sgouritsa, Eleni<BR>
&lt;esgourit@central.uh.edu&gt;wrote:<BR>
<BR>
&gt;&nbsp; Hello,<BR>
&gt; I m sorry if my point was not clear. The problems that I am facing using<BR>
&gt; this class to segment 3D data are the following:<BR>
&gt;<BR>
&gt; 1. My data has negative origin (for example (-100, -100, -60). Unless I do<BR>
&gt; SetOrigin(0,0,0) (or I guess any origin &gt;=0) my program crashes. I have to<BR>
&gt; set my negative origin as zero for my program to run. Does the class assumes<BR>
&gt; positive origin of the input image?<BR>
&gt;<BR>
&gt; 2. I initialize my level set inside (so the evolution should be outwards).<BR>
&gt; In this case I can not make my level set expand (it can just shrink, or not<BR>
&gt; move). I think I should use something like SetReverseExpansionDirection to<BR>
&gt; control the propagation direction, but I can not find this option in the<BR>
&gt; class public members. (This is currently my main problem where I am stuck,<BR>
&gt; any help would be very useful!)<BR>
&gt;<BR>
&gt; 3. The maximum RMS Error serves (along with the number of iterations) as a<BR>
&gt; stopping criteria inside the Halt function (of<BR>
&gt; itkMultiphaseFiniteDifferenceImageFilter.txx file).<BR>
&gt; Lets say I set maximum RMS Error = 0.01.<BR>
&gt; The Halt function returns the following:<BR>
&gt; return ( (this-&gt;GetElapsedIterations() &gt;= this-&gt;m_NumberOfIterations) ||<BR>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( this-&gt;GetMaximumRMSError() &gt;= m_RMSChange ) );<BR>
&gt; While the Halt function returns false the evolution continues.<BR>
&gt; However, in iteration 0, I get m_RMSChange = 0, so the evolution always<BR>
&gt; stops doing 0 iterations (since the condition this-&gt;GetMaximumRMSError() &gt;=<BR>
&gt; m_RMSChange becomes true)<BR>
&gt; If I modify the Halt function to returm the following:<BR>
&gt; if(this-&gt;GetElapsedIterations()==0)<BR>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return false;<BR>
&gt; else<BR>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return ( (this-&gt;GetElapsedIterations() &gt;=<BR>
&gt; this-&gt;m_NumberOfIterations) || ( this-&gt;GetMaximumRMSError() &gt;= m_RMSChange )<BR>
&gt; );<BR>
&gt; then it works.<BR>
&gt; The other problem is that RMSChange in the next iterations (1,2,3...) has<BR>
&gt; very big values (like 5000) even near convergence. So the only way to stop<BR>
&gt; the evolution is through the iterations number.<BR>
&gt;<BR>
&gt; I am still investigating the class and I am not absolutely sure about the<BR>
&gt; points mentioned but I would like to share with you for more advice.<BR>
&gt;<BR>
&gt;<BR>
&gt; Thank you very much,<BR>
&gt; Lena<BR>
&gt;<BR>
&gt;<BR>
&gt;<BR>
&gt;<BR>
&gt; -----Original Message-----<BR>
&gt; From: Luis Ibanez [<A HREF="mailto:luis.ibanez@kitware.com">mailto:luis.ibanez@kitware.com</A>&lt;luis.ibanez@kitware.com&gt;<BR>
&gt; ]<BR>
&gt; Sent: Mon 6/1/2009 8:39 PM<BR>
&gt; To: Sgouritsa, Eleni<BR>
&gt; Cc: alex gouaillard; insight-users@itk.org<BR>
&gt; Subject: Re: [Insight-users] itkScalarChanAndVeseSparseLevelSetImageFilter<BR>
&gt;<BR>
&gt;<BR>
&gt; Hi Lena,<BR>
&gt;<BR>
&gt; Could you please be more specific about your problem report on<BR>
&gt; item (2) ?<BR>
&gt;<BR>
&gt; Is there a chance that you can report this in a form that<BR>
&gt; we could reproduce ?<BR>
&gt;<BR>
&gt;<BR>
&gt;&nbsp;&nbsp;&nbsp;&nbsp; Thanks<BR>
&gt;<BR>
&gt;<BR>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Luis<BR>
&gt;<BR>
&gt;<BR>
&gt; -----------------------------<BR>
&gt; Sgouritsa, Eleni wrote:<BR>
&gt; &gt; Hi,<BR>
&gt; &gt;<BR>
&gt; &gt; I have two more questions.<BR>
&gt; &gt; 1. In the insight journal code there was an option<BR>
&gt; &gt; SetReverseExpansionDirection to control if the expansion of the level<BR>
&gt; &gt; set will be inwards or outwards. How can I do it using the<BR>
&gt; &gt; itkScalarChanAndVeseSparseLevelSetImageFilter of 3.14 release?<BR>
&gt; &gt;<BR>
&gt; &gt; 2. I still have a problem when my 3D data have negative origin.<BR>
&gt; &gt;<BR>
&gt; &gt; Thanks and regards,<BR>
&gt; &gt; Lena<BR>
&gt; &gt;<BR>
&gt; &gt; -----Original Message-----<BR>
&gt; &gt; From: Sgouritsa, Eleni<BR>
&gt; &gt; Sent: Mon 6/1/2009 4:00 AM<BR>
&gt; &gt; To: alex gouaillard<BR>
&gt; &gt; Cc: insight-users@itk.org<BR>
&gt; &gt; Subject: RE: [Insight-users]<BR>
&gt; itkScalarChanAndVeseSparseLevelSetImageFilter<BR>
&gt; &gt;<BR>
&gt; &gt; Hello,<BR>
&gt; &gt; finally the second issue was my fault (I just had a problem while<BR>
&gt; &gt; including the header files of the review folder).<BR>
&gt; &gt; For the first, the inconsistency is just a matter of name. The code runs<BR>
&gt; &gt; correctly, just the name assignement in the txx is a little misleading.<BR>
&gt; &gt;<BR>
&gt; &gt; Thank you very much,<BR>
&gt; &gt; Lena<BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt; -----Original Message-----<BR>
&gt; &gt; From: alex gouaillard [<A HREF="mailto:alexandre_gouaillard@hms.harvard.edu">mailto:alexandre_gouaillard@hms.harvard.edu</A>&lt;alexandre_gouaillard@hms.harvard.edu&gt;<BR>
&gt; ]<BR>
&gt; &gt; Sent: Sun 5/31/2009 11:34 PM<BR>
&gt; &gt; To: Sgouritsa, Eleni<BR>
&gt; &gt; Cc: insight-users@itk.org<BR>
&gt; &gt; Subject: Re: [Insight-users]<BR>
&gt; itkScalarChanAndVeseSparseLevelSetImageFilter<BR>
&gt; &gt;<BR>
&gt; &gt; hi eleni,<BR>
&gt; &gt;<BR>
&gt; &gt; we'll take a look tomorrow.<BR>
&gt; &gt;<BR>
&gt; &gt; thanks for the feedback.<BR>
&gt; &gt;<BR>
&gt; &gt; alex.<BR>
&gt; &gt;<BR>
&gt; &gt; On May 31, 2009, at 11:52 PM, Sgouritsa, Eleni wrote:<BR>
&gt; &gt;<BR>
&gt; &gt;&nbsp; &gt; Hello,<BR>
&gt; &gt;&nbsp; &gt; I know that the review folder of the new release (3.14) is not<BR>
&gt; &gt;&nbsp; &gt; stable yet but I would just like to give two comments:<BR>
&gt; &gt;&nbsp; &gt;<BR>
&gt; &gt;&nbsp; &gt; 1. I found an inconsistency between the template arguments of the<BR>
&gt; &gt;&nbsp; &gt; itkScalarChanAndVeseSparseLevelSetImageFilter.h and the<BR>
&gt; &gt;&nbsp; &gt; itkScalarChanAndVeseSparseLevelSetImageFilter.txx.<BR>
&gt; &gt;&nbsp; &gt; Specifically, in the header file, the template is :<BR>
&gt; &gt;&nbsp; &gt; template &lt; class TInputImage, class TFeatureImage, class<BR>
&gt; &gt;&nbsp; &gt; TOutputImage, class TFunction,<BR>
&gt; &gt;&nbsp; &gt;&nbsp;&nbsp; class TSharedData, typename TIdCell = unsigned int &gt;<BR>
&gt; &gt;&nbsp; &gt; whereas in the txx it is:<BR>
&gt; &gt;&nbsp; &gt; template &lt; class TInput, class TFeature, class TFunction,<BR>
&gt; &gt;&nbsp; &gt; class TOutputPixel, class TSharedData, typename TIdCell &gt;<BR>
&gt; &gt;&nbsp; &gt; There is a change in the order of the 3rd and 4th argument.<BR>
&gt; &gt;&nbsp; &gt;<BR>
&gt; &gt;&nbsp; &gt; 2. I can not find this class<BR>
&gt; &gt;&nbsp; &gt; (itkScalarChanAndVeseSparseLevelSetImageFilter)&nbsp; incuded in the<BR>
&gt; &gt;&nbsp; &gt; CMakeLists of the review folder<BR>
&gt; &gt;&nbsp; &gt;<BR>
&gt; &gt;&nbsp; &gt; Regards,<BR>
&gt; &gt;&nbsp; &gt; Lena<BR>
&gt; &gt;&nbsp; &gt;<BR>
&gt; &gt;&nbsp; &gt; &lt;ATT00001.txt&gt;<BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt; ------------------------------------------------------------------------<BR>
&gt; &gt;<BR>
&gt; &gt; _____________________________________<BR>
&gt; &gt; Powered by www.kitware.com<BR>
&gt; &gt;<BR>
&gt; &gt; Visit other Kitware open-source projects at<BR>
&gt; &gt; <A HREF="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</A><BR>
&gt; &gt;<BR>
&gt; &gt; Please keep messages on-topic and check the ITK FAQ at:<BR>
&gt; <A HREF="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</A><BR>
&gt; &gt;<BR>
&gt; &gt; Follow this link to subscribe/unsubscribe:<BR>
&gt; &gt; <A HREF="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</A><BR>
&gt;<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">http://www.kitware.com/opensource/opensource.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">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">http://www.itk.org/mailman/listinfo/insight-users</A><BR>
&gt;<BR>
&gt;<BR>
<BR>
</FONT>
</P>

</BODY>
</HTML>