def _channel_from_stats(stats):
if stats.standard.units in UNITS:
units = UNITS[stats.standard.units]
else:
units = ''
instrument = stats.standard.instrument
serialnum = stats.standard.sensor_serial_number
if len(instrument) or len(serialnum):
equipment = Equipment(type=instrument,
serial_number=serialnum)
else:
equipment = None
depth = 0.0
azimuth = None
c1 = 'horizontal_orientation' in stats.standard
c2 = c1 and not np.isnan(stats.standard.horizontal_orientation)
if c2:
azimuth = stats.standard.horizontal_orientation
else:
azimuth = 0
if not (azimuth >= 0 and azimuth <= 360):
azimuth = 0
response = None
if 'response' in stats:
response = stats['response']
else:
# we may have instrument sensitivity...
frequency = 1 / stats['standard']['instrument_period']
units = stats.standard.units
if not np.isnan(stats['standard']['instrument_sensitivity']):
sens = stats['standard']['instrument_sensitivity']
else:
sens = 1.0
sensitivity = InstrumentSensitivity(sens,
frequency=frequency,
input_units=units,
output_units='COUNTS')
response = Response(instrument_sensitivity=sensitivity)
comments = Comment(stats.standard.comments)
logging.debug('channel: %s' % stats.channel)
channel = Channel(stats.channel,
stats.location,
stats.coordinates['latitude'],
stats.coordinates['longitude'],
stats.coordinates['elevation'],
depth,
azimuth=azimuth,
sample_rate=stats.sampling_rate,
storage_format=stats.standard.source_format,
calibration_units=units,
comments=[comments],
response=response,
sensor=equipment)
return channel
def _channel_from_stats(stats):
if stats.standard.units in UNITS:
units = UNITS[stats.standard.units]
else:
units = ""
instrument = stats.standard.instrument
serialnum = stats.standard.sensor_serial_number
if len(instrument) or len(serialnum):
equipment = Equipment(type=instrument, serial_number=serialnum)
else:
equipment = None
depth = 0.0
azimuth = None
c1 = "horizontal_orientation" in stats.standard
c2 = c1 and not np.isnan(stats.standard.horizontal_orientation)
if c2:
azimuth = stats.standard.horizontal_orientation
else:
azimuth = 0
if not (azimuth >= 0 and azimuth <= 360):
azimuth = 0
response = None
if "response" in stats:
response = stats["response"]
else:
# we may have instrument sensitivity...
frequency = 1 / stats["standard"]["instrument_period"]
units = stats.standard.units
if not np.isnan(stats["standard"]["instrument_sensitivity"]):
sens = stats["standard"]["instrument_sensitivity"]
else:
sens = 1.0
sensitivity = InstrumentSensitivity(sens,
frequency=frequency,
input_units=units,
output_units="COUNTS")
response = Response(instrument_sensitivity=sensitivity)
comments = Comment(stats.standard.comments)
logging.debug(f"channel: {stats.channel}")
channel = Channel(
stats.channel,
stats.location,
stats.coordinates["latitude"],
stats.coordinates["longitude"],
stats.coordinates["elevation"],
depth,
azimuth=azimuth,
sample_rate=stats.sampling_rate,
calibration_units=units,
comments=[comments],
response=response,
sensor=equipment,
)
return channel
def test_comment_str(self):
"""
Tests the __str__ method of the Comment object.
"""
c = Comment(value='test_comment', id=9,
begin_effective_time=UTCDateTime(1240561632),
end_effective_time=UTCDateTime(1584561632), authors=[])
self.assertEqual(str(c), "Comment:\ttest_comment\n"
"\tBegin Effective Time:\t2009-04-24T08:27:12"
".000000Z\n"
"\tEnd Effective Time:\t2020-03-18T20:00:32.000000Z\n"
"\tAuthors:\t\t[]\n"
"\tId:\t\t\t9")
def _channel_from_stats(stats):
if stats.standard.units in UNITS:
units = UNITS[stats.standard.units]
else:
units = ''
instrument = stats.standard.instrument
serialnum = stats.standard.sensor_serial_number
if len(instrument) or len(serialnum):
equipment = Equipment(type=instrument, serial_number=serialnum)
else:
equipment = None
depth = 0.0
azimuth = None
c1 = 'horizontal_orientation' in stats.standard
c2 = c1 and not np.isnan(stats.standard.horizontal_orientation)
if c2:
azimuth = stats.standard.horizontal_orientation
else:
azimuth = 0
response = None
if 'response' in stats:
response = stats['response']
comments = Comment(stats.standard.comments)
logging.debug('channel: %s' % stats.channel)
channel = Channel(stats.channel,
stats.location,
stats.coordinates['latitude'],
stats.coordinates['longitude'],
stats.coordinates['elevation'],
depth,
azimuth=azimuth,
sample_rate=stats.sampling_rate,
storage_format=stats.standard.source_format,
calibration_units=units,
comments=[comments],
response=response,
sensor=equipment)
return channel