Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

error on delayed STA/LTA #25

Open
fontiela opened this issue Mar 7, 2023 · 4 comments
Open

error on delayed STA/LTA #25

fontiela opened this issue Mar 7, 2023 · 4 comments

Comments

@fontiela
Copy link

fontiela commented Mar 7, 2023

When I select ```delayedstalta´´´ (I also adjusted the remaining triggering settings), I get the following error:

Using config file: settings_saojorge2023.cfg
Opening hdf5 table: SaoJorge_2023.h5
2022-03-26T00:08:01.000000Z
Couldn't find SJ01.EHZ.3k.
Couldn't find SJ02B.EHZ.3k.
Couldn't find SJ08.EHZ.3k.
Couldn't find SJ08B.EHZ.3k.
Couldn't find SJ09.EHZ.3k.
Couldn't find SJ11.EHZ.3k.
Couldn't find SJ12.EHZ.3k.
Couldn't find SJ13.EHZ.3k.
Couldn't find SJ15.EHZ.3k.
Couldn't find SJ16.EHZ.3k.
Couldn't find SJ18.EHZ.3k.
Length of Orphan table: 0
2022-03-26T00:26:50.000000Z
Couldn't find SJ01.EHZ.3k.
Couldn't find SJ02B.EHZ.3k.
Couldn't find SJ08.EHZ.3k.
Couldn't find SJ08B.EHZ.3k.
Couldn't find SJ09.EHZ.3k.
Couldn't find SJ11.EHZ.3k.
Couldn't find SJ12.EHZ.3k.
Couldn't find SJ13.EHZ.3k.
Couldn't find SJ15.EHZ.3k.
Couldn't find SJ16.EHZ.3k.
Couldn't find SJ18.EHZ.3k.
Length of Orphan table: 1
2022-03-26T01:35:39.000000Z
Couldn't find SJ01.EHZ.3k.
Couldn't find SJ02B.EHZ.3k.
Couldn't find SJ08.EHZ.3k.
Couldn't find SJ08B.EHZ.3k.
Couldn't find SJ09.EHZ.3k.
Couldn't find SJ11.EHZ.3k.
Couldn't find SJ12.EHZ.3k.
Couldn't find SJ13.EHZ.3k.
Couldn't find SJ15.EHZ.3k.
Couldn't find SJ16.EHZ.3k.
Couldn't find SJ18.EHZ.3k.
Traceback (most recent call last):
  File "catfill.py", line 74, in <module>
    alltrigs = redpy.trigger.trigger(st, stC, rtable, opt)
  File "/Users/joaofontiela/geophysics/REDPy/redpy/trigger.py", line 207, in trigger
    sta=opt.swin, lta=opt.lwin, details=True)
  File "/Users/joaofontiela/opt/anaconda3/envs/redpy/lib/python3.7/site-packages/obspy/signal/trigger.py", line 697, in coincidence_trigger
    tmp_triggers = trigger_onset(tr.data, thr_on, thr_off, **kwargs)
  File "/Users/joaofontiela/opt/anaconda3/envs/redpy/lib/python3.7/site-packages/obspy/signal/trigger.py", line 354, in trigger_onset
    ind2_[-1] = True
IndexError: index -1 is out of bounds for axis 0 with size 0

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "catfill.py", line 77, in <module>
    except (TypeError, obspy.fdsn.header.FDSNException, Exception):
AttributeError: module 'obspy' has no attribute 'fdsn'
Closing remaining open files:SaoJorge_2023.h5...done

This error happened before, in a much more advanced phase of the analysis.
Path to the files are correct on the settings file and I'm using 22 stations.

@ahotovec
Copy link
Owner

ahotovec commented Mar 7, 2023

Just to confirm a few more details here:

  • Does this work if you use the default STA/LTA parameters, or is it unique to the delayed STA/LTA?
  • Which version of the code are you using? I've been actively developing in a parallel branch and it would be useful to know if the problem persists there.

@fontiela
Copy link
Author

fontiela commented Mar 7, 2023

classic and rec STA/LTA work fine with default and personalised parameters.

@ahotovec
Copy link
Owner

ahotovec commented Mar 7, 2023

Ok, that helps. I will look into it. Default parameters for delayed?

@ahotovec
Copy link
Owner

ahotovec commented Mar 8, 2023

Confirmed that I can reproduce the problem. I'm not sure if the recommended values (lwin, swin, trigon, trigoff) are correct. I'm not having any luck getting delayedstalata to trigger with those values during times where there should be many triggers using classicstalta. I'll circle back around and see if there are better values to use, but in the meantime I'd say avoid using delayedstalta unless you've tested values outside of REDPy and are happy with the triggers they give you.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants