obspy.io.sh - Q and ASC read and write, EVT read support (Seismic Handler)

The obspy.io.sh package contains methods in order to read and write seismogram files in the Q and ASC format used by the Seismic Handler software package (https://www.seismic-handler.org). EVT event files can also be read.

copyright:

The ObsPy Development Team (devs@obspy.org)

license:

GNU Lesser General Public License, Version 3 (https://www.gnu.org/copyleft/lesser.html)

Reading

Importing Q or ASC files is done similar to reading any other waveform data format within ObsPy by using the read() method of the obspy.core module. Examples seismograms files may be found at https://examples.obspy.org.

>>> from obspy import read
>>> st = read("/path/to/QFILE-TEST-ASC.ASC")
>>> st 
<obspy.core.stream.Stream object at 0x...>
>>> print(st)  
3 Trace(s) in Stream:
.TEST..BHN | 2009-10-01T12:46:01.000000Z - ... | 20.0 Hz, 801 samples
.TEST..BHE | 2009-10-01T12:46:01.000000Z - ... | 20.0 Hz, 801 samples
.WET..HHZ  | 2010-01-01T01:01:05.999000Z - ... | 100.0 Hz, 4001 samples

The file format will be determined automatically. Each trace (multiple channels are mapped to multiple traces) will have a stats attribute containing the usual information.

>>> print(st[0].stats) 
             network:
             station: TEST
            location:
             channel: BHN
           starttime: 2009-10-01T12:46:01.000000Z
             endtime: 2009-10-01T12:46:41.000000Z
       sampling_rate: 20.0
               delta: 0.05
                npts: 801
               calib: 1.5
             _format: SH_ASC
                  sh: AttribDict({...})
>>> print(st[0].stats.sh['COMMENT'])
TEST TRACE IN QFILE #1

The actual data is stored as numpy.ndarray in the data attribute of each trace.

>>> st[0].data 
array([ 0.00000000e+00,   0.00000000e+00,   0.00000000e+00, ...
       -4.03387604e+01,  -3.99515305e+01,  -3.95423012e+01], dtype=float32)

Writing

Writing is also done in the usual way:

>>> st.write('file.q', format = 'Q') 

or

>>> st.write('file.asc', format = 'SH_ASC') 

Modules

core

SH bindings to ObsPy core module.

evt

SeismicHandler evt file bindings to ObsPy core module.