[ITK] Fwd: FEM registration, minimal working example

[Luigi Riba]

Dear Matt,

thank you very much for your reply.

I have done some experiments with ITK FEM registration and, following your
advice, I have solved the crashing issue. Nevertheless I have still some
questions.

In fact, I am interested in developing a small piece of code using itk FEM
module such that:
1) takes an image as an input;
2) defines a grid on the image;
3) applies some forces on the boundary;
4) deforms the image accordingly.

 Since ITK has a FEM module I was trying to use it. So, I have started
reading itk::fem::Solver
<http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html>
documentation which states:
The standard way to setup a FE problem in ITK is to use the following
approach.
* typedef itk::fem::FEMObject<3>
<http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1FEMObject.html>
FEMObjectType
<http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html#aa7503c1b694280d0f47bdfb3abc98ebc>
;
* FEMObjectObjectType::Pointer fem = FEMObjectObjectType::New();
* ...
* typedef itk::fem::Solver<3>
<http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html>
FEMSolverType;
* FEMSolverType::Pointer solver = FEMSolverType::New();
*
* solver->SetInput( fem );
* solver->Update( );
* FEMSolverType::Pointer defem = solver->GetOutput( );
* ...

Unfortunately, apart fro the typo (defem should be FEMObjectType), this is
rather succinct. Do you know some places in which to look for doing
something like what I have described?

Best regards and thank you for your work,

Luigi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/community/attachments/20150704/e1df4b49/attachment.html>