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')