altaircode/primeLoop.m at main · psiwex/altaircode · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
clear;
clc;
tic;

subst='.bdf';

dirName = 'C:\Users\John\Documents\MATLAB\soarData\';

numErn=1;
numLst=1;
numNpu=1;

raw_bdf_loadpath = dirName;

raw_dataset_savepath = dirName;
ongoing_dataset_savepath = dirName;

processed_dataset_savepath = strcat(dirName,'\Completed\');

erpset_savepath= strcat(dirName,'\erpLabs\');
raw_eventlist_savepath= strcat(dirName,'\eventList\');

processed_eventlist_savepath= strcat(dirName,'\procEventList\');

binlisterPathErn = 'C:\Users\John\Documents\MATLAB\soarEtl\currentSoar\GorkaBinlister.txt';
binListerPathLst = 'C:\Users\John\Documents\MATLAB\soarEtl\currentSoar\LstBinlister.txt';
binListerPathNpu = 'C:\Users\John\Documents\MATLAB\soarEtl\currentSoar\NpuBinlister.txt';

% channelLocationFile: this path is found in your eeglab dipfit plugin folder: eeglab\plugins\dipfit\standard_BESA\standard-10-5-cap385.elp
channelLocationFile = 'C:\Users\John\Documents\MATLAB\eeglab2021.1\plugins\dipfit\standard_BESA\standard-10-5-cap385.elp';

xx=readlocs(channelLocationFile);

namesErn={};
namesLst={};
namesNpu={};

fName='OSU-00001-04B-01';

fName='UC0004';
nName=strcat(fName,'_NPU.bdf');
lName=strcat(fName,'_LST.bdf');
eName=strcat(fName,'_ERN.bdf');

outEx2='ProcessedEpochs.mat';

%% ern
try
[EEGe,ERPe,erne] = altairErnPreproc(eName,raw_dataset_savepath,ongoing_dataset_savepath,erpset_savepath,binlisterPathErn,channelLocationFile);
[epochsAcceptedPercente] = altairAcceptedEpochs(EEGe);
EEGe.epochsAcceptedPercent=epochsAcceptedPercente;
eOut=strcat(eName,outEx2);
save(eOut,'EEGe')
namesErn{numErn}=eOut;
numErn=numErn+1;

catch
end

%% lst
try
[EEGl,ERPl,ernl] = altairDoorsPreproc(lName,raw_dataset_savepath,ongoing_dataset_savepath,erpset_savepath,binListerPathLst,channelLocationFile);
[epochsAcceptedPercentl] = altairAcceptedEpochs(EEGl);
EEGl.epochsAcceptedPercent=epochsAcceptedPercentl;
lOut=strcat(lName,outEx2);
save(lOut,'EEGl')
namesLst{numLst}=lOut;
numLst=numLst+1;
catch

end

%% npu
try
[EEGn,ERPn,ernn] = altairNpuPreproc(nName,raw_dataset_savepath,ongoing_dataset_savepath,erpset_savepath,binListerPathNpu,channelLocationFile);
[epochsAcceptedPercentn] = altairAcceptedEpochs(EEGn);
EEGn.epochsAcceptedPercent=epochsAcceptedPercentn;
nOut=strcat(nName,outEx2);
save(nOut,'EEGn')
namesNpu{numNpu}=nOut;
numNpu=numNpu+1;
catch
end