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

[Bugfix] Wrong processing flow for checking picks.

Browse Source
This commit is contained in:
Ludger Küperkoch 2017-07-12 12:09:08 +02:00
parent 9ee32ee95e
commit da372b435c
1 changed files with 31 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -189,40 +189,6 @@ class Data(object):
else: else:
evtdata_copy = self.get_evt_data().copy() evtdata_copy = self.get_evt_data().copy()
evtdata_org = self.get_evt_data() evtdata_org = self.get_evt_data()
if upperErrors:
# check for pick uncertainties exceeding adjusted upper errors
# Picks with larger uncertainties will not be saved in output file!
for j in range(len(evtdata_org.picks)):
for i in range(len(evtdata_copy.picks)):
if evtdata_copy.picks[i].phase_hint[0] == 'P':
if evtdata_copy.picks[i].time_errors['lower_uncertainty'] >= upperErrors[0] or \
evtdata_copy.picks[i].time_errors['upper_uncertainty'] >= upperErrors[0] or \
evtdata_copy.picks[i].time_errors['uncertainty'] == None:
print("Uncertainty exceeds or equal adjusted upper time error!")
print("Adjusted uncertainty: {}".format(upperErrors[0]))
print("Pick uncertainty: {}".format(
max([evtdata_copy.picks[i].time_errors['lower_uncertainty'],
evtdata_copy.picks[i].time_errors['upper_uncertainty']])))
print("{1} P-Pick of station {0} will not be saved in outputfile".format(
evtdata_copy.picks[i].waveform_id.station_code,
evtdata_copy.picks[i].method_id))
del evtdata_copy.picks[i]
break
if evtdata_copy.picks[i].phase_hint[0] == 'S':
if evtdata_copy.picks[i].time_errors['lower_uncertainty'] >= upperErrors[1] or \
evtdata_copy.picks[i].time_errors['upper_uncertainty'] >= upperErrors[1] or \
evtdata_copy.picks[i].time_errors['uncertainty'] == None:
print("Uncertainty exceeds or equal adjusted upper time error!")
print("Adjusted uncertainty: {}".format(upperErrors[1]))
print("Pick uncertainty: {}".format(
max([evtdata_copy.picks[i].time_errors['lower_uncertainty'],
evtdata_copy.picks[i].time_errors['upper_uncertainty']])))
print("{1} S-Pick of station {0} will not be saved in outputfile".format(
evtdata_copy.picks[i].waveform_id.station_code,
evtdata_copy.picks[i].method_id))
del evtdata_copy.picks[i]
break
# check for stations picked automatically as well as manually # check for stations picked automatically as well as manually
# Prefer manual picks! # Prefer manual picks!
for i in range(len(evtdata_org.picks)): for i in range(len(evtdata_org.picks)):
@ -239,6 +205,37 @@ class Data(object):
if lendiff is not 0: if lendiff is not 0:
print("Manual as well as automatic picks available. Prefered the {} manual ones!".format(lendiff)) print("Manual as well as automatic picks available. Prefered the {} manual ones!".format(lendiff))
if upperErrors:
# check for pick uncertainties exceeding adjusted upper errors
# Picks with larger uncertainties will not be saved in output file!
for j in range(len(evtdata_org.picks)):
for i in range(len(evtdata_copy.picks)):
if evtdata_copy.picks[i].phase_hint[0] == 'P':
if evtdata_copy.picks[i].time_errors['upper_uncertainty'] >= upperErrors[0] or \
evtdata_copy.picks[i].time_errors['uncertainty'] == None:
print("Uncertainty exceeds or equal adjusted upper time error!")
print("Adjusted uncertainty: {}".format(upperErrors[0]))
print("Pick uncertainty: {}".format(evtdata_copy.picks[i].time_errors['uncertainty']))
print("{1} P-Pick of station {0} will not be saved in outputfile".format(
evtdata_copy.picks[i].waveform_id.station_code,
evtdata_copy.picks[i].method_id))
print("#")
del evtdata_copy.picks[i]
break
if evtdata_copy.picks[i].phase_hint[0] == 'S':
if evtdata_copy.picks[i].time_errors['upper_uncertainty'] >= upperErrors[1] or \
evtdata_copy.picks[i].time_errors['uncertainty'] == None:
print("Uncertainty exceeds or equal adjusted upper time error!")
print("Adjusted uncertainty: {}".format(upperErrors[1]))
print("Pick uncertainty: {}".format(evtdata_copy.picks[i].time_errors['uncertainty']))
print("{1} S-Pick of station {0} will not be saved in outputfile".format(
evtdata_copy.picks[i].waveform_id.station_code,
evtdata_copy.picks[i].method_id))
print("#")
del evtdata_copy.picks[i]
break
if fnext == '.obs': if fnext == '.obs':
try: try:
evtdata_copy.write(fnout + fnext, format=evtformat) evtdata_copy.write(fnout + fnext, format=evtformat)