seismic.pick_harvester package

Submodules

seismic.pick_harvester.createEnsembleXML module

Description:

This script was initially written for inserting new picks into the ISC catalogue. We now use a unified csv catalogue (that Babak has prepared) and this script merges existing picks with those picked by our parallel picker and creates self-consistent SC3ML files to be ingested into Seiscomp3.

CreationDate: 20/11/18 Developer: rakib.hassan@ga.gov.au

Revision History:

LastUpdate: 20/11/18 RH

class seismic.pick_harvester.createEnsembleXML.Arrival(net, sta, loc, cha, lon, lat, elev, phase, utctime, distance)[source]

Bases: object

cha
distance
elev
lat
loc
lon
net
phase
sta
utctime
class seismic.pick_harvester.createEnsembleXML.Catalog(isc_coords_file, fdsn_inventory, our_picks, event_folder, output_path, discard_old_picks=False)[source]

Bases: object

get_id()[source]
class seismic.pick_harvester.createEnsembleXML.Event[source]

Bases: object

origin_list
preferred_magnitude
preferred_origin
public_id
class seismic.pick_harvester.createEnsembleXML.FDSNInv(fn, host=None, port=None)[source]

Bases: object

getClosestStation(lon, lat, maxdist=1000.0)[source]
rtp2xyz(r, theta, phi)[source]
class seismic.pick_harvester.createEnsembleXML.Magnitude(mag, mag_type)[source]

Bases: object

magnitude_type
magnitude_value
class seismic.pick_harvester.createEnsembleXML.Origin(utctime, lat, lon, depthkm)[source]

Bases: object

arrival_list
depthkm
lat
lon
magnitude_list
utctime
class seismic.pick_harvester.createEnsembleXML.OurPicks(fnList, phaseList)[source]

Bases: object

seismic.pick_harvester.pick module

seismic.pick_harvester.quality module

Description:

Computes quality measures of picks using various methods

References:

CreationDate: 24/01/19

Developer: rakib.hassan@ga.gov.au

Revision History:

LastUpdate: 24/01/19 RH LastUpdate: dd/mm/yyyy Who Optional description

seismic.pick_harvester.quality.compute_quality_measures(trc, trc_filtered, scales, plotinfo=None)[source]

Computes quality measures for a given pick based on:

  1. wavelet transforms

  2. waveform complexity analysis (similar to Higuchi fractal dimensions)

Parameters:

  • trc – raw obspy trace centred on pick-time

  • trc_filtered – filtered obspy trace centred on pick-time

  • scales – scales for computing continuous wavelet transforms

  • plotinfo – dictionary containing required plotting information (eventid, origintime, mag, net, sta, phase, ppsnr, pickid, outputfolder)

Returns:

  1. cwtsnr: quality measure based on wavelet analysis

  2. dom_freq: dominant frequency of arrival energy

  3. slope_ratio: quality measure based on waveform complexity analysis

seismic.pick_harvester.utils module

class seismic.pick_harvester.utils.Arrival(net, sta, loc, cha, lon, lat, elev, phase, utctime, distance)[source]

Bases: object

cha
distance
elev
lat
loc
lon
net
phase
sta
utctime
class seismic.pick_harvester.utils.Catalog(event_folder)[source]

Bases: object

get_events()[source]
get_preferred_origin_timestamps()[source]
class seismic.pick_harvester.utils.CatalogCSV(event_folder, events_only=True)[source]

Bases: object

get_events()[source]
get_preferred_origin_timestamps()[source]
class seismic.pick_harvester.utils.Event[source]

Bases: object

origin_list
preferred_magnitude
preferred_origin
public_id
class seismic.pick_harvester.utils.EventParser(xml_filename)[source]

Bases: object

getEvents()[source]
parseEvent(e)[source]
parseMagnitude(m)[source]
parseOrigin(o)[source]
class seismic.pick_harvester.utils.Magnitude(mag, mag_type)[source]

Bases: object

magnitude_type
magnitude_value
class seismic.pick_harvester.utils.Origin(utctime, lat, lon, depthkm)[source]

Bases: object

arrival_list
depthkm
lat
lon
magnitude_list
utctime

Module contents