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