# Plot the magnitude-time density
magnitude_bin = 0.2
time_bin = 5.0 # in Decimal Years
plot_magnitude_time_density(catalogue, magnitude_bin, time_bin)
# In[ ]:
# Depth histogram
plot_depth_histogram(catalogue, 10.)
# In[ ]:
# Time-varying completeness
completeness = np.array([[1980., 3.0], [1985., 4.0], [1964., 5.0],
[1910., 6.5]])
plot_observed_recurrence(catalogue, completeness, 0.2, catalogue.end_year)
# In[ ]:
# Limit the catalogue to depths less than 50 km
valid_depth = catalogue.data['depth'] <= 50.
catalogue.select_catalogue_events(valid_depth)
plot_depth_histogram(catalogue, 2.0)
# In[ ]:
# Set-up the file writer
output_file_name = 'data_output/basic_demo_catalogue_1.csv'
writer = CsvCatalogueWriter(output_file_name)
# Write the catalogue to file
#plot_magnitude_time_scatter(catalogue, fmt_string='o', alpha=0.3, linewidth=0.0)
# Depth histogram
# plot_depth_histogram(catalogue, 10.)
filename = "/Users/pirchiner/Desktop/tmp_plot.png"
# Time-varying completeness
completeness = np.array([[1980., 3.0], [1985., 4.0], [1964., 5.0],
[1910., 6.5], [1900., 9.0]])
plot_observed_recurrence(
catalogue,
completeness,
0.2,
catalogue.end_year,
title="Recurrence [#Eq / Time]",
overlay=False,
markersize=10,
color=['#FEF2D8', '#F18C79'],
#linewidth=3,
alpha=0.6,
)
plt.show()
exit()
input_catalogue_file = 'data_input/hmtk_sa3'
#input_catalogue_file = 'data_input/hmtk_bsb2013'
###
### Catalogue cache or read/cache
get_completeness_adjusted_table,
_get_catalogue_bin_limits)
magnitude_bin = 0.2
end_year = 2014
# earthquake_count = get_completeness_adjusted_table(catalogue,
# completeness_table,
# magnitude_bin,
# catalogue.end_year)
#print 'Magnitude N(OBS) N(CUM) Log10(Nc)'
#for row in earthquake_count:
# print '%6.2f %10.3f %10.3f %10.3f' %(row[0], row[1], row[2], row[3])
plot_observed_recurrence(catalogue, completeness_table, magnitude_bin, catalogue.end_year)
# In[ ]:
# Set up the configuration parameters
recurrence_config = {'reference_magnitude': 0.0,
'magnitude_interval': 0.2,
'Average Type': 'Weighted'}
bml_recurrence = BMaxLikelihood()
bval, sigmab, seismicity_rate, sigma_rate = bml_recurrence.calculate(catalogue,
recurrence_config,
completeness_table)
print 'B-value = %9.4f +/- %9.4f' %(bval, sigmab)
from hmtk.plotting.seismicity.catalogue_plots import (
plot_observed_recurrence, get_completeness_adjusted_table,
_get_catalogue_bin_limits)
magnitude_bin = 0.2
end_year = 2014
# earthquake_count = get_completeness_adjusted_table(catalogue,
# completeness_table,
# magnitude_bin,
# catalogue.end_year)
#print 'Magnitude N(OBS) N(CUM) Log10(Nc)'
#for row in earthquake_count:
# print '%6.2f %10.3f %10.3f %10.3f' %(row[0], row[1], row[2], row[3])
plot_observed_recurrence(catalogue, completeness_table, magnitude_bin,
catalogue.end_year)
# In[ ]:
# Set up the configuration parameters
recurrence_config = {
'reference_magnitude': 0.0,
'magnitude_interval': 0.2,
'Average Type': 'Weighted'
}
bml_recurrence = BMaxLikelihood()
bval, sigmab, seismicity_rate, sigma_rate = bml_recurrence.calculate(
catalogue, recurrence_config, completeness_table)
# plot_magnitude_time_scatter(catalogue, fmt_string='o', alpha=0.3, linewidth=0.0)
# Depth histogram
# plot_depth_histogram(catalogue, 10.)
filename = "/Users/pirchiner/Desktop/tmp_plot.png"
# Time-varying completeness
completeness = np.array([[1980.0, 3.0], [1985.0, 4.0], [1964.0, 5.0], [1910.0, 6.5], [1900.0, 9.0]])
plot_observed_recurrence(
catalogue,
completeness,
0.2,
catalogue.end_year,
title="Recurrence [#Eq / Time]",
overlay=False,
markersize=10,
color=["#FEF2D8", "#F18C79"],
# linewidth=3,
alpha=0.6,
)
plt.show()
exit()
input_catalogue_file = "data_input/hmtk_sa3"
# input_catalogue_file = 'data_input/hmtk_bsb2013'
###
# In[ ]:
# Depth histogram
plot_depth_histogram(catalogue, 10.)
# In[ ]:
# Time-varying completeness
completeness = np.array([[1980., 3.0],
[1985., 4.0],
[1964., 5.0],
[1910., 6.5]])
plot_observed_recurrence(catalogue, completeness, 0.2, catalogue.end_year)
# In[ ]:
# Limit the catalogue to depths less than 50 km
valid_depth = catalogue.data['depth'] <= 50.
catalogue.select_catalogue_events(valid_depth)
plot_depth_histogram(catalogue, 2.0)
# In[ ]:
# Set-up the file writer
output_file_name = 'data_output/basic_demo_catalogue_1.csv'
writer = CsvCatalogueWriter(output_file_name)
