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,