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

mseedpredictor systematically do not picks anything during the last 40 minutes of dayly .mseed file #56

Closed
Luchorigami opened this issue Mar 19, 2021 · 6 comments

Comments

@Luchorigami
Copy link

Hello,

I'm trying to run EQT on several 1day-long .mseed files with mseed_predictor and it seems that i have an issue:

After the picking, when I plot the P-S picks over the raw data it seems that there is systematically a gap of picks during the last 40 min of the trace. (i.e from 23H20min). I'm running EQT over a lot of 1-day long .mseed file and all present the same final gap of picks...

Moreover, it seems that changing the overlap value in the input of the mseed_predictor function make the length of this gap change...

Can you help me on that ? I tried to check deeper on the code but i did not find anything resolving this issue.
Tell me if you need more infos.

Best,
Luc

Note: I checked that there were effectively actually events 'potentially pickable' during the last 40 min of .mseed files.

@Luchorigami Luchorigami changed the title mseedpredictor systematically do not picks anything during the last 40 minutes of a dayly .mseed mseedpredictor systematically do not picks anything during the last 40 minutes of dayly .mseed file Mar 19, 2021
@smousavi05
Copy link
Owner

@Luchorigami could you send me one the plot and mseed file for one of those days that you think there are some events on the last 40 min ?

@Luchorigami
Copy link
Author

Hi thank you for your response,

Here is a plot for a day (23-04-2017) where we can easily see the final gap of picks. Vertical line are picks, red for P and blue for S.
EQTpicksentireday
EQTpickszoomend

(For the plot only:) Seismic traces amplitudes are not consistent over time in the plot because I have normalized each segment of 15 min by it's maximum amplitude to help the visualization of small event when large event occurs during the day.

The related 1 day long .mseed file is:
C1.MT01..HH.zip

I'm running EQT over 4 month of data, here is also the EQT output phase file that include this particular day (The station is MT01). There is no picks later than 23:20 over all the days tested.
X_prediction_results.txt

@smousavi05
Copy link
Owner

smousavi05 commented Mar 23, 2021 via email

@Luchorigami
Copy link
Author

Hi and thank you,

I changed the batch value and it seems to work fine. However, I needed to choose a very low value (ie. batch_size<=5) to make the gap disappear.
In the help of the function mseed_predictor you suggest to use a value between 200 and 1000. Do you think it can cause problems to use 5 or 1 for the batch_size (I'm not really sure that i understand how this variable works) ? Is it only for speed ?

Also I notice that using a strong overlap value increase the number of detection more than changing the batch_size.

Anyway, thanks for your help and for sharing this code,
Luc

@smousavi05
Copy link
Owner

smousavi05 commented Mar 23, 2021 via email

@smousavi05
Copy link
Owner

smousavi05 commented Mar 23, 2021 via email

@smousavi05 smousavi05 pinned this issue Apr 14, 2022
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