response_params.water_level = 600
# --- Create new Archive and set path structure ---
archive = Archive(archive_path=data_in, stations=stations,
archive_format="YEAR/JD/STATION", response_inv=response_inv,
response_removal_params=response_params)
# --- Specify parameters for amplitude measurement ---
amp_params = AttribDict()
amp_params.signal_window = 5.0
amp_params.highpass_filter = True
amp_params.highpass_freq = 2.0
# --- Specify parameters for magnitude calculation ---
mag_params = AttribDict()
mag_params.A0 = "Greenfield2018_bardarbunga"
mag_params.amp_feature = "S_amp"
mags = LocalMag(amp_params=amp_params, mag_params=mag_params,
plot_amplitudes=True)
# --- Load the LUT ---
lut = read_lut(lut_file=lut_file)
# --- Create new Onset ---
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 = {
picker = GaussianPicker(onset=onset)
picker.plot_picks = True
# --- Create new LocalMag object ---
# All parameters are optional: see the documentation for a complete guide.
amp_params = AttribDict()
amp_params.signal_window = 5.
amp_params.bandpass_filter = True
amp_params.bandpass_lowcut = 2.
amp_params.bandpass_highcut = 20.
# A0 attenuation function is required: see the documentation for several
# built-in options, or specify your own function. All other parameters are
# optional - see the documentation for a complete guide.
mag_params = AttribDict()
mag_params.A0 = "Hutton-Boore" # NOTE: REQUIRED PARAMETER!
mag_params.use_hyp_dist = False
mag_params.amp_feature = "S_amp"
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,
