<div dir="ltr"><span style="font-size:12.8000001907349px">Dear all,</span><div><span style="font-size:12.8000001907349px"><br></span><div style="font-size:12.8000001907349px">I have done some experiments with ITK FEM registration and, following preceding advice, I have solved the crashing issue I've had. Nevertheless I have still some questions. </div><div style="font-size:12.8000001907349px"><br>In fact, I am interested in developing a small piece of code using itk FEM module such that:</div><div style="font-size:12.8000001907349px">1) takes an image as an input;</div><div style="font-size:12.8000001907349px">2) defines a grid on the image;</div><div style="font-size:12.8000001907349px">3) applies some forces on the boundary;</div><div style="font-size:12.8000001907349px">4) deforms the image accordingly.</div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px"> Since ITK has a FEM module I was trying to use it. So, I have started reading <a href="http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html" target="_blank">itk::fem::Solver</a> documentation which states:<br></div><div style="font-size:12.8000001907349px"><dl style="font-stretch:normal;font-size:12px;line-height:22px;font-family:'Lucida Grande',Verdana,Geneva,Arial,sans-serif;padding:0px;margin-left:0px;color:rgb(0,0,0)"><dd style="margin-bottom:6px">The standard way to setup a FE problem in ITK is to use the following approach.</dd></dl><div style="font-stretch:normal;font-size:12px;line-height:22px;font-family:'Lucida Grande',Verdana,Geneva,Arial,sans-serif;padding:0px;margin:0px;border:1px solid rgb(196,207,229);color:rgb(0,0,0);background-color:rgb(251,252,253)"><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       <span style="color:rgb(0,128,0)">typedef</span> <a href="http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1FEMObject.html" target="_blank" style="color:rgb(70,101,162);text-decoration:none">itk::fem::FEMObject<3></a>    <a href="http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html#aa7503c1b694280d0f47bdfb3abc98ebc" target="_blank" style="color:rgb(70,101,162);text-decoration:none">FEMObjectType</a>;</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       FEMObjectObjectType::Pointer fem = FEMObjectObjectType::New();</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       ...</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       <span style="color:rgb(0,128,0)">typedef</span> <a href="http://www.itk.org/Doxygen/html/classitk_1_1fem_1_1Solver.html" target="_blank" style="color:rgb(70,101,162);text-decoration:none">itk::fem::Solver<3></a>    FEMSolverType;</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       FEMSolverType::Pointer solver = FEMSolverType::New();</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       solver->SetInput( fem );</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       solver->Update( );</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*       FEMSolverType::Pointer defem = solver->GetOutput( );</div><div style="font-stretch:normal;font-size:13px;line-height:1;font-family:monospace,fixed;min-height:13px;white-space:pre-wrap;word-wrap:break-word;padding-left:53px;padding-bottom:0px;margin:0px">*   ...</div></div></div><div class="gmail_extra" style="font-size:12.8000001907349px"><br><div class="gmail_quote"><br></div><div class="gmail_quote">Unfortunately, apart fro the typo (defem should be FEMObjectType), this is rather succinct. <span style="font-size:12.8000001907349px">Do you know some places in which to look for doing something like what I have described? </span></div><div class="gmail_quote"><br></div><div class="gmail_quote">Best regards,</div><div class="gmail_quote"><br></div><div class="gmail_quote">Luigi</div></div></div></div>