def to_sac_stream(self, st, waveform): """ Return a SACTrace based on the given stream, containing metadata headers from the waveform """ tr = SACTrace.from_obspy_trace(st[0]) tr.kevnm = waveform.event_description[:16] event = waveform.event_ref() if not event: LOGGER.warn("Lost reference to event %s", waveform.event_description) else: origin = get_preferred_origin(waveform.event_ref()) if origin: tr.evla = origin.latitude tr.evlo = origin.longitude tr.evdp = origin.depth / 1000 tr.o = origin.time - waveform.start_time magnitude = get_preferred_magnitude(waveform.event_ref()) if magnitude: tr.mag = magnitude.mag channel = waveform.channel_ref() if not channel: LOGGER.warn("Lost reference to channel %s", waveform.sncl) else: tr.stla = channel.latitude tr.stlo = channel.longitude tr.stdp = channel.depth tr.stel = channel.elevation tr.cmpaz = channel.azimuth tr.cmpinc = channel.dip + 90 tr.kinst = channel.sensor.description[:8] return Stream([tr.to_obspy_trace()])
def to_sac_stream(self, st, waveform): """ Return a SACTrace based on the given stream, containing metadata headers from the waveform """ tr = SACTrace.from_obspy_trace(st[0]) tr.kevnm = waveform.event_description[:16] event = waveform.event_ref() if not event: LOGGER.warn("Lost reference to event %s", waveform.event_description) else: origin = get_preferred_origin(waveform.event_ref()) if origin: tr.evla = origin.latitude tr.evlo = origin.longitude tr.evdp = origin.depth / 1000 tr.o = origin.time - waveform.start_time magnitude = get_preferred_magnitude(waveform.event_ref()) if magnitude: tr.mag = magnitude.mag channel = waveform.channel_ref() if not channel: LOGGER.warn("Lost reference to channel %s", waveform.sncl) else: tr.stla = channel.latitude tr.stlo = channel.longitude tr.stdp = channel.depth tr.stel = channel.elevation tr.cmpaz = channel.azimuth tr.cmpinc = channel.dip + 90 tr.kinst = channel.sensor.description[:8] return Stream([tr.to_obspy_trace()])
def rows(self, data): """ Turn the data into rows (an iterable of lists) of QTableWidgetItems """ for i, event in enumerate(data): origin = get_preferred_origin(event) if not origin: continue magnitude = get_preferred_magnitude(event) if not magnitude: continue time = origin.time.strftime("%Y-%m-%d %H:%M:%S") # use strftime to remove milliseconds event_description = "no description" if len(event.event_descriptions) > 0: event_description = event.event_descriptions[0].text yield [ self.numericWidget(i), self.stringWidget(time), self.numericWidget(magnitude.mag, "%s %s" % (magnitude.mag, magnitude.magnitude_type)), self.numericWidget(origin.longitude, "%.03f°"), self.numericWidget(origin.latitude, "%.03f°"), self.numericWidget(origin.depth / 1000, "%.02f km"), # we wish to report in km self.stringWidget(event_description.title()), ]
def __init__(self, event, network, station, channel, *args, **kwargs): super(WaveformEntry, self).__init__(*args, **kwargs) self.event_ref = weakref.ref(event) self.network_ref = weakref.ref(network) self.station_ref = weakref.ref(station) self.channel_ref = weakref.ref(channel) self.sncl = get_sncl(network, station, channel) self.event_description = get_event_description(event) origin = get_preferred_origin(event) self.event_time = origin.time self.event_time_str = format_time_str(origin.time) self.event_depth = origin.depth / 1000 mag = get_preferred_magnitude(event) self.event_mag = "%s%s" % (mag.mag, mag.magnitude_type) self.event_mag_value = mag.mag self.waveform_id = '%s_%s' % (self.sncl, get_event_id(event)) self.distance = get_distance(origin.latitude, origin.longitude, station.latitude, station.longitude)
def __init__(self, event, network, station, channel, *args, **kwargs): super(WaveformEntry, self).__init__(*args, **kwargs) self.event_ref = weakref.ref(event) self.network_ref = weakref.ref(network) self.station_ref = weakref.ref(station) self.channel_ref = weakref.ref(channel) self.sncl = get_sncl(network, station, channel) self.event_description = get_event_description(event) origin = get_preferred_origin(event) self.event_time = origin.time self.event_time_str = format_time_str(origin.time) self.event_depth = origin.depth / 1000 mag = get_preferred_magnitude(event) self.event_mag = "%s%s" % (mag.mag, mag.magnitude_type) self.event_mag_value = mag.mag self.waveform_id = '%s_%s' % (self.sncl, get_event_id(event)) self.distance = get_distance( origin.latitude, origin.longitude, station.latitude, station.longitude)