def setup(self):
self.catalogs = { 'Geofon': catalog.Geofon(),
'USGS/NEIC PDE': catalog.USGS('pde'),
'USGS/NEIC US': catalog.USGS('us'),
'Global-CMT': catalog.GlobalCMT(),
'Kinherd': catalog.Kinherd(), }
catkeys = sorted(self.catalogs.keys())
self.set_name('Catalog Search')
self.add_parameter(Choice('Catalog', 'catalog', catkeys[0], catkeys))
self.add_parameter(Param('Min Magnitude', 'magmin', 0, 0, 10))
self.set_live_update(False)
def testUSGS(self):
def is_the_haiti_event(ev):
assert near(ev.magnitude, 7.0, 0.1)
assert near(ev.lat, 18.443, 0.01)
assert near(ev.lon, -72.571, 0.01)
assert near(ev.depth, 13000., 1.)
cat = catalog.USGS()
tmin = util.ctimegm('2010-01-12 21:50:00')
tmax = util.ctimegm('2010-01-13 03:17:00')
names = cat.get_event_names(time_range=(tmin, tmax), magmin=5.)
assert len(names) == 13
for name in names:
ev = cat.get_event(name)
if ev.magnitude >= 7.:
is_the_haiti_event(ev)
ident = ev.name
assert ident is not None
cat.flush()
ev = cat.get_event(ident)
is_the_haiti_event(ev)
import numpy as num
from pyrocko import trace, util, io, cake, catalog, automap, pile, model
from pyrocko import orthodrome, weeding
from pyrocko.client import fdsn
from pyrocko.io import resp, enhanced_sacpz as epz, stationxml
from beat import utility
from beat import heart
km = 1000.
g_sites_available = sorted(fdsn.g_site_abbr.keys())
geofon = catalog.Geofon()
usgs = catalog.USGS(catalog=None)
gcmt = catalog.GlobalCMT()
tfade_factor = 1.0
ffade_factors = 0.5, 1.5
fdsn.g_timeout = 60.
class starfill(object):
def __getitem__(self, k):
return '*'
def nice_seconds_floor(s):
nice = [
def USGS(searchparameter):
cat = catalog.USGS()
getevents(searchparameter, cat)
