Hi, I have my own stats files that I would like to incorporate into the freesurfer stream, if possible. Ie. after running FreeSurfer, I create a ?h.mystats (containing something like image intensities for example), where the format of ?h.mystats is just an ascii file, where each line is the value at that vertex.
What I ultimately want to do, is to transform ?h.mystats into fsaverage space, that way I can compare ?h.mystats across subjects. I was hoping I could use the given sphere.reg which has already been calculated to do this.
I renamed ?h.mystats to ?h.thickness and tried the following which I think will work: mris_preproc --s mysubject --hemi lh --meas thickness --target fsaverage --out lh.thickness.fsaverage.mgh
But it complains that:
Reading curvature file /lab1/FreeSurfer/4.0.5/ms-mni/testing//tester/surf/lh.thickness MRISreadBinaryCurvature: incompatible vertex number in file /lab1/FreeSurfer/4.0.5/ms-mni/testing//tester/surf/lh.thickness
I checked and I have the same number of lines in my file as vertices on my object so my questions are:
1. Do I have to convert my ascii ?h.mystats to a curvature (binary) file like ?h.thickness, and if so how does one do this? 2. Is this the right way to get what I ultimately want to do, or is there another method?
thanks, mishkin
Mishkin,
Try appending .asc to your file (ie, lh.mystats.asc) and try again. freesurfer interprets .asc as an ascii file (otherwise it tries to read it as binary).
Nick
On Thu, 2008-10-09 at 17:23 -0400, Mishkin Derakhshan wrote:
Hi, I have my own stats files that I would like to incorporate into the freesurfer stream, if possible. Ie. after running FreeSurfer, I create a ?h.mystats (containing something like image intensities for example), where the format of ?h.mystats is just an ascii file, where each line is the value at that vertex.
What I ultimately want to do, is to transform ?h.mystats into fsaverage space, that way I can compare ?h.mystats across subjects. I was hoping I could use the given sphere.reg which has already been calculated to do this.
I renamed ?h.mystats to ?h.thickness and tried the following which I think will work: mris_preproc --s mysubject --hemi lh --meas thickness --target fsaverage --out lh.thickness.fsaverage.mgh
But it complains that:
Reading curvature file /lab1/FreeSurfer/4.0.5/ms-mni/testing//tester/surf/lh.thickness MRISreadBinaryCurvature: incompatible vertex number in file /lab1/FreeSurfer/4.0.5/ms-mni/testing//tester/surf/lh.thickness
I checked and I have the same number of lines in my file as vertices on my object so my questions are:
- Do I have to convert my ascii ?h.mystats to a curvature (binary)
file like ?h.thickness, and if so how does one do this? 2. Is this the right way to get what I ultimately want to do, or is there another method?
thanks, mishkin _______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu