[Insight-users] Problems with understanding the parameters in the amoeba optimizer

[Arne Hansen]

Hello.
I am using an amoeba optimizer to minimize overlap (itkKappaStatistic)
between 2 binary volumes.

My problem is that it seems that when i set parameters
optimizer->SetMaximumNumberOfIterations(1);
optimizer->SetParametersConvergenceTolerance(10);
optimizer->SetFunctionConvergenceTolerance(10);

they dont seem to have any affect at all. Even though i tell it to maximal
do 1 iteration, it still does 17 iterations in the first layer of my
itkMultiResolutionImageRegistrationMethod.

Furthermore, as i am doing kappastatistic, the maximum function value is of
course 1 for perfectly
gemetric aligned images. So the function value is -1 all the time(Because i
told the optimizer to maximize and not nminimize. But still, it does first
convergeafter the 17 iterationss in the top pyramid. And the same for the
rest of the pyramid.
Here is the result from a 2-layered pyramid.......


Iteration count : 0

///////////////////////////////////////////////////////////

Pyramid level0

///////////////////////////////////////////////////////////

-1:

//[]

[0, 0, 0, 0, 0, 0, 1]

-1:

[0.00025, 0, 0, 0, 0, 0, 1]

-1:

[0, 0.00025, 0, 0, 0, 0, 1]

-1:

[0, 0, 0.00025, 0, 0, 0, 1]

-1:

[0, 0, 0, 0.025, 0, 0, 1]

-1:

[0, 0, 0, 0, 0.025, 0, 1]

-1:

[0, 0, 0, 0, 0, 0.025, 1]

-1:

[0, 0, 0, 0, 0, 0, 1.05]

-1:

[7.1428571429e-005, 7.1428571429e-005, 7.1428571429e-005, 0.0071428571429,
0.0071428571429, 0.0071428571429, 1.0142857143]

-1:

[1.7857142857e-005, 1.7857142857e-005, 1.7857142857e-005, 0.0017857142857,
0.0017857142857, 0.0017857142857, 1.0035714286]

-1:

[0.000125, 0.000125, 0, 0, 0, 0, 1]

-1:

[0.000125, 0, 0.000125, 0, 0, 0, 1]

-1:

[0.000125, 0, 0, 0.0125, 0, 0, 1]

-1:

[0.000125, 0, 0, 0, 0.0125, 0, 1]

-1:

[0.000125, 0, 0, 0, 0, 0.0125, 1]

-1:

[0.000125, 0, 0, 0, 0, 0, 1.025]

-1:

[0.000125, 0, 0, 0, 0, 0, 1]

Iteration count : 17

///////////////////////////////////////////////////////////

Pyramid level1

///////////////////////////////////////////////////////////

-1:

[0.000125, 0.000125, 0, 0, 0, 0, 1]

-1:

[0.00025, 0.000125, 0, 0, 0, 0, 1]

-1:

[0.000125, 0.00025, 0, 0, 0, 0, 1]

-1:

[0.000125, 0.000125, 0.00025, 0, 0, 0, 1]

-1:

[0.000125, 0.000125, 0, 0.025, 0, 0, 1]

-1:

[0.000125, 0.000125, 0, 0, 0.025, 0, 1]

-1:

[0.000125, 0.000125, 0, 0, 0, 0.025, 1]

-0.97433615191:

[0.000125, 0.000125, 0, 0, 0, 0, 1.05]

-0.96254003621:

[0.00016071428571, 0.00016071428571, 7.1428571429e-005, 0.0071428571429,
0.0071428571429, 0.0071428571429, 0.95]

-1:

[0.00013392857143, 0.00013392857143, 1.7857142857e-005, 0.0017857142857,
0.0017857142857, 0.0017857142857, 1.025]

-1:

[0.0001875, 0.0001875, 0, 0, 0, 0, 1]

-1:

[0.0001875, 0.000125, 0.000125, 0, 0, 0, 1]

-1:

[0.0001875, 0.000125, 0, 0.0125, 0, 0, 1]

-1:

[0.0001875, 0.000125, 0, 0, 0.0125, 0, 1]

-1:

[0.0001875, 0.000125, 0, 0, 0, 0.0125, 1]

-1:

[0.0001875, 0.000125, 0, 0, 0, 0, 1]

-1:

[0.0001875, 0.000125, 0, 0, 0, 0, 1.025]

----------------------------Registration completed-------------------



Result :

Translation X = 0

Translation Y = 0

Translation Z = 0

Rotation X = 0.70710678119

Rotation Y = 0.70710678119

Rotation Z = 0

GetAngle = 0.0005303300921

Offset = -0.021839980215 0.021839980215 -0.0025979353519

Matrix =

0.99999992969 7.03125e-008 0.00037499998682

7.03125e-008 0.99999992969 -0.00037499998682

-0.00037499998682 0.00037499998682 0.99999985938

Offset =

[-0.021839980215, 0.021839980215, -0.0025979353519]

Culd someone help my. Than you very much
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20061129/ba727014/attachment.html