Hi again,
OK. I managed to reverse engineer things back from mri_tessalate.c Here is the code snippet that should be added to read_surf to parse these new files: NEW_QUAD_FILE_MAGIC_NUMBER =16777213;
elseif (magic == NEW_QUAD_FILE_MAGIC_NUMBER) vnum = fread3(fid); fnum = fread3(fid); vertex_coords = fread(fid, vnum*3, 'float32') ; tmp = fread(fid, fnum*3*3, 'uchar') ; tmp = reshape(tmp, 3, 3, fnum) ; faces = reshape(bitshift(tmp(1,:,:), 16) + bitshift(tmp(2,:,:),8) + tmp(3,:,:),3, fnum)'; end
-- Shay
Sorry, I guess it got out of date. We don't it much Bruce
On May 29, 2012, at 9:12 PM, Shay Ohayon shay.ohayon@gmail.com wrote:
Hi again,
OK. I managed to reverse engineer things back from mri_tessalate.c Here is the code snippet that should be added to read_surf to parse these new files: NEW_QUAD_FILE_MAGIC_NUMBER =16777213;
elseif (magic == NEW_QUAD_FILE_MAGIC_NUMBER) vnum = fread3(fid); fnum = fread3(fid); vertex_coords = fread(fid, vnum*3, 'float32') ; tmp = fread(fid, fnum*3*3, 'uchar') ; tmp = reshape(tmp, 3, 3, fnum) ; faces = reshape(bitshift(tmp(1,:,:), 16) + bitshift(tmp(2,:,:),8) + tmp(3,:,:),3, fnum)'; end
-- Shay
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu