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

[bugfix] small workaround, will still cause error in case Z-component ist tilted by more than 5deg

Browse Source 
before: wfdat[0] was using HHE component of previously rotated stream, which could not be found in metadata (dataless fname still HH2)
This commit is contained in:
Marcel Paffrath 2019-04-03 15:00:08 +02:00
parent 0995350697
commit dfc9de69de
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
pylot/core/pick/autopick.py
View File

@ -251,6 +251,7 @@ def autopickstation(wfstream, pickparam, verbose=False,
Ao = None # Wood-Anderson peak-to-peak amplitude Ao = None # Wood-Anderson peak-to-peak amplitude
picker = 'auto' # type of picks picker = 'auto' # type of picks
# TODO: REMOVE HARDCODED CHANNEL MAPPING
# split components # split components
zdat = wfstream.select(component="Z") zdat = wfstream.select(component="Z")
if len(zdat) == 0: # check for other components if len(zdat) == 0: # check for other components
@ -298,7 +299,7 @@ def autopickstation(wfstream, pickparam, verbose=False,
if not metadata: if not metadata:
print('Warning: Could not use TauPy to estimate onsets as there are no metadata given.') print('Warning: Could not use TauPy to estimate onsets as there are no metadata given.')
else: else:
station_id = wfstream[0].get_id() station_id = zdat[0].get_id()
station_coords = metadata.get_coordinates(station_id, time=wfstream[0].stats.starttime) station_coords = metadata.get_coordinates(station_id, time=wfstream[0].stats.starttime)
if station_coords and origin: if station_coords and origin:
source_origin = origin[0] source_origin = origin[0]