def test_is_zmap_binary_files(self):
"""
Test zmap format detection on non-ZMAP (e.g. binary) files, see #1022.
"""
# Non-ZMAP file, binary
for filename in ["test.mseed", "test.sac"]:
file_ = get_example_file(filename)
self.assertFalse(zmap._is_zmap(file_))
def test_is_zmap_binary_files(self):
"""
Test zmap format detection on non-ZMAP (e.g. binary) files, see #1022.
"""
# Non-ZMAP file, binary
for filename in ["test.mseed", "test.sac"]:
file_ = get_example_file(filename)
self.assertFalse(zmap._is_zmap(file_))
def test_is_zmap(self):
"""
Test zmap format detection
"""
# Regular ZMAP
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertTrue(zmap._is_zmap(f.name))
# Pre-opened file
f.seek(0)
self.assertTrue(zmap._is_zmap(f))
# Extended ZMAP (13 columns)
self.zmap_fields += _EXT_ZMAP_FIELDS
self.test_data.update({
'h_err': '0.138679',
'z_err': '0.000000',
'm_err': '0.000000'
})
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertTrue(zmap._is_zmap(f.name))
# ZMAP string
test_string = self._serialize(test_events)
self.assertTrue(zmap._is_zmap(test_string))
# Non-ZMAP string
test_string = '0.000000\t' + test_string
self.assertFalse(zmap._is_zmap(test_string + '\n'))
# Non-ZMAP file (14 columns)
self.zmap_fields += ('dummy', )
self.test_data.update({'dummy': '0'})
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertFalse(zmap._is_zmap(f.name))
# Non-ZMAP file (non-numeric columns)
self.zmap_fields = _STD_ZMAP_FIELDS + _EXT_ZMAP_FIELDS
self.test_data.update({'mag': 'bad'})
test_events = [self.test_data]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertFalse(zmap._is_zmap(f.name))
def test_is_zmap(self):
"""
Test zmap format detection
"""
# Regular ZMAP
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertTrue(zmap._is_zmap(f.name))
# Pre-opened file
f.seek(0)
self.assertTrue(zmap._is_zmap(f))
# Extended ZMAP (13 columns)
self.zmap_fields += _EXT_ZMAP_FIELDS
self.test_data.update({'h_err': '0.138679', 'z_err': '0.000000',
'm_err': '0.000000'})
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertTrue(zmap._is_zmap(f.name))
# ZMAP string
test_string = self._serialize(test_events)
self.assertTrue(zmap._is_zmap(test_string))
# Non-ZMAP string
test_string = '0.000000\t' + test_string
self.assertFalse(zmap._is_zmap(test_string + '\n'))
# Non-ZMAP file (14 columns)
self.zmap_fields += ('dummy',)
self.test_data.update({'dummy': '0'})
test_events = [self.test_data, dict(self.test_data, mag='5.1')]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertFalse(zmap._is_zmap(f.name))
# Non-ZMAP file (non-numeric columns)
self.zmap_fields = _STD_ZMAP_FIELDS + _EXT_ZMAP_FIELDS
self.test_data.update({'mag': 'bad'})
test_events = [self.test_data]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode('utf-8'))
self.assertFalse(zmap._is_zmap(f.name))
def test_is_zmap(self):
"""
Test zmap format detection
"""
# Regular ZMAP
test_events = [self.test_data, dict(self.test_data, mag="5.1")]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode("utf-8"))
self.assertTrue(zmap._is_zmap(f.name))
# Pre-opened file
f.seek(0)
self.assertTrue(zmap._is_zmap(f))
# Extended ZMAP (13 columns)
self.zmap_fields += _EXT_ZMAP_FIELDS
self.test_data.update({"h_err": "0.138679", "z_err": "0.000000", "m_err": "0.000000"})
test_events = [self.test_data, dict(self.test_data, mag="5.1")]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode("utf-8"))
self.assertTrue(zmap._is_zmap(f.name))
# ZMAP string
test_string = self._serialize(test_events)
self.assertTrue(zmap._is_zmap(test_string))
# Non-ZMAP string
test_string = "0.000000\t" + test_string
self.assertFalse(zmap._is_zmap(test_string + "\n"))
# Non-ZMAP file (14 columns)
self.zmap_fields += ("dummy",)
self.test_data.update({"dummy": "0"})
test_events = [self.test_data, dict(self.test_data, mag="5.1")]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode("utf-8"))
self.assertFalse(zmap._is_zmap(f.name))
# Non-ZMAP file (non-numeric columns)
self.zmap_fields = _STD_ZMAP_FIELDS + _EXT_ZMAP_FIELDS
self.test_data.update({"mag": "bad"})
test_events = [self.test_data]
with NamedTemporaryFile() as f:
f.write(self._serialize(test_events).encode("utf-8"))
self.assertFalse(zmap._is_zmap(f.name))
tempout_dir,
ev_catalog,
start_itemp,
stop_itemp,
taup_model,
] = read_input_par(filepar)
st = Stream()
tr = Trace()
# half time length of templates is used to calculate origin_time
half_time = tlen_bef
# read event coordinates from catalog
print("event catalog should be ZMAP or QUAKEML")
if _is_zmap(ev_catalog):
print("reading ZMAP catalog")
elif _is_quakeml(ev_catalog):
print("reading QUAKEML catalog")
else:
print("warning error in reading ZMAP or QUAKEML")
cat = read_events(ev_catalog)
ncat = len(cat)
print(cat.__str__(print_all=True))
for iev in range(start_itemp, stop_itemp):
inplist = tempinp_dir + str(
