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.createEnsembleXML.setup_logger(name, log_file, level=20)[source]

Function to setup a logger; adapted from stackoverflow

seismic.pick_harvester.pick module

Description:

Harvests picks from ASDF data-sets in parallel

References:

CreationDate: 13/09/18 Developer: rakib.hassan@ga.gov.au

Revision History:

LastUpdate: 13/09/18 RH LastUpdate: dd/mm/yyyy Who Optional description

seismic.pick_harvester.pick.dropBogusTraces(st, sampling_rate_cutoff=5)[source]
seismic.pick_harvester.pick.extract_p(taupy_model, pickerlist, event, station_longitude, station_latitude, st, win_start=- 50, win_end=50, resample_hz=20, bp_freqmins=[0.5, 2.0, 5.0], bp_freqmaxs=[5.0, 10.0, 10.0], margin=None, max_amplitude=100000000.0, plot_output_folder=None)[source]
seismic.pick_harvester.pick.extract_s(taupy_model, pickerlist, event, station_longitude, station_latitude, stn, ste, ba, win_start=- 50, win_end=50, resample_hz=20, bp_freqmins=[0.01, 0.01, 0.5], bp_freqmaxs=[1, 2.0, 5.0], margin=None, max_amplitude=100000000.0, plot_output_folder=None)[source]
seismic.pick_harvester.pick.getWorkloadEstimate(fds, originTimestamps)[source]
seismic.pick_harvester.pick.merge_results(output_path)[source]

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.Catalog(event_folder)[source]

Bases: object

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

Bases: object

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

Bases: object

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

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

Bases: object

class seismic.pick_harvester.utils.ProgressTracker(output_folder, restart_mode=False)[source]

Bases: object

increment()[source]
seismic.pick_harvester.utils.recursive_glob(treeroot, pattern)[source]
seismic.pick_harvester.utils.split_list(lst, npartitions)[source]

Module contents