I think this should work with FSL data in 2x2x2mm standard space, but someone correct me if I'm wrong. I tested this on a few voxels in MRIcron and it seems to work fine...
function voxelspace = mni2voxel(inputvoxel) %inputvoxel should be [x y z] MNI coordinates
%my assumptions:
%voxel size == 2mm iso
%datasize=(91,109,91)
%mniorigin=[46 64 37] %identified as MNI coordinate [0 0 0] in MRIcron for fsl standard brain
mniorigin=[46 64 37];
yourvoxel=inputvoxel/2;
voxelspace = yourvoxel + mniorigin;
Obviously this could be made more flexible by loosening the hard coded assumptions...
David