예제 #1
0
    def setUp(self):
        import obspy.realtime
        rt_dict= obspy.realtime.rttrace.REALTIME_PROCESS_FUNCTIONS
        rt_dict['neg_to_zero']=(am_rt_signal.neg_to_zero,0)
        rt_dict['convolve']=(am_rt_signal.convolve,1)
        
        self.wo = RtWavelocOptions()
        self.wo.opdict['base_path'] = 'test_data'
        self.wo.opdict['outdir'] = 'Test'
        self.wo.opdict['time_grid'] = 'Slow_len.100m.P'
        self.wo.opdict['max_length'] = 120
        self.wo.opdict['safety_margin'] = 20
        self.wo.opdict['dt'] = 0.01
        self.wo.opdict['syn'] = True

        self.wo.verifyDirectories()

        # make synthetic data
        self.obs_list, self.ot, (x0,y0,z0) = make_synthetic_data(self.wo)

        self.starttime=self.obs_list[0].stats.starttime
        self.dt=self.obs_list[0].stats.delta

        # split data files to simulate packets of real-time data
        self.obs_split=[]
        for obs in self.obs_list:
            obs_split = obs / 3
            self.obs_split.append(obs_split)

        # generate ttimes_files for test
        n_test=150
        generate_random_test_points(self.wo, n_test, (x0, y0, z0))
예제 #2
0
    def setUp(self):
        import obspy.realtime
        rt_dict = obspy.realtime.rttrace.REALTIME_PROCESS_FUNCTIONS
        rt_dict['neg_to_zero'] = (am_rt_signal.neg_to_zero, 0)
        rt_dict['convolve'] = (am_rt_signal.convolve, 1)

        self.wo = RtWavelocOptions()
        self.wo.opdict['base_path'] = 'test_data'
        self.wo.opdict['outdir'] = 'Test'
        self.wo.opdict['time_grid'] = 'Slow_len.100m.P'
        self.wo.opdict['max_length'] = 120
        self.wo.opdict['safety_margin'] = 20
        self.wo.opdict['dt'] = 0.01
        self.wo.opdict['syn'] = True

        self.wo.verifyDirectories()

        # make synthetic data
        self.obs_list, self.ot, (x0, y0, z0) = make_synthetic_data(self.wo)

        self.starttime = self.obs_list[0].stats.starttime
        self.dt = self.obs_list[0].stats.delta

        # split data files to simulate packets of real-time data
        self.obs_split = []
        for obs in self.obs_list:
            obs_split = obs / 3
            self.obs_split.append(obs_split)

        # generate ttimes_files for test
        n_test = 150
        generate_random_test_points(self.wo, n_test, (x0, y0, z0))
예제 #3
0
    st = read(name)
    obs_list.append(st[0])

starttime = obs_list[0].stats.starttime
dt = obs_list[0].stats.delta
wo.opdict['dt'] = dt

# split data files to simulate packets of real-time data
obs_split = []
for obs in obs_list:
    split = obs / 3
    obs_split.append(split)

# generate ttimes_files for test
#generate_random_test_points(wo, n_test, (x0, y0, z0))
generate_random_test_points(wo, n_test)

tic = time.time()

# set up migrator
migrator = RtMigrator(wo)
nsta = migrator.nsta

ntr = len(obs_split[0])
#########################
# start loops
#########################
# loop over segments (simulate real-time data)
t_update_data = []
t_update_stacks = []
t_update_max = []
예제 #4
0
    st=read(name)
    obs_list.append(st[0])

starttime=obs_list[0].stats.starttime
dt=obs_list[0].stats.delta
wo.opdict['dt'] = dt

# split data files to simulate packets of real-time data
obs_split=[]
for obs in obs_list:
    split = obs / 3
    obs_split.append(split)

# generate ttimes_files for test
#generate_random_test_points(wo, n_test, (x0, y0, z0))
generate_random_test_points(wo, n_test)

tic=time.time()

# set up migrator
migrator = RtMigrator(wo)
nsta = migrator.nsta


ntr=len(obs_split[0])
#########################
# start loops
#########################
# loop over segments (simulate real-time data)
t_update_data=[]
t_update_stacks=[]