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