\n','plotly',650,800); fprintf(fid,' \n'); fprintf(fid,' \n','plotly'); fprintf(fid,' \n*/\n'); fclose(fid); else % JSON if exist('jsonencode','builtin') fid = fopen(filename,'wt'); if fid == -1 error('Unable to write file %s: permission denied.',filename); end fprintf(fid,'%s',jsonencode(struct('data',{data},'layout',struct()))); fclose(fid); else spm_jsonwrite(filename,struct('data',{data},'layout',struct()),struct('indent',' ')); end end %========================================================================== % function save_dae(s,filename) %========================================================================== function save_dae(s,filename) o = @(x) blanks(x*3); % Split the mesh into connected components %-------------------------------------------------------------------------- try C = spm_mesh_label(s.faces); d = []; for i=1:numel(unique(C)) d(i).faces = s.faces(C==i,:); u = unique(d(i).faces); d(i).vertices = s.vertices(u,:); a = 1:max(d(i).faces(:)); a(u) = 1:size(d(i).vertices,1); %a = sparse(1,double(u),1:1:size(d(i).vertices,1)); d(i).faces = a(d(i).faces); end s = d; end % Open file for writing %-------------------------------------------------------------------------- fid = fopen(filename,'wt'); if fid == -1 error('Unable to write file %s: permission denied.',filename); end % Prolog & root of the Collada XML file %-------------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'\n'); % Assets %-------------------------------------------------------------------------- fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(2)); fprintf(fid,'%s%s\n',o(3),... 'http://www.fil.ion.ucl.ac.uk/spm/'); fprintf(fid,'%s%s\n',o(3),'SPM'); fprintf(fid,'%s\n',o(2)); fprintf(fid,'%s%s\n',o(2),datestr(now,'yyyy-mm-ddTHH:MM:SSZ')); fprintf(fid,'%s%s\n',o(2),datestr(now,'yyyy-mm-ddTHH:MM:SSZ')); fprintf(fid,'%s\n',o(2)); fprintf(fid,'%sZ_UP\n',o(2)); fprintf(fid,'%s\n',o(1)); % Image, Materials, Effects %-------------------------------------------------------------------------- %fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(1)); for i=1:numel(s) fprintf(fid,'%s\n',o(2),i,i); fprintf(fid,'%s\n',o(3),i); fprintf(fid,'%s\n',o(2)); end fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(1)); for i=1:numel(s) fprintf(fid,'%s\n',o(2),i,i); fprintf(fid,'%s\n',o(3)); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s%f %f %f %d\n',o(7),[0 0 0 1]); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s%f %f %f %d\n',o(7),[0 0 0 1]); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s%f %f %f %d\n',o(7),[0.5 0.5 0.5 1]); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s%d %d %d %d\n',o(7),[1 1 1 1]); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s%f\n',o(7),0); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(3)); fprintf(fid,'%s\n',o(2)); end fprintf(fid,'%s\n',o(1)); % Geometry %-------------------------------------------------------------------------- fprintf(fid,'%s\n',o(1)); for i=1:numel(s) fprintf(fid,'%s\n',o(2),i,i); fprintf(fid,'%s\n',o(3)); fprintf(fid,'%s\n',o(4),i); fprintf(fid,'%s',o(5),i,numel(s(i).vertices)); fprintf(fid,'%f ',reshape(s(i).vertices',1,[])); fprintf(fid,'\n'); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(6),size(s(i).vertices,1),i); fprintf(fid,'%s\n',o(7)); fprintf(fid,'%s\n',o(7)); fprintf(fid,'%s\n',o(7)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(4),i); fprintf(fid,'%s\n',o(5),i); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(4),i,size(s(i).faces,1)); fprintf(fid,'%s\n',o(5),i); fprintf(fid,'%s

',o(5)); fprintf(fid,'%d ',reshape(s(i).faces',1,[])-1); fprintf(fid,'

\n'); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(3)); fprintf(fid,'%s\n',o(2)); end fprintf(fid,'%s\n',o(1)); % Scene %-------------------------------------------------------------------------- fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(2)); for i=1:numel(s) fprintf(fid,'%s\n',o(3),i); fprintf(fid,'%s\n',o(4),i); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(7),i,i); fprintf(fid,'%s\n',o(6)); fprintf(fid,'%s\n',o(5)); fprintf(fid,'%s\n',o(4)); fprintf(fid,'%s\n',o(3)); end fprintf(fid,'%s\n',o(2)); fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(1)); fprintf(fid,'%s\n',o(2)); fprintf(fid,'%s\n',o(1)); % End of XML %-------------------------------------------------------------------------- fprintf(fid,'\n'); % Close file %-------------------------------------------------------------------------- fclose(fid); %========================================================================== % function save_idtf(s,filename) %========================================================================== function save_idtf(s,filename) o = @(x) blanks(x*3); % Compute normals %-------------------------------------------------------------------------- if ~isfield(s,'normals') try s.normals = spm_mesh_normals(... struct('vertices',s.vertices,'faces',s.faces),true); catch s.normals = []; end end % Split the mesh into connected components %-------------------------------------------------------------------------- try C = spm_mesh_label(s.faces); d = []; try if size(s.cdata,2) == 1 && (any(s.cdata>1) || any(s.cdata<0)) mi = min(s.cdata); ma = max(s.cdata); s.cdata = (s.cdata-mi)/ (ma-mi); else end end for i=1:numel(unique(C)) d(i).faces = s.faces(C==i,:); u = unique(d(i).faces); d(i).vertices = s.vertices(u,:); d(i).normals = s.normals(u,:); a = 1:max(d(i).faces(:)); a(u) = 1:size(d(i).vertices,1); %a = sparse(1,double(u),1:1:size(d(i).vertices,1)); d(i).faces = a(d(i).faces); d(i).mat = s.mat; try d(i).cdata = s.cdata(u,:); if size(d(i).cdata,2) == 1 d(i).cdata = repmat(d(i).cdata,1,3); end end end s = d; end % Open file for writing %-------------------------------------------------------------------------- fid = fopen(filename,'wt'); if fid == -1 error('Unable to write file %s: permission denied.',filename); end % FILE_HEADER %-------------------------------------------------------------------------- fprintf(fid,'FILE_FORMAT "IDTF"\n'); fprintf(fid,'FORMAT_VERSION 100\n\n'); % NODES %-------------------------------------------------------------------------- for i=1:numel(s) fprintf(fid,'NODE "MODEL" {\n'); fprintf(fid,'%sNODE_NAME "%s"\n',o(1),sprintf('Mesh%04d',i)); fprintf(fid,'%sPARENT_LIST {\n',o(1)); fprintf(fid,'%sPARENT_COUNT %d\n',o(2),1); fprintf(fid,'%sPARENT %d {\n',o(2),0); fprintf(fid,'%sPARENT_NAME "%s"\n',o(3),''); fprintf(fid,'%sPARENT_TM {\n',o(3)); I = s(i).mat; % eye(4); for j=1:size(I,2) fprintf(fid,'%s',o(4)); fprintf(fid,'%f ',I(:,j)'); fprintf(fid,'\n'); end fprintf(fid,'%s}\n',o(3)); fprintf(fid,'%s}\n',o(2)); fprintf(fid,'%s}\n',o(1)); fprintf(fid,'%sRESOURCE_NAME "%s"\n',o(1),sprintf('Mesh%04d',i)); %fprintf(fid,'%sMODEL_VISIBILITY "BOTH"\n',o(1)); fprintf(fid,'}\n\n'); end % NODE_RESOURCES %-------------------------------------------------------------------------- for i=1:numel(s) fprintf(fid,'RESOURCE_LIST "MODEL" {\n'); fprintf(fid,'%sRESOURCE_COUNT %d\n',o(1),1); fprintf(fid,'%sRESOURCE %d {\n',o(1),0); fprintf(fid,'%sRESOURCE_NAME "%s"\n',o(2),sprintf('Mesh%04d',i)); fprintf(fid,'%sMODEL_TYPE "MESH"\n',o(2)); fprintf(fid,'%sMESH {\n',o(2)); fprintf(fid,'%sFACE_COUNT %d\n',o(3),size(s(i).faces,1)); fprintf(fid,'%sMODEL_POSITION_COUNT %d\n',o(3),size(s(i).vertices,1)); fprintf(fid,'%sMODEL_NORMAL_COUNT %d\n',o(3),size(s(i).normals,1)); if ~isfield(s(i),'cdata') || isempty(s(i).cdata) c = 0; else c = size(s(i).cdata,1); end fprintf(fid,'%sMODEL_DIFFUSE_COLOR_COUNT %d\n',o(3),c); fprintf(fid,'%sMODEL_SPECULAR_COLOR_COUNT %d\n',o(3),0); fprintf(fid,'%sMODEL_TEXTURE_COORD_COUNT %d\n',o(3),0); fprintf(fid,'%sMODEL_BONE_COUNT %d\n',o(3),0); fprintf(fid,'%sMODEL_SHADING_COUNT %d\n',o(3),1); fprintf(fid,'%sMODEL_SHADING_DESCRIPTION_LIST {\n',o(3)); fprintf(fid,'%sSHADING_DESCRIPTION %d {\n',o(4),0); fprintf(fid,'%sTEXTURE_LAYER_COUNT %d\n',o(5),0); fprintf(fid,'%sSHADER_ID %d\n',o(5),0); fprintf(fid,'%s}\n',o(4)); fprintf(fid,'%s}\n',o(3)); fprintf(fid,'%sMESH_FACE_POSITION_LIST {\n',o(3)); fprintf(fid,'%d %d %d\n',s(i).faces'-1); fprintf(fid,'%s}\n',o(3)); fprintf(fid,'%sMESH_FACE_NORMAL_LIST {\n',o(3)); fprintf(fid,'%d %d %d\n',s(i).faces'-1); fprintf(fid,'%s}\n',o(3)); fprintf(fid,'%sMESH_FACE_SHADING_LIST {\n',o(3)); fprintf(fid,'%d\n',zeros(size(s(i).faces,1),1)); fprintf(fid,'%s}\n',o(3)); if c fprintf(fid,'%sMESH_FACE_DIFFUSE_COLOR_LIST {\n',o(3)); fprintf(fid,'%d %d %d\n',s(i).faces'-1); fprintf(fid,'%s}\n',o(3)); end fprintf(fid,'%sMODEL_POSITION_LIST {\n',o(3)); fprintf(fid,'%f %f %f\n',s(i).vertices'); fprintf(fid,'%s}\n',o(3)); fprintf(fid,'%sMODEL_NORMAL_LIST {\n',o(3)); fprintf(fid,'%f %f %f\n',s(i).normals'); fprintf(fid,'%s}\n',o(3)); if c fprintf(fid,'%sMODEL_DIFFUSE_COLOR_LIST {\n',o(3)); fprintf(fid,'%f %f %f\n',s(i).cdata'); fprintf(fid,'%s}\n',o(3)); end fprintf(fid,'%s}\n',o(2)); fprintf(fid,'%s}\n',o(1)); fprintf(fid,'}\n'); end % Close file %-------------------------------------------------------------------------- fclose(fid); %========================================================================== % function save_obj(s,filename) %========================================================================== function save_obj(s,filename) % Open file for writing %-------------------------------------------------------------------------- fid = fopen(filename,'wt'); if fid == -1 error('Unable to write file %s: permission denied.',filename); end % Vertices & faces %-------------------------------------------------------------------------- fprintf(fid,'# Wavefront OBJ file saved by %s\n','SPM'); % spm('Version') fprintf(fid,'v %f %f %f\n',s.vertices'); fprintf(fid,'f %d %d %d\n',s.faces'); % Close file %-------------------------------------------------------------------------- fclose(fid); %========================================================================== % function write_freesurfer_file(s,filename,format) %========================================================================== function write_freesurfer_file(s,filename,format) switch lower(format) case 'surf' fid = fopen(filename,'wb','ieee-be'); if fid == -1, error('Cannot open "%s".',filename); end fwrite(fid,[255 255 254],'uchar'); fprintf(fid,'created by @gifti on %s\n\n',... datestr(now,'ddd mmm dd HH:MM:SS yyyy')); fwrite(fid,size(s.vertices,1),'int32'); % nv fwrite(fid,size(s.faces,1),'int32'); % nf fwrite(fid,s.vertices','float32'); fwrite(fid,s.faces'-1,'int32'); fclose(fid); case 'curv' fid = fopen(filename,'wb','ieee-be'); if fid == -1, error('Cannot open "%s".',filename); end fwrite(fid,[255 255 255],'uchar'); fwrite(fid,size(s.cdata,1),'int32'); % nv if isfield(s,'faces') nf = size(s.faces,1); else nf = 0; end fwrite(fid,nf,'int32'); % nf fwrite(fid,size(s.cdata,2),'int32'); % n fwrite(fid,s.cdata','float32'); fclose(fid); otherwise error('File format not supported.'); end