mesh_write_fs_curv

 mesh_write_fs_curv - Save mesh CData to FreeSurfer curvature file

 USEAGE: mesh_write_fs_curv(p)
 
 Write a binary curvature file for each mesh in p.mesh.data.  
 If any cells of p.mesh.data.meshtype are empty, they are skipped.
 
 This function outputs an *.curv curvature file from the
 vertex scalar data contained in p.mesh.data.Cdata - if it 
 contains a timeseries, the output values are taken from 
 the column specified in p.mesh.samplePoint.
 
 All output is handled by fs_write_curv.
 
 See the FreeSurfer website at http://surfer.nmr.mgh.harvard.edu/
 for more information.

0001 function mesh_write_fs_curv(p)
0002 
0003 % mesh_write_fs_curv - Save mesh CData to FreeSurfer curvature file
0004 %
0005 % USEAGE: mesh_write_fs_curv(p)
0006 %
0007 % Write a binary curvature file for each mesh in p.mesh.data.
0008 % If any cells of p.mesh.data.meshtype are empty, they are skipped.
0009 %
0010 % This function outputs an *.curv curvature file from the
0011 % vertex scalar data contained in p.mesh.data.Cdata - if it
0012 % contains a timeseries, the output values are taken from
0013 % the column specified in p.mesh.samplePoint.
0014 %
0015 % All output is handled by fs_write_curv.
0016 %
0017 % See the FreeSurfer website at http://surfer.nmr.mgh.harvard.edu/
0018 % for more information.
0019 %
0020 
0021 % $Revision: 1.1 $ $Date: 2004/11/12 01:32:35 $
0022 
0023 % Licence:  GNU GPL, no implied or express warranties
0024 % History:  10/2002, Darren.Weber_at_radiology.ucsf.edu
0025 %           01/2003, Darren.Weber_at_radiology.ucsf.edu
0026 %                    now using fs_write_surf and fs_write_wfile
0027 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0028 
0029 fprintf('\nMESH_WRITE_FS_CURV...\n');
0030 
0031 if ~exist('p','var'),
0032     error('...no input p struct.\n');
0033 elseif isempty(p),
0034     error('...input p struct is empty.\n');
0035 elseif isempty(p.mesh.data),
0036     error('...input p struct has no mesh data.\n');
0037 end
0038 
0039 [path,name,ext] = fileparts(strcat(p.mesh.path,filesep,p.mesh.file));
0040 file = fullfile(path,[name ext]);
0041 
0042 fprintf('...writing FreeSurfer curvatures to:\n\t%s\n',path);
0043 
0044 tic;
0045 
0046 Meshes = p.mesh.data.meshtype;
0047 
0048 [path,name,ext] = fileparts(file);
0049 
0050 for i=1:size(Meshes,2),
0051     
0052     if Meshes{i},
0053         
0054         % check for mesh Cdata
0055         if size(p.mesh.data.Cdata,2) >= i,
0056             if ~isempty(p.mesh.data.Cdata{1}),
0057                 if size(p.mesh.data.Cdata{i},2) > 1,
0058                     % Obtain the Cdata at the selected time point
0059                     curv = p.mesh.data.Cdata{i}(:,p.mesh.samplePoint);
0060                 else
0061                     curv = p.mesh.data.Cdata{i};
0062                 end
0063                 
0064                 Nfaces = size(p.mesh.data.faces{i},1);
0065                 
0066                 outputExt = [ext,'.curv'];
0067                 outputName = strcat(name,'.',Meshes{i});
0068                 outputFile = fullfile(path,[outputName outputExt]);
0069                 
0070                 fprintf('...writing curvature: %s\n',[outputName outputExt]);
0071                 fs_write_curv(outputFile,curv,Nfaces);
0072             else
0073                 msg = sprintf('\np.mesh.data.Cdata{%d} is empty\n',i);
0074                 warning(msg);
0075             end
0076         else
0077             msg = sprintf('\nno p.mesh.data.Cdata{%d}\n',i);
0078             warning(msg);
0079         end
0080     end
0081 end
0082 
0083 t=toc; fprintf('...done (%5.2f sec).\n\n',t);
0084 
0085 return