<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Hui,<div><br></div><div>The principle of Chan and Vese formulation is to assume there are two parts in your image: background and foreground and each of them are represented by a constant mean value.</div><div>In the implementation, we did follow the same principle, e.g. the mean is a scalar if you are working with scalar images, a vector if you are working with vector images, etc.</div><div><br></div><div>So, now you want to create a new term / formulation (there might be a specific name for that, but I am not aware of it).</div><div>Once you provide the mean image, you need to replace m_Mean by the pixel value for the given input location iP. In your special case, most likely you would only need to change m_Mean by m_MeanImage->GetPixel( iP )</div><div><br></div><div>Please, let us know if it works.</div><div><br></div><div>Best,</div><div>Arnaud</div><div><br></div><div> </div><div><div><div>On Oct 2, 2011, at 1:55 PM, Hui Tang wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Arnaud,</div> <div> </div> <div>That works!!Thanks!</div> <div>Now may I ask if my mean intensity is a constant over time but varies in the image, so I will have a mean intensity image, how can I add that to interal term? I know how to get and set image to get the mean intensity image, but in update(), it seems that this->mean can only be a scalar.......</div> <div> </div> <div>Do you have any idea on solve this?</div> <div> </div> <div>thanks,</div> <div> </div> <div>Best,</div> <div> </div> <div>Hui</div> <div><br><br> </div> <div class="gmail_quote">On Fri, Sep 30, 2011 at 3:14 PM, Arnaud Gelas <span dir="ltr"><<a href="mailto:arnaud_gelas@hms.harvard.edu">arnaud_gelas@hms.harvard.edu</a>></span> wrote:<br> <blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><u></u> <div text="#000000" bgcolor="#ffffff">Hi Hui,<br><br>Sorry, there was a typo in my previous email :-/...<br><br>Here you go:<br><br> $ git checkout -t origin/v4.0b01<br><br>Here are the *.h you are looking for<br><br> $ find . | grep itkBinaryImageToWhitakerSparseLevelSetAdaptor.h<br>./Modules/Segmentation/LevelSetsv4/include/itkBinaryImageToWhitakerSparseLevelSetAdaptor.hxx<br>./Modules/Segmentation/LevelSetsv4/include/itkBinaryImageToWhitakerSparseLevelSetAdaptor.h<br> <br>To be able to use the LevelSetsv4, you should either build ITK with all modules, or explicitly turn on the LevelSetsv4 one (which is part of the Segmentation group).<br><br>Please keep the conversation on the mailing list, other users may have the same problem.<br> <br>Best,<br><font color="#888888">Arnaud</font> <div> <div></div> <div class="h5"><br><br>On 09/30/2011 09:06 AM, Hui Tang wrote: <blockquote type="cite"> <div>Dear Arnaud, </div> <div> </div> <div>Thanks for your help. I still can not find itkBinaryImageToWhitakerSparseLevelSetAdaptor.h in v4.0b1. </div> <div>git shows the version number, so it should be right.....</div> <div>Best,</div> <div> </div> <div>Hui</div> <div class="gmail_quote">On Thu, Sep 29, 2011 at 7:05 PM, Hui Tang <span dir="ltr"><<a href="mailto:tanghui.seu@gmail.com" target="_blank">tanghui.seu@gmail.com</a>></span> wrote:<br> <blockquote style="BORDER-LEFT: rgb(204,204,204) 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote"><p>Hi arnaud thanks a loooooot for this! I will try, I will let you know how it works:) best hui </p> <div class="gmail_quote">Op 29 sep. 2011 18:55 schreef "Gelas, Arnaud Joel Florent" <<a href="mailto:Arnaud_Gelas@hms.harvard.edu" target="_blank">Arnaud_Gelas@hms.harvard.edu</a>> het volgende: <div> <div><br type="attribution">> Hi Hui Tang,<br>> <br>> I am really happy to see you are interested and are testing the new level set infrastructure!<br>> <br>> If you get this error, I can only assume that you are now using the git master version of ITK, this file has been remove recently from the repository...<br> > I would strongly recommend to stick to ITK tags, since API keeps on changing all around the place in the toolkit.<br>> <br>> Here is how you'd get the last tagged version:<br>> <br>> $ git fetch --all<br> > $ git checkout -t origin/v4.0b1<br>> <br>> With this version, you should be able to find itkBinaryImageToWhitakerSparseLevelSetAdaptor.h<br>> <br>> If you are more adventurous, and want to be part of the development, you can use master, but expect API / files / class name to change.<br> > You can also create a gerrit account to review patches (instructions are on the wiki). Note that all contributions / inputs on this effort are welcome!!<br>> <br>> Now regarding your question, this is really easy to make it happened in the v4 infrastructure. You have two possibilities:<br> > <br>> 1- grab the code from iternal or external term (depending on your case); copy the code; you modify and remove the part in which the mean gets updated.<br>> and you add itkSetMacro to set the mean at the beginning.<br> > <br>> 2- you inherit, might need to change some method to virtual (if needed, let me know that'd be easy to fix).<br>> <br>> If you need more details on how to proceed for 1 or 2, just let me know, we can definitively give you some more details on how to proceed.<br> > <br>> Thanks,<br>> <br>> Best,<br>> Arnaud<br>> <br>> ________________________________________<br>> From: Hui Tang [<a href="mailto:tanghui.seu@gmail.com" target="_blank">tanghui.seu@gmail.com</a>]<br> > Sent: Thursday, September 29, 2011 12:20 PM<br>> To: Gelas, Arnaud Joel Florent<br>> Subject: Questions about ITK4<br>> <br>> Dear Arnaud,<br>> <br>> How are you?<br>> I am using ITK4, but this headfile "itkBinaryImageToWhitakerSparseLevelSetAdaptor.h" can not be found in ITK4. Did I do anything wrong?<br> > And maybe you could help me with another question:) in chanvase model, if my I would like to have a fixed average intensity ("c1") for the foreground how can I implement that using ITK4?<br>> <br>> Thanks so much.<br> > <br>> Looking forward to your reply.<br>> <br>> Best,<br>> <br>> Hui Tang<br></div></div></div></blockquote></div><br></blockquote><br></div></div></div></blockquote></div><br></blockquote></div><br></div></body></html>