Questions tagged [obspy]

ObsPy is a seismology toolbox for Python

ObsPy is a Python library that provides routines useful in seismology. Routines include parsers for common file formats and signal processing for data analysis and manipulation.

42 questions
1
vote
1 answer

Need more precision to save to SEG-Y with ObsPy

I am trying to save data in SEG-Y format using Python's ObsPy. I am having trouble with the data types. I need higher precision than float32, because I have a very high sampling rate (10 MSamples/sec), so a time interval of only 0.1 microseconds. …
user3403779
  • 157
  • 2
  • 10
1
vote
1 answer

Obspy.read() not working when called from mod_python apache module

I'm trying to make a simple web service which reads an MSEED file and outputs some information on that. The most reasonable choice seemed to be using python with the excellent obspy module (A Python Toolbox for seismology/seismological…
Daniele B
  • 3,117
  • 2
  • 23
  • 46
0
votes
0 answers

Obspy installation error: 'ERROR: Dependency "OpenBLAS" not found'

I am trying to install obspy on my raspberry pi but I keep getting this error: user@user:~ $ pip3 install obspy Defaulting to user installation because normal site-packages is not writeable Looking in indexes: https://pypi.org/simple,…
RNG
  • 101
  • 5
0
votes
0 answers

How to get name of all stations in a network in the form of list using obspy?

I know that get_station method of a client returns an inventory object, but I can't find any attributes or method related to the name of a particular station. Could anyone help me with it? Thank you. I need to handle data from a series of station in…
0
votes
0 answers

The way to stack files in obspy

I got 10 data files which have 24 traces (MASW results). I need to stack 10 data files and made the code but the append part in the code to stack. Could you please edit the append part to stack in the code? Or if you know the way to stack output…
Joe
  • 1
  • 1
0
votes
0 answers

Key character reading segy format with variable names

I have several earthquakes in segy format. I am trying to create individual Streams(ObsPy) for each event, but this is not being possible as the name of each event varies according to the date. from obspy.io.segy.core import _read_segy from…
Pahan
  • 5
  • 2
0
votes
0 answers

Geographiclib installed but not working properly

(re-posted question!) I am trying to use the Obspy get_ray_paths_geo tool. I already installed geographiclib Requirement already satisfied: geographiclib in c:\users\pahan\anaconda3\lib\site-packages (2.0) , but when plotting what I need I get this…
Pahan
  • 5
  • 2
0
votes
1 answer

how can i plot a specrtogram in python using Obspy?

First, I display the waveform by SAC file and I want to plot a spectrogram. st = read('SAC file HERE', header=0, index_col=0, parse_dates=True, squeeze=True) pyplot.show() st.plot() That will show the waveform, so how can I put it in a spectrogram?
0
votes
1 answer

Copy traces in a loop OBSPY

I have several traces which I imported from a folder, and I need to copy them for further processing. I was trying to copy them into a loop with the original procedure for copying one trace for i in range(550): tr[i] = st[i].copy() but it…
Pahan
  • 5
  • 2
0
votes
0 answers

Obspy Issue: Acceleration Data from GeoNet

I am trying to create ASCII files of acceleration data from various New Zealand earthquakes, code attached below. I was previously able to generate ASCII files of the acceleration data, however, now running the same code I get the…
b_ang
  • 1
0
votes
0 answers

How could I interpolate the data in obspy?

I was using obspy to look some Deep-ocean Assessment and Reporting of Tsunamis (DART) data. But the data is uncontinuous, I find the gap in the time series such as this. enter image description here And then, how can i fill the gap ? The raw data…
Frank
  • 1
  • 1
0
votes
0 answers

how to create .npy file from SAC or mseed files

how to create .npy files for both single trace and multiple traces of SAC or mseed files of seismographs.
Pragnath
  • 1
  • 1
0
votes
0 answers

Could I get the original data in the Obspy spectrograms?

I was useing Obspy spectrogram to look some Deep-ocean Assessment and Reporting of Tsunamis (DART) data. The raw data is seafloor pressure[dbar].(The variation of pressure with Time) The spectrogram x scale is time series, y scale is frequency. I…
Frank
  • 1
  • 1
0
votes
1 answer

Adding more than one station to xml using obspy

I'm a new python user. I'm trying to use obspy to create xml for a seismic array. I downloaded the template found at https://docs.obspy.org/tutorial/code_snippets/stationxml_file_from_scratch.html. import obspy from obspy.core.inventory import…
Ric
  • 3
  • 2
0
votes
1 answer

obspy plot streams as one plot with different color

Hi I am new to use obspy. I want to plot two streams to one plot. I made code as below. st1=read('/path/1.SAC') st1+=read('/path/2.SAC') st1.plot() I succeed to plot two plots but what I want to do is plotting them as two colors. When I put the…
Joe
  • 51
  • 7