[Insight-users] slice to volume registration: optimization is correct?
stan
stan.oubel at gmail.com
Fri Jan 8 10:18:18 EST 2010
Hi all,
I'm doing some tests on slice to volume affine registration.
To test my code, I register a specific slice of the
volume (fixed image) with the whole volume (moving image) by
providing an initial transformation (an affine transformation
which only rotates the slice 10 deg around the x-axis).
Even when the result is fine, I've found that some parameters
are not optimized at all ... (see the 6th and the 9th parameters):
Initial parameters = [1, 0, 0, 0, 0.984808, 0.173649, 0, -0.173649,
0.984808, 0, 0, 0]
0 0.2 -0.557382 [1.00186, -0.000366774, 0, 0.00025731, 0.986009,
0.173649, -0.000127247, -0.167419, 0.984808, -0.0105518, 0.196942,
-0.032534]
1 0.2 -0.577705 [0.998905, 0.000928042, 0, 0.000672959, 0.986261,
0.173649, -0.000406448, -0.159618, 0.984808, -0.0102843, 0.391734,
-0.077084]
2 0.2 -0.595902 [1.0119, 0.000218071, 0, 0.00082287, 0.989483,
0.173649, -8.3977e-05, -0.15157, 0.984808, 0.00914714, 0.589726,
-0.0903846]
...
77 0.000195313 -1.64248 [0.999986, -2.18366e-05, 0, -2.51637e-06,
0.999994, 0.173649, -4.839e-05, -3.68753e-05, 0.984808, 0.00165371,
0.00172779, 0.00317447]
Final parameters = [0.999986, -2.18366e-05, 0, -2.51637e-06, 0.999994,
0.173649, -4.839e-05, -3.68753e-05, 0.984808, 0.00165371, 0.00172779,
0.00317447]
I wonder if this is an expected possible behaviour of the optimizer, and
if so, how to explain this ... I would have expected to obtain a result
equal to the identity, as I obtain when the initial rotation is around
the z-axis ...
I've attached my code (I've removed the irrelevant parts of it).
I run it with rx=-10*3.1416/180 and ry=rz=tx=ty=tz=0.
Thanks !
Stan
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CodeForITK.cxx
Type: text/x-c++src
Size: 6306 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20100108/da0ba0b4/attachment.cxx>
More information about the Insight-users
mailing list