<HTML><BODY style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><DIV><DIV>Alpha is used solely by the solver routine (itkFEMSolverCrankNicolson) and determines the relative weighting of previous solutions in calculating the current solution.  Perhaps looking  at the header of the specified class will give you more insight.  Gamma is used in a single place in the registration filter for updating the current solution.  I'd have to ask my colleague (who, unfortunately, is not here at the moment) regarding it's use.  However, I do know that the default value should be 1.0 and that should be adequate for most applications.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>You're right about the elasticity and rho parameters.  Try altering the time step for that could have a similar effect.   Regarding RhoC --- it is used only in the material class.  I'll have to look closer to see if we can dispense with its specification and simply derive it from Rho and E.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Nick</DIV><DIV><BR class="khtml-block-placeholder"></DIV><BR><DIV><DIV>On Apr 6, 2006, at 10:15 AM, Markus Weigert wrote:</DIV><BLOCKQUOTE type="cite"><DIV><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">Hi Nick,</SPAN></FONT></DIV><DIV> </DIV><DIV><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">thanks a lot for your help. </SPAN></FONT><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">I will try it with your values again.</SPAN></FONT></DIV><DIV><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">Can you perhaps explain the parameters alpha and gamma</SPAN></FONT></DIV><DIV><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">a little bit more detailed for me, since I didn't find information about the </SPAN></FONT><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">real meaning of these parameters.</SPAN></FONT></DIV><DIV><FONT class="Apple-style-span" face="Arial" size="2"><SPAN class="Apple-style-span" style="font-size: 10px;">Nearly all the posts only deal with elasticity and rho.</SPAN></FONT></DIV><DIV> </DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV></DIV><BR><BLOCKQUOTE type="cite"><SPAN class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">Another question is:</SPAN></FONT></DIV><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">If I change </SPAN><FONT face="Times New Roman" size="3"><SPAN class="Apple-style-span" style="font-family: Times New Roman; ">elasticity to 1.0e4 and Rho to 1.0e4, than the deformations should</SPAN></FONT></FONT></DIV><DIV>even become bigger (I used 1.0e5 and the deformation is already much to high).</DIV><DIV><FONT face="Arial" size="2"></FONT> </DIV><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">Also one has to specify rho and rhoc in the filter, but as far as I remember, rhoC should always</SPAN></FONT></DIV><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">be approx. 1.0 and equals elasticity divided by rho (or vice versa), so why is it specified separately?</SPAN></FONT></DIV><DIV><FONT face="Arial" size="2"></FONT> </DIV><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">Best regards,</SPAN></FONT></DIV><DIV><FONT face="Arial" size="2"><SPAN class="Apple-style-span" style="font-family: Arial; font-size: 10px; ">Markus</SPAN></FONT></DIV><DIV><FONT face="Arial" size="2"></FONT> </DIV><DIV><FONT face="Arial" size="2"></FONT> </DIV><BLOCKQUOTE style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px"><DIV style="FONT: 10pt arial; font-family: arial; font-size: 13.3333px; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; ">----- Original Message -----</SPAN></DIV><DIV style="BACKGROUND: #e4e4e4; FONT: 10pt arial; font-color: black; font-family: arial; font-size: 13.3333px; "><B style="font-family: arial; font-size: 13.3333px; font-weight: bold; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; font-weight: bold; ">From:</SPAN></B><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "> </SPAN><A title="ntustison@gmail.com" href="mailto:ntustison@gmail.com"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); font-family: arial; font-size: 13.3333px; -khtml-text-decorations-in-effect: underline; ">Nicholas Tustison</SPAN></A><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "></SPAN></DIV><DIV style="FONT: 10pt arial; font-family: arial; font-size: 13.3333px; "><B style="font-family: arial; font-size: 13.3333px; font-weight: bold; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; font-weight: bold; ">To:</SPAN></B><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "> </SPAN><A title="m.weigert@fz-juelich.de" href="mailto:m.weigert@fz-juelich.de"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); font-family: arial; font-size: 13.3333px; -khtml-text-decorations-in-effect: underline; ">Markus Weigert</SPAN></A><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "></SPAN></DIV><DIV style="FONT: 10pt arial; font-family: arial; font-size: 13.3333px; "><B style="font-family: arial; font-size: 13.3333px; font-weight: bold; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; font-weight: bold; ">Cc:</SPAN></B><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "> </SPAN><A title="insight-users@itk.org" href="mailto:insight-users@itk.org"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); font-family: arial; font-size: 13.3333px; -khtml-text-decorations-in-effect: underline; ">insight-users@itk.org</SPAN></A><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "></SPAN></DIV><DIV style="FONT: 10pt arial; font-family: arial; font-size: 13.3333px; "><B style="font-family: arial; font-size: 13.3333px; font-weight: bold; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; font-weight: bold; ">Sent:</SPAN></B><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "> Thursday, April 06, 2006 3:2PM</SPAN></DIV><DIV style="FONT: 10pt arial; font-family: arial; font-size: 13.3333px; "><B style="font-family: arial; font-size: 13.3333px; font-weight: bold; "><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; font-weight: bold; ">Subject:</SPAN></B><SPAN class="Apple-style-span" style="font-family: arial; font-size: 13.3333px; "> Re: [Insight-users] Setting up FEM - Registration parameters</SPAN></DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR></DIV><DIV>Hi Markus,</DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR class="khtml-block-placeholder"></DIV><DIV>The registration scenario you describe seems reasonably straightforward.  However, there are a couple of items I noticed from the images that you might want to change.  Generally, you want to keep alpha = 0.5 (as opposed to 1.0) so that the solution is stable for the time step that you choose.  The chosen material governs the regularization of your registration results.  The default values that I have are Elasticity = 1.0e4 and Rho = 1.0e4 which give me reasonable registration results.  You should be able to achieve reasonable results using the specified elastic material.  Unfortunately, without running the registration myself, it is difficult to know exactly what the problem is.  </DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR class="khtml-block-placeholder"></DIV><DIV>On a related note, it is interesting that you chose to forgo the .config file.  I actually rewrote the itkFEMRegistrationFilter and plan to include it in an upcoming release along with additional modifications and extensions to the library (including B-spline elements).  The first thing I got rid of was the .config file along with some parameters that weren't being used.</DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR class="khtml-block-placeholder"></DIV><DIV>Nick Tustison </DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR class="khtml-block-placeholder"></DIV><DIV><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><FONT face="Arial" size="2"></FONT><BR class="khtml-block-placeholder"></DIV><BR><DIV><DIV>On Apr 6, 2006, at 7:36 AM, Markus Weigert wrote:</DIV><BR class="Apple-interchange-newline"><BLOCKQUOTE type="cite"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Hi,</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">after scanning the list on FEM - registration, I still have huge problems </SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">to configure the paramters for my specific registration problem with the FEM - method.</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">I want to register two 3D MR 256^3 Thorax images with only small deformation.</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">The problem is, that however I set the parameters </SPAN><B style="color: rgb(0, 0, 221); font-weight: bold; "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); font-weight: bold; ">elasticity</SPAN></B><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); "> and </SPAN><B style="color: rgb(0, 0, 221); font-weight: bold; "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); font-weight: bold; ">rho / rhoc</SPAN></B><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">,</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">it doesn't deliver any useful results.</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">I started with elasticity 100000 (even higher than in the examples), but the image became completely</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">deformed. Also, it spends only 2 iterations in the process (has this something to do with</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">the energyreductionfactor???) and after the registration, the metric (I use SSD) is much bader then before.</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">I think this should not be the case, after what I understood about the theory behind this method.</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Also with elasticity 1000.000, it still gets completely deformed.</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">With even higher the deformation decreases of course, but the missregistration of the lungs is never</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">caught (see middle of the images at the bottom of this message).</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Please note, that I don't use a config file but set all the parameters manually.</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Also I wonder what's the meaning of the parameters alpha and gamma (changing  these doesn't</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">seem to have a great effect on the result).</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Indeed, I call the following methods to setup the FEMRegistrationFilter (public functions),</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">while I didn't find a function to set the scaling, which is also used in the config file (is this intended???): </SPAN><BR style="color: rgb(0, 0, 221); "><FONT color="#ff0000"><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetFixedImage(m_fixedImage);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetMovingImage(m_floatingImage);</SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;DoMultiRes(true);    </SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetMeshPixelsPerElementAtEachResolution(this-&gt;m_meshPixPerElement);</SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetNumberOfIntegrationPoints(this-&gt;m_numberOfIntegrationPoints);</SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetRho(this-&gt;m_rho);</SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetElasticity(this-&gt;m_elasticity);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetGamma(this-&gt;m_gamma);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetWidthOfMetricRegion(this-&gt;m_widthOfMetricRegion);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;ChooseMetric((int) this-&gt;m_metric);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetMaximumIterations(this-&gt;m_maxNumberOfIterations,0);</SPAN><BR style="color: rgb(255, 0, 0); "><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetTimeStep(this-&gt;m_timestep);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetLineSearchMaximumIterations(this-&gt;m_maxLineSearchIterations);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetEnergyReductionFactor(this-&gt;m_energyReductionFactor);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  if(m_metric==3) m_registrationFilter-&gt;SetDescentDirectionMaximize();</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  else{ m_registrationFilter-&gt;SetDescentDirectionMinimize();}</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetAlpha(this-&gt;m_alpha);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;DoLineSearch(2);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;EmployRegridding(1);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;UseLandmarks(false);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetWriteDisplacements(false);</SPAN></FONT><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Also what is about the material properties:</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Does it make sense to change something here??? </SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); "> </SPAN><FONT color="#ff0000"><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">itk::fem::MaterialLinearElasticity::Pointer m;</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m = itk::fem::MaterialLinearElasticity::New();</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;GN = 0;                  // Global number of the material</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;E = m_registrationFilter-&gt;GetElasticity();  // Young's modulus -- used in the membrane</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;A = 1.0;                 // Cross-sectional area</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;h = 1.0;                 // Thickness</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;I = 1.0;                 // Moment of inertia</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;nu = 0.;                 // Poisson's ratio -- DONT CHOOSE 1.0!!</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m-&gt;RhoC = this-&gt;m_density;  // Density</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  // Create the element type </SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  Element3DType::Pointer e1=Element3DType::New();</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  e1-&gt;m_mat=dynamic_cast&lt;itk::fem::MaterialLinearElasticity*&gt;( m );</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetElement(e1);</SPAN><BR style="color: rgb(255, 0, 0); "><SPAN class="Apple-style-span" style="color: rgb(255, 0, 0); ">  m_registrationFilter-&gt;SetMaterial(m);</SPAN></FONT><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Additionally, you can find some screenshots which show the problem and the RegistrationFilter</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">outputs here:</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter2.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter2.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter1.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter1.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/before_registration.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/before_registration.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter2.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/E10%5E8_FEMRegOutput_iter2.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/result_E100000.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/result_E100000.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/result_E1000000.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/result_E1000000.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><A class="moz-txt-link-freetext" href="http://www.markusweigert.de/images/result_E10000000.jpg"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.markusweigert.de/images/result_E10000000.jpg</SPAN></A><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">I hope somebody will have the time to answer at least some of my questions.</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Best regards,</SPAN><BR style="color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Markus</SPAN><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><BR style="color: rgb(0, 0, 221); "><DIV style="MARGIN: 0px; color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">_______________________________________________</SPAN></DIV><DIV style="MARGIN: 0px; color: rgb(0, 0, 221); "><SPAN class="Apple-style-span" style="color: rgb(0, 0, 221); ">Insight-users mailing list</SPAN></DIV><DIV style="MARGIN: 0px; color: rgb(0, 0, 221); "><A href="mailto:Insight-users@itk.org"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">Insight-users@itk.org</SPAN></A></DIV><DIV style="MARGIN: 0px; color: rgb(0, 0, 221); "><A href="http://www.itk.org/mailman/listinfo/insight-users"><SPAN class="Apple-style-span" style="color: rgb(0, 0, 238); -khtml-text-decorations-in-effect: underline; ">http://www.itk.org/mailman/listinfo/insight-users</SPAN></A></DIV></BLOCKQUOTE></DIV><BR></BLOCKQUOTE><BR class="Apple-interchange-newline"></SPAN></BLOCKQUOTE></DIV><BR></BODY></HTML>