Sorry, my previous post had a mistake...
I actually then this for the first step:
1) obtain the transformation matrix:
pointflirt -i involcoords -r refvolcoords -o pointflirtoutputmatrix_4pt
and then applied step 2..
After 1), the transformation matrix (pointflirtoutputmatrix_4pt) was:
0.9694 -0.1187 0.0639 -8.1211
0.1035 0.9663 0.2204 38.2453
-0.0483 -0.4508 0.5423 -44.1636
0 0 0 1.0000
My involcoords, were:
-2.5000 -12.5000 15.5000 -3.5000
-76.5000 -17.5000 -20.5000 -54.5000
19.0000 69.0000 69.0000 75.0000
I padded 1's to the 4th row to make it 4x4 and named it involcoords_wzeropad.
When I multiply the two matrices, i.e. when I apply the transformation to the input coordinates, I obtain:
(pointflirtoutputmatrix_4pt*involcoords_wzeropad in Matlab)
-0.2500 -13.7500 13.7500 -0.2500
-31.7500 35.2500 35.2500 1.7500
0.7499 1.7500 1.7500 21.2500
1.0000 1.0000 1.0000 1.0000
The first 3 rows of this is exactly the same as my refvolcoords, which tells me that the pointflirt output, the transfomation matrix, is correct. The pointflirt worked perfectly. But, it is just that I am not able to apply it to my input volume with "flirt -applyxfm" as I wrote in my previous post. Is there another way to apply this transformation?
..
Unal
Unal
|