How to import my own SAC data and calculate SRF

[trichter]

1. Please have look here:

RF, reading SAC downloaded via SOD

Fairly new to ObsPy; sorry a likely naive question, I have SAC files (network=US; station name=GLMI; components; BH1, BH2, BHZ) downloaded using SOD and want to calculate a receiver function using RF. Following the manual, I succeed to read: from rf import read_rf, rfstats, RFStream data = read_rf(‘US.GLMI.00.BH?.sac’) stream = RFStream() print('aa ', data[0].stats.sac.kstnm) returns: aa GLMI but the command: stats = rfstats(station=station, event=event, phase=‘P’, dist_range=(30,90)) r…

2. Probably not, rf only supports waveforms which can be read by ObsPy.