[rf] needs to cut teleseismic data in 30-90 degree range with events info

[Tom Eulenfeld]

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.
>