def test_ichinose_file():
"""Test building verified QuakeML from Ichinose solution file"""
from qmlutil import (dumps, Rounder, ResourceURIGenerator, Root,
timestamp2isostr)
from qmlutil.plugins.xml import validate
from qmlutil.ichinose import IchinoseToQmlConverter
# Preprocessor for XML serialization
my_preproc = Rounder()
# Test QuakeML file from text MT solution
MT_FILE = os.path.join(PWD, 'data', 'mt_509589.txt')
with open(MT_FILE) as f:
ichicnv = IchinoseToQmlConverter(
f,
rid_factory=ResourceURIGenerator("quakeml", CFG['authority_id']),
utc_factory=timestamp2isostr,
agency=CFG['agency_id'],
)
event = ichicnv.get_event(anss=True)
qmlroot = ichicnv.event2root(event)
assert isinstance(qmlroot, dict)
assert 'q:quakeml' in qmlroot
# Generate QuakeML and validate
qmls = dumps(qmlroot, indent=" ", pretty=True, preprocessor=my_preproc)
if pytest.config.getoption("--writefiles"):
with open('/tmp/qmlutil-test-ichinose-mt.xml', 'w') as f:
f.write(qmls)
assert validate(qmls)
def test_nullphase():
"""Test case of no phase in arrival"""
import logging
from qmlutil import dumps, Rounder
from qmlutil.plugins.xml import validate
from qmlutil.plugins.antelope import Db2Quakeml
#from qmlutil.ichinose import IchinoseToQmlConverter
# Preprocessor for XML serialization
my_preproc = Rounder()
#-- From args
dsn = PWD + "/data/db/dbnullphase/db"
orid = 1414472
evid = 540344
# Convert and check for everything we asked for (ANSS, phases, etc)
logging.basicConfig() # config root logger TODO: change to mudule logger:
# 'qmlutil.plugins.antelope'
conv = Db2Quakeml(**CFG)
# ------------------------------------------------------------------------
# Test event generation
event = conv.get_event(dsn,
orid,
pick=True,
focalMechanism=True,
anss=True)
# Check event stuff, like anss...
assert event['type'] == "earthquake"
if isinstance(event['description'], dict):
assert event['description'].get('type') == "nearest cities"
assert 'origin' in event and len(event['origin']) > 0
assert 'magnitude' in event and len(event['magnitude']) > 0
assert 'pick' in event and len(event['pick']) > 0
assert event['@catalog:eventid'] == "00540344"
assert event['@catalog:dataid'] == "nn00540344"
assert event['@catalog:datasource'] == "nn"
assert event['@catalog:eventsource'] == "nn"
assert event['@publicID'] == "quakeml:edu.unr.seismo/event/540344"
qmlroot = conv.event2root(event)
assert isinstance(qmlroot, dict)
assert 'q:quakeml' in qmlroot
# Generate QuakeML and validate
qmls = dumps(qmlroot, indent=" ", pretty=True, preprocessor=my_preproc)
if pytest.config.getoption("--writefiles"):
with open('/tmp/qmlutil-test-nullphase.xml', 'w') as f:
f.write(qmls)
assert validate(qmls)
def test_db2qml():
"""Test the whole integrated shebang"""
import logging
from qmlutil import dumps, Rounder
from qmlutil.plugins.xml import validate
from qmlutil.plugins.antelope import Db2Quakeml
from qmlutil.ichinose import IchinoseToQmlConverter
# Preprocessor for XML serialization
my_preproc = Rounder()
#-- From args
dsn = PWD + "/data/db/dbfull/reno"
orid = 1371545
evid = 524465
# Convert and check for everything we asked for (ANSS, phases, etc)
logging.basicConfig() # config root logger TODO: change to mudule logger:
# 'qmlutil.plugins.antelope'
conv = Db2Quakeml(**CFG)
# ------------------------------------------------------------------------
# Test event generation
event = conv.get_event(dsn,
orid,
pick=True,
focalMechanism=True,
anss=True)
# Check event stuff, like anss...
assert event['type'] == "earthquake"
if isinstance(event['description'], dict):
assert event['description'].get('type') == "nearest cities"
assert 'origin' in event and len(event['origin']) > 0
assert 'magnitude' in event and len(event['magnitude']) > 0
assert 'pick' in event and len(event['pick']) > 0
assert 'focalMechanism' in event and len(event['focalMechanism']) > 0
assert event['@catalog:eventid'] == "00524465"
assert event['@catalog:dataid'] == "nn00524465"
assert event['@catalog:datasource'] == "nn"
assert event['@catalog:eventsource'] == "nn"
assert event['@publicID'] == "quakeml:edu.unr.seismo/event/524465"
qmlroot = conv.event2root(event)
assert isinstance(qmlroot, dict)
assert 'q:quakeml' in qmlroot
# Generate QuakeML and validate
qmls = dumps(qmlroot, indent=" ", pretty=True, preprocessor=my_preproc)
if pytest.config.getoption("--writefiles"):
with open('/tmp/qmlutil-test.xml', 'w') as f:
f.write(qmls)
assert validate(qmls)
# ------------------------------------------------------------------------
# Test delete from ORID
event = conv.get_deleted_event(dsn, orid=orid, anss=True)
assert event['type'] == "not existing"
assert event['@catalog:eventid'] == "00524465"
assert event['@catalog:dataid'] == "nn00524465"
assert event['@catalog:datasource'] == "nn"
assert event['@catalog:eventsource'] == "nn"
assert event['@publicID'] == "quakeml:edu.unr.seismo/event/524465"
qmlroot = conv.event2root(event)
assert isinstance(qmlroot, dict)
assert 'q:quakeml' in qmlroot
# Generate QuakeML and validate
qmls = dumps(qmlroot, indent=" ", pretty=True, preprocessor=my_preproc)
if pytest.config.getoption("--writefiles"):
with open('/tmp/qmlutil-test-delete-orid.xml', 'w') as f:
f.write(qmls)
assert validate(qmls)
# ------------------------------------------------------------------------
# Test delete from EVID -- non-existing
event = conv.get_deleted_event("/tmp/4c267822-3f72-4501-91fc-651851fd50a5",
evid=evid,
anss=True)
assert event['type'] == "not existing"
assert event['@catalog:eventid'] == "00524465"
assert event['@catalog:dataid'] == "nn00524465"
assert event['@catalog:datasource'] == "nn"
assert event['@catalog:eventsource'] == "nn"
assert event['@publicID'] == "quakeml:edu.unr.seismo/event/524465"
qmlroot = conv.event2root(event)
assert isinstance(qmlroot, dict)
assert 'q:quakeml' in qmlroot
# Generate QuakeML and validate
qmls = dumps(qmlroot, indent=" ", pretty=True, preprocessor=my_preproc)
if pytest.config.getoption("--writefiles"):
with open('/tmp/qmlutil-test-delete-evid.xml', 'w') as f:
f.write(qmls)
assert validate(qmls)
