def test_event(self):
"""
Test Event object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
event = catalog[0]
self.assertEqual(event.resource_id, 'smi:local/event/280435')
self.assertEqual(event.event_type, 'earthquake')
self.assertEqual(event.event_type_certainty, "known")
self.assertNotEqual(event.creation_info, None)
self.assertEqual(event.preferred_origin_id, 'smi:local/origin/282672')
self.assertEqual(event.preferred_magnitude_id,
'smi:local/origin/282672/magnitude/0')
# event descriptions
self.assertEqual(len(event.event_descriptions), 1)
event_description = event.event_descriptions[0]
self.assertEqual(
event_description.text,
'GREECE-ALBANIA BORDER REGION')
self.assertEqual(event_description.type, 'region name')
# comments
self.assertEqual(len(event.comments), 1)
comment = event.comments[0]
self.assertEqual(comment.text, 'GSE2.0:evtype=ke')
self.assertEqual(comment.resource_id.id[:10], 'smi:local/')
def test_station_magnitude(self):
"""
Test StationMagnitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
station_magnitudes = catalog[0].station_magnitudes
self.assertEqual(len(station_magnitudes), 4)
# Test first StationMagnitude
sta_mag_1 = station_magnitudes[0]
self.assertEqual(sta_mag_1.resource_id,
'smi:local/magnitude/station/3586432/0')
self.assertEqual(sta_mag_1.origin_id, 'smi:local/origin/282672')
self.assertEqual(sta_mag_1.mag, 4.0)
self.assertEqual(sta_mag_1.station_magnitude_type, 'ML')
self.assertEqual(sta_mag_1.amplitude_id, 'smi:local/amplitude/3586432')
self.assertEqual(sta_mag_1.method_id, None)
self.assertEqual(sta_mag_1.waveform_id, None)
self.assertNotEqual(sta_mag_1.creation_info, None)
self.assertEqual(len(sta_mag_1.comments), 0)
# Test second StationMagnitude
sta_mag_2 = station_magnitudes[1]
self.assertEqual(sta_mag_2.resource_id,
'smi:local/magnitude/station/3586555/0')
self.assertEqual(sta_mag_2.origin_id, 'smi:local/origin/282672')
self.assertEqual(sta_mag_2.mag, 3.7)
self.assertEqual(sta_mag_2.station_magnitude_type, 'mb')
self.assertEqual(sta_mag_2.amplitude_id, 'smi:local/amplitude/3586555')
self.assertEqual(sta_mag_2.method_id, None)
self.assertEqual(sta_mag_2.waveform_id, None)
self.assertNotEqual(sta_mag_2.creation_info, None)
self.assertEqual(len(sta_mag_2.comments), 0)
def test_event(self):
"""
Test Event object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
event = catalog[0]
self.assertEqual(event.resource_id, 'smi:local/event/280435')
self.assertEqual(event.event_type, 'earthquake')
self.assertEqual(event.event_type_certainty, "known")
self.assertNotEqual(event.creation_info, None)
self.assertEqual(event.preferred_origin_id, 'smi:local/origin/282672')
self.assertEqual(event.preferred_magnitude_id,
'smi:local/origin/282672/magnitude/0')
# event descriptions
self.assertEqual(len(event.event_descriptions), 1)
event_description = event.event_descriptions[0]
self.assertEqual(event_description.text,
'GREECE-ALBANIA BORDER REGION')
self.assertEqual(event_description.type, 'region name')
# comments
self.assertEqual(len(event.comments), 1)
comment = event.comments[0]
self.assertEqual(comment.text, 'GSE2.0:evtype=ke')
self.assertEqual(comment.resource_id.id[:10], 'smi:local/')
def test_station_magnitude(self):
"""
Test StationMagnitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
station_magnitudes = catalog[0].station_magnitudes
self.assertEqual(len(station_magnitudes), 4)
# Test first StationMagnitude
sta_mag_1 = station_magnitudes[0]
self.assertEqual(
sta_mag_1.resource_id, 'smi:local/magnitude/station/3586432/0')
self.assertEqual(sta_mag_1.origin_id, 'smi:local/origin/282672')
self.assertEqual(sta_mag_1.mag, 4.0)
self.assertEqual(sta_mag_1.station_magnitude_type, 'ML')
self.assertEqual(sta_mag_1.amplitude_id, 'smi:local/amplitude/3586432')
self.assertEqual(sta_mag_1.method_id, None)
self.assertEqual(sta_mag_1.waveform_id, None)
self.assertNotEqual(sta_mag_1.creation_info, None)
self.assertEqual(len(sta_mag_1.comments), 0)
# Test second StationMagnitude
sta_mag_2 = station_magnitudes[1]
self.assertEqual(
sta_mag_2.resource_id, 'smi:local/magnitude/station/3586555/0')
self.assertEqual(sta_mag_2.origin_id, 'smi:local/origin/282672')
self.assertEqual(sta_mag_2.mag, 3.7)
self.assertEqual(sta_mag_2.station_magnitude_type, 'mb')
self.assertEqual(sta_mag_2.amplitude_id, 'smi:local/amplitude/3586555')
self.assertEqual(sta_mag_2.method_id, None)
self.assertEqual(sta_mag_2.waveform_id, None)
self.assertNotEqual(sta_mag_2.creation_info, None)
self.assertEqual(len(sta_mag_2.comments), 0)
def test_several_begin(self):
"""
Test with several events.
"""
filename = os.path.join(self.path, 'gse_2.0_2_begins.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 2)
# Test firt event
event_1 = catalog[0]
self.assertEqual(event_1.resource_id, 'smi:local/event/280435')
self.assertEqual(len(event_1.event_descriptions), 1)
self.assertEqual(event_1.event_descriptions[0].text,
'GREECE-ALBANIA BORDER REGION')
self.assertEqual(len(event_1.comments), 1)
self.assertEqual(len(event_1.picks), 9)
self.assertEqual(len(event_1.amplitudes), 6)
self.assertEqual(len(event_1.origins), 1)
self.assertEqual(len(event_1.magnitudes), 2)
self.assertEqual(len(event_1.station_magnitudes), 4)
# Test second event
event_2 = catalog[1]
self.assertEqual(event_2.resource_id, 'smi:local/event/280436')
self.assertEqual(len(event_2.event_descriptions), 1)
self.assertEqual(event_2.event_descriptions[0].text,
'VANCOUVER ISLAND REGION')
self.assertEqual(len(event_2.comments), 1)
self.assertEqual(len(event_2.picks), 7)
self.assertEqual(len(event_2.amplitudes), 5)
self.assertEqual(len(event_2.origins), 1)
self.assertEqual(len(event_2.magnitudes), 1)
self.assertEqual(len(event_2.station_magnitudes), 2)
def test_several_begin(self):
"""
Test with several events.
"""
filename = os.path.join(self.path, 'gse_2.0_2_begins.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 2)
# Test firt event
event_1 = catalog[0]
self.assertEqual(event_1.resource_id, 'smi:local/event/280435')
self.assertEqual(len(event_1.event_descriptions), 1)
self.assertEqual(
event_1.event_descriptions[0].text,
'GREECE-ALBANIA BORDER REGION')
self.assertEqual(len(event_1.comments), 1)
self.assertEqual(len(event_1.picks), 9)
self.assertEqual(len(event_1.amplitudes), 6)
self.assertEqual(len(event_1.origins), 1)
self.assertEqual(len(event_1.magnitudes), 2)
self.assertEqual(len(event_1.station_magnitudes), 4)
# Test second event
event_2 = catalog[1]
self.assertEqual(event_2.resource_id, 'smi:local/event/280436')
self.assertEqual(len(event_2.event_descriptions), 1)
self.assertEqual(
event_2.event_descriptions[0].text,
'VANCOUVER ISLAND REGION')
self.assertEqual(len(event_2.comments), 1)
self.assertEqual(len(event_2.picks), 7)
self.assertEqual(len(event_2.amplitudes), 5)
self.assertEqual(len(event_2.origins), 1)
self.assertEqual(len(event_2.magnitudes), 1)
self.assertEqual(len(event_2.station_magnitudes), 2)
def test_origin(self):
"""
Test Origin object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
self.assertEqual(len(catalog[0].origins), 1)
origin = catalog[0].origins[0]
self.assertEqual(origin.resource_id, 'smi:local/origin/282672')
self.assertEqual(origin.time, UTCDateTime('1995-01-16T07:26:52.4'))
self.assertEqual(origin.time_errors.uncertainty, 12.69)
self.assertEqual(origin.latitude, 39.45)
self.assertEqual(origin.latitude_errors.uncertainty, None)
self.assertEqual(origin.longitude, 20.44)
self.assertEqual(origin.longitude_errors.uncertainty, None)
self.assertEqual(origin.depth, 66800)
self.assertEqual(origin.depth_errors.uncertainty, 83800)
self.assertEqual(origin.depth_type, OriginDepthType('from location'))
self.assertEqual(origin.time_fixed, False)
self.assertEqual(origin.epicenter_fixed, False)
self.assertEqual(origin.reference_system_id, None)
self.assertEqual(origin.method_id, 'smi:local/method/inversion')
self.assertEqual(origin.earth_model_id, None)
self.assertEqual(origin.origin_type, None)
self.assertEqual(origin.region, None)
self.assertEqual(origin.evaluation_mode, "manual")
self.assertEqual(origin.evaluation_status, None)
# quality
self.assertEqual(origin.quality.associated_phase_count, 9)
self.assertEqual(origin.quality.used_phase_count, 9)
self.assertEqual(origin.quality.associated_station_count, 8)
self.assertEqual(origin.quality.used_station_count, 8)
self.assertEqual(origin.quality.depth_phase_count, None)
self.assertEqual(origin.quality.standard_error, 0.53)
self.assertEqual(origin.quality.azimuthal_gap, 322)
self.assertEqual(origin.quality.secondary_azimuthal_gap, None)
self.assertEqual(origin.quality.ground_truth_level, None)
self.assertEqual(origin.quality.minimum_distance, 10.56)
self.assertEqual(origin.quality.maximum_distance, 78.21)
self.assertEqual(origin.quality.median_distance, None)
# origin uncertainty
u = origin.origin_uncertainty
self.assertEqual(u.horizontal_uncertainty, None)
self.assertEqual(u.min_horizontal_uncertainty, 83700)
self.assertEqual(u.max_horizontal_uncertainty, 93600)
self.assertEqual(u.azimuth_max_horizontal_uncertainty, 27)
self.assertEqual(u.confidence_ellipsoid, None)
self.assertEqual(u.preferred_description, 'uncertainty ellipse')
self.assertEqual(u.confidence_level, None)
# creation info
self.assertEqual(origin.creation_info.author, 'GSE_IDC')
# comments
self.assertEqual(len(origin.comments), 0)
# composite times
self.assertEqual(len(origin.composite_times), 0)
def test_origin(self):
"""
Test Origin object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
self.assertEqual(len(catalog[0].origins), 1)
origin = catalog[0].origins[0]
self.assertEqual(origin.resource_id, 'smi:local/origin/282672')
self.assertEqual(origin.time, UTCDateTime('1995-01-16T07:26:52.4'))
self.assertEqual(origin.time_errors.uncertainty, 12.69)
self.assertEqual(origin.latitude, 39.45)
self.assertEqual(origin.latitude_errors.uncertainty, None)
self.assertEqual(origin.longitude, 20.44)
self.assertEqual(origin.longitude_errors.uncertainty, None)
self.assertEqual(origin.depth, 66800)
self.assertEqual(origin.depth_errors.uncertainty, 83800)
self.assertEqual(origin.depth_type, OriginDepthType('from location'))
self.assertEqual(origin.time_fixed, False)
self.assertEqual(origin.epicenter_fixed, False)
self.assertEqual(origin.reference_system_id, None)
self.assertEqual(origin.method_id, 'smi:local/method/inversion')
self.assertEqual(origin.earth_model_id, None)
self.assertEqual(origin.origin_type, None)
self.assertEqual(origin.region, None)
self.assertEqual(origin.evaluation_mode, "manual")
self.assertEqual(origin.evaluation_status, None)
# quality
self.assertEqual(origin.quality.associated_phase_count, 9)
self.assertEqual(origin.quality.used_phase_count, 9)
self.assertEqual(origin.quality.associated_station_count, 8)
self.assertEqual(origin.quality.used_station_count, 8)
self.assertEqual(origin.quality.depth_phase_count, None)
self.assertEqual(origin.quality.standard_error, 0.53)
self.assertEqual(origin.quality.azimuthal_gap, 322)
self.assertEqual(origin.quality.secondary_azimuthal_gap, None)
self.assertEqual(origin.quality.ground_truth_level, None)
self.assertEqual(origin.quality.minimum_distance, 10.56)
self.assertEqual(origin.quality.maximum_distance, 78.21)
self.assertEqual(origin.quality.median_distance, None)
# origin uncertainty
u = origin.origin_uncertainty
self.assertEqual(u.horizontal_uncertainty, None)
self.assertEqual(u.min_horizontal_uncertainty, 83700)
self.assertEqual(u.max_horizontal_uncertainty, 93600)
self.assertEqual(u.azimuth_max_horizontal_uncertainty, 27)
self.assertEqual(u.confidence_ellipsoid, None)
self.assertEqual(u.preferred_description, 'uncertainty ellipse')
self.assertEqual(u.confidence_level, None)
# creation info
self.assertEqual(origin.creation_info.author, 'GSE_IDC')
# comments
self.assertEqual(len(origin.comments), 0)
# composite times
self.assertEqual(len(origin.composite_times), 0)
def test_pick(self):
"""
Test Pick object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
picks = catalog[0].picks
self.assertEqual(len(picks), 9)
# Test first Pick
pick_1 = picks[0]
self.assertEqual(pick_1.resource_id, 'smi:local/pick/3586432')
self.assertEqual(pick_1.time, UTCDateTime('1995-01-16T07:29:20.7'))
# WaveformStreamId
waveform_1 = pick_1.waveform_id
self.assertEqual(waveform_1.network_code, 'XX')
self.assertEqual(waveform_1.station_code, 'GERES')
self.assertEqual(waveform_1.channel_code, None)
self.assertEqual(waveform_1.location_code, None)
self.assertEqual(waveform_1.resource_uri, None)
self.assertEqual(pick_1.filter_id, None)
self.assertEqual(pick_1.method_id, None)
self.assertEqual(pick_1.horizontal_slowness, 13.8)
self.assertEqual(pick_1.backazimuth, 163.7)
self.assertEqual(pick_1.slowness_method_id, None)
self.assertEqual(pick_1.onset, 'emergent')
self.assertEqual(pick_1.phase_hint, 'P')
self.assertEqual(pick_1.polarity, None)
self.assertEqual(pick_1.evaluation_mode, 'manual')
self.assertEqual(pick_1.evaluation_status, None)
self.assertNotEqual(pick_1.creation_info, None)
self.assertEqual(len(pick_1.comments), 0)
# Test second Pick
pick_2 = picks[1]
self.assertEqual(pick_2.resource_id, 'smi:local/pick/3586513')
self.assertEqual(pick_2.time, UTCDateTime('1995-01-16T07:31:17.5'))
# WaveformStreamId
waveform_2 = pick_2.waveform_id
self.assertEqual(waveform_2.network_code, 'XX')
self.assertEqual(waveform_2.station_code, 'GERES')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
self.assertEqual(waveform_2.resource_uri, None)
self.assertEqual(pick_2.filter_id, None)
self.assertEqual(pick_2.method_id, None)
self.assertEqual(pick_2.horizontal_slowness, 23.4)
self.assertEqual(pick_2.backazimuth, 153.4)
self.assertEqual(pick_2.slowness_method_id, None)
self.assertEqual(pick_2.onset, None)
self.assertEqual(pick_2.phase_hint, 'S')
self.assertEqual(pick_2.polarity, PickPolarity.POSITIVE)
self.assertEqual(pick_2.evaluation_mode, None)
self.assertEqual(pick_2.evaluation_status, None)
self.assertNotEqual(pick_2.creation_info, None)
self.assertEqual(len(pick_2.comments), 0)
def test_pick(self):
"""
Test Pick object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
picks = catalog[0].picks
self.assertEqual(len(picks), 9)
# Test first Pick
pick_1 = picks[0]
self.assertEqual(pick_1.resource_id, 'smi:local/pick/3586432')
self.assertEqual(pick_1.time, UTCDateTime('1995-01-16T07:29:20.7'))
# WaveformStreamId
waveform_1 = pick_1.waveform_id
self.assertEqual(waveform_1.network_code, 'XX')
self.assertEqual(waveform_1.station_code, 'GERES')
self.assertEqual(waveform_1.channel_code, None)
self.assertEqual(waveform_1.location_code, None)
self.assertEqual(waveform_1.resource_uri, None)
self.assertEqual(pick_1.filter_id, None)
self.assertEqual(pick_1.method_id, None)
self.assertEqual(pick_1.horizontal_slowness, 13.8)
self.assertEqual(pick_1.backazimuth, 163.7)
self.assertEqual(pick_1.slowness_method_id, None)
self.assertEqual(pick_1.onset, 'emergent')
self.assertEqual(pick_1.phase_hint, 'P')
self.assertEqual(pick_1.polarity, None)
self.assertEqual(pick_1.evaluation_mode, 'manual')
self.assertEqual(pick_1.evaluation_status, None)
self.assertNotEqual(pick_1.creation_info, None)
self.assertEqual(len(pick_1.comments), 0)
# Test second Pick
pick_2 = picks[1]
self.assertEqual(pick_2.resource_id, 'smi:local/pick/3586513')
self.assertEqual(pick_2.time, UTCDateTime('1995-01-16T07:31:17.5'))
# WaveformStreamId
waveform_2 = pick_2.waveform_id
self.assertEqual(waveform_2.network_code, 'XX')
self.assertEqual(waveform_2.station_code, 'GERES')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
self.assertEqual(waveform_2.resource_uri, None)
self.assertEqual(pick_2.filter_id, None)
self.assertEqual(pick_2.method_id, None)
self.assertEqual(pick_2.horizontal_slowness, 23.4)
self.assertEqual(pick_2.backazimuth, 153.4)
self.assertEqual(pick_2.slowness_method_id, None)
self.assertEqual(pick_2.onset, None)
self.assertEqual(pick_2.phase_hint, 'S')
self.assertEqual(pick_2.polarity, PickPolarity.POSITIVE)
self.assertEqual(pick_2.evaluation_mode, None)
self.assertEqual(pick_2.evaluation_status, None)
self.assertNotEqual(pick_2.creation_info, None)
self.assertEqual(len(pick_2.comments), 0)
def test_catalog(self):
"""
Test Catalog object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(catalog.resource_id, 'smi:local/event/evid')
self.assertEqual(len(catalog.comments), 1)
comment = catalog.comments[0]
text = ("Reviewed Event Bulletin (REB) of the GSE_IDC for January 16, "
"1995")
self.assertEqual(comment.text, text)
self.assertEqual(comment.resource_id.id[:10], 'smi:local/')
def test_catalog(self):
"""
Test Catalog object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(catalog.resource_id, 'smi:local/event/evid')
self.assertEqual(len(catalog.comments), 1)
comment = catalog.comments[0]
text = ("Reviewed Event Bulletin (REB) of the GSE_IDC for January 16, "
"1995")
self.assertEqual(comment.text, text)
self.assertEqual(comment.resource_id.id[:10], 'smi:local/')
def test_arrival(self):
"""
Test Arrival object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
self.assertEqual(len(catalog[0].origins), 1)
arrivals = catalog[0].origins[0].arrivals
self.assertEqual(len(arrivals), 9)
# Test first Arrival
arrival_1 = arrivals[0]
self.assertEqual(arrival_1.resource_id,
'smi:local/origin/282672/arrival/3586432')
self.assertEqual(arrival_1.pick_id, 'smi:local/pick/3586432')
self.assertEqual(arrival_1.phase, 'P')
self.assertEqual(arrival_1.time_correction, None)
self.assertEqual(arrival_1.azimuth, 150.3)
self.assertEqual(arrival_1.distance, 10.56)
self.assertEqual(arrival_1.takeoff_angle, None)
self.assertEqual(arrival_1.time_residual, -0.2)
self.assertEqual(arrival_1.horizontal_slowness_residual, 0.1)
self.assertEqual(arrival_1.backazimuth_residual, 13.4)
self.assertEqual(arrival_1.time_weight, 1)
self.assertEqual(arrival_1.backazimuth_weight, None)
self.assertEqual(arrival_1.horizontal_slowness_weight, None)
self.assertEqual(arrival_1.earth_model_id, None)
self.assertNotEqual(arrival_1.creation_info, None)
self.assertEqual(len(arrival_1.comments), 0)
# Test second Arrival
arrival_2 = arrivals[1]
self.assertEqual(arrival_2.resource_id,
'smi:local/origin/282672/arrival/3586513')
self.assertEqual(arrival_2.pick_id, 'smi:local/pick/3586513')
self.assertEqual(arrival_2.phase, 'S')
self.assertEqual(arrival_2.time_correction, None)
self.assertEqual(arrival_2.azimuth, 150.3)
self.assertEqual(arrival_2.distance, 10.56)
self.assertEqual(arrival_2.takeoff_angle, None)
self.assertEqual(arrival_2.time_residual, -0.6)
self.assertEqual(arrival_2.horizontal_slowness_residual, -1.0)
self.assertEqual(arrival_2.backazimuth_residual, 3.1)
self.assertEqual(arrival_2.time_weight, None)
self.assertEqual(arrival_2.backazimuth_weight, 1)
self.assertEqual(arrival_2.horizontal_slowness_weight, None)
self.assertEqual(arrival_2.earth_model_id, None)
self.assertNotEqual(arrival_2.creation_info, None)
self.assertEqual(len(arrival_2.comments), 0)
def test_arrival(self):
"""
Test Arrival object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
self.assertEqual(len(catalog[0].origins), 1)
arrivals = catalog[0].origins[0].arrivals
self.assertEqual(len(arrivals), 9)
# Test first Arrival
arrival_1 = arrivals[0]
self.assertEqual(
arrival_1.resource_id, 'smi:local/origin/282672/arrival/3586432')
self.assertEqual(arrival_1.pick_id, 'smi:local/pick/3586432')
self.assertEqual(arrival_1.phase, 'P')
self.assertEqual(arrival_1.time_correction, None)
self.assertEqual(arrival_1.azimuth, 150.3)
self.assertEqual(arrival_1.distance, 10.56)
self.assertEqual(arrival_1.takeoff_angle, None)
self.assertEqual(arrival_1.time_residual, -0.2)
self.assertEqual(arrival_1.horizontal_slowness_residual, 0.1)
self.assertEqual(arrival_1.backazimuth_residual, 13.4)
self.assertEqual(arrival_1.time_weight, 1)
self.assertEqual(arrival_1.backazimuth_weight, None)
self.assertEqual(arrival_1.horizontal_slowness_weight, None)
self.assertEqual(arrival_1.earth_model_id, None)
self.assertNotEqual(arrival_1.creation_info, None)
self.assertEqual(len(arrival_1.comments), 0)
# Test second Arrival
arrival_2 = arrivals[1]
self.assertEqual(
arrival_2.resource_id, 'smi:local/origin/282672/arrival/3586513')
self.assertEqual(arrival_2.pick_id, 'smi:local/pick/3586513')
self.assertEqual(arrival_2.phase, 'S')
self.assertEqual(arrival_2.time_correction, None)
self.assertEqual(arrival_2.azimuth, 150.3)
self.assertEqual(arrival_2.distance, 10.56)
self.assertEqual(arrival_2.takeoff_angle, None)
self.assertEqual(arrival_2.time_residual, -0.6)
self.assertEqual(arrival_2.horizontal_slowness_residual, -1.0)
self.assertEqual(arrival_2.backazimuth_residual, 3.1)
self.assertEqual(arrival_2.time_weight, None)
self.assertEqual(arrival_2.backazimuth_weight, 1)
self.assertEqual(arrival_2.horizontal_slowness_weight, None)
self.assertEqual(arrival_2.earth_model_id, None)
self.assertNotEqual(arrival_2.creation_info, None)
self.assertEqual(len(arrival_2.comments), 0)
def test_inventory_with_multiple_channels(self):
filename = os.path.join(self.path, 'gse_2.0_non_standard.txt')
inventory_filename = os.path.join(self.path,
'inventory_multiple_channels.xml')
inventory = read_inventory(inventory_filename)
fields = {
'line_1': {
'author': slice(105, 113),
'id': slice(114, 123),
},
'line_2': {
'az': slice(40, 46),
'antype': slice(105, 106),
'loctype': slice(107, 108),
'evtype': slice(109, 111),
},
'arrival': {
'amp': slice(94, 104),
},
}
warnings.simplefilter("ignore", UserWarning)
catalog = _read_gse2(filename,
inventory,
fields=fields,
event_point_separator=True)
warnings.filters.pop(0)
self.assertEqual(len(catalog), 2)
# Test a station present in the inventory
event = catalog[0]
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
waveform = pick.waveform_id
self.assertEqual(waveform.network_code, 'ZU')
self.assertEqual(waveform.channel_code, 'HHZ')
self.assertEqual(waveform.location_code, '1')
# Test a station with several channels
pick_2 = event.picks[2]
waveform_2 = pick_2.waveform_id
self.assertEqual(waveform_2.network_code, 'ZU')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
# Test a station not present in the inventory
pick_3 = event.picks[3]
waveform_3 = pick_3.waveform_id
self.assertEqual(waveform_3.network_code, 'XX')
self.assertEqual(waveform_3.channel_code, None)
self.assertEqual(waveform_3.location_code, None)
def test_non_standard_format(self):
"""
Test non-standard GSE2 format which can normally be parsed too.
"""
filename = os.path.join(self.path, 'gse_2.0_non_standard.txt')
fields = {
'line_1': {
'author': slice(105, 113),
'id': slice(114, 123),
},
'line_2': {
'az': slice(40, 46),
'antype': slice(105, 106),
'loctype': slice(107, 108),
'evtype': slice(109, 111),
},
'arrival': {
'amp': slice(94, 104),
},
}
warnings.simplefilter("ignore", UserWarning)
catalog = _read_gse2(filename,
fields=fields,
event_point_separator=True)
warnings.filters.pop(0)
self.assertEqual(len(catalog), 2)
event = catalog[0]
self.assertEqual(len(event.origins), 1)
origin = event.origins[0]
# Test fields with non-standard indexes
self.assertEqual(origin.creation_info.author, 'bulletin_test')
self.assertEqual(origin.resource_id, 'smi:local/origin/282672')
self.assertEqual(origin.quality.standard_error, 0.55)
uncertainty = origin.origin_uncertainty
self.assertEqual(uncertainty.azimuth_max_horizontal_uncertainty, 27.05)
self.assertEqual(origin.evaluation_mode, "manual")
self.assertEqual(origin.method_id, 'smi:local/method/inversion')
self.assertEqual(event.event_type, 'earthquake')
self.assertEqual(event.event_type_certainty, "known")
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
self.assertEqual(pick.phase_hint, 'Pg')
self.assertEqual(len(event.amplitudes), 4)
amplitude = event.amplitudes[0]
self.assertEqual(amplitude.generic_amplitude, 2.9)
def test_inventory_with_multiple_channels(self):
filename = os.path.join(self.path, 'gse_2.0_non_standard.txt')
inventory_filename = os.path.join(self.path,
'inventory_multiple_channels.xml')
inventory = read_inventory(inventory_filename)
fields = {
'line_1': {
'author': slice(105, 113),
'id': slice(114, 123),
},
'line_2': {
'az': slice(40, 46),
'antype': slice(105, 106),
'loctype': slice(107, 108),
'evtype': slice(109, 111),
},
'arrival': {
'amp': slice(94, 104),
},
}
warnings.simplefilter("ignore", UserWarning)
catalog = _read_gse2(filename, inventory, fields=fields,
event_point_separator=True)
warnings.filters.pop(0)
self.assertEqual(len(catalog), 2)
# Test a station present in the inventory
event = catalog[0]
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
waveform = pick.waveform_id
self.assertEqual(waveform.network_code, 'ZU')
self.assertEqual(waveform.channel_code, 'HHZ')
self.assertEqual(waveform.location_code, '1')
# Test a station with several channels
pick_2 = event.picks[2]
waveform_2 = pick_2.waveform_id
self.assertEqual(waveform_2.network_code, 'ZU')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
# Test a station not present in the inventory
pick_3 = event.picks[3]
waveform_3 = pick_3.waveform_id
self.assertEqual(waveform_3.network_code, 'XX')
self.assertEqual(waveform_3.channel_code, None)
self.assertEqual(waveform_3.location_code, None)
def test_non_standard_format(self):
"""
Test non-standard GSE2 format which can normally be parsed too.
"""
filename = os.path.join(self.path, 'gse_2.0_non_standard.txt')
fields = {
'line_1': {
'author': slice(105, 113),
'id': slice(114, 123),
},
'line_2': {
'az': slice(40, 46),
'antype': slice(105, 106),
'loctype': slice(107, 108),
'evtype': slice(109, 111),
},
'arrival': {
'amp': slice(94, 104),
},
}
warnings.simplefilter("ignore", UserWarning)
catalog = _read_gse2(filename, fields=fields,
event_point_separator=True)
warnings.filters.pop(0)
self.assertEqual(len(catalog), 2)
event = catalog[0]
self.assertEqual(len(event.origins), 1)
origin = event.origins[0]
# Test fields with non-standard indexes
self.assertEqual(origin.creation_info.author, 'bulletin_test')
self.assertEqual(origin.resource_id, 'smi:local/origin/282672')
self.assertEqual(origin.quality.standard_error, 0.55)
uncertainty = origin.origin_uncertainty
self.assertEqual(uncertainty.azimuth_max_horizontal_uncertainty, 27.05)
self.assertEqual(origin.evaluation_mode, "manual")
self.assertEqual(origin.method_id, 'smi:local/method/inversion')
self.assertEqual(event.event_type, 'earthquake')
self.assertEqual(event.event_type_certainty, "known")
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
self.assertEqual(pick.phase_hint, 'Pg')
self.assertEqual(len(event.amplitudes), 4)
amplitude = event.amplitudes[0]
self.assertEqual(amplitude.generic_amplitude, 2.9)
def test_magnitude(self):
"""
Test Magnitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
magnitudes = catalog[0].magnitudes
self.assertEqual(len(magnitudes), 2)
# Test first Magnitude
mag_1 = magnitudes[0]
self.assertEqual(
mag_1.resource_id, 'smi:local/origin/282672/magnitude/0')
self.assertEqual(mag_1.mag, 3.6)
self.assertEqual(mag_1.mag_errors.uncertainty, 0.2)
self.assertEqual(mag_1.magnitude_type, 'mb')
self.assertEqual(mag_1.origin_id, 'smi:local/origin/282672')
self.assertEqual(mag_1.method_id, None)
self.assertEqual(mag_1.station_count, 3)
self.assertEqual(mag_1.azimuthal_gap, None)
self.assertEqual(mag_1.evaluation_mode, None)
self.assertEqual(mag_1.evaluation_status, None)
self.assertNotEqual(mag_1.creation_info, None)
self.assertEqual(len(mag_1.comments), 0)
self.assertEqual(len(mag_1.station_magnitude_contributions), 0)
# Test second Magnitude
mag_2 = magnitudes[1]
self.assertEqual(
mag_2.resource_id, 'smi:local/origin/282672/magnitude/1')
self.assertEqual(mag_2.mag, 4.0)
self.assertEqual(mag_2.mag_errors.uncertainty, None)
self.assertEqual(mag_2.magnitude_type, 'ML')
self.assertEqual(mag_2.origin_id, 'smi:local/origin/282672')
self.assertEqual(mag_2.method_id, None)
self.assertEqual(mag_2.station_count, 1)
self.assertEqual(mag_2.azimuthal_gap, None)
self.assertEqual(mag_2.evaluation_mode, None)
self.assertEqual(mag_2.evaluation_status, None)
self.assertNotEqual(mag_2.creation_info, None)
self.assertEqual(len(mag_1.comments), 0)
self.assertEqual(len(mag_1.station_magnitude_contributions), 0)
def test_magnitude(self):
"""
Test Magnitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
magnitudes = catalog[0].magnitudes
self.assertEqual(len(magnitudes), 2)
# Test first Magnitude
mag_1 = magnitudes[0]
self.assertEqual(mag_1.resource_id,
'smi:local/origin/282672/magnitude/0')
self.assertEqual(mag_1.mag, 3.6)
self.assertEqual(mag_1.mag_errors.uncertainty, 0.2)
self.assertEqual(mag_1.magnitude_type, 'mb')
self.assertEqual(mag_1.origin_id, 'smi:local/origin/282672')
self.assertEqual(mag_1.method_id, None)
self.assertEqual(mag_1.station_count, 3)
self.assertEqual(mag_1.azimuthal_gap, None)
self.assertEqual(mag_1.evaluation_mode, None)
self.assertEqual(mag_1.evaluation_status, None)
self.assertNotEqual(mag_1.creation_info, None)
self.assertEqual(len(mag_1.comments), 0)
self.assertEqual(len(mag_1.station_magnitude_contributions), 0)
# Test second Magnitude
mag_2 = magnitudes[1]
self.assertEqual(mag_2.resource_id,
'smi:local/origin/282672/magnitude/1')
self.assertEqual(mag_2.mag, 4.0)
self.assertEqual(mag_2.mag_errors.uncertainty, None)
self.assertEqual(mag_2.magnitude_type, 'ML')
self.assertEqual(mag_2.origin_id, 'smi:local/origin/282672')
self.assertEqual(mag_2.method_id, None)
self.assertEqual(mag_2.station_count, 1)
self.assertEqual(mag_2.azimuthal_gap, None)
self.assertEqual(mag_2.evaluation_mode, None)
self.assertEqual(mag_2.evaluation_status, None)
self.assertNotEqual(mag_2.creation_info, None)
self.assertEqual(len(mag_1.comments), 0)
self.assertEqual(len(mag_1.station_magnitude_contributions), 0)
def test_parameters(self):
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
fields = {
'line_1': {
'time': slice(0, 21),
'time_fixf': slice(22, 23),
'lat': slice(25, 33),
'lon': slice(34, 43),
'epicenter_fixf': slice(44, 45),
'depth': slice(47, 52),
'depth_fixf': slice(53, 54),
'n_def': slice(56, 60),
'n_sta': slice(61, 65),
'gap': slice(66, 69),
'mag_type_1': slice(71, 73),
'mag_1': slice(73, 77),
'mag_n_sta_1': slice(78, 80),
'mag_type_2': slice(82, 84),
'mag_2': slice(84, 88),
'mag_n_sta_2': slice(89, 91),
'mag_type_3': slice(93, 95),
'mag_3': slice(95, 99),
'mag_n_sta_3': slice(100, 102),
'author': slice(104, 112),
'id': slice(114, 122),
},
'line_2': {
'rms': slice(5, 10),
'ot_error': slice(15, 21),
's_major': slice(25, 31),
's_minor': slice(32, 38),
'az': slice(40, 43),
'depth_err': slice(49, 54),
'min_dist': slice(56, 62),
'max_dist': slice(63, 69),
'mag_err_1': slice(74, 77),
'mag_err_2': slice(85, 88),
'mag_err_3': slice(96, 99),
'antype': slice(104, 105),
'loctype': slice(106, 107),
'evtype': slice(108, 110),
},
'arrival': {
'sta': slice(0, 5),
'dist': slice(6, 12),
'ev_az': slice(13, 18),
'picktype': slice(19, 20),
'direction': slice(20, 21),
'detchar': slice(21, 22),
'phase': slice(23, 30),
'time': slice(31, 52),
't_res': slice(53, 58),
'azim': slice(59, 64),
'az_res': slice(65, 71),
'slow': slice(72, 77),
's_res': slice(78, 83),
't_def': slice(84, 85),
'a_def': slice(85, 86),
's_def': slice(86, 87),
'snr': slice(88, 93),
'amp': slice(94, 103),
'per': slice(104, 109),
'mag_type_1': slice(110, 112),
'mag_1': slice(112, 116),
'mag_type_2': slice(117, 119),
'mag_2': slice(119, 123),
'id': slice(124, 132),
},
}
catalog = _read_gse2(filename, None, 'TZ', '00', 'SHZ', 'quakeml:idc',
fields, False, 'AGE')
self.assertEqual(len(catalog), 1)
# Test Comment ResourceIdentifier
self.assertEqual(len(catalog.comments), 1)
comment = catalog.comments[0]
self.assertEqual(comment.resource_id.id[:12], 'quakeml:idc/')
# Test Event ResourceIdentifier
event = catalog[0]
self.assertEqual(event.resource_id, 'quakeml:idc/event/280435')
self.assertEqual(event.creation_info.agency_id, 'AGE')
# Test Origin ResourceIdentifier
self.assertEqual(len(event.origins), 1)
origin = event.origins[0]
self.assertEqual(origin.resource_id, 'quakeml:idc/origin/282672')
self.assertEqual(origin.method_id, 'quakeml:idc/method/inversion')
# Test Pick ResourceIdentifier
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
self.assertEqual(pick.resource_id, 'quakeml:idc/pick/3586432')
# Test Arrival ResourceIdentifier
self.assertEqual(len(origin.arrivals), 9)
arrival = origin.arrivals[0]
self.assertEqual(arrival.resource_id,
'quakeml:idc/origin/282672/arrival/3586432')
# Test Magnitude ResourceIdentifier
self.assertEqual(len(event.magnitudes), 2)
magnitude = event.magnitudes[0]
self.assertEqual(magnitude.resource_id,
'quakeml:idc/origin/282672/magnitude/0')
# Test StationMagnitude ResourceIdentifier
self.assertEqual(len(event.station_magnitudes), 4)
sta_mag = event.station_magnitudes[0]
self.assertEqual(sta_mag.resource_id,
'quakeml:idc/magnitude/station/3586432/0')
# Test Amplitude ResourceIdentifier
self.assertEqual(len(event.amplitudes), 6)
amplitude = event.amplitudes[0]
self.assertEqual(amplitude.resource_id,
'quakeml:idc/amplitude/3586432')
# Test network code
waveform = pick.waveform_id
self.assertEqual(waveform.network_code, 'TZ')
self.assertEqual(waveform.channel_code, 'SHZ')
self.assertEqual(waveform.location_code, '00')
def test_amplitude(self):
"""
Test Amplitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
amplitudes = catalog[0].amplitudes
self.assertEqual(len(amplitudes), 6)
# Test first amplitude
amplitude_1 = amplitudes[0]
self.assertEqual(amplitude_1.resource_id,
'smi:local/amplitude/3586432')
self.assertEqual(amplitude_1.generic_amplitude, 0.6)
self.assertEqual(amplitude_1.type, None)
self.assertEqual(amplitude_1.category, None)
self.assertEqual(amplitude_1.unit, None)
self.assertEqual(amplitude_1.method_id, None)
self.assertEqual(amplitude_1.period, 0.3)
self.assertEqual(amplitude_1.snr, 6.8)
self.assertEqual(amplitude_1.time_window, None)
self.assertEqual(amplitude_1.pick_id, 'smi:local/pick/3586432')
# WaveformStreamId
waveform_1 = amplitude_1.waveform_id
self.assertEqual(waveform_1.network_code, 'XX')
self.assertEqual(waveform_1.station_code, 'GERES')
self.assertEqual(waveform_1.channel_code, None)
self.assertEqual(waveform_1.location_code, None)
self.assertEqual(waveform_1.resource_uri, None)
self.assertEqual(amplitude_1.filter_id, None)
self.assertEqual(amplitude_1.scaling_time, None)
self.assertEqual(amplitude_1.magnitude_hint, 'ML')
self.assertEqual(amplitude_1.evaluation_mode, None)
self.assertEqual(amplitude_1.evaluation_status, None)
self.assertNotEqual(amplitude_1.creation_info, None)
self.assertEqual(len(amplitude_1.comments), 0)
# Test second amplitude
amplitude_2 = amplitudes[1]
self.assertEqual(amplitude_2.resource_id,
'smi:local/amplitude/3586513')
self.assertEqual(amplitude_2.generic_amplitude, 2.9)
self.assertEqual(amplitude_2.type, None)
self.assertEqual(amplitude_2.category, None)
self.assertEqual(amplitude_2.unit, None)
self.assertEqual(amplitude_2.method_id, None)
self.assertEqual(amplitude_2.period, 0.6)
self.assertEqual(amplitude_2.snr, 4.9)
self.assertEqual(amplitude_2.time_window, None)
self.assertEqual(amplitude_2.pick_id, 'smi:local/pick/3586513')
# WaveformStreamId
waveform_2 = amplitude_2.waveform_id
self.assertEqual(waveform_2.network_code, 'XX')
self.assertEqual(waveform_2.station_code, 'GERES')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
self.assertEqual(waveform_2.resource_uri, None)
self.assertEqual(amplitude_2.filter_id, None)
self.assertEqual(amplitude_2.scaling_time, None)
self.assertEqual(amplitude_2.magnitude_hint, None)
self.assertEqual(amplitude_2.evaluation_mode, None)
self.assertEqual(amplitude_2.evaluation_status, None)
self.assertNotEqual(amplitude_2.creation_info, None)
self.assertEqual(len(amplitude_2.comments), 0)
def test_amplitude(self):
"""
Test Amplitude object.
"""
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
catalog = _read_gse2(filename)
self.assertEqual(len(catalog), 1)
amplitudes = catalog[0].amplitudes
self.assertEqual(len(amplitudes), 6)
# Test first amplitude
amplitude_1 = amplitudes[0]
self.assertEqual(
amplitude_1.resource_id, 'smi:local/amplitude/3586432')
self.assertEqual(amplitude_1.generic_amplitude, 0.6)
self.assertEqual(amplitude_1.type, None)
self.assertEqual(amplitude_1.category, None)
self.assertEqual(amplitude_1.unit, None)
self.assertEqual(amplitude_1.method_id, None)
self.assertEqual(amplitude_1.period, 0.3)
self.assertEqual(amplitude_1.snr, 6.8)
self.assertEqual(amplitude_1.time_window, None)
self.assertEqual(amplitude_1.pick_id, 'smi:local/pick/3586432')
# WaveformStreamId
waveform_1 = amplitude_1.waveform_id
self.assertEqual(waveform_1.network_code, 'XX')
self.assertEqual(waveform_1.station_code, 'GERES')
self.assertEqual(waveform_1.channel_code, None)
self.assertEqual(waveform_1.location_code, None)
self.assertEqual(waveform_1.resource_uri, None)
self.assertEqual(amplitude_1.filter_id, None)
self.assertEqual(amplitude_1.scaling_time, None)
self.assertEqual(amplitude_1.magnitude_hint, 'ML')
self.assertEqual(amplitude_1.evaluation_mode, None)
self.assertEqual(amplitude_1.evaluation_status, None)
self.assertNotEqual(amplitude_1.creation_info, None)
self.assertEqual(len(amplitude_1.comments), 0)
# Test second amplitude
amplitude_2 = amplitudes[1]
self.assertEqual(
amplitude_2.resource_id, 'smi:local/amplitude/3586513')
self.assertEqual(amplitude_2.generic_amplitude, 2.9)
self.assertEqual(amplitude_2.type, None)
self.assertEqual(amplitude_2.category, None)
self.assertEqual(amplitude_2.unit, None)
self.assertEqual(amplitude_2.method_id, None)
self.assertEqual(amplitude_2.period, 0.6)
self.assertEqual(amplitude_2.snr, 4.9)
self.assertEqual(amplitude_2.time_window, None)
self.assertEqual(amplitude_2.pick_id, 'smi:local/pick/3586513')
# WaveformStreamId
waveform_2 = amplitude_2.waveform_id
self.assertEqual(waveform_2.network_code, 'XX')
self.assertEqual(waveform_2.station_code, 'GERES')
self.assertEqual(waveform_2.channel_code, None)
self.assertEqual(waveform_2.location_code, None)
self.assertEqual(waveform_2.resource_uri, None)
self.assertEqual(amplitude_2.filter_id, None)
self.assertEqual(amplitude_2.scaling_time, None)
self.assertEqual(amplitude_2.magnitude_hint, None)
self.assertEqual(amplitude_2.evaluation_mode, None)
self.assertEqual(amplitude_2.evaluation_status, None)
self.assertNotEqual(amplitude_2.creation_info, None)
self.assertEqual(len(amplitude_2.comments), 0)
def test_parameters(self):
filename = os.path.join(self.path, 'gse_2.0_standard.txt')
fields = {
'line_1': {
'time': slice(0, 21),
'time_fixf': slice(22, 23),
'lat': slice(25, 33),
'lon': slice(34, 43),
'epicenter_fixf': slice(44, 45),
'depth': slice(47, 52),
'depth_fixf': slice(53, 54),
'n_def': slice(56, 60),
'n_sta': slice(61, 65),
'gap': slice(66, 69),
'mag_type_1': slice(71, 73),
'mag_1': slice(73, 77),
'mag_n_sta_1': slice(78, 80),
'mag_type_2': slice(82, 84),
'mag_2': slice(84, 88),
'mag_n_sta_2': slice(89, 91),
'mag_type_3': slice(93, 95),
'mag_3': slice(95, 99),
'mag_n_sta_3': slice(100, 102),
'author': slice(104, 112),
'id': slice(114, 122),
},
'line_2': {
'rms': slice(5, 10),
'ot_error': slice(15, 21),
's_major': slice(25, 31),
's_minor': slice(32, 38),
'az': slice(40, 43),
'depth_err': slice(49, 54),
'min_dist': slice(56, 62),
'max_dist': slice(63, 69),
'mag_err_1': slice(74, 77),
'mag_err_2': slice(85, 88),
'mag_err_3': slice(96, 99),
'antype': slice(104, 105),
'loctype': slice(106, 107),
'evtype': slice(108, 110),
},
'arrival': {
'sta': slice(0, 5),
'dist': slice(6, 12),
'ev_az': slice(13, 18),
'picktype': slice(19, 20),
'direction': slice(20, 21),
'detchar': slice(21, 22),
'phase': slice(23, 30),
'time': slice(31, 52),
't_res': slice(53, 58),
'azim': slice(59, 64),
'az_res': slice(65, 71),
'slow': slice(72, 77),
's_res': slice(78, 83),
't_def': slice(84, 85),
'a_def': slice(85, 86),
's_def': slice(86, 87),
'snr': slice(88, 93),
'amp': slice(94, 103),
'per': slice(104, 109),
'mag_type_1': slice(110, 112),
'mag_1': slice(112, 116),
'mag_type_2': slice(117, 119),
'mag_2': slice(119, 123),
'id': slice(124, 132),
},
}
catalog = _read_gse2(filename, None, 'TZ', '00', 'SHZ', 'quakeml:idc',
fields, False, 'AGE')
self.assertEqual(len(catalog), 1)
# Test Comment ResourceIdentifier
self.assertEqual(len(catalog.comments), 1)
comment = catalog.comments[0]
self.assertEqual(comment.resource_id.id[:12], 'quakeml:idc/')
# Test Event ResourceIdentifier
event = catalog[0]
self.assertEqual(event.resource_id, 'quakeml:idc/event/280435')
self.assertEqual(event.creation_info.agency_id, 'AGE')
# Test Origin ResourceIdentifier
self.assertEqual(len(event.origins), 1)
origin = event.origins[0]
self.assertEqual(origin.resource_id, 'quakeml:idc/origin/282672')
self.assertEqual(origin.method_id, 'quakeml:idc/method/inversion')
# Test Pick ResourceIdentifier
self.assertEqual(len(event.picks), 9)
pick = event.picks[0]
self.assertEqual(pick.resource_id, 'quakeml:idc/pick/3586432')
# Test Arrival ResourceIdentifier
self.assertEqual(len(origin.arrivals), 9)
arrival = origin.arrivals[0]
self.assertEqual(
arrival.resource_id, 'quakeml:idc/origin/282672/arrival/3586432')
# Test Magnitude ResourceIdentifier
self.assertEqual(len(event.magnitudes), 2)
magnitude = event.magnitudes[0]
self.assertEqual(
magnitude.resource_id, 'quakeml:idc/origin/282672/magnitude/0')
# Test StationMagnitude ResourceIdentifier
self.assertEqual(len(event.station_magnitudes), 4)
sta_mag = event.station_magnitudes[0]
self.assertEqual(
sta_mag.resource_id, 'quakeml:idc/magnitude/station/3586432/0')
# Test Amplitude ResourceIdentifier
self.assertEqual(len(event.amplitudes), 6)
amplitude = event.amplitudes[0]
self.assertEqual(
amplitude.resource_id, 'quakeml:idc/amplitude/3586432')
# Test network code
waveform = pick.waveform_id
self.assertEqual(waveform.network_code, 'TZ')
self.assertEqual(waveform.channel_code, 'SHZ')
self.assertEqual(waveform.location_code, '00')
