<DIV>
<DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">
<DIV>Hi Edi,</DIV>
<DIV> </DIV>
<DIV>I actually went ahead and installed ITK 1.8.0 and tried to integrate the new optimizer. However the program just runs through without doing anything (not even going through any iterations) if I have the following code:</DIV>
<DIV> </DIV>
<DIV> typedef itk::LBFGSBOptimizer OptimizerType;<BR> <BR> typedef itk::MeanSquaresImageToImageMetric< <BR> FixedImageType, <BR> MovingImageType > MetricType;</DIV>
<DIV> </DIV>
<DIV> MetricType::Pointer metric = MetricType::New();<BR> OptimizerType::Pointer optimizer = OptimizerType::New();</DIV>
<DIV> </DIV>
<DIV> // optimizer->SetCostFunction( metric.GetPointer() );</DIV>
<DIV> </DIV>
<DIV> const double F_Convergence_Factor = 1e+7; // Function value tolerance<BR> const double Projected_G_Tolerance = 1e-5; // Proj gradient tolerance <BR> const int Max_Iterations = 100; // Maximum number of iterations</DIV>
<DIV> optimizer->SetCostFunctionConvergenceFactor( F_Convergence_Factor );<BR> optimizer->SetProjectedGradientTolerance( Projected_G_Tolerance );<BR> optimizer->SetMaximumNumberOfIterations( Max_Iterations );<BR> optimizer->SetMaximumNumberOfEvaluations( Max_Iterations );</DIV>
<DIV><BR> const unsigned int SpaceDimension = 5;</DIV>
<DIV><BR> OptimizerType::BoundValueType lower(SpaceDimension);<BR> OptimizerType::BoundValueType upper(SpaceDimension);<BR> OptimizerType::BoundSelectionType select(SpaceDimension);</DIV>
<DIV><BR> lower.Fill( -5 );<BR> upper.Fill( 5 );<BR> select.Fill( SpaceDimension );</DIV>
<DIV> optimizer->SetLowerBound( lower );<BR> optimizer->SetUpperBound( upper );<BR> optimizer->SetBoundSelection( select );</DIV>
<DIV><BR>However if I uncomment the line: optimizer->SetCostFunction( metric.GetPointer() );</DIV>
<DIV>then it crashes when it tries to access m_transform within itk::MeanSquaresImageToImageMetric. I am wondering how I can include the cost function (since it doesn't seem to do anything otherwise) and am I doing something else wrong that is causing the optimizer to have no effect. </DIV>
<DIV> </DIV>
<DIV>Thanks again and sorry for the repeated questions but I have been trying to get this thing to run for a long time!<BR></DIV>
<DIV>Robert</DIV>
<DIV><BR><B><I>Eduard Schreibmann <eduards@stanford.edu></I></B> wrote:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">
<DIV>
<DIV>Hi Robert.</DIV>
<DIV> </DIV>
<DIV>I think it is a fairlly new adition to ITK, even newer than 1.8 maybe. it appeared suddenly in the latest manual documentation.</DIV>
<DIV> </DIV>
<DIV>Attached is the code from my computer, maybe you can get along without compiling, by just adding itkLBFGSB files to your project (that is compiling it in your oproject rather than in the ITK library). It is not an elegant solution but I think it may work, and anyway it is not a big deal to try this before recompiling ITK. There is also an example of how to set limits. With my images, I had problems with LBFGSB because it was not making many iterations. it took days to find out it is important to use a good BSpline interpolator, so if you have problems with LBFGSB not converging, it may be that you need a good interpolator </DIV>
<DIV> </DIV>
<DIV>Cheers and good luck with the code.</DIV>
<DIV>Edi</DIV>
<DIV> </DIV>
<DIV>Edi<BR><BR></DIV></DIV></BLOCKQUOTE></BLOCKQUOTE></DIV></DIV><p>
<hr size=1>Do you Yahoo!?<br>
<a href="http://us.rd.yahoo.com/mail_us/taglines/aac/*http://promotions.yahoo.com/new_mail/static/ease.html">Yahoo! Mail Address AutoComplete</a> - You start. We finish.