@ThomasLecocq I wrote a small code after following the manual of obspy for slicing.But my code only slice first 20 seconds of the given trace, but i need to slice whole 1 hr length continuous trace by shifting 20 seconds....can you please suggest whats wrong with my code ....my code is here `
from obspy import read
for wav in list:
st = read(wav)
tr=st.copy()
tr1=tr[0]
start = tr1.stats.starttime
stop = tr1.stats.starttime+20
tr2 = tr1.slice(start,stop)
print(tr2)
tr2.write('sample_data.SAC',format='SAC')`
@Kay089 it looks like this error is raised by this line which suggests that you should reduce the width of your dayplot to get around this.
1 sample-per minite is quite a low sampling-rate, so I doubt this plotting routine is designed with this in mind, hence the slightly opaque error message. Try passing the argument st.plot(..., size=(200, 600)) - in the size argument, width is the first part of the tuple. The default is 600, so reduce this until you get a plot!
3 replies
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
'''
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?
@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.
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?
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 ?
MassDownloader indexes the "already downloaded data", @krischer ?
Please see this link
@79seismo yeah i am using " rotate to GCP"