function test_suite = test_meeg_layout_collection
% tests for cosmo_meeg_layout_collection
%
% # For CoSMoMVPA's copyright information and license terms, #
% # see the COPYING file distributed with CoSMoMVPA. #
try % assignment of 'localfunctions' is necessary in Matlab >= 2016
test_functions = localfunctions();
catch % no problem; early Matlab versions can use initTestSuite fine
end
initTestSuite;
function test_fieldtrip_correspondence
if cosmo_skip_test_if_no_external('fieldtrip')
return
end
lc = cosmo_meeg_layout_collection();
rp_lay = randperm(numel(lc));
ntest_layout = 2;
ntest_chan_max = 10;
% may or may not be present
ignore_labels = {'COMNT', 'SCALE'};
% temporarily switch off warnings (emitted by fieldtrip)
warning_state = warning('query', 'all');
cleaner = onCleanup(@()warning(warning_state));
warning('off', 'all');
% shorter form
tolerance = 1e-4;
aeae = @(x, y)max(abs(x(:) - y(:))) < tolerance && isequal(size(x), size(y));
for ii = 1:ntest_layout
layout = lc{rp_lay(ii)};
cfg = struct();
cfg.layout = layout.name;
layout_ft = ft_prepare_layout(cfg);
% test outlnie
for j = 1:max(numel(layout_ft.outline), numel(layout.outline))
aeae(layout_ft.outline{j}, layout.outline{j});
end
% test mask
for j = 1:max(numel(layout_ft.mask), numel(layout.mask))
aeae(layout_ft.outline{j}, layout.outline{j});
end
% test labels
keep_idxs = find(~cosmo_match(layout.label, ignore_labels));
keep_ft_idxs = find(~cosmo_match(layout_ft.label, ignore_labels));
assertEqual(sort(layout.label(keep_idxs)), ...
sort(layout_ft.label(keep_ft_idxs)));
% test positions
ntest_chan = min(numel(keep_idxs), ntest_chan_max);
rp_chan = randperm(numel(keep_idxs));
for jj = 1:ntest_chan
test_label_idx = keep_idxs(rp_chan(jj));
label = layout.label{test_label_idx};
if strcmp(layout_ft.label{test_label_idx}, label)
% little optimization
pos = test_label_idx;
else
pos = find(cosmo_match(layout_ft.label, {label}));
assert(numel(pos) == 1);
end
aeae(layout_ft.width(pos, :), layout.width(test_label_idx, :));
aeae(layout_ft.height(pos, :), layout.height(test_label_idx, :));
aeae(layout_ft.pos(pos, :), layout.pos(test_label_idx, :));
end
end
function test_meeg_layout_collection_
if cosmo_skip_test_if_no_external('fieldtrip')
return
end
% get layout properties
% order is:
% - layout name
% - number of channels
% - some channel positions
% - labels of channel positions
lay_props = get_layout_properties();
% test each one
n = numel(lay_props);
clear('cosmo_meeg_layout_collection');
lc = cosmo_meeg_layout_collection();
lc_cached = cosmo_meeg_layout_collection();
assertEqual(lc, lc_cached);
lc_names = cellfun(@(x) x.name, lc, 'UniformOutput', false);
for k = 1:n
lay_prop = lay_props{k};
name = lay_prop{1};
nchan = lay_prop{2};
pos = lay_prop{3};
label = lay_prop{4};
i = find(cosmo_match(lc_names, name));
assertEqual(numel(i), 1);
lay = lc{i};
keep = find(~cosmo_match(lay.label, {'COMNT', 'SCALE'}));
assertEqual(numel(keep), nchan);
for j = 1:numel(label)
lab_i = find(cosmo_match(lay.label, label{j}));
assert(numel(lab_i) == 1);
assertElementsAlmostEqual(lay.pos(lab_i, :), pos(j, :), ...
'absolute', 1e-4);
end
end
function lay_props = get_layout_properties()
% the layout properties below were generated with the following code:
% names={'yokogawa440.lay', 'neuromag306planar.lay',...
% 'neuromag306mag.lay', 'neuromag306cmb.lay',...
% 'neuromag306all.lay', 'elec1020.lay',...
% 'elec1010.lay', 'elec1005.lay', 'biosemi32.lay',...
% 'biosemi128.lay', 'biosemi256.lay', 'EEG1020.lay',...
% 'EEG1010.lay', 'EEG1005.lay', 'CTF151.lay',...
% 'CTF275.lay', '4D148.lay', '4D248.lay'};
% lc=cosmo_meeg_layout_collection();
% lcn=cellfun(@(x) x.name,lc,'UniformOutput',false);
% for k=1:numel(names)
% fn=[names{k}];
% m=find(cosmo_match(lcn,fn));
% lay=lc{m};
% keep=find(~cosmo_match(lay.label,{'COMNT','SCALE'}));
% nch=numel(keep);
% desc=sprintf(['{''%s'',%d,...\n\t\t\t[%.4f %.4f;%.4f %.4f],'...
% '...\n\t\t\t{''%s'',''%s''}},...'],...
% lay.name,nch,lay.pos(keep([1 end]),:)',...
% lay.label{keep(1)},lay.label{keep(end)});
% fprintf('%s\n',desc);
% end
lay_props = { {'4D248.lay', 248, ...
[0.0038 0.0232; 0.3897 0.3453], ...
{'A1', 'A248'}}, ...
{'yokogawa440.lay', 440, ...
[-0.3337 0.3729; -0.0108 0.0726], ...
{'AG001', 'AG440'}}, ...
{'neuromag306planar.lay', 204, ...
[-0.4084 0.2532; 0.3733 -0.0820], ...
{'MEG0113', 'MEG2643'}}, ...
{'neuromag306mag.lay', 102, ...
[-0.4084 0.2732; 0.3733 -0.1036], ...
{'MEG0111', 'MEG2641'}}, ...
{'neuromag306cmb.lay', 102, ...
[-0.4084 0.2732; 0.3733 -0.1036], ...
{'MEG0112+0113', 'MEG2642+2643'}}, ...
{'neuromag306all.lay', 306, ...
[-0.4099 0.2532; 0.3761 -0.0976], ...
{'MEG0113', 'MEG2641'}}, ...
{'elec1020.lay', 21, ...
[-0.1710 0.4010; 0.1710 -0.4010], ...
{'Fp1', 'O2'}}, ...
{'elec1010.lay', 86, ...
[-0.1112 0.4261; 0.1391 -0.4257], ...
{'Fp1', 'I2'}}, ...
{'elec1005.lay', 335, ...
[-0.1112 0.3881; 0.1252 -0.3814], ...
{'Fp1', 'OI2'}}, ...
{'biosemi32.lay', 32, ...
[-0.1391 0.4500; 0.0000 0.0113], ...
{'Fp1', 'Cz'}}, ...
{'biosemi128.lay', 128, ...
[0.0000 0.0500; -0.4096 -0.2145], ...
{'A1', 'D32'}}, ...
{'biosemi256.lay', 256, ...
[0.0000 0.0587; -0.2290 -0.4060], ...
{'A1', 'H32'}}, ...
{'EEG1020.lay', 21, ...
[-0.1390 0.4280; 0.1390 -0.4280], ...
{'Fp1', 'O2'}}, ...
{'EEG1010.lay', 86, ...
[-0.1112 0.4260; 0.1390 -0.4256], ...
{'Fp1', 'I2'}}, ...
{'EEG1005.lay', 335, ...
[-0.1112 0.3880; 0.1252 -0.3814], ...
{'Fp1', 'OI2'}}, ...
{'CTF151.lay', 151, ...
[-0.0344 0.1732; 0.0008 -0.2668], ...
{'MLC11', 'MZP02'}}, ...
{'CTF275.lay', 275, ...
[-0.0171 0.1815; 0.0002 -0.2056], ...
{'MLC11', 'MZP01'}}, ...
{'4D148.lay', 148, ...
[-0.0109 0.0939; 0.3709 0.3364], ...
{'A1', 'A148'}}, ...
{'4D248.lay', 248, ...
[0.0038 0.0232; 0.3897 0.3453], ...
{'A1', 'A248'}}};
