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