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