pick_path = path.join('record_picks', recfile)
# check if pick file exists
if not path.isfile(pick_path):
# fix stream channels
fix_stream_channels(path.join('iris_dump', mseedfile))
# read mseed
#st = read(path.join('mseed_dump', mseedfile))
st = read(path.join('iris_dump', mseedfile))
# remove junk channels
cannotMerge = False
try:
st = remove_low_sample_data(st)
except:
cannotMerge = True
###############################################################################
# associate event and get fe
###############################################################################
evFound = False
for evnum, ev in enumerate(evdict):
#ev['datetime'] = UTCDateTime(1996,9,25,7,49,56)
if st[0].stats.starttime > UTCDateTime(ev['datetime']-timedelta(seconds=601)) \
and st[0].stats.starttime < UTCDateTime(ev['datetime']+timedelta(seconds=300)):
evFound = True
eqlo = ev['lon']
eqla = ev['lat']
if review_wave == True:
# reopen for appending
f = open('mseed_qa.csv', 'a')
st = read(path.join('mseed_dump', msf))
st_filt = st.copy()
print(msf)
#st.plot()
fig = plt.figure(1, figsize=(13,9))
i = 0
jj = []
plt.suptitle(msf)
st_filt = remove_low_sample_data(st_filt)
freqmax = st_filt[0].stats.sampling_rate * 0.45
st_filt = st_filt.filter('bandpass', freqmin=0.1, freqmax=freqmax)
for j, tr in enumerate(st_filt):
if tr.stats.channel[1] != 'N' and len(st_filt) > 3:
i += 1
ax = plt.subplot(3, 1, i)
idxs = range(0, len(tr.data))
plt.plot(idxs, tr.data, '-', lw=0.5)
plt.title(tr.stats.channel)
jj.append(j)
# select
idx1, idx2 = plt.ginput(2)
plt.close()
