def setup(self):
self.data_pile = pile.make_pile(
self.data_paths, fileformat=self.data_format)
if self.data_pile.is_empty():
sys.exit('Data pile is empty!')
self.deltat_want = self.config.deltat_want or \
min(self.data_pile.deltats.keys())
self.n_samples = int(
(self.config.sample_length + self.config.tpad) / self.deltat_want)
logger.debug('loading marker file %s' % self.fn_markers)
# loads just plain markers:
markers = marker.load_markers(self.fn_markers)
if self.fn_events:
markers.extend(
[marker.EventMarker(e) for e in
load_events(self.fn_events)])
if self.sort_markers:
logger.info('sorting markers!')
markers.sort(key=lambda x: x.tmin)
marker.associate_phases_to_events(markers)
markers_by_nsl = {}
for m in markers:
if not m.match_nsl(self.config.reference_target.codes[:3]):
continue
if m.get_phasename().upper() != self.align_phase:
continue
markers_by_nsl.setdefault(m.one_nslc()[:3], []).append(m)
assert(len(markers_by_nsl) == 1)
# filter markers that do not have an event assigned:
self.markers = list(markers_by_nsl.values())[0]
if not self.labeled:
dummy_event = Event(lat=0., lon=0., depth=0.)
for m in self.markers:
if not m.get_event():
m.set_event(dummy_event)
self.markers = [m for m in self.markers if m.get_event() is not None]
if not len(self.markers):
raise Exception('No markers left in dataset')
self.config.channels = list(self.data_pile.nslc_ids.keys())
self.config.channels.sort()
def setup(self):
self.data_pile = pile.make_pile(
self.data_paths, fileformat=self.data_format)
if self.data_pile.is_empty():
sys.exit('Data pile is empty!')
self.deltat_want = self.config.deltat_want or \
min(self.data_pile.deltats.keys())
self.n_samples = int(
(self.config.sample_length + self.config.tpad) / self.deltat_want)
logger.debug('loading marker file %s' % self.fn_markers)
# loads just plain markers:
markers = marker.load_markers(self.fn_markers)
if self.fn_events:
markers.extend(
[marker.EventMarker(e) for e in
load_events(self.fn_events)])
marker.associate_phases_to_events(markers)
markers = [m for m in markers if isinstance(m, marker.PhaseMarker)]
markers_dict = defaultdict(list)
for m in markers:
if m.get_phasename().upper() != self.align_phase:
continue
markers_dict[m.get_event()].append(m)
self.markers = []
for e, _markers in markers_dict.items():
first = min(_markers, key=lambda x: x.tmin)
self.markers.append(first)
if not self.labeled:
dummy_event = Event(lat=0., lon=0., depth=0.)
for m in self.markers:
if not m.get_event():
m.set_event(dummy_event)
self.markers = [m for m in self.markers if m.get_event() is not None]
if not len(self.markers):
raise Exception('No markers left in dataset')
self.config.channels = list(self.data_pile.nslc_ids.keys())
self.config.channels.sort()
event_id = int(line[:4])
date = line[5:27]
try:
t = datetime.strptime(date, '%d-%b-%Y %H:%M:%S.%f')
except ValueError:
continue
t = (t - datetime(1970, 1, 1)).total_seconds()
lat = float(line[28:34])
lon = float(line[35:41])
depth = float(line[43:47])
print line[56:60]
if year == 'temp':
mag = float(line[56:60])
else:
print line[60:64]
mag = float(line[60:64])
name = line[69::]
e = Event(lat=lat,
lon=lon,
depth=depth,
magnitude=mag,
name=name,
time=t,
catalog='Colm')
print e
m = EventMarker(e)
markers.append(m)
EventMarker.save_markers(markers, 'event_markers_Colm%s.txt' % year)
