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
- 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
- class seismic.pick_harvester.createEnsembleXML.Magnitude(mag, mag_type)[source]¶
Bases:
object
- magnitude_type¶
- magnitude_value¶
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.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.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:
wavelet transforms
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
cwtsnr: quality measure based on wavelet analysis
dom_freq: dominant frequency of arrival energy
slope_ratio: quality measure based on waveform complexity analysis