Does anyone have code available for reading the vox2ras-tkr transform from .mgz files? I am working on implementing the transform outlined on the Freesurfer wiki:
MNI305RAS = TalXFM*Norig*inv(Torig)*[tkrR tkrA tkrS 1]'
TalXFM: subject/orig/transforms/talairach.xfm Norig: mri_info --vox2ras orig.mgz Torig: mri_info --vox2ras-tkr orig.mgz
I can get vox2ras-tkr (a.k.a., Torig) using mri_info, but it would be better if they were available in the python package nibabel for use in mne-python (for easy conversion from surface to MNI coordinates). nibabel already supports reading the --vox2ras transform, but I need to code the other one in order to avoid a dependecy in our code. Does anyone have code to extract that information? MATLAB, Python, C code, or heck, even pseudocode would help. If not, if anyone has helpful information on the specs on the .mgz specification, maybe it won't be too difficult to code from scratch...
Eric