sure, if you save it in the right format. The easiest thing to do would
be to save it as a volume, where the 1st spatial dimension is # of
vertices on the surface, and the last (4th one) is the # of time points
(the 2nd and 3rd would both be 1). Or there are other formats you can save
in if you like, include our "w-file" format. I've attached a matlab file
for writing it.
cheers,
Bruce
On Fri, 2 Feb 2007, T. Song wrote:
Hi Bruce:
Is there a way that I can load the MEG timecourse data on the freesurfer
surface? Is it possible that we could get the timecourse waveform for a
specific point on the surface?
Thanks!
Tao
function err = write_wfile(fname, w, v)
% err = write_wfile(fname, w, <v>)
%
% writes a vector into a binary 'w' file
% fname - name of file to write to
% w - vector of values to be written
% v - 0-based vertex numbers
% (assumes 0 to N-1 if not present or empty).
%
% See also read_wfile.
%
%
% write_wfile.m
%
% Original Author: Doug Greve
% CVS Revision Info:
% $Author: nicks $
% $Date: 2007/01/10 22:55:10 $
% $Revision: 1.2 $
%
% Copyright (C) 2002-2007,
% The General Hospital Corporation (Boston, MA).
% All rights reserved.
%
% Distribution, usage and copying of this software is covered under the
% terms found in the License Agreement file named 'COPYING' found in the
% FreeSurfer source code root directory, and duplicated here:
% https://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferOpenSourceLicense
%
% General inquiries: [email protected]
% Bug reports: [EMAIL PROTECTED]
%
err = 1;
if(nargin ~= 2 & nargin ~= 3)
fprintf('USAGE: err = write_wfile(fname, w, <v>) \n');
return;
end
vnum = length(w) ;
% Handle when v not given or is empty %
if(exist('v') ~= 1) v = []; end
if(isempty(v)) v = [0:vnum-1]; end
% open it as a big-endian file
fid = fopen(fname, 'wb', 'b') ;
if(fid == -1)
fprintf('ERROR: could not open %s\n',fname);
return;
end
fwrite(fid, 0, 'int16') ;
fwrite3(fid, vnum) ;
for i=1:vnum
fwrite3(fid, v(i)) ; % vertex number (0-based)
fwrite(fid, w(i), 'float') ; % vertex value
end
fclose(fid) ;
err = 0;
return
_______________________________________________
Freesurfer mailing list
[email protected]
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
