-
Notifications
You must be signed in to change notification settings - Fork 0
/
mkTPWS_perDir.m
153 lines (131 loc) · 6.35 KB
/
mkTPWS_perDir.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
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
% mkTPWS_perDir.m
% Script takes output from de_detector.m and put it into a format for use in
% detEdit.
% Output:
% A TPWS.m file containing 4 variables:
% MTT: An Nx2 vector of detection start and end times, where N is the
% number of detections
% MPP: An Nx1 vector of recieved level (RL) amplitudes.
% MSP: An NxF vector of detection spectra, where F is dictated by the
% parameters of the fft used to generate the spectra and any
% normalization preferences.
% f = An Fx1 frequency vector associated with MSP
clearvars
% Setup variables:
baseDir = 'D:\DTmetadata\DT_postDecompFix\DT08'; % directory containing de_detector output
outDir = 'D:\DTmetadata\DT_postDecompFix\LTSAs'; % directory where you want to save your TPWS file
siteName = 'GofMX_DT08'; % site name, used to name the output file
ppThresh = 120; % minimum RL in dBpp. If detections have RL below this
% threshold, they will be excluded from the output file. Useful if you have
% an unmanageable number of detections.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% check if the output file exists, if not, make it
if ~exist(outDir,'dir')
fprintf('Creating output directory %s\n',outDir)
mkdir(outDir)
end
letterCode = 97:122;
dirSet = dir(fullfile(baseDir,[siteName,'*']));
for itr0 = 1:length(dirSet)
fSave = [];
if dirSet(itr0).isdir &&~strcmp(dirSet(itr0).name,'.')&&...
~strcmp(dirSet(itr0).name,'..')
letterFlag = 0; % flag for knowing if a letter should be appended to disk name
inDir = fullfile(baseDir,dirSet(itr0).name);
fileSet = what(inDir);
lfs = length(fileSet.mat);
clickTimesVec = [];
ppSignalVec = [];
specClickTfVec = [];
tsVecStore = [];
subTP = 1;
for itr2 = 1:lfs
thisFile = fileSet.mat(itr2);
load(char(fullfile(inDir,thisFile)),'-mat','clickTimes','hdr',...
'ppSignal','specClickTf','yFiltBuff','f','durClick')
if ~isempty(clickTimes)&& size(specClickTf,2)>1
% specClickTf = specClickTfHR;
keepers = find(ppSignal >= ppThresh);
ppSignal = ppSignal(keepers);
clickTimes = clickTimes(keepers,:);
[~,keepers2] = unique(clickTimes(:,1));
clickTimes = clickTimes(keepers2,:);
ppSignal = ppSignal(keepers2);
fileStart = datenum(hdr.start.dvec);
posDnum = (clickTimes(:,1)/(60*60*24)) + fileStart +...
datenum([2000,0,0,0,0,0]);
clickTimesVec = [clickTimesVec; posDnum];
ppSignalVec = [ppSignalVec; ppSignal];
tsWin = 200;
tsVec = zeros(length(keepers),tsWin);
for iTS = 1:length(keepers(keepers2))
thisClick = yFiltBuff{keepers(keepers2(iTS))};
[~,maxIdx] = max(thisClick);
% want to align clicks by max cycle
% f = fHR;
if isempty(fSave)
fSave = f;
end
dTs = (tsWin/2) - maxIdx; % which is bigger, the time series or the window?
dTe = (tsWin/2)- (length(thisClick)-maxIdx); % is the length after the peak bigger than the window?
if dTs<=0 % if the signal starts more than N samples ahead of the peak
% the start position in the TS vector has to be 1
posStart = 1;
% the start of the click ts has to be peak - 1/N
sigStart = maxIdx - (tsWin/2)+1;
else
% if it's smaller
posStart = dTs+1; % the start has to make up the difference
sigStart = 1; % and use the first index of the signal
end
if dTe<=0 % if it ends after the cut off of N samples
posEnd = tsWin; % the end in the TS vector has to be at N
sigEnd = maxIdx + (tsWin/2); % and last index is N/2 points after the peak.
else % if it ends before the end of the TS vector
posEnd = tsWin-dTe; % the end point has to be the
% difference between the window length and the click length
sigEnd = length(thisClick);
end
tsVec(iTS,posStart:posEnd) = thisClick(sigStart:sigEnd);
end
tsVecStore = [tsVecStore;tsVec];
if iscell(specClickTf)
spv = cell2mat(specClickTf');
specClickTfVec = [specClickTfVec; spv(:,keepers(keepers2))'];
else
specClickTfVec = [specClickTfVec; specClickTf(keepers(keepers2),:)];
end
clickTimes = [];
hdr = [];
specClickTf = [];
ppSignal = [];
end
fprintf('Done with file %d of %d \n',itr2,lfs)
if (size(clickTimesVec,1)>= 1800000 && (lfs-itr2>=10))|| itr2 == lfs
MSN = tsVecStore;
MTT = clickTimesVec;
MPP = ppSignalVec;
MSP = specClickTfVec;
if itr2 == lfs && letterFlag == 0
ttppOutName = [fullfile(outDir,dirSet(itr0).name),'_Delphin_TPWS1','.mat'];
fprintf('Done with directory %d of %d \n',itr0,length(dirSet))
subTP = 1;
else
ttppOutName = [fullfile(outDir,dirSet(itr0).name),char(letterCode(subTP)),'_Delphin_TPWS1','.mat'];
subTP = subTP+1;
letterFlag = 1;
end
f = fSave;
save(ttppOutName,'MTT','MPP','MSP','MSN','f','-v7.3')
MTT = [];
MPP = [];
MSP = [];
MSN = [];
clickTimesVec = [];
ppSignalVec = [];
specClickTfVec = [];
tsVecStore = [];
end
end
end
end