Manage catalog
With this tool, the user can create an Obspy catalog merging the information from the location, source spectrum and moment tensor inversion output. Finally, the user can filer the catalog by time spam or geographically. We have added a method to write in human language the catalog. See an example:
Event 46: Date 01/02/2022 02:03:00.598999 rms 0.47 s Lat 42.5250 Lon 1.4252 Depth -2.0 km +- 1.3 min_dist 0.089 max_dist 0.621 smin 0.6 km smax 0.7 km ell_azimuth 161.3 gap 65.3 conf_lev 90.0 %
Magnitudes: Mw 3.58 +- 0.11
Moment Tensor Solution:
Mw 3.77 Mo 5.07e+14 Nm DC 44.64 % CLVD 22.70 % iso 32.65 % variance_red 34.54
Nodal Plane: Strike 319.9 Dip 40.0 Rake -74.7
Moment Tensor: mrr -1.05e+14 mtt 2.64e+14 mpp 3.38e+14 mrp 1.09e+13 mrt -7.33e+13 mrp -7.33e+13
station phase polarity date time time_residual distance_degrees distance_km azimuth takeoff_angle
PAND P ? 01/02/2022 02:03:01.900000 -0.32 0.1 9.9 90.8 89.4
PAND S ? 01/02/2022 02:03:03.330000 0.00 0.1 9.9 90.8 89.4
ARBS P ? 01/02/2022 02:03:02.340000 -0.48 0.1 13.6 137.9 90.8
ARBS S ? 01/02/2022 02:03:04.290000 -0.04 0.1 13.6 137.9 90.8
CEST P ? 01/02/2022 02:03:03.200000 -0.08 0.1 16.3 300.1 87.5
CEST S ? 01/02/2022 02:03:05.210000 0.10 0.1 16.3 300.1 87.5
CSOR P ? 01/02/2022 02:03:05.170000 -0.19 0.3 29.1 234.8 87.9
CSOR S ? 01/02/2022 02:03:09.050000 0.44 0.3 29.1 234.8 87.9
SALF S ? 01/02/2022 02:03:09.580000 -0.39 0.3 32.7 322.8 87.4
CORG P ? 01/02/2022 02:03:05.770000 -0.40 0.3 34.0 194.3 87.1
CORG S ? 01/02/2022 02:03:10.930000 0.97 0.3 34.0 194.3 87.1
GENF P ? 01/02/2022 02:03:06.150000 -0.21 0.3 35.1 19.0 87.6
GENF S ? 01/02/2022 02:03:10.530000 0.26 0.3 35.1 19.0 87.6
VALC P ? 01/02/2022 02:03:08.590000 -0.70 0.5 53.0 106.7 89.4
CARF P ? 01/02/2022 02:03:09.900000 -0.38 0.5 59.1 69.0 88.7
CARF S ? 01/02/2022 02:03:17.300000 0.43 0.5 59.1 69.0 88.7
FNEB P ? 01/02/2022 02:03:11.830000 -0.09 0.6 69.1 52.5 88.3
FNEB S ? 01/02/2022 02:03:20.730000 1.11 0.6 69.1 52.5 88.3
Manage Catalog from CLI
Usage
>> surfquake buildcatalog -e [path_event_files_folder] -s [path_source_summary_file] -m [path_mti_summary_file] -t [type_of_catalog] -o [path_to_ouput_folder]
Interactive help
Run Create Metadatafrom CLI
>> surfquake buildcatalog -e ./locations -s source_summary.txt -m mti_source_summary.txt -t "QUAKEXML" -o ./catalog
Manage Catalog from Library
Classes
class BuildCatalog:
def __init__(self, loc_folder, output_path, format="QUAKEML", source_summary_file=None, mti_summary_file=None):
"""
BuildCatalog class helps to join information from all surfquake outputs and create a catalog
Attributes:
- loc_folder (str): Path to the folder where the user have the locations files *hyp
- output_path (str): Output folder path where catalog object and file will be saved
- format (str): https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html
- source_summary_file (str): Path to the output file from source module
- mti_summary_file (str): Path to the output file from mti module
Methods:
- __init__(root_path): Initialize a new instance of BuildCatalog.
- __merge_info(catalog: Catalog): Merges the information from loc, source and mti
- build_catalog(): Starts the process to create the catalog from loc files, then calls __merge_info
"""
class WriteCatalog:
def __init__(self, path_catalog):
"""
WriteCatalog class helps to filter the catalog obj write in the most Readable
Attributes:
- path_catalog (str): Path to the pickle file saved when run build_catalog_loc from BuildCatalog class
Methods:
- __init__(root_path): Initialize a new instance of BuildCatalog.
- filter_time_catalog(verbose=True, **kwargs): filter the catalog in time span
- filter_geographic_catalog(verbose=True, **kwargs): filter the catalog in geographic contrain and magnitude
- write_catalog(catalog: Catalog, format, output_path)
- write_catalog_surf(catalog: Union[Catalog, None], output_path)
"""
self.path_catalog = path_catalog
self.catalog = []
self.__test_catalog()
Methods
build_catalog_loc
Instance method from BuildCatalog Class
filter_time_catalog
Instance method from WriteCatalog Class
def filter_time_catalog(self, verbose=True, **kwargs):
"""
Filter the catalog readed in the class instantiation
verbose: bool:
**kwargs
starttime :str: starttime to filter the catalog in format %d/%m/%Y, %H:%M:%S.%f
endtime :str: endtime to filter the catalog in format %d/%m/%Y, %H:%M:%S.%f
example:
wc = WriteCatalog(catalog_path)
catalog_filtered = wc.filter_time_catalog(starttime="30/01/2022, 00:00:00.0",
endtime="20/02/2022, 00:00:00.0")
return :catalog obj:
"""
filter_geographic_catalog
Instance method from WriteCatalog Class
def filter_geographic_catalog(self, catalog: Union[Catalog, None], verbose= True, **kwargs):
"""
Filter the catalog readed in the class instantiation or the catalog provided in bu the user when the method
is called
verbose: bool:
catalog obj (optional), if not found it uses catalog from the catalog attribute
**kwargs --> All keys must be used to filter success.
lat_min:float
lat_max:float
lon_min:float
lon_max:float
depth_min:float: km
depth_max:float: km
mag_min:float
mag_max:float
return :catalog obj:
"""
write_catalog_surf
Instance method from WriteCatalog Class
def write_catalog_surf(self, catalog: Union[Catalog, None], output_path):
"""
Writes in human language the catalog instantiated with the class
verbose: bool:
catalog obj (optional), if not found it uses catalog from the catalog attribute
"""
Examples using library
First, We can join all infrmation from located events in your loc folder, with the information of magnitudes and source parameters plus the ourput from Moment Tensor.
Secondly, We can play with the catalog, filtering it (filter_time_catalog and filter_geographic_catalog) or writing a human readable version using write_catalog_surf method.