Dear All,

I ran recon-all and generated a brain.mgz and say, lh.pial.

I converted brain.mgz to brain.nii, imported the latter to matlab/spm and realigned it to rbrain.nii, generating a transform matrix M=rbrain.reorient.mat.

if MM is the header for rbrain.nii, the header O for the original brain.nii can be calculated with

O=inv(M) *MM.

How do I reorient the lh.pial to the rbrain.nii? I suppose I need a reader for the pial header in matlab.

Many thanks,

Octavian