error on delayed STA/LTA

fontiela commented 1 year ago

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 commented 1 year ago

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 commented 1 year ago

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

ahotovec commented 1 year ago

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

ahotovec commented 1 year ago

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.

ahotovec / REDPy

error on delayed STA/LTA #25