In matlab you can use the load_mgh.m save_mgh.m routines found in FREESURFER_HOME/matlab.

[T1, M, mr_parms] = load_mgh('t1.mgh'); R1 = 1./T1; save_mgh(R1,'r1.mgh', M, mr_parms);

M and mr_params are header info that should stay with the data.

doug

Irina Sigalovsky wrote:

Hi,

I was wondering if there is a way to invert *.mgh files (i.e., for each voxel, newvalue = 1/value). For example, is there a way to load *.mgh files into matlab?

I need this to calculate R1 from T1.mgh.

Thanks, Irina _______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer