<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2800.1491" name=GENERATOR></HEAD>
<BODY>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2>Stefan, Marius, Mathieu,</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff size=2>
Excellent. I liked the try / catch blocks, and the stopping criteria in
Mathieu's implementation. Essentially, the code is the same across all 3
implementations... As you would expect from the nicely written
paper.</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff size=2>- what
should be the default values for a, A, C?</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2> - I used a=A=c=0.1 (Don't trust my
implementation...)</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2> - Mathieu uses a=01.6, A = 60(!),
C=0.01</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2> - Stefan uses a=A=c=1.0</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff size=2>- Are
we OK with Alpha = 0.602, Gamma = 0.101?</FONT></SPAN></DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=287370114-17032005><FONT face=Arial color=#0000ff
size=2>-dan</FONT></SPAN></DIV>
<BLOCKQUOTE>
<DIV class=OutlookMessageHeader dir=ltr align=left><FONT face=Tahoma
size=2>-----Original Message-----<BR><B>From:</B> Stefan Klein
[mailto:stefan@isi.uu.nl]<BR><B>Sent:</B> Thursday, March 17, 2005 8:55
AM<BR><B>To:</B> Blezek, Daniel J (Research); Mathieu De Craene<BR><B>Cc:</B>
Insight-users@itk.org; insight-developers@itk.org<BR><B>Subject:</B> RE:
[Insight-users] some contributions<BR><BR></FONT></DIV><FONT size=3>Eh, but
then we'll have 2 SPSA-optimizers....<BR>Wouldn't it be the best to have one
that combines the good of all 3 versions?<BR><BR>plan:<BR>- we include the
scales thing and the test and some elements from Mathieu's implementation<BR>-
we'll mail it to the list.<BR>- you check it (if we combined all features
correctly) and correct it<BR>- Daniel commits it to the
itk.<BR><BR>Marius&Stefan.<BR><BR><BR><BR>At 08:39 17/03/05, Blezek,
Daniel J (Research) wrote:<BR>
<BLOCKQUOTE class=cite cite="" type="cite">Oddly enough, I also have an
implementation of the SPSA algorithm! I offered to put Stefan's
implementation in, provided he account for m_Scale and write a test.
I'll make the same offer for you, Mathieu. I've put my code on the itk
Wiki @ <A href="http://itk.org/Wiki/SPSAOptimizer%A0"
eudora="autourl">http://itk.org/Wiki/SPSAOptimizer </A>however, I haven't
tested my code very well, as it's new.<BR><BR>-dan<BR><BR>-----Original
Message-----<BR>From: insight-users-bounces@itk.org<BR>[<A
href="mailto:insight-users-bounces@itk.org"
eudora="autourl">mailto:insight-users-bounces@itk.org</A>]On Behalf Of
Mathieu De Craene<BR>Sent: Thursday, March 17, 2005 8:20 AM<BR>To: Stefan
Klein<BR>Cc: Insight-users@itk.org; insight-developers@itk.org<BR>Subject:
Re: [Insight-users] some contributions<BR><BR><BR>Le mercredi 16 mars 2005 à
21:34 +0100, Stefan Klein a écrit :<BR>> Hi all,<BR>> <BR>> We have
created some classes that are useful for us. We hope also for other <BR>>
people, so we would like to contribute them. The code can be found for
<BR>> download on:<BR>>
<X-TAB> </X-TAB><A
href="http://www.isi.uu.nl/People/Stefan/"
eudora="autourl">http://www.isi.uu.nl/People/Stefan/</A><BR>> <BR>> -
itkEulerTransform.h: This class is a wrap around the Euler2DTransform and
<BR>> the Euler3DTransform. In this way it is templated over the
dimension, like <BR>> other Transform-classes.<BR>> <BR>> -
SimultaneousPerturbation.zip: This class implements the so-called <BR>>
Simultaneous Perturbation method as an ITK optimizer. For more info on this
<BR>> optimization method, go to:<BR>>
<X-TAB> </X-TAB><A
href="http://www.jhuapl.edu/SPSA"
eudora="autourl">http://www.jhuapl.edu/SPSA</A><BR>> <BR><BR>I have also
have some code for SPSA optimization. I found it quite<BR>useful for rigid
registration (it seems to be pretty robust). The source<BR>code for a rigid
registration application based on ITK is available on <BR><BR><A
href="http://euterpe.tele.ucl.ac.be/Waleo2/vesale/"
eudora="autourl">http://euterpe.tele.ucl.ac.be/Waleo2/vesale/</A><BR><BR>I
would like to compare our implementations of this
optimizer.<BR><BR>Regards,<BR><BR>Mathieu.<BR><BR>> - FullSearch.zip:
This class implements a (semi-)Full Search optimization <BR>> routine.
This can be useful for research, to scan the optimization surface, <BR>>
or to evaluate the result of other optimizers. It allows the user to<BR>>
specify which parameters should be evaluated, in which range.<BR>>
<BR>> Since we have contributed 3 things now, we also have a feature
request; as <BR>> a reward :-)<BR>> <BR>> Would it be possible to
add a public function to the <BR>>
MattesMutualInformationImageToImageMetric that forces the metric to select
<BR>> new spatial samples? This would be useful if you would want new
spatial <BR>> samples every iteration for example.<BR>> <BR>> We
hope that the contributions are useful. Any comments are welcome!!<BR>>
<BR>> Marius and Stefan.<BR>> <BR>> <BR>> Marius Staring, Stefan
Klein<BR>> Image Sciences Institute<BR>> University Medical Centre
Utrecht<BR>> Heidelberglaan 100, 3584 CX Utrecht, The Netherlands<BR>>
phone: +31 (0)30 250 3186, fax: +31 (0)30 251 3399<BR>>
{marius,stefan}@isi.uu.nl,<BR>> <A
href="http://www.isi.uu.nl/People/Marius"
eudora="autourl">http://www.isi.uu.nl/People/Marius</A><BR>> <A
href="http://www.isi.uu.nl/People/?stefan"
eudora="autourl">http://www.isi.uu.nl/People/?stefan</A><BR>> <BR>>
<BR>> <BR>> <BR>> <BR>> <BR>> <BR>> <BR>>
_______________________________________________<BR>> Insight-users
mailing list<BR>> Insight-users@itk.org<BR>> <A
href="http://www.itk.org/mailman/listinfo/insight-users"
eudora="autourl">http://www.itk.org/mailman/listinfo/insight-users</A><BR>>
<BR><BR>_______________________________________________<BR>Insight-users
mailing list<BR>Insight-users@itk.org<BR><A
href="http://www.itk.org/mailman/listinfo/insight-users"
eudora="autourl">http://www.itk.org/mailman/listinfo/insight-users</A>
</FONT></BLOCKQUOTE><BR></BLOCKQUOTE></BODY></HTML>