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))
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))
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 = []
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=[]