ParseFSSegments.m

function segdetails = parseFSSegments(fname) %%Given a filename fname, pull out the names of each segment. fname is the %%name of a file produced by mri_segstats

if (~exist('fname'))

   fname=uigetfile('*.sum.txt');

end

%open the file fid=fopen(fname, 'r'); if(fid<=0)

   error(['Specified segment file  fname  not found']);

end

%First, count number of entries %read in file, line by line, until we reach end of file linecount=0; while ~feof(fid)

   line=fgets(fid);
   if(line(1)~='#')
       %only interested in lines not starting with hash
       linecount=linecount+1;
   end

end %rewind to beginning of file now that we know how many lines of interest frewind(fid) %preallocate variables segname=cell(linecount,1); colidx=nan(linecount,1); segid=nan(linecount,1); nvertices=nan(linecount,1); segarea=nan(linecount,1); segmean=nan(linecount,1); segstd=nan(linecount,1); segmin=nan(linecount,1); segmax=nan(linecount,1); segrange=nan(linecount,1); linecount=0; %recycle linecount variable to track where to store interesting lines while ~feof(fid)

   line=fgets(fid);
   if(line(1)~='#')
       %STILL only interested in lines not starting with hash!
       linecount=linecount+1;
       A=sscanf(line, '%d %d %d %f %s %f %f %f %f %f');
       colidx(linecount)=A(1);
       segid(linecount)=A(2);
       nvertices(linecount)=A(3);
       segarea(linecount)=A(4);
       segname{linecount}=char(A(5:length(A)-5))';
       segmean(linecount)=A(length(A)-4);
       segstd(linecount)=A(length(A)-3);
       segmin(linecount)=A(length(A)-2);
       segmax(linecount)=A(length(A)-1);
       segrange(linecount)=A(length(A));
   end

end

segdetails.colidx=colidx; segdetails.segid=segid; segdetails.nvertices=nvertices; segdetails.segarea=segarea; segdetails.segname=segname; segdetails.segmean=segmean; segdetails.segstd=segstd; segdetails.segmin=segmin; segdetails.segmax=segmax; segdetails.segrange=segrange; end