def test_evalresp(self):
"""
Tests evaluating instrument response information.
"""
client = Client()
dt = UTCDateTime("2005-01-01")
# plot as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="plot", filename=tempfile
)
self.assertEqual(open(tempfile, "rb").read(4)[1:4], "PNG")
os.remove(tempfile)
# plot-amp as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="plot-amp", filename=tempfile
)
self.assertEqual(open(tempfile, "rb").read(4)[1:4], "PNG")
os.remove(tempfile)
# plot-phase as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="plot-phase", filename=tempfile
)
self.assertEqual(open(tempfile, "rb").read(4)[1:4], "PNG")
os.remove(tempfile)
# fap as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="fap", filename=tempfile
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 1.202802E+04 1.792007E+02\n")
os.remove(tempfile)
# cs as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="cs", filename=tempfile
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 -1.202685E+04 1.677835E+02\n")
os.remove(tempfile)
# fap & def as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output="fap",
units="def",
filename=tempfile,
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 1.202802E+04 1.792007E+02\n")
os.remove(tempfile)
# fap & dis as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output="fap",
units="dis",
filename=tempfile,
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 7.557425E-01 2.692007E+02\n")
os.remove(tempfile)
# fap & vel as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output="fap",
units="vel",
filename=tempfile,
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 1.202802E+04 1.792007E+02\n")
os.remove(tempfile)
# fap & acc as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(
network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output="fap",
units="acc",
filename=tempfile,
)
self.assertEqual(open(tempfile, "rt").readline(), "1.000000E-05 1.914318E+08 8.920073E+01\n")
os.remove(tempfile)
# fap as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="fap")
np.testing.assert_array_equal(data[0], [1.00000000e-05, 1.20280200e04, 1.79200700e02])
# cs as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00", channel="BHZ", time=dt, output="cs")
np.testing.assert_array_equal(data[0], [1.00000000e-05, -1.20268500e04, 1.67783500e02])
def test_evalresp(self):
"""
Tests evaluating instrument response information.
"""
client = Client()
dt = UTCDateTime("2005-01-01")
# plot as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# plot-amp as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot-amp',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# plot-phase as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot-phase',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# fap as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# cs as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='cs',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 -1.202685E+04 1.677835E+02\n')
os.remove(tempfile)
# fap & def as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='def',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# fap & dis as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='dis',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 7.557425E-01 2.692007E+02\n')
os.remove(tempfile)
# fap & vel as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='vel',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# fap & acc as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='acc',
filename=tempfile)
self.assertEqual(open(tempfile, 'rt').readline(),
'1.000000E-05 1.914318E+08 8.920073E+01\n')
os.remove(tempfile)
# fap as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap')
np.testing.assert_array_equal(data[0],
[1.00000000e-05, 1.20280200e+04, 1.79200700e+02])
# cs as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='cs')
np.testing.assert_array_equal(data[0],
[1.00000000e-05, -1.20268500e+04, 1.67783500e+02])
def test_evalresp(self):
"""
Tests evaluating instrument response information.
"""
client = Client()
dt = UTCDateTime("2005-01-01")
# plot as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='plot',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# plot-amp as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='plot-amp',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# plot-phase as PNG file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='plot-phase',
filename=tempfile)
self.assertEqual(open(tempfile, 'rb').read(4)[1:4], 'PNG')
os.remove(tempfile)
# fap as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# cs as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='cs',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 -1.202685E+04 1.677835E+02\n')
os.remove(tempfile)
# fap & def as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap',
units='def',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# fap & dis as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap',
units='dis',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 7.557425E-01 2.692007E+02\n')
os.remove(tempfile)
# fap & vel as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap',
units='vel',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 1.202802E+04 1.792007E+02\n')
os.remove(tempfile)
# fap & acc as ASCII file
tempfile = NamedTemporaryFile().name
client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap',
units='acc',
filename=tempfile)
self.assertEqual(
open(tempfile, 'rt').readline(),
'1.000000E-05 1.914318E+08 8.920073E+01\n')
os.remove(tempfile)
# fap as NumPy ndarray
data = client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='fap')
np.testing.assert_array_equal(
data[0], [1.00000000e-05, 1.20280200e+04, 1.79200700e+02])
# cs as NumPy ndarray
data = client.evalresp(network="IU",
station="ANMO",
location="00",
channel="BHZ",
time=dt,
output='cs')
np.testing.assert_array_equal(
data[0], [1.00000000e-05, -1.20268500e+04, 1.67783500e+02])
from obspy.iris import Client
import numpy as np
import matplotlib.pyplot as plt
from obspy.core import UTCDateTime
from obspy.signal import cornFreq2Paz, seisSim
client = Client()
# dt = UTCDateTime(2005, 4, 26, 21, 0, 0, 31225480)
dt = UTCDateTime("2006-5-26")
data = client.evalresp("TA", "M03C", "--", "BHE", dt)
# data[0]
def test_evalresp(self):
"""
Tests evaluating instrument response information.
"""
client = Client()
dt = UTCDateTime("2005-01-01")
# plot as PNG file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot',
filename=tempfile)
with open(tempfile, 'rb') as fp:
self.assertEqual(fp.read(4)[1:4], b'PNG')
# plot-amp as PNG file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot-amp',
filename=tempfile)
with open(tempfile, 'rb') as fp:
self.assertEqual(fp.read(4)[1:4], b'PNG')
# plot-phase as PNG file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='plot-phase',
filename=tempfile)
with open(tempfile, 'rb') as fp:
self.assertEqual(fp.read(4)[1:4], b'PNG')
# fap as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 1.055999E+04 1.792007E+02\n')
# cs as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='cs',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 -1.055896E+04 1.473054E+02\n')
# fap & def as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='def',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 1.055999E+04 1.792007E+02\n')
# fap & dis as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='dis',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 6.635035E-01 2.692007E+02\n')
# fap & vel as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='vel',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 1.055999E+04 1.792007E+02\n')
# fap & acc as ASCII file
with NamedTemporaryFile() as tf:
tempfile = tf.name
client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap', units='acc',
filename=tempfile)
with open(tempfile, 'rt') as fp:
self.assertEqual(fp.readline(),
'1.000000E-05 1.680674E+08 8.920073E+01\n')
# fap as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='fap')
np.testing.assert_array_equal(
data[0], [1.00000000e-05, 1.05599900e+04, 1.79200700e+02])
# cs as NumPy ndarray
data = client.evalresp(network="IU", station="ANMO", location="00",
channel="BHZ", time=dt, output='cs')
np.testing.assert_array_equal(
data[0], [1.00000000e-05, -1.05589600e+04, 1.47305400e+02])