def test_rotate_one_station_stream_synt():
syn = testsyn.copy()
inv = deepcopy(teststaxml)
syn1 = syn.copy()
syn_r1 = rotate.rotate_one_station_stream(syn1,
0.0,
0.0,
inventory=inv,
mode="NE->RT")
syn2 = syn.copy()
syn_r2 = rotate.rotate_one_station_stream(syn2,
0.0,
0.0,
station_latitude=34.5408,
station_longitude=69.0432,
mode="NE->RT")
assert syn_r1 == syn_r2
syn3 = syn.copy()
syn_r3 = rotate.rotate_one_station_stream(syn3,
0.0,
0.0,
inventory=inv,
mode="NE->RT")
assert syn_r3 == syn_r1
def test_rotate_one_station_stream_obsd():
return
obs = testobs.copy()
inv = deepcopy(teststaxml)
obs1 = obs.copy()
obs_r1 = rotate.rotate_one_station_stream(obs1,
0.0,
0.0,
inventory=inv,
mode="NE->RT",
sanity_check=False)
obs2 = obs.copy()
obs_r2 = rotate.rotate_one_station_stream(obs2,
0.0,
0.0,
station_latitude=34.5408,
station_longitude=69.0432,
mode="NE->RT",
sanity_check=False)
assert obs_r1 == obs_r2
obs3 = obs.copy()
obs_r3 = rotate.rotate_one_station_stream(obs3,
0.0,
0.0,
inventory=inv,
mode="NE->RT",
sanity_check=True)
assert len(obs_r3) == 1
assert obs_r3[0].stats.channel == "BHZ"
assert obs_r3[0] == obs_r1.select(channel="BHZ")[0]
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_rotate_one_station_stream_synt():
syn = testsyn.copy()
inv = deepcopy(teststaxml)
syn1 = syn.copy()
syn_r1 = rotate.rotate_one_station_stream(
syn1, 0.0, 0.0, inventory=inv, mode="NE->RT")
syn2 = syn.copy()
syn_r2 = rotate.rotate_one_station_stream(
syn2, 0.0, 0.0, station_latitude=34.5408,
station_longitude=69.0432, mode="NE->RT")
assert syn_r1 == syn_r2
syn3 = syn.copy()
syn_r3 = rotate.rotate_one_station_stream(
syn3, 0.0, 0.0, inventory=inv, mode="NE->RT")
assert syn_r3 == syn_r1
def test_rotate_one_station_stream_obsd():
return
obs = testobs.copy()
inv = deepcopy(teststaxml)
obs1 = obs.copy()
obs_r1 = rotate.rotate_one_station_stream(
obs1, 0.0, 0.0, inventory=inv, mode="NE->RT",
sanity_check=False)
obs2 = obs.copy()
obs_r2 = rotate.rotate_one_station_stream(
obs2, 0.0, 0.0, station_latitude=34.5408,
station_longitude=69.0432, mode="NE->RT",
sanity_check=False)
assert obs_r1 == obs_r2
obs3 = obs.copy()
obs_r3 = rotate.rotate_one_station_stream(
obs3, 0.0, 0.0, inventory=inv, mode="NE->RT",
sanity_check=True)
assert len(obs_r3) == 1
assert obs_r3[0].stats.channel == "BHZ"
assert obs_r3[0] == obs_r1.select(channel="BHZ")[0]