<div>Hi,</div>
<div>&nbsp;</div>
<div>When I used Regular Step Gradient Descent optimizer to register two images I saw an abnormal behavior. The metric value first decreased but after some iterations it became larger. See the following output:</div>
<div>&nbsp;</div>
<div>0&nbsp;&nbsp; -0.835581&nbsp;&nbsp; [-0.00321978, -0.0466639, 0.0883861]<br>1&nbsp;&nbsp; -0.836504&nbsp;&nbsp; [-0.00422488, -0.0933436, 0.176817]<br>2&nbsp;&nbsp; -0.83737&nbsp;&nbsp; [-0.00293453, -0.139983, 0.265265]<br>3&nbsp;&nbsp; -0.838174&nbsp;&nbsp; [0.000731481, -0.18652, 0.353701]<br>
4&nbsp;&nbsp; -0.838908&nbsp;&nbsp; [0.00685175, -0.232884, 0.442092]<br>5&nbsp;&nbsp; -0.839563&nbsp;&nbsp; [0.015502, -0.278993, 0.530404]<br>6&nbsp;&nbsp; -0.840129&nbsp;&nbsp; [0.0267537, -0.32476, 0.618602]<br>7&nbsp;&nbsp; -0.840598&nbsp;&nbsp; [0.0406717, -0.370081, 0.706649]<br>8&nbsp;&nbsp; -0.84096&nbsp;&nbsp; [
0.0573129, -0.414844, 0.794509]<br>9&nbsp;&nbsp; -0.84121&nbsp;&nbsp; [0.0767229, -0.458919, 0.882148]<br>10&nbsp;&nbsp; -0.841343&nbsp;&nbsp; [0.0989338, -0.50216, 0.969537]<br>11&nbsp;&nbsp; -0.841358&nbsp;&nbsp; [0.123961, -0.544401, 1.05665]<br>12&nbsp;&nbsp; -0.841256&nbsp;&nbsp; [0.151799, -0.58545
, 1.14349]<br>13&nbsp;&nbsp; -0.841041&nbsp;&nbsp; [0.182419, -0.625084, 1.23004]<br>14&nbsp;&nbsp; -0.840719&nbsp;&nbsp; [0.215761, -0.663034, 1.31634]<br>15&nbsp;&nbsp; -0.840298&nbsp;&nbsp; [0.214669, -0.680096, 1.41487]<br>16&nbsp;&nbsp; -0.839375&nbsp;&nbsp; [0.216594, -0.696938, 1.51342]<br>17&nbsp;&nbsp; -
0.838339&nbsp;&nbsp; [0.221686, -0.713544, 1.6119]<br>18&nbsp;&nbsp; -0.837215&nbsp;&nbsp; [0.230089, -0.729893, 1.7102]<br>19&nbsp;&nbsp; -0.836032&nbsp;&nbsp; [0.241937, -0.74596, 1.80819]<br>20&nbsp;&nbsp; -0.83482&nbsp;&nbsp; [0.257347, -0.76171, 1.90573]<br>21&nbsp;&nbsp; -0.833613&nbsp;&nbsp; [0.276409, -
0.7771, 2.00268]<br>22&nbsp;&nbsp; -0.832445&nbsp;&nbsp; [0.299184, -0.792076, 2.09889]<br>23&nbsp;&nbsp; -0.831349&nbsp;&nbsp; [0.325688, -0.806571, 2.19422]<br>24&nbsp;&nbsp; -0.830357&nbsp;&nbsp; [0.35589, -0.820502, 2.28853]<br>25&nbsp;&nbsp; -0.829494&nbsp;&nbsp; [0.389705, -0.833768, 2.3817]<br>
26&nbsp;&nbsp; -0.828783&nbsp;&nbsp; [0.426981, -0.84624, 2.47365]<br>27&nbsp;&nbsp; -0.828236&nbsp;&nbsp; [0.467495, -0.857749, 2.56435]<br>28&nbsp;&nbsp; -0.827859&nbsp;&nbsp; [0.510938, -0.868055, 2.65383]<br>29&nbsp;&nbsp; -0.827633&nbsp;&nbsp; [0.556896, -0.876774, 2.74221]<br>30&nbsp;&nbsp; -0.827525&nbsp;&nbsp; [
0.60475, -0.883034, 2.8298]<br>31&nbsp;&nbsp; -0.827502&nbsp;&nbsp; [0.650812, -0.875321, 2.91822]<br>32&nbsp;&nbsp; -0.827512&nbsp;&nbsp; [0.629235, -0.881114, 2.87349]<br>33&nbsp;&nbsp; -0.827519&nbsp;&nbsp; [0.603546, -0.883645, 2.83067]<br>34&nbsp;&nbsp; -0.827483&nbsp;&nbsp; [0.620407, -0.880267
, 2.84881]<br>35&nbsp;&nbsp; -0.827564&nbsp;&nbsp; [0.608674, -0.882388, 2.84506]<br>36&nbsp;&nbsp; -0.827459&nbsp;&nbsp; [0.613038, -0.88246, 2.84059]<br>37&nbsp;&nbsp; -0.827529&nbsp;&nbsp; [0.610107, -0.882339, 2.84166]<br>38&nbsp;&nbsp; -0.827493&nbsp;&nbsp; [0.611332, -0.882386, 2.8407]</div>
<div>&nbsp;</div>
<div>However, if I change the optimizer to the LBFGSB optimizer, I got the following output:</div>
<div>&nbsp;</div>
<div>0&nbsp;&nbsp; -0.840282&nbsp;&nbsp; [-0.0321978, -0.466639, 0.883861]<br>1&nbsp;&nbsp; -0.840732&nbsp;&nbsp; [0.0493198, -0.53362, 1.01199]<br>2&nbsp;&nbsp; -0.840735&nbsp;&nbsp; [0.0592661, -0.540774, 1.03432]<br>3&nbsp;&nbsp; -0.840751&nbsp;&nbsp; [0.100359, -0.570829, 1.11035]</div>
<div>&nbsp;</div>
<div>This looks reasonable. The setup for the registration is the following:</div>
<div>&nbsp;</div>
<div>
<p>&nbsp;&nbsp; typedef itk::BSplineInterpolateImageFunction&lt;ImageType, double &gt;&nbsp;&nbsp; InterpolatorType;<br>&nbsp; typedef itk::ImageRegistrationMethod&lt;ImageType, ImageType &gt;&nbsp;&nbsp;&nbsp; RegistrationType;</p>
<p>&nbsp; typedef itk::NormalizedCorrelationImageToImageMetric&lt;ImageType,ImageType &gt;&nbsp;&nbsp;&nbsp; MetricType;<br>&nbsp;&nbsp;&nbsp;&nbsp;<br>&nbsp;&nbsp;&nbsp; typedef itk::TranslationTransform&lt; double, Dimension &gt;&nbsp;&nbsp;&nbsp;&nbsp; TransformType;</p>
<p>Can anybody tell me why I have so different results with different optimizers? </p>
<p>Thanks,</p>
<p>Ming</p>
<p>&nbsp;</p></div>