Dear Freesurfer community,
I have processed a subject (name:P061_pre, 3D MPRAGE T1) with Freesurfer
5.1. A resection area has been defiened on this volume and exists as a
binary volume (1:voxels of interest, 0:others) lesionMask.mgz. I now
would like to project this volume on the surface of the subject, e.g.
lh.white.
I have tried using mri_vol2surf :
mri_vol2surf --src lesionMask.mgz --out lesionMask_surf.mgh --out_type
mgh --hemi lh --regheader P061_pre
The command seems to work:
srcvol = lesionMask.mgz
srcreg unspecified
srcregold = 0
srcwarp unspecified
surf = white
hemi = lh
reshape = 0
interp = nearest
float2int = round
GetProjMax = 0
INFO: float2int code = 0
INFO: changing type to float
Done loading volume
Computing registration from header.
Using /root/data/FS51/freesurfer/subjects/P061_pre/mri/orig.mgz as
target reference.
Reading surface /root/data/FS51/freesurfer/subjects/P061_pre/surf/lh.white
Done reading source surface
Mapping Source Volume onto Source Subject Surface
1 0 0 0
using old
Done mapping volume to surface
Number of source voxels hit = 79054
Writing to ./lesionMask_surf.mgh
Dim: 128233 1 1
But when trying to open the surface with freeview, it returns an error
message :
freeview -f lesionMask_surf.mgh
ERROR: MRISread: cannot read surface data from file
/home/epilepsy/sorted/P061/preMR/lesionMask_surf.mgh!
Segmentation fault
What did I miss ?
Thanks for your help,
--
Christopher Coello, Ph.D.
Preclinical PET/CT unit
Institute of Basic Medical Sciences
Faculty of Medicine, University of Oslo, Norway
tel. +4745154516