<div dir="ltr">Hi Matt, <div><br></div><div>Thanks the quick answer and welcome :) . That is indeed the ppt. There they define two energies and then write an integral in terms of some "similarity" which wasn't really defined. Maybe it is the E_D term? Hopefully the same applies for smoothness (E_S), but I'm not sure. In the code I had trouble because RunRegistration uses MultiResSolve, and there they instantiate a solver which uses Update(). I couldn't find Update in either .h or .hxx files, and apparently update is redirected (somehow) to AssembleKandM. The Doxygen page for the CrankNicolson solver shows a formulation of the method but not for what problem, and that leaves other questions without answer, such as the boundary conditions, elements used for the FEM formulation, etc. Thanks so much for your time. </div><div><br></div><div>Best regards</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 23 March 2017 at 18:34, Matt McCormick <span dir="ltr"><<a href="mailto:matt.mccormick@kitware.com" target="_blank">matt.mccormick@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Nicolás,<br>
<br>
Welcome to ITK!<br>
<br>
<br>
Is this the PowerPoint?<br>
<br>
  <a href="https://itk.org/CourseWare/Training/NonRigidRegistrationMethods" rel="noreferrer" target="_blank">https://itk.org/CourseWare/<wbr>Training/<wbr>NonRigidRegistrationMethods</a><br>
<br>
Note that DeformableRegistration1.cxx (FEM registration) is discussed<br>
after the optical flow example.<br>
<br>
<br>
Some more information on the FEMRegistrationFilter can be found in its<br>
class documentation:<br>
<br>
  <a href="https://itk.org/Doxygen/html/classitk_1_1fem_1_1FEMRegistrationFilter.html" rel="noreferrer" target="_blank">https://itk.org/Doxygen/html/<wbr>classitk_1_1fem_1_<wbr>1FEMRegistrationFilter.html</a><br>
<br>
Multiple metrics are available:<br>
<br>
  <a href="https://itk.org/Doxygen/html/classitk_1_1fem_1_1FEMRegistrationFilter.html#acc636f4752e592cd780503a5fbfeba82" rel="noreferrer" target="_blank">https://itk.org/Doxygen/html/<wbr>classitk_1_1fem_1_<wbr>1FEMRegistrationFilter.html#<wbr>acc636f4752e592cd780503a5fbfeb<wbr>a82</a><br>
<br>
<br>
In general, "the code reveals all" and should be given the highest<br>
amount of trust.<br>
<br>
<br>
Hope this helps,<br>
Matt<br>
<div><div class="h5"><br>
On Thu, Mar 23, 2017 at 4:32 PM, Nicolás Barnafi <<a href="mailto:nabw91@gmail.com">nabw91@gmail.com</a>> wrote:<br>
> Hi everyone,<br>
><br>
> I have been trying to understand exactly what is happening in the<br>
> DeformableRegistration1.h file without much success. The problem is:<br>
> according to the software guide, the problem being solved comes from the<br>
> variational problem given by<br>
><br>
> min D[image1, image2; u] + S[u]<br>
><br>
> where D is just the SSD (or the L2 norm of (Im1 - Im2 o phi), with phi the<br>
> unknown displacement field) and S is a linear elastic potential. From here<br>
> you get  the euler lagrange equations (asuming some unspecified boundary<br>
> ocndition) and solve it using some semi implicit newton-raphson scheme. This<br>
> is where it starts getting blurry, because the ITK ppt on deformable<br>
> registration first shows an optical flow formulation, which would mean that<br>
> the SSD metric isn't really what is being used, and also if I dig deeper in<br>
> the code, I find actually a Crank-Nicolson scheme being used, which really<br>
> implies some kind of temporality that really does not exist in the<br>
> variational formulation. The only hint I have found was in Modersitzki's<br>
> book where a fixed point scheme is artificially stabilized:<br>
><br>
> A(u[k+1]) = f_u[k]<br>
><br>
> => u[k+1] + t A(u[k+1]) = t f_u[k] + u[k].<br>
><br>
> I would want to know what is exactly happening in that example to be able to<br>
> validate an example I implemented in python. Thanks for your time.<br>
><br>
> Best regards<br>
><br>
><br>
><br>
> --<br>
> Nicolás Alejandro Barnafi Wittwer<br>
><br>
</div></div>> ______________________________<wbr>_________________<br>
> Community mailing list<br>
> <a href="mailto:Community@itk.org">Community@itk.org</a><br>
> <a href="http://public.kitware.com/mailman/listinfo/community" rel="noreferrer" target="_blank">http://public.kitware.com/<wbr>mailman/listinfo/community</a><br>
><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">Nicolás Alejandro Barnafi Wittwer</div></div></div></div>
</div>