-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathrun_workflow-bp0520.py
More file actions
48 lines (38 loc) · 1.14 KB
/
Copy pathrun_workflow-bp0520.py
File metadata and controls
48 lines (38 loc) · 1.14 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
import datetime
import json
import logging
from datetime import timedelta
import alaska_utils
import dask
import eqt_utils
import matplotlib.pyplot as plt
import numpy as np
import obspy
import pandas as pd
import postprocess
import requests
import seisbench
import seisbench.models as sbm
# Progress bar for dask operations:
from dask.diagnostics import ProgressBar
from obspy.clients.fdsn import Client
pbar = dask.diagnostics.ProgressBar()
pbar.register()
import warnings
warnings.filterwarnings("ignore")
waveform_length = 60
waveform_overlap = 30
starttime = datetime.datetime(2019, 5, 27)
endtime = datetime.datetime(2019, 5, 28)
# Pre=saved station list in pandas dataframe format
dfS = pd.read_parquet(
"https://github.com/zoekrauss/alaska_catalog/raw/main/data_acquisition/alaska_stations.parquet"
)
dfS = alaska_utils.station_list(dfS,starttime,endtime,elevation=False,network=False)
filt_type = 1
f1 = 5
f2 = 20
pick_info,gamma_picks = eqt_utils.ml_pick(
dfS, starttime, endtime, waveform_length, waveform_overlap, filt_type, f1=f1, f2=f2
)
pick_info.to_parquet("/home/adminuser/lynx/picks_20190527_bp0520_overlapped.parquet", version="2.6")