def dummy_cat(self):
scnls1 = ["UK.STA1..HHZ", "UK.STA2..HHZ"]
cat = ev.Catalog()
eve1 = ev.Event()
eve1.origins.append(ev.Origin(time=UTCDateTime()))
eve1.amplitudes = make_amplitudes(scnls1)
eve1.station_magnitudes = sm_generator(amplitudes=eve1.amplitudes)
scnls2 = ["UK.STA3..HHZ", "UK.STA4..HHZ", "UK.STA5..HHZ"]
eve2 = ev.Event()
eve2.origins.append(ev.Origin(time=UTCDateTime()))
eve2.station_magnitudes = sm_generator(scnls=scnls2)
cat.events = [eve1, eve2]
return cat
def dummy_cat(self):
scnls1 = ["UK.STA1..HHZ", "UK.STA2..HHZ"]
cat = ev.Catalog()
eve1 = ev.Event()
eve1.origins.append(ev.Origin(time=UTCDateTime()))
eve1.picks = pick_generator(scnls1)
eve1.amplitudes = amp_generator(picks=eve1.picks)
scnls2 = ["UK.STA3..HHZ", "UK.STA4..HHZ", "UK.STA5..HHZ"]
eve2 = ev.Event()
eve2.origins.append(ev.Origin(time=UTCDateTime()))
eve2.amplitudes = amp_generator(scnls=scnls2)
cat.events = [eve1, eve2]
return cat
def _create_origins():
ori = ev.Origin(
resource_id=ev.ResourceIdentifier('smi:local/First'),
time=UTCDateTime('2016-05-04T12:00:00'),
time_errors={'uncertainty': .01},
longitude=-111.12525,
longitude_errors={'uncertainty': .020},
latitude=47.48589325,
latitude_errors={'uncertainty': .021},
depth=2.123,
depth_errors={'uncertainty': 1.22},
depth_type='from location',
time_fixed=False,
epicenter_fixed=False,
reference_system_id=ev.ResourceIdentifier(),
method_id=ev.ResourceIdentifier(),
earth_model_id=ev.ResourceIdentifier(),
arrivals=[_get_arrival()],
composite_times=[_get_composite_times()],
quality=_get_origin_quality(),
origin_type='hypocenter',
origin_uncertainty=_get_origin_uncertainty(),
region='US',
evaluation_mode='manual',
evaluation_status='final',
)
state['origin_id'] = ori.resource_id
return ori
def catalog_no_magnitude(self):
""" get a events with no magnitudes (should just fill with NaN) """
t1 = obspy.UTCDateTime("2099-04-01T00-01-00")
ori = ev.Origin(time=t1, latitude=47.1, longitude=-100.22)
event = ev.Event(origins=[ori])
cat = ev.Catalog(events=[event])
return events_to_df(cat)
def create_origin(origintime, cinfo, latitude, longitude, depth):
'''
create_origin - function to create an ObsPy Origin
:param origintime: the origins time of occurence
:type origintime: :class: `~obspy.core.utcdatetime.UTCDateTime` object
:param cinfo:
:type cinfo:
:param latitude: latitude in decimal degree of the origins location
:type latitude: float
:param longitude: longitude in decimal degree of the origins location
:type longitude: float
:param depth: hypocentral depth of the origin
:type depth: float
:return: An ObsPy :class: `~obspy.core.event.Origin` object
'''
assert isinstance(origintime, UTCDateTime), "origintime has to be " \
"a UTCDateTime object, but " \
"actually is of type " \
"'%s'" % type(origintime)
origin = ope.Origin()
origin.time = origintime
origin.creation_info = cinfo
origin.latitude = latitude
origin.longitude = longitude
origin.depth = depth
return origin
def dummy_cat(self):
cat = ev.Catalog()
eve = ev.Event()
eve.origins.append(ev.Origin(time=UTCDateTime()))
eve.magnitudes = mag_generator(["ML", "Md", "MW"])
cat.append(eve)
return cat
def test_put_new_events(self, bing_ebank):
""" ensure a new event can be put into the bank """
ori = ev.Origin(time=obspy.UTCDateTime("2016-01-01"))
event = ev.Event(origins=[ori])
event.origins[0].depth_errors = None # see obspy 2173
bing_ebank.put_events(event)
event_out = bing_ebank.get_events(event_id=event.resource_id)
assert len(event_out) == 1
assert event_out[0] == event
def dummy_cat(self):
scnls1 = ["UK.STA1..HHZ", "UK.STA2..HHZ"]
cat = ev.Catalog()
eve1 = ev.Event()
eve1.origins.append(ev.Origin(time=UTCDateTime()))
eve1.preferred_origin_id = eve1.origins[0].resource_id
picks = pick_generator(scnls1)
eve1.picks = picks
eve1.preferred_origin().arrivals = make_arrivals(picks)
scnls2 = ["UK.STA3..HHZ", "UK.STA4..HHZ", "UK.STA5..HHZ"]
eve2 = ev.Event()
eve2.origins.append(ev.Origin(time=UTCDateTime()))
eve2.preferred_origin_id = eve2.origins[0].resource_id
picks = pick_generator(scnls2)
eve2.picks = picks
eve2.preferred_origin().arrivals = make_arrivals(picks)
cat.events = [eve1, eve2]
return cat
def _create_first_pick_origin(first_pick, channel_ser, depth):
""" Create an origin based on first pick and a channel series. """
msg = ("origin fixed to location and time of earliest pick by "
f"obsplus version {obsplus.__last_version__}")
comment = ev.Comment(text=msg)
odict = dict(
time=to_utc(first_pick["time"]),
latitude=channel_ser["latitude"],
longitude=channel_ser["longitude"],
depth=depth,
time_fixed=True,
comments=[comment],
)
return ev.Origin(**odict)
def event_non_orphaned_rejected_pick(self, event_rejected_pick):
""" Change both picks to rejected but reference one from arrival. """
eve = event_rejected_pick.copy()
picks = eve.picks
for pick in picks:
pick.evaluation_status = "rejected"
first_pick_rid = str(picks[0].resource_id)
# now create origin and arrival
arrival = ev.Arrival(pick_id=first_pick_rid)
origin = ev.Origin(
latitude=0, longitude=0, depth=10, time=picks[0].time, arrivals=[arrival]
)
eve.origins.append(origin)
return eve
def ebank_with_event_no_time(self, tmp_path_factory):
""" Create an event bank which has one file with no time. """
tmp_path = Path(tmp_path_factory.mktemp("basic"))
cat = obspy.read_events()
# clear origin from first event and add an empty one
cat[0].origins.clear()
new_origin = ev.Origin()
cat[0].origins.append(new_origin)
cat[0].preferred_origin_id = new_origin.resource_id
# now save the events
for num, event in enumerate(cat):
path = tmp_path / f"{num}.xml"
event.write(path, format="quakeml")
# init, update, return bank
return obsplus.EventBank(tmp_path).update_index()
def _create_origin(ser):
""" create an origin and attach to event """
event = get_object(ser.evid)
origin = oe.Origin(
resource_id=rid(ser.orid),
time=UTC(ser.time),
latitude=ser.lat,
longitude=ser.lon,
depth=ser.depth * 1000, # convert to m
)
# temporarily attach event reference to origin
origin.__dict__['event'] = event
event.origins.append(origin)
def event_origins(self, event_description):
"""Add an origin to the event"""
event_description.origins.append(
ev.Origin(time=obspy.UTCDateTime(), longitude=-111, latitude=37))
return event_description