Coverage for /usr/local/lib/python3.7/dist-packages/pyrocko/client/isc.py : 17%

# http://pyrocko.org - GPLv3 

# 

# The Pyrocko Developers, 21st Century 

# ---|P------/S----------~Lg---------- 

from __future__ import absolute_import, division 

import logging 

from pyrocko import util 

from pyrocko.util import urlopen 

from pyrocko.io import quakeml 

from .base_catalog import EarthquakeCatalog 

logger = logging.getLogger('pyrocko.client.isc') 

km = 1000. 

class ISC(EarthquakeCatalog): 

'''Interfacing the catalog of the Internation Seismological Centre (ISC). 

''' 

def __init__(self, catalog=None): 

self.events = {} 

def flush(self): 

self.events = {} 

def append_time_params(self, a, time_range): 

date_start_s, tstart_s = util.time_to_str( 

time_range[0], format='%Y-%m-%d %H:%M:%S').split() 

date_end_s, tend_s = util.time_to_str( 

time_range[1], format='%Y-%m-%d %H:%M:%S').split() 

date_start_s = date_start_s.split('-') 

date_end_s = date_end_s.split('-') 

a('start_year=%s' % date_start_s[0]) 

a('start_month=%s' % date_start_s[1]) 

a('start_day=%s' % date_start_s[2]) 

a('start_time=%s' % tstart_s) 

a('end_year=%s' % date_end_s[0]) 

a('end_month=%s' % date_end_s[1]) 

a('end_day=%s' % date_end_s[2]) 

a('end_time=%s' % tend_s) 

def iter_event_names( 

self, 

time_range=None, 

magmin=None, 

magmax=None, 

latmin=-90., 

latmax=90., 

lonmin=-180., 

lonmax=180.): 

p = [] 

a = p.append 

a('out_format=CATQuakeML') 

a('request=REVIEWED') 

a('searchshape=RECT') 

self.append_time_params(a, time_range) 

if magmin: 

a('min_mag=%g' % magmin) 

if magmax: 

a('max_mag=%g' % magmax) 

a('bot_lat=%g' % latmin) 

a('top_lat=%g' % latmax) 

a('left_lon=%g' % lonmin) 

a('right_lon=%g' % lonmax) 

url = 'http://www.isc.ac.uk/cgi-bin/web-db-v4?' + '&'.join(p) 

logger.debug('Opening URL: %s' % url) 

page = urlopen(url).read().decode() 

logger.debug('Received page (%i bytes)' % len(page)) 

if 'The search could not be run due to problems' in page: 

logger.warn('%s\nurl: %s' % (page, url)) 

return 

elif 'No events were found.' in page: 

logger.info('No events were found.') 

events = [] 

else: 

data = quakeml.QuakeML.load_xml(string=page) 

events = data.get_pyrocko_events() 

for ev in events: 

self.events[ev.name] = ev 

for ev in events: 

if time_range[0] <= ev.time and ev.time <= time_range[1]: 

yield ev.name 

def get_event(self, name): 

if name not in self.events: 

t = self._name_to_date(name) 

for name2 in self.iter_event_names( 

time_range=(t-24*60*60, t+24*60*60)): 

if name2 == name: 

break 

return self.events[name] 

def get_phase_markers(self, time_range, station_codes, phases): 

'''Download phase picks from ISC catalog and return them as a list 

of `pyrocko.gui.PhaseMarker` instances. 

:param time_range: Tuple with (tmin tmax) 

:param station_codes: List with ISC station codes 

(see http://www.isc.ac.uk/cgi-bin/stations?lista). 

If `station_codes` is 'global', query all ISC stations. 

:param phases: List of seismic phases. (e.g. ['P', 'PcP'] 

''' 

p = [] 

a = p.append 

a('out_format=QuakeML') 

a('request=STNARRIVALS') 

if station_codes == 'global': 

a('stnsearch=GLOBAL') 

else: 

a('stnsearch=STN') 

a('sta_list=%s' % ','.join(station_codes)) 

a('phaselist=%s' % ','.join(phases)) 

self.append_time_params(a, time_range) 

url = 'http://www.isc.ac.uk/cgi-bin/web-db-v4?' + '&'.join(p) 

logger.debug('Opening URL: %s' % url) 

page = urlopen(url) 

page = page.read().decode() 

if 'No stations were found.' in page: 

logger.info('No stations were found.') 

return [] 

logger.debug('Received page (%i bytes)' % len(page)) 

data = quakeml.QuakeML.load_xml(string=page) 

markers = data.get_pyrocko_phase_markers() 

markers = self.replace_isc_codes(markers) 

return markers 

def replace_isc_codes(self, markers): 

for m in markers: 

new_nslc_ids = [] 

for (n, s, l_, c) in m.get_nslc_ids(): 

l_ = l_.replace('--', '') 

c = c.replace('???', '*') 

new_nslc_ids.append((n, s, l_, c)) 

m.nslc_ids = new_nslc_ids 

return markers 

def _name_to_date(self, name): 

ds = name[-23:] 

t = util.str_to_time(ds, format='%Y-%m-%d_%H-%M-%S.3FRAC') 

return t