Note: The section between the two 3-line segments has been added specifically for Obspy-seismo. The rest of the instructions pertain to the main Obspy page.
Custom functions for working with big (more than 2GB) SEG Y files.
The list of functions are:
- seismo_read_segy_
- seismo_segy_read_textual_header
- seismo_segy_read_trace_headers
- seismo_segy_get_number_of_traces
- seismo_segy_remove_extended_headers
- seismo_segy_read_extended_textual_headers
- seismo_segy_read_num_samples_from_textual_header
- seismo_segy_read_binary_header
from obspy.io.segy.core import seismo_read_segy_
file_path = "home/user/segyfile.sgy"
streams = seismo_read_segy_(file_path, trace_duration=60000)
print(streams) # will return an array of streams for different components
streams[0].plot()Obspy-Seismo Installation
before Installing, you should have Git already installed on your machine.
If ObsPy is not yet installed on your machine:
pip install git+https://github.com/alemkhodadadi/obspy_seismo.git
If ObsPy is already installed:
pip install --upgrade git+https://github.com/alemkhodadadi/obspy_seismo.git
ObsPy is an open-source project dedicated to provide a Python framework for processing seismological data. It provides parsers for common file formats, clients to access data centers and seismological signal processing routines which allow the manipulation of seismological time series (see Beyreuther et al. 2010, Megies et al. 2011, Krischer et al. 2015).
The goal of the ObsPy project is to facilitate rapid application development for seismology.
ObsPy is licensed under the GNU Lesser General Public License (LGPL) v3.0.
A one-hour introduction to ObsPy is available at YouTube.
Installation instructions can be found in the wiki.
Read about how to get started in the wiki and in our Tutorial section in the documentation.
ObsPy Tutorial notebooks -- and much more on specific seismology topics -- can also be found on Seismo-Live, both as a static preview and as interactively runnable version.
from obspy import read
st = read() # load example seismogram
st.filter(type='highpass', freq=3.0)
st = st.select(component='Z')
st.plot()
The detailed changelog is available here, our docs can be found at docs.obspy.org.
Please see details on how to contribute to the project here.
- Moritz Beyreuther, Robert Barsch, Lion Krischer, Tobias Megies, Yannik Behr and Joachim Wassermann (2010), ObsPy: A Python Toolbox for Seismology, SRL, 81(3), 530-533, doi:
10.1785/gssrl.81.3.530. - Tobias Megies, Moritz Beyreuther, Robert Barsch, Lion Krischer, Joachim Wassermann (2011), ObsPy – What can it do for data centers and observatories? Annals Of Geophysics, 54(1), 47-58, doi:
10.4401/ag-4838. - Lion Krischer, Tobias Megies, Robert Barsch, Moritz Beyreuther, Thomas Lecocq, Corentin Caudron and Joachim Wassermann (2015), ObsPy: a bridge for seismology into the scientific Python ecosystem, Computational Science & Discovery, 8(1), 014003, doi:
10.1088/1749-4699/8/1/014003 - concept DOI, that always resolves to latest ObsPy versions:
10.5281/zenodo.1040769(see Zenodo FAQ) - specific version DOIs for all released ObsPy versions
