Source code for obspy.sac.core

# -*- coding: utf-8 -*-
"""
SAC bindings to ObsPy core module.

:copyright:
    The ObsPy Development Team (devs@obspy.org) & C. J. Ammon
:license:
    GNU Lesser General Public License, Version 3
    (http://www.gnu.org/copyleft/lesser.html)
"""

from obspy import Trace, Stream
from obspy.sac.sacio import SacIO, _isText
import os
import struct


[docs]def isSAC(filename): """ Checks whether a file is a SAC file or not. :type filename: string :param filename: SAC file to be checked. :rtype: bool :return: ``True`` if a SAC file. .. rubric:: Example >>> isSAC('/path/to/test.sac') #doctest: +SKIP """ try: f = open(filename, 'rb') # 70 header floats, 9 position in header integers f.seek(4 * 70 + 4 * 9) data = f.read(4) f.close() npts = struct.unpack('<i', data)[0] except: return False # check file size st = os.stat(filename) sizecheck = st.st_size - (632 + 4 * npts) if sizecheck != 0: # check if file is big-endian npts = struct.unpack('>i', data)[0] sizecheck = st.st_size - (632 + 4 * npts) if sizecheck != 0: # File-size and theoretical size inconsistent! return False return True
[docs]def isSACXY(filename): """ Checks whether a file is alphanumeric SAC file or not. :type filename: string :param filename: Alphanumeric SAC file to be checked. :rtype: bool :return: ``True`` if a alphanumeric SAC file. .. rubric:: Example >>> isSACXY('/path/to/testxy.sac') #doctest: +SKIP """ ### First find out if it is a text or a binary file. This should ### always be true if a file is a text-file and only true for a ### binary file in rare occasions (Recipe 173220 found on ### http://code.activestate.com/ if not _isText(filename, blocksize=512): return False try: f = open(filename) hdcards = [] # read in the header cards for _i in xrange(30): hdcards.append(f.readline()) npts = int(hdcards[15].split()[-1]) # read in the seismogram seis = f.read(-1).split() except: return False # check that npts header value and seismogram length are consistent if npts != len(seis): return False return True
[docs]def readSACXY(filename, headonly=False, debug_headers=False, **kwargs): # @UnusedVariable """ Reads an alphanumeric SAC file and returns an ObsPy Stream object. .. warning:: This function should NOT be called directly, it registers via the ObsPy :func:`~obspy.core.stream.read` function, call this instead. :type filename: str :param filename: Alphanumeric SAC file to be read. :type headonly: bool, optional :param headonly: If set to True, read only the head. This is most useful for scanning available data in huge (temporary) data sets. :type debug_headers: bool, optional :param debug_headers: Extracts also the SAC headers ``'nzyear', 'nzjday', 'nzhour', 'nzmin', 'nzsec', 'nzmsec', 'delta', 'scale', 'npts', 'knetwk', 'kstnm', 'kcmpnm'`` which are usually directly mapped to the :class:`~obspy.core.stream.Stream` object if set to ``True``. Those values are not synchronized with the Stream object itself and won't be used during writing of a SAC file! Defaults to ``False``. :rtype: :class:`~obspy.core.stream.Stream` :return: A ObsPy Stream object. .. rubric:: Example >>> from obspy import read # doctest: +SKIP >>> st = read("/path/to/testxy.sac") # doctest: +SKIP """ t = SacIO(debug_headers=debug_headers) if headonly: t.ReadSacXYHeader(filename) else: t.ReadSacXY(filename) # assign all header entries to a new dictionary compatible with ObsPy header = t.get_obspy_header() if headonly: tr = Trace(header=header) else: tr = Trace(header=header, data=t.seis) return Stream([tr])
[docs]def writeSACXY(stream, filename, **kwargs): # @UnusedVariable """ Writes a alphanumeric SAC file. .. warning:: This function should NOT be called directly, it registers via the the :meth:`~obspy.core.stream.Stream.write` method of an ObsPy :class:`~obspy.core.stream.Stream` object, call this instead. :type stream: :class:`~obspy.core.stream.Stream` :param stream: The ObsPy Stream object to write. :type filename: str :param filename: Name of file to write. .. rubric:: Example >>> from obspy import read >>> st = read() >>> st.write("testxy.sac", format="SACXY") #doctest: +SKIP """ # Translate the common (renamed) entries base, ext = os.path.splitext(filename) for i, trace in enumerate(stream): t = SacIO(trace) if len(stream) != 1: filename = "%s%02d%s" % (base, i + 1, ext) t.WriteSacXY(filename) return
[docs]def readSAC(filename, headonly=False, debug_headers=False, **kwargs): # @UnusedVariable """ Reads an SAC file and returns an ObsPy Stream object. .. warning:: This function should NOT be called directly, it registers via the ObsPy :func:`~obspy.core.stream.read` function, call this instead. :type filename: str :param filename: SAC file to be read. :type headonly: bool, optional :param headonly: If set to True, read only the head. This is most useful for scanning available data in huge (temporary) data sets. :type debug_headers: bool, optional :param debug_headers: Extracts also the SAC headers ``'nzyear', 'nzjday', 'nzhour', 'nzmin', 'nzsec', 'nzmsec', 'delta', 'scale', 'npts', 'knetwk', 'kstnm', 'kcmpnm'`` which are usually directly mapped to the :class:`~obspy.core.stream.Stream` object if set to ``True``. Those values are not synchronized with the Stream object itself and won't be used during writing of a SAC file! Defaults to ``False``. :rtype: :class:`~obspy.core.stream.Stream` :return: A ObsPy Stream object. .. rubric:: Example >>> from obspy import read # doctest: +SKIP >>> st = read("/path/to/test.sac") # doctest: +SKIP """ # read SAC file t = SacIO(debug_headers=debug_headers) if headonly: t.ReadSacHeader(filename) else: t.ReadSacFile(filename) # assign all header entries to a new dictionary compatible with an ObsPy header = t.get_obspy_header() if headonly: tr = Trace(header=header) else: tr = Trace(header=header, data=t.seis) return Stream([tr])
[docs]def writeSAC(stream, filename, **kwargs): # @UnusedVariable """ Writes a SAC file. .. warning:: This function should NOT be called directly, it registers via the the :meth:`~obspy.core.stream.Stream.write` method of an ObsPy :class:`~obspy.core.stream.Stream` object, call this instead. :type stream: :class:`~obspy.core.stream.Stream` :param stream: The ObsPy Stream object to write. :type filename: str :param filename: Name of file to write. .. rubric:: Example >>> from obspy import read >>> st = read() >>> st.write("test.sac", format="SAC") #doctest: +SKIP """ # Translate the common (renamed) entries base, ext = os.path.splitext(filename) for i, trace in enumerate(stream): t = SacIO(trace) if len(stream) != 1: filename = "%s%02d%s" % (base, i + 1, ext) t.WriteSacBinary(filename) return