onset = STALTAOnset(position="centred", sampling_rate=50)
onset.phases = ["P", "S"]
onset.bandpass_filters = {
"P": [2, 16, 2],
"S": [2, 16, 2]}
onset.sta_lta_windows = {
"P": [0.2, 1.0],
"S": [0.2, 1.0]}
# --- Create new PhasePicker ---
picker = GaussianPicker(onset=onset)
picker.plot_picks = True
# --- Create new QuakeScan ---
scan = QuakeScan(archive, lut, onset=onset, picker=picker, mags=mags,
run_path=run_path, run_name=run_name, log=True,
loglevel="info")
# --- Set locate parameters ---
# For a complete list of parameters and guidance on how to choose them, please
# see the manual and read the docs.
scan.marginal_window = 1.0
scan.threads = 4 # NOTE: increase as your system allows to increase speed!
# --- Toggle plotting options ---
scan.plot_event_summary = True
scan.xy_files = "./inputs/XY_FILES/dike_xyfiles.csv"
# --- Toggle writing of waveforms ---
scan.write_cut_waveforms = True
# --- Read in station file ---
stations = read_stations(station_file)
# --- Create new Archive and set path structure ---
archive = Archive(archive_path=data_in, stations=stations,
archive_format="YEAR/JD/*_STATION_*")
# --- Load the LUT ---
lut = read_lut(lut_file=lut_out)
# --- Create new Onset ---
onset = STALTAOnset(position="classic", sampling_rate=500)
onset.phases = ["P", "S"]
onset.bandpass_filters = {
"P": [20, 200, 4],
"S": [10, 125, 4]}
onset.sta_lta_windows = {
"P": [0.01, 0.25],
"S": [0.05, 0.5]}
# --- Create new QuakeScan ---
scan = QuakeScan(archive, lut, onset=onset, run_path=run_path,
run_name=run_name, log=True, loglevel="info")
# --- Set detect parameters ---
scan.timestep = 0.75
scan.threads = 4 # NOTE: increase as your system allows to increase speed!
# --- Run detect ---
scan.detect(starttime, endtime)
mag_params.station_corrections = {}
mag_params.trace_filter = ".[BH]H[NE]$"
mag_params.noise_filter = 3.
mag_params.station_filter = ["KVE", "LIND"] # List of stations to exclude.
mag_params.dist_filter = False
mags = LocalMag(amp_params=amp_params, mag_params=mag_params)
mags.plot_amplitudes = True
# --- Create new QuakeScan ---
# If you do not want to calculate local magnitudes, specify `mags=None`
scan = QuakeScan(archive,
lut,
onset=onset,
picker=picker,
mags=mags,
run_path=run_path,
run_name=run_name,
log=True,
loglevel="info")
# --- Set locate parameters ---
# For a complete list of parameters and guidance on how to choose them, please
# see the manual and read the docs.
scan.marginal_window = 1
# NOTE: increase the thread-count as your system allows. The core migration
# routines are compiled against OpenMP, so the compute time (particularly for
# detect), will decrease roughly linearly with the number of threads used.
scan.threads = 4
# --- Toggle plotting options ---
# --- Create new Onset ---
onset = STALTAOnset(position="centred", sampling_rate=1000)
onset.phases = ["P", "S"]
onset.bandpass_filters = {"P": [20, 200, 4], "S": [10, 125, 4]}
onset.sta_lta_windows = {"P": [0.01, 0.25], "S": [0.05, 0.5]}
# --- Create new PhasePicker ---
picker = GaussianPicker(onset=onset)
picker.plot_picks = True
# --- Create new QuakeScan ---
scan = QuakeScan(archive,
lut,
onset=onset,
picker=picker,
run_path=run_path,
run_name=run_name,
log=True,
loglevel="info")
# --- Set locate parameters ---
# For a complete list of parameters and guidance on how to choose them, please
# see the manual and read the docs.
scan.marginal_window = 0.1
scan.threads = 4 # NOTE: increase as your system allows to increase speed!
# --- Toggle plotting options ---
scan.plot_event_summary = True
# --- Toggle writing of waveforms ---
scan.write_cut_waveforms = True
lut = read_lut(lut_file=lut_file)
# --- Decimate the lookup table ---
lut = lut.decimate([2, 2, 2])
# --- Create new Onset ---
onset = STALTAOnset(position="classic", sampling_rate=50)
onset.phases = ["P", "S"]
onset.bandpass_filters = {"P": [2, 16, 2], "S": [2, 14, 2]}
onset.sta_lta_windows = {"P": [0.2, 1.0], "S": [0.2, 1.0]}
# --- Create new QuakeScan ---
scan = QuakeScan(archive,
lut,
onset=onset,
run_path=run_path,
run_name=run_name,
log=True,
loglevel="info")
# --- Set detect parameters ---
# For a complete list of parameters and guidance on how to choose them, please
# see the manual and read the docs.
scan.timestep = 120.
# NOTE: increase the thread-count as your system allows. The core migration
# routines are compiled against OpenMP, so the compute time (particularly for
# detect), will decrease roughly linearly with the number of threads used.
scan.threads = 4
# --- Run detect ---
scan.detect(starttime, endtime)