def test_process_adjoint_2(): st, meta = prepare_real_adj_data() inv = obspy.read_inventory() event = obspy.read_events()[0] adjs = pa.convert_stream_to_adjs(st, meta) starttime = adjs[0].starttime - 20 delta = adjs[0].dt / 2.0 npts = 40 + 2 * len(adjs[0].adjoint_source) new_adj = pa.process_adjoint(adjs, interp_flag=True, interp_starttime=starttime, interp_delta=delta, interp_npts=npts, sum_over_comp_flag=True, weight_flag=False, filter_flag=True, pre_filt=[0.02, 0.025, 0.059, 0.073], taper_percentage=0.05, taper_type="hann", add_missing_comp_flag=True, rotate_flag=True, inventory=inv, event=event) assert len(new_adj) == 3
def test_convert_stream_to_adjs(): array = np.array([1., 2., 3., 4., 5.]) starttime = UTCDateTime(1990, 1, 1) adjsrcs = get_sample_adjsrcs(array, starttime) st, meta = pa.convert_adjs_to_stream(adjsrcs) adjsrcs_new = pa.convert_stream_to_adjs(st, meta) for adj, adj_new in zip(adjsrcs, adjsrcs_new): assert_adj_same(adj, adj_new)
def rotate_one_station_adjsrcs(sta_adjs, slat, slon, elat, elon): adj_stream, meta_info = convert_adjs_to_stream(sta_adjs) add_missing_components(adj_stream) rotate_one_station_stream( adj_stream, elat, elon, station_latitude=slat, station_longitude=slon, mode="RT->NE") new_adjs = convert_stream_to_adjs(adj_stream, meta_info) adj_dict = {} for _adj in new_adjs: adj_id = "%s_%s_%s" % (_adj.network, _adj.station, _adj.component) adj_dict[adj_id] = _adj return adj_dict
def test_process_adjoint_2(): st, meta = prepare_real_adj_data() inv = obspy.read_inventory() event = obspy.read_events()[0] adjs = pa.convert_stream_to_adjs(st, meta) starttime = adjs[0].starttime - 20 delta = adjs[0].dt / 2.0 npts = 40 + 2 * len(adjs[0].adjoint_source) new_adj = pa.process_adjoint( adjs, interp_flag=True, interp_starttime=starttime, interp_delta=delta, interp_npts=npts, sum_over_comp_flag=True, weight_flag=False, filter_flag=True, pre_filt=[0.02, 0.025, 0.059, 0.073], taper_percentage=0.05, taper_type="hann", add_missing_comp_flag=True, rotate_flag=True, inventory=inv, event=event ) assert len(new_adj) == 3