def test_process_synt():
staxmlfile = os.path.join(DATA_DIR, "stationxml", "IU.KBL.syn.xml")
inv = obspy.read_inventory(staxmlfile)
st = testsyn.copy()
event = obspy.read_events(testquakeml)[0]
origin = event.preferred_origin() or event.origins[0]
event_lat = origin.latitude
event_lon = origin.longitude
event_time = origin.time
pre_filt = [1 / 90., 1 / 60., 1 / 27.0, 1 / 22.5]
t1 = event_time
t2 = event_time + 6000.0
st_new = proc.process_stream(st,
remove_response_flag=False,
inventory=inv,
filter_flag=True,
pre_filt=pre_filt,
starttime=t1,
endtime=t2,
resample_flag=True,
sampling_rate=2.0,
taper_type="hann",
taper_percentage=0.05,
rotate_flag=True,
event_latitude=event_lat,
event_longitude=event_lon)
bmfile = os.path.join(DATA_DIR, "proc", "IU.KBL.syn.proc.mseed")
st_compare = obspy.read(bmfile)
assert compare_stream_kernel(st_new, st_compare)
def test_process_obsd_2():
st = testobs.copy()
inv = deepcopy(teststaxml)
event = obspy.read_events(testquakeml)[0]
origin = event.preferred_origin() or event.origins[0]
event_lat = origin.latitude
event_lon = origin.longitude
event_time = origin.time
pre_filt = [1 / 90., 1 / 60., 1 / 27.0, 1 / 22.5]
t1 = event_time
t2 = event_time + 6000.0
st_new = proc.process_stream(st,
remove_response_flag=True,
water_level=60,
inventory=inv,
filter_flag=True,
pre_filt=pre_filt,
starttime=t1,
endtime=t2,
resample_flag=True,
sampling_rate=2.0,
taper_type="hann",
taper_percentage=0.05,
rotate_flag=True,
event_latitude=event_lat,
event_longitude=event_lon,
sanity_check=True)
bmfile = os.path.join(DATA_DIR, "proc", "IU.KBL.obs.proc.mseed")
st_compare = obspy.read(bmfile)
assert len(st_new) == 3
assert compare_trace_kernel(
st_new.select(channel="BHZ")[0],
st_compare.select(channel="BHZ")[0])