GetSegCoords.m

The current code for this custom MATLAB function is provided below for convenient dissemination.

GetSegCoords.m

function [xyz, hemi, segnames] = getSegCoords(varargin)
% Given a set of label file names, parse the label files and return the mean XYZ coordinate of each label
% Can be called with a list of parameters for easy scripting, or interactively without parameters, which will 
% generate a series of ui dialog boxes.
%
% returns:
% xyz, a nx3 array of XYZ coordinates, where n is the number of labels parsed
% hemi, a nx1 binary vector indicating the corresponding label is left (1) or right (0) hemisphere
% segnames: a nx1 cell array of label strings.
%
% (C) 2018 Chris McNorgan
%  Use at your own risk.
%
 
%% Part 1: Parse arguments
options = struct( ...
   'lhsegdir',, ...
   'rhsegdir',, ...
   'lsegnames', , ...
   'rsegnames', , ...
   'automatch', true, ...
   'hemi', 'lhrh', ...
   'directory', pwd() );

% read the acceptable names
optionNames = fieldnames(options);

% count arguments
nArgs = length(varargin);
if round(nArgs/2)~=nArgs/2
   error('TSExtractor needs propertyName/propertyValue pairs')
end

for pair = reshape(varargin,2,[]) % pair is {propName;propValue}
   inpName = lower(pair{1}); % make case insensitive
   
   if any(strcmp(inpName,optionNames))
       % overwrite options. If you want you can test for the right class here
       % Also, if you find out that there is an option you keep getting wrong,
       % you can use "if strcmp(inpName,'problemOption'),testMore,end"-statements
       options.(inpName) = pair{2};
   else
       error('%s is not a recognized parameter name',inpName)
   end
end

lh=~isempty(strfind(options.hemi, 'lh'));
rh=~isempty(strfind(options.hemi, 'rh'));
 
%% Part 2: Determine folders
%If we want lh coords and we haven't already provided a directory, open
%dialog box
if (lh && isempty(options.lhsegdir))
   options.lhsegdir = uigetdir(options.directory,'Select LH labels directory');
end
if (rh && isempty(options.rhsegdir))
   if(options.automatch)
       options.rhsegdir = strrep(options.lhsegdir, 'lh', 'rh');
   else
       options.rhsegdir = uigetdir(options.directory,'Select RH labels directory');
   end
end

%% Part 3: Verify segment names
% If segment names were not specified, call to parseFSSegments
if (lh && isempty(options.lsegnames))
   %If you've indicated you want lh segment info but haven't provided the
   %names, prompt for them
   warntxt=sprintf('You have not specified the names of the lh segments.\nYou will now be prompted to locate an appropriate file (LEFT HEMISPHERE).');
   d=warndlg(warntxt,'!! Warning !!');
   uiwait(d);
   options.lsegnames=parseFSSegments();
end
if (rh && isempty(options.rsegnames))
   %If you've indicated you want rh segment info but haven't provided the
   %names, prompt for them
   warntxt=sprintf('You have not specified the names of the rh segments.\nYou will now be prompted to locate an appropriate file (RIGHT HEMISPHERE).');
   d=warndlg(warntxt,'!! Warning !!');
   uiwait(d);
   options.rsegnames=parseFSSegments();
end

%% Part 4: Iterate through files associated with segments
xyz=nan(length(options.lsegnames)+length(options.rsegnames),3);
hemi=[repmat('lh',length(options.lsegnames),1);repmat('rh',length(options.rsegnames),1)];
segnames=[(options.lsegnames);(options.rsegnames)];
i=0;
if(lh)
   for i=1:length(options.lsegnames)
       fname=fullfile(options.lhsegdir, ['lh.' options.lsegnames{i} '.label']);
       xyz(i,:)=getMeanLabelCoords(fname);
   end
end
rstart=i;
if(rh)
   for i=1:length(options.rsegnames)
       fname=fullfile(options.rhsegdir, ['rh.' options.rsegnames{i} '.label']);
       xyz(rstart+i,:)=getMeanLabelCoords(fname);
   end
end

%% Helper function: Find mean XYZ of named file
   function [xyz]=getMeanLabelCoords(labelfile)
        
       %First line is a header, second line is nvertices. Skip those two
       %and read the rest of the file: 5 columns, x y z are columns 2, 3,
       %and 4 (except when dlmread interprets it, the 2-space delimiters
       %between the columns is read as an extra column of zeros)
       M=dlmread(labelfile, ' ', 2, 0);
       x=mean(M(:,3));
       y=mean(M(:,5));
       z=mean(M(:,7));
       xyz=[x y z];
   end
end