-
Notifications
You must be signed in to change notification settings - Fork 0
/
aggregate_kogia.m
42 lines (40 loc) · 1.72 KB
/
aggregate_kogia.m
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
baseDir = 'I:\JAX13D_broad_metadata'; % directory containing de_detector output
outDir = 'I:\JAX13D_broad_metadata\TPWS_116'; % directory where you want to save your TPWS file
siteName = 'JAX_D_13_disk'; % site name, used to name the output file
tsWin = 300;
dirSet = dir(fullfile(baseDir,[siteName,'*']));
MSP = [];
MPP = [];
MSN = [];
MTT = [];
for itr0 = 1:length(dirSet)
if dirSet(itr0).isdir &&~strcmp(dirSet(itr0).name,'.')&&...
~strcmp(dirSet(itr0).name,'..')
fMatch = fullfile(baseDir,[dirSet(itr0).name,'\*tsMat.mat']);
fileSet = dir(fMatch);
lfs = length(fileSet);
for itr2 = 1:lfs
thisFullFileTs = fullfile(fileSet(itr2).folder,fileSet(itr2).name);
thisFullFilePred = strrep(thisFullFileTs,'tsMat.mat','predLab.mat');
if exist(thisFullFilePred)==2
load(thisFullFilePred)
ksppLabels = predLabels==2;
if sum(ksppLabels)>0
load(strrep(thisFullFileTs,'_tsMat.mat','.mat'),...
'specClickTf','ppSignal','f','clickTimes','hdr');
load(thisFullFileTs)
fileStart = datenum(hdr.start.dvec);
posDnum = (clickTimes(:,1)/(60*60*24)) + fileStart +...
datenum([2000,0,0,0,0,0]);
MSP = [MSP;specClickTf(ksppLabels,:)];
MPP = [MPP;ppSignal(ksppLabels)];
MSN = [MSN;tsMat(ksppLabels,:)];
MTT = [MTT;posDnum(ksppLabels,:)];
end
end
fprintf('done with file %0.0f of %0.0f\n',itr2,lfs)
end
end
1;
end
save(fullfile(baseDir,'Kogia_TPWS33.mat'),'MTT','MSN','MSP','MPP','f','-v7.3')