test meeg layout collectionΒΆ

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'}}};