From 2212bijaya at gmail.com Mon Dec 2 16:58:51 2019 From: 2212bijaya at gmail.com (bijaya dalei) Date: Mon, 2 Dec 2019 21:28:51 +0530 Subject: [rf] needs to cut teleseismic data in 30-90 degree range with events info Message-ID: Good evening sir.I am Bijay from India,currently works as ph.d research scholar.I have a very small query regarding sac data manipulations(cutting).However, my raw data sets are of sac format that consists of day file(one raw sac file/day/station) and i donot know more about sac files in which portion of the waveform the earthquake is recorded. i just want to apply receiver functions(rf package python) on my raw data sets.For this, data needs to be in teleseismic range(30-90 degree) with event magnitude > 4/5 and cutting of the waveforms some seconds before p onset and some seconds after p-onset and so on....can you please suggest how to make my raw sac file data to that mentioned form(by using something like ndksearch i.e. i want to make raw data like what is given in example data in rf package).Thank you. -------------- next part -------------- An HTML attachment was scrubbed... URL: From tom.eulenfeld at uni-jena.de Tue Dec 3 16:26:07 2019 From: tom.eulenfeld at uni-jena.de (Tom Eulenfeld) Date: Tue, 3 Dec 2019 16:26:07 +0100 Subject: [rf] needs to cut teleseismic data in 30-90 degree range with events info In-Reply-To: References: Message-ID: <2b7f109b-eb2b-7fb0-1fbe-f52d6898b1f6@uni-jena.de> Hi Bijay, welcome to the mailing list. Firstly, you need to download an earthquake catalog with events in the specific distance range. Either via a web page or with a FDSN request or with an obspy script (see obspy tutorial). After that you need to calculate the travel time to your stations, load the corresponding day file(s) and trim it accordingly. For this task you could use your own script or use the rf.util.iter_event_data iterator (see rf documentation). It can be used like from rf.util import iter_event_data events = read_events(...) inventory = read_inventory(...) def get_waveforms(...): .... for stream in iter_event_data(events, inventory, get_waveforms, ...): stream.write('raw...') stream.rf(...) stream.write('rf...') Best regards, Tom On 02.12.19 16:58, bijaya dalei wrote: > Good evening sir.I am Bijay from India,currently works as ph.d research > scholar.I have a very small query regarding sac data > manipulations(cutting).However, my raw data sets are of sac format that > consists of day file(one raw sac file/day/station) and i donot know more > about sac files in which portion of the waveform the earthquake is > recorded. i just want to apply receiver functions(rf package python) on > my raw data sets.For this, data needs to be in teleseismic range(30-90 > degree) with event magnitude > 4/5 and cutting of the waveforms some > seconds before p onset and some seconds after p-onset and so on....can > you please suggest how to make my raw sac file data to that mentioned > form(by using something like ndksearch i.e. i want to make raw data like > what is given in example data in rf package).Thank you. >