Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • 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
  • Nov 27 14:20
    megies edited #2759
  • Nov 27 14:20
    megies edited #2759
  • Nov 27 14:19
    megies edited #2759
  • Nov 27 14:19
    megies labeled #2759
  • Nov 27 14:19
    megies opened #2759
  • Nov 27 14:19
    megies milestoned #2759
  • Nov 27 14:10

    megies on read_inventory_level

    read_inventory/stationxml: add … (compare)

DJMO3
@DJMO3
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
the install is likely correct, but obspy doens't find the 'entry points' (=all format readers are defined via entry points) if the kernels are not restarted
Emmanuel D. Castillo Taborda
@ecastillot
@ThomasLecocq have you tried it with conda? Just like I told to @calum-chamberlain , this solution is ok for pip. However, when I restart the kernel after the conda installation, nothing happens.
Thomas Lecocq
@ThomasLecocq
I never used google collab. but it shouldn't be different (and I didn't understand from your answer that it actually worked for you)
and are you working in the right env ?
Emmanuel D. Castillo Taborda
@ecastillot
Yes, I'm agree with you, it shouldn't be different. For this reason I don't understand what is the mistake... I guess that I am working in the right env because
as you can see in the link that I posted, obspy is well imported.
Thomas Lecocq
@ThomasLecocq
So when you restart the kernel, do you actually reinit the conda env? Could you deactivate and activate, evne if it is the base env?