<div>Dear Monica<br><br>&nbsp;Thanks very much.<br></div>&nbsp;This help me a lot.<br><br>Thanks.<br>Kun&nbsp;  <br><br><br>From: <b class="gmail_sendername">Monica Hernandez</b> <span dir="ltr">&lt;<a href="mailto:mhg@unizar.es">mhg@unizar.es</a>&gt;</span><br><blockquote id="isReplyContent" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;"><div dir="ltr"><div class="gmail_quote">
Date: Thu, Aug 28, 2008 at 12:25 PM<br>Subject: Re: [Insight-users] Insight-users Digest, Vol 52, Issue 49 needs help for "Gaussian Curvature" (Kun)<br>To: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
<br><br>Dear Khun,<br>
<br>
the formula for computing the Gaussian curvature of a surface given its<br>
implicit representation<br>
\phi is given by (I use LaTex notation)<br>
<br>
\begin{eqnarray}<br>
K_\textnormal{\scriptsize Gauss} &amp; = &amp; \frac{\phi_x^2(\phi_{yy}\phi_{zz}<br>
- \phi_{yz}^2) +<br>
\phi_y^2(\phi_{xx}\phi_{zz} - \phi_{xz}^2) +<br>
\phi_z^2(\phi_{xx}\phi_{yy} - \phi_{xy}^2)}<br>
{(\phi_x^2 + \phi_y^2 + \phi_z^2)^2} &nbsp;\\<br>
&amp; + &amp; 2 \phantom{c} \frac{\phi_x\phi_y(\phi_{xz} \phi_{yz} - \phi_{xy}<br>
\phi_{zz}) +<br>
\phi_y\phi_z(\phi_{xy} \phi_{xz} - \phi_{yz} \phi_{xx}) +<br>
\phi_x\phi_z(\phi_{xy} \phi_{yz} - \phi_{xz} \phi_{yy})}<br>
{(\phi_x^2 + \phi_y^2 + \phi_z^2)^2} \nonumber<br>
\end{eqnarray}<br>
<br>
The computation of Gaussian curvature is used in the computation of<br>
minimum curvature.As far<br>
as I know, either Gaussian or minimum curvature are not currently<br>
implemented at ITK. Anyway,<br>
it is quite easy to extend the mean curvature filter at ITK in order to<br>
compute Gaussian and minimum<br>
curvatures using this formula.<br>
<br>
The level set method uses the implicit representation \phi for all the<br>
computations. Therefore, the<br>
Gaussian curvature given in the formula is not just the curvature at the<br>
surface but at all the level<br>
sets of the surface in the whole computation domain. In order to recover<br>
the Gaussian curvature<br>
at the points of the surface you have to seek for the values of K_Gauss<br>
at \phi = 0.<br>
<br>
Best,<br>
Monica.<br>
<br>
<a href="mailto:insight-users-request@itk.org">insight-users-request@itk.org</a> wrote:<br>
&gt; Send Insight-users mailing list submissions to<br>
&gt; &nbsp; &nbsp; &nbsp; <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
&gt;<br>
&gt; To subscribe or unsubscribe via the World Wide Web, visit<br>
&gt; &nbsp; &nbsp; &nbsp; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt; or, via email, send a message with subject or body 'help' to<br>
&gt; &nbsp; &nbsp; &nbsp; <a href="mailto:insight-users-request@itk.org">insight-users-request@itk.org</a><br>
&gt;<br>
&gt; You can reach the person managing the list at<br>
&gt; &nbsp; &nbsp; &nbsp; <a href="mailto:insight-users-owner@itk.org">insight-users-owner@itk.org</a><br>
&gt;<br>
&gt; When replying, please edit your Subject line so it is more specific<br>
&gt; than "Re: Contents of Insight-users digest..."<br>
&gt;<br>
&gt;<br>
&gt; Today's Topics:<br>
&gt;<br>
&gt; &nbsp; &nbsp;1. needs help for "Gaussian Curvature" (Kun)<br>
&gt; &nbsp; &nbsp;2. need help with my Code (Nadia Alvi)<br>
&gt; &nbsp; &nbsp;3. which version of FLTK must be uesd with Itk ? (Martine Lefevre)<br>
&gt; &nbsp; &nbsp;4. ITK -TCON 2.0 : Second Life (Luis Ibanez)<br>
&gt; &nbsp; &nbsp;5. How to concatenate itk transformations (Jian Yang)<br>
&gt;<br>
&gt;<br>
&gt; ----------------------------------------------------------------------<br>
&gt;<br>
&gt; Message: 1<br>
&gt; Date: Thu, 28 Aug 2008 06:01:26 +0800 (CST)<br>
&gt; From: Kun &lt;<a href="mailto:kun_ch@163.com">kun_ch@163.com</a>&gt;<br>
&gt; Subject: [Insight-users] needs help for "Gaussian Curvature"<br>
&gt; To: "<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>" &lt;<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>&gt;<br>
&gt; Message-ID:<br>
&gt; &nbsp; &nbsp; &nbsp; &lt;<a href="mailto:27132818.509911219874486411.JavaMail.coremail@bj163app121.163.com">27132818.509911219874486411.JavaMail.coremail@bj163app121.163.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset="gbk"<br>
&gt;<br>
&gt; &nbsp;Hi, all<br>
&gt;<br>
&gt; &nbsp; &nbsp; I meet a problem these days, and need help.<br>
&gt; &nbsp; &nbsp; I try to use Gaussian Curvature in the Level Set Method. Suppose the Level Set Function phi = f(x,y,z), and the zero-level-set-surface is phi=0, then how can I calculate the Gaussian Curvature of the surface?<br>

&gt;<br>
&gt; And I think the curvature should be calculated using phi, is it right?<br>
&gt;<br>
&gt; Or is there any existing code to do this? Currently, I only find the "Mean Curvature" code in ITK, but can't &nbsp;find the "Gaussian Curvature" Code.<br>
&gt;<br>
&gt; &nbsp; &nbsp; Could anybody help me?<br>
&gt;<br>
&gt; Thanks all.<br>
&gt;<br>
&gt; Kun<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL: &lt;<a href="http://www.itk.org/pipermail/insight-users/attachments/20080828/df5ecea6/attachment-0001.htm" target="_blank">http://www.itk.org/pipermail/insight-users/attachments/20080828/df5ecea6/attachment-0001.htm</a>&gt;<br>

&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 2<br>
&gt; Date: Thu, 28 Aug 2008 01:14:49 -0700 (PDT)<br>
&gt; From: Nadia Alvi &lt;<a href="mailto:nadiaalvi@yahoo.com">nadiaalvi@yahoo.com</a>&gt;<br>
&gt; Subject: [Insight-users] need help with my Code<br>
&gt; To: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
&gt; Message-ID: &lt;<a href="mailto:410155.11425.qm@web31006.mail.mud.yahoo.com">410155.11425.qm@web31006.mail.mud.yahoo.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset="iso-8859-1"<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; Hi All,<br>
&gt; ?<br>
&gt; Please see the code attached.<br>
&gt; It should read the 3 slices and write them as a tiff file but it only reads the 2nd slice and write it 3 times in a tiff file. Apparently it is incrementing the start index by 1 and reading that slice 3 times.<br>
&gt; Can some body please help me in finding what i m doing wrong?<br>
&gt; Many many thanks<br>
&gt; Nadia<br>
&gt; ?<br>
&gt; ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////<br>
&gt;<br>
&gt; typedef unsigned char?????????????????????? PixelType;<br>
&gt; ? const unsigned int Dimension = 3;<br>
&gt; ?<br>
&gt; ? typedef itk::Image&lt; PixelType, Dimension &gt;? ImageType;<br>
&gt; ?<br>
&gt; ? typedef itk::ImageSeriesReader&lt; ImageType &gt;? ReaderType;<br>
&gt; ? typedef itk::ImageFileWriter&lt;ImageType&gt; WriterType;<br>
&gt;<br>
&gt; ? WriterType::Pointer writer = WriterType::New();<br>
&gt; ? ReaderType::Pointer reader1 = ReaderType::New();<br>
&gt; ? typedef itk::NumericSeriesFileNames??? NameGeneratorType;<br>
&gt; ? NameGeneratorType::Pointer nameGenerator = NameGeneratorType::New();<br>
&gt; ? nameGenerator-&gt;SetSeriesFormat( "clint_binary_080320_noDescription%04d.tiff" );<br>
&gt; ? nameGenerator-&gt;SetStartIndex(1);<br>
&gt; ? nameGenerator-&gt;SetEndIndex(3);<br>
&gt; ? nameGenerator-&gt;SetIncrementIndex(1);<br>
&gt;<br>
&gt; ? reader1-&gt;SetImageIO( itk::TIFFImageIO::New() );<br>
&gt; ? reader1-&gt;SetFileNames( nameGenerator-&gt;GetFileNames()? );<br>
&gt;<br>
&gt; ? writer-&gt;SetFileName("labeledImage.tif");<br>
&gt; ? writer-&gt;SetInput(reader1-&gt;GetOutput());<br>
&gt; ? try<br>
&gt; ??? {<br>
&gt; ??<br>
&gt; ??writer-&gt;Update();<br>
&gt; ??? }<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL: &lt;<a href="http://www.itk.org/pipermail/insight-users/attachments/20080828/3f5fe742/attachment-0001.htm" target="_blank">http://www.itk.org/pipermail/insight-users/attachments/20080828/3f5fe742/attachment-0001.htm</a>&gt;<br>

&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 3<br>
&gt; Date: Thu, 28 Aug 2008 10:30:59 +0000 (GMT)<br>
&gt; From: Martine Lefevre &lt;<a href="mailto:martine_lef@yahoo.fr">martine_lef@yahoo.fr</a>&gt;<br>
&gt; Subject: [Insight-users] which version of FLTK must be uesd with Itk ?<br>
&gt; To: insight itk &lt;<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>&gt;<br>
&gt; Message-ID: &lt;<a href="mailto:246756.49418.qm@web28311.mail.ukl.yahoo.com">246756.49418.qm@web28311.mail.ukl.yahoo.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset="utf-8"<br>
&gt;<br>
&gt; Hi<br>
&gt; I?try to build InsightApplications-3.8.0 and as I have read there is some application which need FLTK.<br>
&gt; For this reason I have downloded FLTK version 1.1.6 since based on the mailing list its the best to use with itk. I configure it using Cmake and when I try to build it using VC Net, I get a lot of errors and one fatal error related to some library.<br>

&gt; So my questions are:<br>
&gt; Is it relay the?version to be used with itk ? if yes, am I wrong to configure it using Cmake ?<br>
&gt; What's the way to configure and built such package? I have read in?readme file of this package that we need?just? to open the solution FLTK (in the visualc folder) and built it using the compiler ....please could you give me more suggestions on how to built and integrate FLTK especialy wich version should be used with itk ?<br>

&gt; Waiting for your reply !<br>
&gt; Many thanks<br>
&gt; Martine<br>
&gt;<br>
&gt;<br>
&gt; &nbsp; &nbsp; &nbsp; _____________________________________________________________________________<br>
&gt; Envoyez avec Yahoo! Mail. Une boite mail plus intelligente <a href="http://mail.yahoo.fr" target="_blank">http://mail.yahoo.fr</a><br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL: &lt;<a href="http://www.itk.org/pipermail/insight-users/attachments/20080828/1eeadbb0/attachment-0001.htm" target="_blank">http://www.itk.org/pipermail/insight-users/attachments/20080828/1eeadbb0/attachment-0001.htm</a>&gt;<br>

&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 4<br>
&gt; Date: Thu, 28 Aug 2008 11:47:53 -0400<br>
&gt; From: Luis Ibanez &lt;<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>&gt;<br>
&gt; Subject: [Insight-users] ITK -TCON 2.0 : Second Life<br>
&gt; To: Insight Users &lt;<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>&gt;<br>
&gt; Message-ID: &lt;<a href="mailto:48B6C8A9.9030406@kitware.com">48B6C8A9.9030406@kitware.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
&gt;<br>
&gt;<br>
&gt; Just a reminder that tomorrow we will hold again the<br>
&gt; ITK Tcon 2.0 in Second Life at 1pm EST.<br>
&gt;<br>
&gt; You can teleport to the Tcon site (the "Linden Open Source Project HQ")<br>
&gt; by using the following link:<br>
&gt; <a href="http://slurl.com/secondlife/Hippotropolis/238/14/24/?img=http%3A//public.kitware.com/Insight/Doxygen/html/itkLogo.jpg&amp;title=ITK%20Tcon&amp;msg=Insight%20Toolkit%20%28ITK%29%20Tcon" target="_blank">http://slurl.com/secondlife/Hippotropolis/238/14/24/?img=http%3A//public.kitware.com/Insight/Doxygen/html/itkLogo.jpg&amp;title=ITK%20Tcon&amp;msg=Insight%20Toolkit%20%28ITK%29%20Tcon</a><br>

&gt;<br>
&gt;<br>
&gt; &nbsp; &nbsp; Luis<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; Message: 5<br>
&gt; Date: Thu, 28 Aug 2008 11:51:55 -0400<br>
&gt; From: "Jian Yang" &lt;<a href="mailto:yaland1977@gmail.com">yaland1977@gmail.com</a>&gt;<br>
&gt; Subject: [Insight-users] How to concatenate itk transformations<br>
&gt; To: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
&gt; Message-ID:<br>
&gt; &nbsp; &nbsp; &nbsp; &lt;<a href="mailto:a17c5eac0808280851l175da73fp32e7c09dc94c7b6@mail.gmail.com">a17c5eac0808280851l175da73fp32e7c09dc94c7b6@mail.gmail.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset=ISO-8859-1<br>
&gt;<br>
&gt; Hi all, I have serials of itk transformations, which includes Affine<br>
&gt; and BSpline transformations, I am wondering how to concatenate them<br>
&gt; together and get a combined transformation?<br>
&gt;<br>
&gt; I checked the insight journal and find two classes related to my<br>
&gt; topic: Combining<br>
&gt; Transforms(<a href="http://www.insight-journal.org/dspace/view_item.php?itemid=438" target="_blank">http://www.insight-journal.org/dspace/view_item.php?itemid=438</a>)<br>
&gt; and Serial Transformations<br>
&gt; (<a href="http://www.insight-journal.org/dspace/view_item.php?itemid=1080" target="_blank">http://www.insight-journal.org/dspace/view_item.php?itemid=1080</a>).<br>
&gt; However it seems both of them doing the transformations in sequence<br>
&gt; instead of combining them together to get a final form of<br>
&gt; transformation. If the number of the transformations is very big, the<br>
&gt; sequence transformation of image will no doubt downgrade the image<br>
&gt; qualities.<br>
&gt;<br>
&gt; Your help will greatly appreciated.<br>
&gt;<br>
&gt; Jian<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Insight-users mailing list<br>
&gt; <a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><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>
&gt; End of Insight-users Digest, Vol 52, Issue 49<br>
&gt; *********************************************<br>
&gt;<br>
&gt;<br>
<br>
_______________________________________________<br>
Insight-users mailing list<br>
<a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</div><br>

</div>
</blockquote>