% readneurolocs() - read neuroscan polar location file (.asc)
%
% Usage:
% >> CHANLOCS = readneurolocs( filename );
% >> CHANLOCS = readneurolocs( filename, method, 'key1', val1, ...);
%
% Inputs:
% filename - file name or matlab cell array { names x_coord y_coord }
% method - [1 2, 3, 4 or 5] different import methods
%
% Optional inputs:
% same as caliblocs()
% note that if no optional input are provided, re-centering will be
% performed automatically and re-scaling of coordinates will be
% performed for '.asc' files (not '.dat' files).
%
% Outputs:
% CHANLOCS - EEGLAB channel location data structure. See
% help readlocs()
%
% Author: Arnaud Delorme, CNL / Salk Institute, 4 March 2003
%
% See also: readlocs()
% Copyright (C) 2003 Arnaud Delorme, Salk Institute, arno@salk.edu
%
% This program is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
function chanlocs = readneurolocs( filename, method, varargin)
if nargin < 1
help readneurolocs;
return;
end
if nargin < 2
plottag = 0;
end
if nargin < 2
disp('There are 5 methods to import Neuroscan files, if one does not work, try others')
try
chanlocs = readneurolocs( filename, 1);
return
catch
disp('Import method 1 failed, trying method 2');
try
chanlocs = readneurolocs( filename, 2);
return
catch
disp('Import method 2 failed, trying method 3');
try
chanlocs = readneurolocs( filename, 3);
return
catch
disp('Import method 3 failed, trying method 4');
try
chanlocs = readneurolocs( filename, 4);
return
catch
disp('Import method 3 failed, trying method 5');
try
chanlocs = readneurolocs( filename, 5);
return
catch
error('No import method worked');
end
end
end
end
end
end
% try new method
if method == 1
locs = readtable(filename, 'filetype', 'text', 'Delimiter', { ';' ' ' '\t' }, 'ConsecutiveDelimitersRule', 'join');
locs = table2cell(locs);
if mod(size(locs,1),2) == 1
error('Issue with file format')
end
halfHeight = size(locs,1)/2;
loc1 = locs(1:halfHeight,:);
loc2 = locs(halfHeight+1:end,:);
x = [loc2{:,3}]; x = x-0.5;
y = [loc2{:,4}]; y = y-0.5;
radius = sqrt(x.^2 + y.^2);
theta = atan2d(x, y);
for iChan = 1:halfHeight
chanlocs(iChan).label = loc1{iChan,2};
chanlocs(iChan).theta = theta(iChan);
chanlocs(iChan).radius = radius(iChan);
end
chanlocs = convertlocs( chanlocs, 'topo2all');
return
elseif method == 2
% read location file
% ------------------
if ischar(filename)
locs = loadtxt( filename, 'delim', 9 );
end
if ~ischar(filename) || locs{1,1}(1) == ';' || size(locs,2) < 5
if ~ischar(filename)
names = filename{1};
x = filename{2};
y = filename{3};
else
if locs{1,1}(1) == ';'
% remove trailing control channels
% --------------------------------
while isnumeric( locs{end,1} ) & locs{end,1} ~= 0
locs = locs(1:end-1,:);
end
% find first numerical index
% --------------------------
index = 1;
while ischar( locs{index,1} ) && index < size(locs,1)
index = index + 1;
end
% extract location array
% ----------------------
nchans = size( locs, 1 ) - index +1;
chans = [locs{end-nchans+1:end, 1:5}];
chans = reshape(chans,nchans,5); %% Added this line in order to get x = chans(:,3)
names = locs(end-nchans*2+1: end-nchans, 2);
for index = 1:length(names)
if ~ischar(names{index})
names{index} = int2str(names{index});
end
end
x = chans(:,3);
y = -chans(:,4);
else
[tmp2 tmpind] = sort( [ locs{:,1} ]);
locs = locs(tmpind,:);
y = [ locs{:,end} ];
x = [ locs{:,end-1} ];
x = x/513.1617*44;
y = y/513.1617*44;
names = locs(:,end-2);
end
end
% second solution using angle
% ---------------------------
[phi,theta] = cart2pol(x, y);
phi = phi/pi*180;
% convert to other types of coordinates
% -------------------------------------
labels = names';
labels = cellfun(@num2str, labels, 'UniformOutput', false);
chanlocs = struct('labels', labels', 'sph_theta_besa', mattocell(theta)', 'sph_phi_besa', mattocell(phi)'); %% labels instead of labels(:)
chanlocs = convertlocs( chanlocs, 'sphbesa2all');
for index = 1:length(chanlocs)
chanlocs(index).labels = num2str(chanlocs(index).labels);
end
% re-calibration
% --------------
chanlocs = adjustlocs(chanlocs, 'autoscale', 'on', 'autorotate', 'off', varargin{:});
end
elseif method == 3
chanlocs = readneurodat( filename);
elseif method == 4
% read location file
% ------------------
if ischar(filename)
locs = loadtxt( filename, 'delim', [9 ' ']);
end
if size(locs,2) == 5
% 5 rows, xyz positions
for index = 1:size(locs,1)
locs{index,3} = - locs{index,3};
end
chanlocs = struct('labels', locs(:,1), 'type', locs(:,2), 'X', locs(:,4), 'Y', locs(:,3), 'Z', locs(:,5));
chanlocs = convertlocs( chanlocs, 'cart2all');
elseif size(locs,2) == 4
chanlocs = readlocs(filename, 'filetype', 'custom', 'format', { 'labels' '-Y' 'X' 'Z' });
end
elseif method == 5
chanlocs = readlocs(filename, 'filetype', 'custom', 'format', { 'labels' 'ignore' '-Y' 'X' 'Z' });
end
Datasets
Models
