So I am doing some matlab based analysis where I have to align several fmri volumes.

To load the data from Freesurfer I used f = MRIread('filename.nii') where I have f.vol is the voxel intensity series.

How is it possible to extract the [x,y,z] coordinates of the voxels from that structure to match the registered volume space?

Thank you,

Best,

Mohamed =)