fs_read_curv.m

Matlab下的EEG处理程序库

function [curv, fnum] = fs_read_curv(fname)

% FS_READ_CURV - FREESURFER I/O function to read a curvature file
%
% [curv, fnum] = fs_read_curv(fname)
% 
% reads a binary curvature file into a vector
%
% After reading an associated surface, with fs_read_surf, try:
% patch('vertices',vert,'faces',face,...
%       'facecolor','interp','edgecolor','none',...
%       'FaceVertexCData',curv); light
% 
% See also FS_WRITE_CURV, FS_READ_SURF, FS_READ_WFILE


% open it as a big-endian file
fid = fopen(fname, 'rb', 'b') ;
if (fid < 0),
    str = sprintf('could not open curvature file %s.', fname) ;
    error(str) ;
end

vnum = fs_fread3(fid);
NEW_VERSION_MAGIC_NUMBER = 16777215;
if (vnum == NEW_VERSION_MAGIC_NUMBER),
    fprintf('...reading new version (float)\n');
    vnum = fread(fid, 1, 'int32') ;
    fnum = fread(fid, 1, 'int32') ;
    vals_per_vertex = fread(fid, 1, 'int32') ;
    curv = fread(fid, vnum, 'float') ; 
else
    fprintf('...reading old version (int16)\n');
    fnum = fs_fread3(fid) ;
    curv = fread(fid, vnum, 'int16') ./ 100 ; 
end

fclose(fid) ;

return