> I have two data sets which I can use SPM to register one to the other.
> Here comes the question: can I possibly convert (register) say a voxel
> of coordinate (x,y,z) at floating volume to (x',y',z') on the space of the
> reference space?
You can. How easily depends on how you registered them. If they are
just rigidly or affinely registered (e.g. with realign, or coreg,)
then the matrices in the vol.mat fields (after mapping the images with
spm_vol) relate the coordinate systems quite simply.
E.g. something like:
xyz1 = [x y z 1]'; % -- voxel (homogeneous) coords, image 1
XYZ1 = vol1.mat*[xyz1 1];
% -- world coords (in both *registered* images)
xyz2 = vol2.mat\XYZ1; % voxel coords, image 2
If you've used non-linear normalisation, or SPM5's unified
segmentation and normalisation, then the easiest way appears to be to
convert the sn.mat (or inv_sn.mat, check the first post below
carefully to be sure which you want) into a deformation field, and
then read the voxels from this, see the following:
and the parts of this which are relevant for reading from the
deformation field, after you have read the above,