def test_issue623(self):
"""
obspy.iris bulkdataselect only returns last trace in result
"""
t1 = UTCDateTime("2011-03-11T06:31:30Z")
t2 = UTCDateTime("2011-03-11T06:48:00Z")
client = Client()
st = client.getWaveform("GE", "EIL", "", "BHZ", t1, t2)
self.assertEqual(len(st), 5)
def test_issue623(self):
"""
obspy.iris bulkdataselect only returns last trace in result
"""
t1 = UTCDateTime("2011-03-11T06:31:30Z")
t2 = UTCDateTime("2011-03-11T06:48:00Z")
client = Client()
st = client.getWaveform("GE", "EIL", "", "BHZ", t1, t2)
self.assertEqual(len(st), 5)
def test_getWaveform(self):
"""
Testing simple waveform request method.
"""
# simple example
client = Client()
start = UTCDateTime("2010-02-27T06:30:00.019538Z")
end = start + 20
stream = client.getWaveform("IU", "ANMO", "00", "BHZ", start, end)
self.assertEquals(len(stream), 1)
self.assertEquals(stream[0].stats.starttime, start)
self.assertEquals(stream[0].stats.endtime, end)
self.assertEquals(stream[0].stats.network, "IU")
self.assertEquals(stream[0].stats.station, "ANMO")
self.assertEquals(stream[0].stats.location, "00")
self.assertEquals(stream[0].stats.channel, "BHZ")
# no data raises an exception
self.assertRaises(Exception, client.getWaveform, "YY", "XXXX", "00", "BHZ", start, end)
def test_getWaveform(self):
"""
Testing simple waveform request method.
"""
# simple example
client = Client()
start = UTCDateTime("2010-02-27T06:30:00.019538Z")
end = start + 20
stream = client.getWaveform("IU", "ANMO", "00", "BHZ", start, end)
self.assertEquals(len(stream), 1)
self.assertEquals(stream[0].stats.starttime, start)
self.assertEquals(stream[0].stats.endtime, end)
self.assertEquals(stream[0].stats.network, 'IU')
self.assertEquals(stream[0].stats.station, 'ANMO')
self.assertEquals(stream[0].stats.location, '00')
self.assertEquals(stream[0].stats.channel, 'BHZ')
# no data raises an exception
self.assertRaises(Exception, client.getWaveform, "YY", "XXXX", "00",
"BHZ", start, end)
from obspy.iris import Client
from obspy.core import UTCDateTime
from obspy.core.util import NamedTemporaryFile
import matplotlib.pyplot as plt
import numpy as np
# MW 7.1 Darfield earthquake, New Zealand
t1 = UTCDateTime("2010-09-3T16:30:00.000")
t2 = UTCDateTime("2010-09-3T17:00:00.000")
# Fetch waveform from IRIS web service into a ObsPy stream object
client = Client()
st = client.getWaveform('NZ', 'BFZ', '10', 'HHZ', t1, t2)
# Download and save instrument response file into a temporary file
with NamedTemporaryFile() as tf:
respf = tf.name
client.saveResponse(respf, 'NZ', 'BFZ', '10', 'HHZ', t1, t2, format="RESP")
# make a copy to keep our original data
st_orig = st.copy()
# define a filter band to prevent amplifying noise during the deconvolution
pre_filt = (0.005, 0.006, 30.0, 35.0)
# this can be the date of your raw data or any date for which the
# SEED RESP-file is valid
date = t1
seedresp = {'filename': respf, # RESP filename
# when using Trace/Stream.simulate() the "date" parameter can
print (st[0].stats)
st.write('REF.EHZ.2009:082.mseed', format='MSEED')
#st.write('REF.EHZ.2009:082.wave', format='WAV', framerate=6000)
st.write('REF.EHZ.2009:082.sac', format='SAC')
# This needs basemap and demonstrates the event and catalog classes
from obspy.core.event import *
cat = readEvents(\
"http://www.seismicportal.eu/services/event/search?magMin=8.0")
cat.plot()
# IRIS DMC example
from obspy.iris import Client
from obspy.core import UTCDateTime
client = Client()
t = UTCDateTime("2012-08-05T06:00:00.000")
st = client.getWaveform('IU', 'ANMO', '00', 'BHZ', t, t + 300)
st.plot()
# Earthworm wave server example - connection is refused though
from obspy.earthworm import Client
client = Client("pele.ess.washington.edu", 16017)
response = client.availability("UW", "TUCA", channel="BHZ")
print response
t = response[0][4]
st = client.getWaveform('UW', 'TUCA', '', 'BH*', t + 100, t + 130)
st.plot()
print st
print(st[0].stats)
st.write('REF.EHZ.2009:082.mseed', format='MSEED')
#st.write('REF.EHZ.2009:082.wave', format='WAV', framerate=6000)
st.write('REF.EHZ.2009:082.sac', format='SAC')
# This needs basemap and demonstrates the event and catalog classes
from obspy.core.event import *
cat = readEvents(\
"http://www.seismicportal.eu/services/event/search?magMin=8.0")
cat.plot()
# IRIS DMC example
from obspy.iris import Client
from obspy.core import UTCDateTime
client = Client()
t = UTCDateTime("2012-08-05T06:00:00.000")
st = client.getWaveform('IU', 'ANMO', '00', 'BHZ', t, t + 300)
st.plot()
# Earthworm wave server example - connection is refused though
from obspy.earthworm import Client
client = Client("pele.ess.washington.edu", 16017)
response = client.availability("UW", "TUCA", channel="BHZ")
print response
t = response[0][4]
st = client.getWaveform('UW', 'TUCA', '', 'BH*', t + 100, t + 130)
st.plot()
