marcel/pylot
3
0
Fork 0
Code Issues 7 Pull Requests 1 Releases Wiki Activity

Before exporting phases, check for manual picks and export them instead of automatic picks.

Browse Source
This commit is contained in:
Ludger Küperkoch 2017-06-28 14:23:24 +02:00
parent f52b42f009
commit d31fec8781
1 changed files with 18 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
pylot/core/io/data.py
View File

@ -188,11 +188,26 @@ class Data(object):
# try exporting event via ObsPy # try exporting event via ObsPy
elif fnext == '.obs': elif fnext == '.obs':
# check for stations picked automatically as well as manually
# Prefer manual picks!
evtdata_copy = self.get_evt_data().copy()
evtdata_org = self.get_evt_data()
for i in range(len(evtdata_org.picks)):
if evtdata_org.picks[i].method_id == 'manual':
mstation = evtdata_org.picks[i].waveform_id.station_code
mstation_ext = mstation + '_'
for k in range(len(evtdata_copy.picks)):
if evtdata_copy.picks[k].waveform_id.station_code == mstation or \
evtdata_copy.picks[k].waveform_id.station_code == mstation_ext and \
evtdata_copy.picks[k].method_id == 'auto':
del evtdata_copy.picks[k]
break
try: try:
self.get_evt_data().write(fnout + fnext, format=evtformat) evtdata_copy.write(fnout + fnext, format=evtformat)
# write header afterwards # write header afterwards
evtdata = self.get_evt_data() evid = str(evtdata_org.resource_id).split('/')[1]
evid = str(evtdata.resource_id).split('/')[1]
header = '# EQEVENT: Label: EQ%s Loc: X 0.00 Y 0.00 Z 10.00 OT 0.00 \n' % evid header = '# EQEVENT: Label: EQ%s Loc: X 0.00 Y 0.00 Z 10.00 OT 0.00 \n' % evid
nllocfile = open(fnout + fnext) nllocfile = open(fnout + fnext)
l = nllocfile.readlines() l = nllocfile.readlines()