def select_rfs_from_region():
"""Select rfs from region."""
for file_ in glob(PATH + '*_mout.QHD'):
ms = read(file_).select(expr='not st.mark')
for region in REGIONS:
ms2 = ms.select(around=region)
print len(ms2), 'vs', len(ms), file_
station = ms2[0].stats.station
num = get_region_number(dist=region[2])
ms2.setHI('eventsfrom', region)
ms2.setHI('eventsfromfe', util.feregion(region[0], region[1]))
newfile_ = TIMEVARI_PATH + ('rf_%s_R%s' % (station, num))
ms2.write(newfile_, 'Q')
def select_rfs_from_region():
"""Select rfs from region."""
for file_ in glob(PATH + '*_mout.QHD'):
ms = read(file_).select(expr='not st.mark')
for region in REGIONS:
ms2 = ms.select(around=region)
print len(ms2), 'vs', len(ms), file_
station = ms2[0].stats.station
num = get_region_number(dist=region[2])
ms2.setHI('eventsfrom', region)
ms2.setHI('eventsfromfe', util.feregion(region[0], region[1]))
newfile_ = TIMEVARI_PATH + ('rf_%s_R%s' % (station, num))
ms2.write(newfile_, 'Q')
def __init__(self, event_list):
"""
Init Event instance from event_list of events
"""
event_list = [AttribDict(event) for event in event_list]
for event in event_list:
event.datetime = UTC(event.datetime)
if feregion is not None:
event.flynn_region = feregion(event.latitude, event.longitude)
for item in ['datetime_quality', 'depth_quality', 'magnitude_type', 'author', 'quality', 'information', 'origin_id', 'flynn_region']:
if event.get(item) == None:
event[item] = ''
#if not event.get('magnitude_type'):
# event['magnitude_type'] = 'xx'
if event.get('id') == None:
event['id'] = (str(event['datetime']))[:-4].replace('-', '').replace(':', '').replace('.', '')
super(Events, self).__init__(event_list)
def __init__(self, event_list):
"""
Init Event instance from event_list of events
"""
event_list = [AttribDict(event) for event in event_list]
for event in event_list:
event.datetime = UTC(event.datetime)
if feregion is not None:
event.flynn_region = feregion(event.latitude, event.longitude)
for item in [
'datetime_quality', 'depth_quality', 'magnitude_type',
'author', 'quality', 'information', 'origin_id',
'flynn_region'
]:
if event.get(item) == None:
event[item] = ''
#if not event.get('magnitude_type'):
# event['magnitude_type'] = 'xx'
if event.get('id') == None:
event['id'] = (str(event['datetime']))[:-4].replace(
'-', '').replace(':', '').replace('.', '')
super(Events, self).__init__(event_list)
def test_util_region(self):
reg = util.feregion(10, 10)
self.assertEqual(reg, 'Nigeria')
