Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • 11:42
    megies closed #2363
  • 11:42
    megies commented #2363
  • 11:36
    ThomasLecocq commented #2452
  • 11:21
    ThomasLecocq synchronize #2452
  • 11:21

    ThomasLecocq on feature_io_dmx

    reabased and corrected based on… remaining stuff (compare)

  • 11:16
    megies commented #2456
  • 11:15
    megies labeled #2456
  • 11:15
    megies milestoned #2456
  • 10:10
    megies milestoned #2424
  • 10:05
    megies edited #2457
  • 10:05
    megies edited #2457
  • 10:05
    megies edited #2457
  • 10:05
    megies edited #2457
  • 10:05
    megies edited #2457
  • 10:05
    megies edited #2457
  • 09:49
    megies commented #2457
  • 09:45
    megies closed #2409
  • 09:45
    megies commented #2409
  • 09:41
    megies commented #2454
  • 09:39
    megies edited #2455
weitianli
@weitianli
I try to install by conda..but conda shows me PackagesNotFoundError: The following packages are not available from current channels:..
when I try using pip install. It comes out 'Command "python setup.py egg_info" failed with error code 1 in /private/tmp/pip-build-8UCf_3/obspy/'
Thank you for your advice.
Tobias Megies
@megies

I rebooted my obspy

0_o

Thomas Lecocq
@ThomasLecocq
woow, waveloc :) good package, please go and make it obspy 1.1.1 compatible !
Tobias Megies
@megies
to behonest, switchin from 1.0 to 1.1 should just be a matter of renaming some functions and methods, nothin crazy should need to be done
weitianli
@weitianli
yep waveloc is good package. But it seems like nobody maintain
Thomas Lecocq
@ThomasLecocq
Yep, it's hard to keep maintaining packages, but you're lucky at least Nadege and Alessia published it and made it available, so it's up to you, other users, the community, to perpetuate it!
titi-dev
@titi-dev
Hii guys good afternoon.Today I come with a very important query. Do anybody work on receiver functions in seismology.Anyway generally for teleseismic receiver function computation people usually need the earthquake data having body wave magnitude > 5.5 with epicentral distance range 30-90 degrees. My query is, if i have the raw field data how can i extract the earthquake data of body wave magnitude > 5.5 and cut it into 30-90 degree epicentral range so that i can use it for imaging the subsurface. I think it can be done by NDK search and global earthquake catalog. is there any facility in obspy to do it.please suggest.
Thomas Lecocq
@ThomasLecocq
:yellow_heart:
Mattina Contreras
@mattinacm_twitter
@megies Hello. I just started to work with some REFTEK data, and i found myself in need to convert it to mseed format.
I found on gitter.im https://gitter.im/obspy/obspy/archives/2015/02/17 a feedback on this subject,
If it isn´t too much trouble, could you share withme : rt_mseed to convert to miniseed for archiving, and the python script around rt_mseed (to make daily mseed files per channel from a top-level folder containing RT130 recorded data). Everything else you could share on this subject ... i would be very thankful
Tobias Megies
@megies
obspy 1.1.x has reftek read support @mattinacm_twitter, so its trivial to read and write to a different format
titi-dev
@titi-dev
@megies but how to write the reftek format to mseed .....error usually comes like PermissionError: [Errno 13] Permission denied: 'filename.mseed'
jkmacc-LANL
@jkmacc-LANL
Permission errors don’t sound like an ObsPy problem. You may need to consider your file permissions.
Muhammad Wahyu Putra Indi
@putraindi
how to processing signal with feature extraction polarization in osbpy?
79seismo
@79seismo
Is there a way to compute rms of a stream? couldn't find any documentation for this. Thanks!
Calum Chamberlain
@calum-chamberlain

AFAIK there isn't a method - just loop over the traces and compute using numpy methods?

trace_rms = {}
for tr in st:
    rms = np.sqrt(np.mean(tr.data ** 2))
    trace_rms.update({tr.id: rms})

might do what you want (should generate a dictionary of individual trace RMS values). Not tested though.

79seismo
@79seismo
cool! thanks @calum-chamberlain
79seismo
@79seismo
How do you extract the time series after applying new_tr = trace.remove_response(....)? I tried new_tr.data but that attribute is not defined. Thanks!
Calum Chamberlain
@calum-chamberlain

Seems strange, when I try, new_tr.data is fine:

from obspy import read, read_inventory
st = read()
tr = st[0].copy()
inv = read_inventory()
new_tr = tr.remove_response(inventory=inv)
print(new_tr.data)

works happily and prints some data.

79seismo
@79seismo
I get "AttributeError: 'Stream' object has no attribute 'data'" but this new_tr.plot() shows that the data is there .
Calum Chamberlain
@calum-chamberlain
That suggests that what you have is a stream, not a trace.
79seismo
@79seismo
new_tr[0].data works! Thanks!
Calum Chamberlain
@calum-chamberlain
No worries, if you find yourself getting errors like that in the future, the same thing applies. Generally the error tries to be helpful.
79seismo
@79seismo
Thanks, still learning the ABCs of Obspy
Calum Chamberlain
@calum-chamberlain
Yup, Obspy is massive! I remain impressed. Also, amusingly, ABCs are a thing in Python
79seismo
@79seismo
Thanks!
Muhammad Wahyu Putra Indi
@putraindi
Thanks!
Muhammad Wahyu Putra Indi
@putraindi
I have 300 seismic data, how to read data seismic from the folder?
Tobias Megies
@megies
read('*')?
79seismo
@79seismo
@putraindi Lets say you have your data in a folder named data/. Then you can do the following to read one file at a time;
import obspy
from obspy import read
from obspy.io.xseed import Parser
import os, shutil, glob

for msd_file in glob.glob('data/*'):      # looping through files
       st = read(msd_file)
Tobias Megies
@megies
Any segy users care to have a look at obspy/obspy#2416 There's also a mail on users list
Fred Massin
@FMassin
Hi, is there a way to compute PPSD with signal which have already been response-corrected? Something like:
ppsd = obspy.signal.PPSD(tr.stats, None)
ppsd.add(tr,verbose=True)
Fred Massin
@FMassin
I just checked the code, and now I think there isn't such functionality... :(
Tobias Megies
@megies
you might be able to misuse special_handling="ringlaser" for this
Fred Massin
@FMassin
Yes but there is still a sensitivity correct?
Tobias Megies
@megies
or set up a flat response and feed that in like #2408
well, check the documentation..? but answer is yes
Fred Massin
@FMassin
The flat response might do the job !
Tobias Megies
@megies
in any case, it is definitly not ideal to feed in corrected data, see #931
Fred Massin
@FMassin
I know :) ! But some networks around the world don't have standard practice ...
Tobias Megies
@megies
yeah, I know can't be helped sometimes, I also get some pre-corrected industry data sometimes :(
Fred Massin
@FMassin
anyway thanks a lot for the response !
Tobias Megies
@megies
no problem
Fred Massin
@FMassin
The flat resp works but it got nasty, any simplification would be welcomed :
meta = nrl.get_response(datalogger_keys=['Generic','Unity'],
                            sensor_keys=['Generic','Unity Velocity Sensor'])
f = io.StringIO(meta.get_sacpz())
obspy.io.sac.sacpz.attach_paz(tr,f)
meta=dict(tr.stats.paz)
ppsd = PPSD(tr.stats, meta)
ppsd.add(tr,verbose=False)
Tobias Megies
@megies
How about let's try to push 1.2.0? All hands on deck https://github.com/obspy/obspy/projects/2#card-26310701
Calum Chamberlain
@calum-chamberlain
Looks like #2277 got bumped to 2.0.0 - should that card be moved out from the 1.2.0 project?
Thomas Lecocq
@ThomasLecocq

How about let's try to push 1.2.0? All hands on deck https://github.com/obspy/obspy/projects/2#card-26310701

owowowo yéééééééé