olap = overlap
samp = 0
df = tr.stats.sampling_rate
if trId(
tr.stats
)[1] != last_id or tr.stats.starttime - last_endtime > 1.0 / df:
data_buf = np.array([], dtype='float64')
olap = 0
while samp < tr.stats.npts:
data = tr.data[samp:samp + nfft - olap].astype('float64')
data = np.concatenate((data_buf, data))
data = detrend(data)
# Correct for frequency response of instrument
data = seisSim(data, tr.stats.sampling_rate, paz, inst_sim=inst)
data /= (paz['sensitivity'] / 1e9
) #V/nm/s correct for overall sensitivity
data = recStalta(data, s2p(2.5, tr), s2p(10.0, tr))
picked_values = triggerOnset(data, 3.0, 0.5, max_len=overlap)
#
for i, j in picked_values:
begin = tr.stats.starttime + float(i + samp - olap) / df
end = tr.stats.starttime + float(j + samp - olap) / df
f.write("%s,%s,%s\n" %
(str(begin), str(end), tr.stats.station))
olap = overlap # only needed for first time in loop
samp += nfft - overlap
data_buf = data[-overlap:]
print '.', # Progress Bar
last_endtime, last_id = trId(tr.stats)
f.close()
data_buf = np.array([], dtype='float64')
olap = 0
while samp < tr.stats.npts:
data = tr.data[samp:samp + nfft - olap].astype('float64')
data = np.concatenate((data_buf, data))
data = detrend(data)
# Correct for frequency response of instrument
data = seisSim(data,
df,
paz_remove=tr.stats.paz,
paz_simulate=inst,
remove_sensitivity=True)
# XXX is removed in seisSim... ?!
# XXX data /= (paz['sensitivity'] / 1e9) #V/nm/s correct for overall sensitivity
data = bandpass(data, LOW, HIGH, df)
data = recStalta(data, s2p(STA, tr), s2p(LTA, tr))
picked_values = triggerOnset(data, ON, OFF, max_len=overlap)
#
for i, j in picked_values:
begin = tr.stats.starttime + float(i + samp - olap) / df
end = tr.stats.starttime + float(j + samp - olap) / df
trigger_list.append(
(begin.timestamp, end.timestamp, tr.stats.station))
olap = overlap # only needed for first time in loop
samp += nfft - overlap
data_buf = data[-overlap:]
last_endtime, last_id = trId(tr.stats)
###############################################################################
# start of coincidence part
###############################################################################
