Hi,
I am having several problems with the feat2surf script when using this on my own data. For privacy purposes I renamed my subject to bert, however this refers not to the sample dataset:
The problems that I face with the feat2surf script are as follows:
1) If I do not edit the script manually it stops when not finding a data file (zstat1.img) in analyze format. If I edit the feat2surf script so it only finds a zstat1.nii file it continues.
2) Results are not stored under 4 new subdirectories in the .feat main directory (as suggested in the tutorial). They are stored (and overwritten) in the /root directory. (I know it’s a bad habit of mine to work under the root account, but I expect the outputdir will only change to /home/$username/ if I change user. Ẃhen the script only looks for .nii files (and finds them) it exits nicely (after commenting away some of the images the feat2surf script looks for such as cluster* and rendered*).
3) When opening for example the rh volume in tksurfer, I get an error message that the .mgh volume is incorrect.
This is the output when I run the feat2surf script without editing it.
Mon Jan 30 13:32:00 CET 2006
/media/userdata/experiments/stereo_mtl/pilot/bert/results/fwhm5_hp50/1st_lev el
mri_vol2surf --src mblok1+.feat/stats/*.img --srcreg mblok1+.feat/reg/freesurfer/anat2exf.register.dat --trgsubject bert2 --hemi lh --out /*.mgh --surf white --out_type mgh
srcvol = mblok1+.feat/stats/*.img
srcreg = mblok1+.feat/reg/freesurfer/anat2exf.register.dat
srcregold = 0
srcwarp unspecified
surf = white
hemi = lh
trgsubject = bert2
surfreg = sphere.reg
interp = nearest
float2int = round
INFO: float2int code = 0
ERROR: analyzeRead(): cannot find any files for mblok1+.feat/stats/*.img
ERROR: could not read mblok1+.feat/stats/*.img as type 16
ERROR: with mri_vol2surf --src mblok1+.feat/stats/*.img --srcreg mblok1+.feat/reg/freesurfer/anat2exf.register.dat --trgsubject bert2 --hemi lh --out /*.mgh --surf white --out_type mgh
When I ran mri_vol2surf manually (see command line below; $out1 = the output dir I created manually to store the results in), the computer runs forever displaying a vol2surf_linear message. I left it running overnight and it didn't finish. Any suggestions in what I could alter to get it right? The functional overlay works fine with tkmedit.
Thanks in advance,
Niels van Strien
mri_vol2surf --src $datalocation/mblok1+.feat/stats/zstat1.nii --srcreg $datalocation/mblok1+.feat/reg/freesurfer/anat2exf.register.dat --trgsubject bert2 --hemi lh --out $out1/zstat1.nii.mgh --surf white --out_type mg
srcvol = /media/userdata/experiments/stereo_mtl/pilot/bert/results/fwhm5_hp50/1st_lev el/mblok1+.feat/stats/zstat1.nii
srcreg = /media/userdata/experiments/stereo_mtl/pilot/bert/results/fwhm5_hp50/1st_lev el/mblok1+.feat/reg/freesurfer/anat2exf.register.dat
srcregold = 0
srcwarp unspecified
surf = white
hemi = lh
trgsubject = bert2
surfreg = sphere.reg
interp = nearest
float2int = round
INFO: float2int code = 0
WARNING: missing NIfTI-1 orientation (qform_code = 0)
WARNING: your volume will probably be incorrectly oriented
Done loading volume
Reading surface /media/userdata/subjects/bert2/surf/lh.white
INFO: NOT fixing vertex area
Done reading source surface
Mapping Source Volume onto Source Subject Surface
1 0 0 0
vol2surf_linear: nhits = 74277/126030
2 0 0 0
vol2surf_linear: nhits = 74277/126030
3 0 0 0
vol2surf_linear: nhits = 74277/126030
4 0 0 0
vol2surf_linear: nhits = 74277/126030
5 0 0 0
vol2surf_linear: nhits = 74277/126030
6 0 0 0
vol2surf_linear: nhits = 74277/126030
7 0 0 0
vol2surf_linear: nhits = 74277/126030
8 0 0 0