Hello,
I was successful in implementing the voxel normalisation procedure (As
described in this mailing list, using Utils -> Deformations ->
Composition -> Inverse), in order to find the position of a voxel in
stereotactic space given it's location in native space.
The issue is that, in some volumes, when I visualise the results,
there seems to be a systematic shift (~4 voxels) on the y axis from
where I would expect to find the normalised voxel.
I'm extracting the coordinates from native space ([x y z] in voxel
position), and using the N.dat from the y_* deformation file
(N.dat(x,y,z,:,:) = [x' y' z'] - point in standard space in mm). In
order to convert the result to voxel positions, I'm using:
M = inv(spm_get_space("mwc file"))
output = [x' y' z']*M(1:3,1:3) + M(1:3,4)
Any ideas as to what is wrong, or is it just a display issue?
Thank you for your help!
Joao
|