'figname': output_base+"_pp_decluster.png",
}
declustered_catalogue = decluster_catalogue(catalogue, config)
plot_rate(catalogue, cumulative=True, color='#FF9D84', new_figure=False, overlay=True, linewidth=1.5, label="Gardner-Knopoff/Gardner-Knopoff", alpha=0.8)
catalogue = read_catalog(input_catalogue_file)
valid_months = catalogue.data['day'] != 0
catalogue.select_catalogue_events(valid_months)
valid_years = catalogue.data['year'] >= 1900
catalogue.select_catalogue_events(valid_years)
config = {
'decluster_method': 'gardner_knopoff', # afteran | gardner_knopoff
#'time_distance_window': GardnerKnopoffWindow(),
#'time_distance_window': GruenthalWindow(),
'time_distance_window': UhrhammerWindow(),
# afteran
'time_window': 60.,
# gardner_knopoff
'fs_time_prop': 1.0,
# write declustered
'filename': None,
#'filename': output_base+"_pp_decluster.csv",
# plot cluster
'plot':False,
'figname': output_base+"_pp_decluster.png",
}
declustered_catalogue = decluster_catalogue(catalogue, config)
plot_rate(catalogue, cumulative=True, color='#FF5460', new_figure=False, overlay=True, linewidth=1.5, label="Gardner-Knopoff/Uhrhammer", alpha=0.8)
def setUp(self):
self.gardner_knopoff_window = GardnerKnopoffWindow()
self.gruenthal_window = GruenthalWindow()
self.uhrhammer_window = UhrhammerWindow()
