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 ###############################################################################