Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • 15:02
    qitiandashengsunwukong starred obspy/obspy
  • 10:50
    zacherle synchronize #2736
  • Nov 28 18:50
  • Nov 28 10:38
    ThomasLecocq commented #2758
  • Nov 28 06:41
    kamacpherson commented #2740
  • Nov 28 03:42
    msaezuchile opened #2760
  • Nov 28 01:00
    d-chambers commented #2743
  • Nov 27 16:23
    megies commented #2740
  • Nov 27 16:15
    heavelock commented #2722
  • Nov 27 16:14
    heavelock synchronize #2722
  • Nov 27 15:56
    megies edited #2759
  • Nov 27 15:56
    megies edited #2759
  • Nov 27 15:55
    megies synchronize #2759
  • Nov 27 15:55

    megies on read_inventory_level

    changelog add test for "level" in reading… (compare)

  • Nov 27 15:52
    megies commented #2757
  • Nov 27 15:45
    megies labeled #2757
  • Nov 27 14:39
    megies labeled #2756
  • Nov 27 14:39
    megies commented #2756
  • Nov 27 14:20
    megies edited #2759
  • Nov 27 14:20
    megies edited #2759
uky-jps
@uky-jps
Hi all. I was wondering if anyone knew of a way to make a plot of 3 ppsd subplots, one for each component of a seismometer in a single "summary" type plot for the station. I know you can mute the output of the ppsd.plot(show=False) to then manipulate the active ppsd.plot. Any guidance on this would be extremely helpful!
3 replies
Glenn Nelson
@dr-glenn
I have very large mseed files (12 GB) and am using buffered read as suggested in issue #1419. That works well, but I was surprised to find that if I decreased chunksize from 100K to 10K, I actually had MORE traces returned in each buffer. I expected fewer. I don't know how miniSeed files are organized. I suspect a guru can easily explain this strange observation to me.
DJMO3
@DJMO3
I'm very very new to obspy and can't seem to get the data I need correctly. Anybody willing to help out?
Tobias Megies
@megies
@dr-glenn hmm.. in general this should not influence what you get out of the read operations, although miniseed reading is merging together chunks with sub-sampling rate gaps and if you read smaller parts at once you might not get them merged together automatically if your chunk is right at such a subsample gap
@DJMO3 uhm.. need more details?
DJMO3
@DJMO3
I have to capture the waveform data from all stations in a specific network for signal processing. I figured that I will need the mass_downloader module to do so. I have not understood how the CircularDomain or other Domain method work
Does it return all the stations in the domain or does it return the events that occurred during the specified time?
Thomas Lecocq
@ThomasLecocq
generic question related to "best practice": how far away from each other could you have two sensors sharing net+sta , but not loc ids ? I imagine on a field with up to 100m apart, logged on a single unit? or borehole strings with different sensors at different depths ?
Lula
@LulaLopezM_twitter
Hello! I need help, I try to run isola but I have got this problem:
runfile('/Users/lucianalopez/Sismologia/ISOLA/1_Efpalio_SAC_orig/src/class_isola.py', wdir='/Users/lucianalopez/Sismologia/ISOLA/1_Efpalio_SAC_orig/src')
File "/Users/lucianalopez/Sismologia/ISOLA/1_Efpalio_SAC_orig/src/class_isola.py", line 64
return dict(map(lambda (key, value): (prefix+str(key)+suffix, value), somedict.items()))
^
SyntaxError: invalid syntax
Elliott Sales de Andrade
@QuLogic
what does this have to do with ObsPy?
Your code is invalid Python
uky-jps
@uky-jps
Obspy's PPSD function has the ability or calculates the Mean and Mode as an array that can be plotted, but does not appear to calculate the Median. Is there a way to extract the Median from PPSD?
-> 50
uky-jps
@uky-jps
@ThomasLecocq Thanks!
uky-jps
@uky-jps
I am trying to pull some data from IRIS using obspy's Client.get_waveform function
'''
client = Client('IRIS')
st = client.get_waveforms("IU", "CCM", "00,10", "HH,BHZ", t1, t2)
'''
I pulled both BH
and HH to see if I could get CCM station to pull data, and I can pull for BH components, but for some reason, I cannot get data from the HH channels. I'm able to get HH for other stations but I cannot get these sepcifically. As far as what I can see on IRIS Station Query page for this station, HH* channels should be available.
https://ds.iris.edu/mda/IU/CCM/
I will also note that there are several other stations that I am having this same issue with, specifically stations in the IU network (IU.WCI, IU.WVT). Is this maybe something to do with this network?
Calum Chamberlain
@calum-chamberlain

@uky-jps you are missing the required wildcard character. - HH is not a valid channel code, but HH? will find the valid codes (HH1, HH2 and HHZ in this case).
I tried:

from obspy import UTCDateTime
from obspy.clients.fdsn import Client

client = Client("IRIS")
t1, t2 = UTCDateTime(2019, 1, 1), UTCDateTime(2019, 1, 1) + 600
st = client.get_waveforms("IU", "CCM", "00,10", "HH?,BHZ", t1, t2)

and I get 5 traces including HH1, HH2 and HHZ.

uky-jps
@uky-jps
My apologies, you are correct, my wildcard should have been there. I must have accidentally deleted the "*" symbol I had there originally, although I know "?" does work for single character wildcards. I tried it again but this time for an earlier time as you did (2019) and it does appear to pull the station I was searching for. I was trying to pull data from earlier this month (Oct 2020) to do the comparison to match my dates, but I guess its related to IRIS data availablility... Dang. Thanks for the help!
vivigb
@vivigb
Hello,
I rotated the IRIS data component to R and T components , for the surface tomography, using both Obspy and SAC , but the output seems to be different . Is there anyway to find which one is wrong? or anyone I can consult here?
Emmanuel D. Castillo Taborda
@ecastillot

Good moorning!
I downloaded one year of data with massdownloader. However, at the first time of the downloading, I writed wrong the name of the files.Due to I already had half downloaded. Therefore, I developed a script to change the name of these files automatically.

After that, I wanted to rerunning the download with massdownloader, so, I corrected the filenames in mseed_storage parameter of massdownloader. However, when I ran the downloading with massdownloader again, then it repeats the file with the same name. I hoped that it didn't repeat the downloads

Do you know what happen here?

Thomas Lecocq
@ThomasLecocq

Hello,
I rotated the IRIS data component to R and T components , for the surface tomography, using both Obspy and SAC , but the output seems to be different . Is there anyway to find which one is wrong? or anyone I can consult here?

without a screenshot, an example etc... difficult to say. did you apply the same filters (if any), etc ?

@ecastillot nope, no idea, difficult from your explanation to understand what you expect, but I guess it's linked to the way the MassDownloader indexes the "already downloaded data", @krischer ?
vivigb
@vivigb
@ThomasLecocq I didn't apply the filter at all. I resampled all the files into 1Hz. Is there any way to share screenshots or attach images here?
Thomas Lecocq
@ThomasLecocq
yep, just copy paste
vivigb
@vivigb
New Microsoft Word Document.pdf
these are the result I got from both obspy and SAC.
Thomas Lecocq
@ThomasLecocq
SAC doesn't rename the components ?
I don't know SAC but did you follow all the commadns to set the orientations of the E and N sensors ? http://geophysics.eas.gatech.edu/classes/SAC/ --> Rotate a pair of data.
?
79seismo
@79seismo
@vivigb I have a feeling that you're doing something wrong with SAC. As @ThomasLecocq mentioned, your channel codes haven't changed after the rotation. SAC does change these codes. Are you using "rotate to GCP"?
vivigb
@vivigb
@ThomasLecocq @79seismo I think SAC won't rename the components when using the rotate command.
Please see this link
@79seismo yeah i am using " rotate to GCP"
brioglade
@brioglade_gitlab
from obspy.imaging.radpattern import plot_3drpattern
mt = [0.91, -0.89, -0.02, 1.78, -1.55, 0.47]
plot_3drpattern(mt,kind='p_sphere')

ModuleNotFoundError Traceback (most recent call last)

<ipython-input-2-4718903f829c> in <module>
1 import obspy
----> 2 from obspy.imaging.radpattern import plot_3drpattern
3 mt = [0.91, -0.89, -0.02, 1.78, -1.55, 0.47]
4 plot_3drpattern(mt,kind='p_sphere')

ModuleNotFoundError: No module named 'obspy.imaging.radpattern'

hi,how can I fix it?
Elliott Sales de Andrade
@QuLogic
What is radpattern? There's no such thing in ObsPy
brioglade
@brioglade_gitlab
Elliott Sales de Andrade
@QuLogic
That is very old, and never released like that, which you can see in that PR
brioglade
@brioglade_gitlab
@QuLogic So what's the right way to write to do seismic source ridiation pattern?I am a new hand.Thank you very much.
What's the CMTSOLUTION file format to write?Do you have an example?Thank you.
Calum Chamberlain
@calum-chamberlain
This first thing that comes up when I google "obspy radiation plot" is this which is what appears to have been the merged change in PR #1192. It looks like you might want to check that page of the obspy docs out and see if that is what you want.
Calum Chamberlain
@calum-chamberlain

For your second question, CMTSOLUTION is a standard format that looks like this.

However, the Event.plot() with defaults to kind=('beachball', 'p_sphere', 's_sphere') calls the obspy.imaging.source.plot_radiation_pattern under the hood, which if you have the moment-tensor and no Event might be simpler to call.

Calum Chamberlain
@calum-chamberlain
I didn't know about this though - they look like really cool plots!
brioglade
@brioglade_gitlab
@calum-chamberlain Thank you very much!
79seismo
@79seismo
@vivigb yes, you are right, KCMPNM has to be changed manually. However, you should read in only the horizontal components for rotate to GCP. The fact that your Z is also changing after the rotation suggests that you're doing something additional to rotating? Rotate to GCP is only a horizontal operation.
Emmanuel D. Castillo Taborda
@ecastillot

Good moorning!
I tried to install obspy with miniconda in google collab. Apparently, it was installed ok. However, when I wanted to use ' obspy.read() ' appears this error: 'TypeError: Unknown format for file /usr/local/lib/python3.6/site-packages/obspy/core/data/BW_RJOB.xml'

can you help me please? Maybe obspy can't be installed with conda in google collab? When I installed with pip, it does his job ok. But for certain reasons I need to install it with conda.

In the next link, you can see the collab notebook that I'm using:
https://colab.research.google.com/drive/1eXNROuqsYeWjHqEyCbIDbZFlgSlPKNMS?usp=sharing

Calum Chamberlain
@calum-chamberlain
Kia Ora, this is the same as obspy/obspy#2547 - the solution then was to use pip and restart the run-time. I don't know if simply restarting the run-time is enough with conda installs, but it is worth a try.
Emmanuel D. Castillo Taborda
@ecastillot
@calum-chamberlain Regards,
I usually use this solution for pip. However, with conda installs are not the same. The error persists.
Thomas Lecocq
@ThomasLecocq
you have to restart the runtime/kernel