# Related: rf

**URL:** https://discourse.obspy.org/c/obspy-related-projects/rf/14.md

[Latest](https://discourse.obspy.org/latest.md) · [Categories](https://discourse.obspy.org/categories.md) · [Tags](https://discourse.obspy.org/tags.md)

---

## [About the Related: rf category](https://discourse.obspy.org/t/about-the-related-rf-category/1011)

<div class="topic-metadata">

**Author:** [@megies](https://discourse.obspy.org/u/megies)\
**Replies:** 0

</div>

rf – Receiver function calculation in seismology rf is a Python framework for receiver function analysis. Read and write support of necessary metadata is provided for SAC, SeismicHandler and HDF5 waveform files. Data is…

---

## [Rfstats function returns empty string](https://discourse.obspy.org/t/rfstats-function-returns-empty-string/2019)

<div class="topic-metadata">

**Author:** [@Ali\_Aslanli](https://discourse.obspy.org/u/Ali_Aslanli)\
**Replies:** 11\
**Last updated:** [August 10, 2026, 12:05pm UTC](https://discourse.obspy.org/t/rfstats-function-returns-empty-string/2019 "2026-08-10T12:05:06Z")

</div>

Hey fellas, Fairly new to the seismology domain, trying to calculate receiver functions using RF package (making some experiments with local high frequency events), however, I get an empty array after executing rfstats(…

---

## [Plot\_profile - S RFs](https://discourse.obspy.org/t/plot-profile-s-rfs/2358)

<div class="topic-metadata">

**Author:** [@seismolab-uct](https://discourse.obspy.org/u/seismolab-uct)\
**Replies:** 1\
**Last updated:** [February 24, 2026, 2:41pm UTC](https://discourse.obspy.org/t/plot-profile-s-rfs/2358 "2026-02-24T14:41:15Z")

</div>

Hi everyone, working on Sp RFs and found that plotting the profile with plot\_profile results in no positive thickmarks for both time and depth axes, which is not the case as for Ps RFs (like in the tutorials). See the fi…

---

## [How to identify the Ps converted phase from receiver functions](https://discourse.obspy.org/t/how-to-identify-the-ps-converted-phase-from-receiver-functions/2279)

<div class="topic-metadata">

**Author:** [@Sarmat\_Maqsat](https://discourse.obspy.org/u/Sarmat_Maqsat)\
**Replies:** 1\
**Last updated:** [August 5, 2025, 11:08am UTC](https://discourse.obspy.org/t/how-to-identify-the-ps-converted-phase-from-receiver-functions/2279 "2025-08-05T11:08:25Z")

</div>

I am an undergraduate student working on an array of seismic stations in Antarctica. I am very new to Seismology and receiver functions, and relatively new to python programming. I’ve used the rf package and followed the…

---

## [Plot\_profile trimming traces](https://discourse.obspy.org/t/plot-profile-trimming-traces/2158)

<div class="topic-metadata">

**Author:** [@TobiasRocks](https://discourse.obspy.org/u/TobiasRocks)\
**Replies:** 3\
**Last updated:** [February 20, 2025, 3:48pm UTC](https://discourse.obspy.org/t/plot-profile-trimming-traces/2158 "2025-02-20T15:48:50Z")

</div>

Hi all, in my attempts to learn the rf package, I’ve been following the tutorials provided on the docs, and found that the plot\_profile() function seems to be trimming the traces. Here’s what I mean: Here’s what the pl…

---

## [Gaussian parameter (iterative deconvolution)](https://discourse.obspy.org/t/gaussian-parameter-iterative-deconvolution/1886)

<div class="topic-metadata">

**Author:** [@Ademar](https://discourse.obspy.org/u/Ademar)\
**Replies:** 6\
**Last updated:** [January 24, 2025, 1:53pm UTC](https://discourse.obspy.org/t/gaussian-parameter-iterative-deconvolution/1886 "2025-01-24T13:53:43Z")

</div>

Dear all, I am a user of ITERDECON (from Dr. Ammon) and SACITERD (from Dr. Herrmann), and I realized that the RFs from both programs with a Gaussian parameter of 2.5 resemble the traces of the rf package when the Gaussi…

---

## [Trimming traces with different start times](https://discourse.obspy.org/t/trimming-traces-with-different-start-times/2132)

<div class="topic-metadata">

**Author:** [@Ali\_Aslanli](https://discourse.obspy.org/u/Ali_Aslanli)\
**Replies:** 1\
**Last updated:** [January 24, 2025, 1:13pm UTC](https://discourse.obspy.org/t/trimming-traces-with-different-start-times/2132 "2025-01-24T13:13:53Z")

</div>

Hello everyone, I have traces with start times differing in microseconds with the same onsets calculated by rfstats(). When I cut around the onsets, start times do not align which is problematic for the RF calculation. …

---

## [Rotation: Peak at 0 s on the Q component](https://discourse.obspy.org/t/rotation-peak-at-0-s-on-the-q-component/1745)

<div class="topic-metadata">

**Author:** [@rafa0110](https://discourse.obspy.org/u/rafa0110)\
**Replies:** 1\
**Last updated:** [August 6, 2023, 6:47pm UTC](https://discourse.obspy.org/t/rotation-peak-at-0-s-on-the-q-component/1745 "2023-08-06T18:47:58Z")

</div>

Hi Tom! By chance do you know some articles that describe the non-perfect rotation problem? I am trying to understand a little bit more the peak at 0 on the Q component, but I can not find any literature. I look for the…

---

## [S receiver functions Peaks](https://discourse.obspy.org/t/s-receiver-functions-peaks/1715)

<div class="topic-metadata">

**Author:** [@rafa0110](https://discourse.obspy.org/u/rafa0110)\
**Replies:** 2\
**Last updated:** [May 24, 2023, 11:24am UTC](https://discourse.obspy.org/t/s-receiver-functions-peaks/1715 "2023-05-24T11:24:59Z")

</div>

Hello! I am working with the SRF and I have noticed that the move-out correction causes some peaks for later times(stream3c.moveout(phase=‘Sp’)). I think it is caused specifically by the stretch\_delay\_times() . Are th…

---

## [Energy conversion: two interphases](https://discourse.obspy.org/t/energy-conversion-two-interphases/1711)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 1\
**Last updated:** [May 22, 2023, 9:43am UTC](https://discourse.obspy.org/t/energy-conversion-two-interphases/1711 "2023-05-22T09:43:57Z")

</div>

Hi Tom! I have a theoretical question. I sketched the situation in the attached plot. With my data, I am able to reproduce the sketched situation, which was expected. Where we expect double signals is where I have th…

---

## [Plotting receiver functions for thesis/publications](https://discourse.obspy.org/t/plotting-receiver-functions-for-thesis-publications/1668)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 1\
**Last updated:** [March 22, 2023, 7:39am UTC](https://discourse.obspy.org/t/plotting-receiver-functions-for-thesis-publications/1668 "2023-03-22T07:39:42Z")

</div>

Hi Tom! Just a quick question: For plotting your receiver functions for your thesis did you output them directly from Python or did you used something like GMT for making the plots? If you used GMT, do you know where I …

---

## [QC receiver functions](https://discourse.obspy.org/t/qc-receiver-functions/1635)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 4\
**Last updated:** [February 20, 2023, 5:38pm UTC](https://discourse.obspy.org/t/qc-receiver-functions/1635 "2023-02-20T17:38:01Z")

</div>

Hi Tom! As I see it is more convenient to do the QC right before deconvolution. For processing the RF, some of the lines of my code are: Script 1 example stream3c.rf() stream3c.moveout() stream.extend(stream3c) Sin…

---

## [Can't rotate seismograms for some events, works for example events though](https://discourse.obspy.org/t/cant-rotate-seismograms-for-some-events-works-for-example-events-though/1616)

<div class="topic-metadata">

**Author:** [@csogden](https://discourse.obspy.org/u/csogden)\
**Replies:** 7\
**Last updated:** [February 6, 2023, 5:13pm UTC](https://discourse.obspy.org/t/cant-rotate-seismograms-for-some-events-works-for-example-events-though/1616 "2023-02-06T17:13:07Z")

</div>

Hi Tom, I’ve been working on other projects since we chatted via email but to pick that thread up from November. I am using your example template for calculating basic RFs for a chosen event. I can run your example data…

---

## [S receiver functions](https://discourse.obspy.org/t/s-receiver-functions/1576)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 5\
**Last updated:** [December 5, 2022, 10:24pm UTC](https://discourse.obspy.org/t/s-receiver-functions/1576 "2022-12-05T22:24:04Z")

</div>

Hi Tom!! Now, I am trying to produce a SRF stacked secction. When computing the stats, instead of phase=‘P’ I used phase=‘S’. For example: for s in iter\_event\_data(catalog,inventory,client.get\_waveforms, phase=‘S’,pbar…

---

## [Variable piercing point (ppoints) depth](https://discourse.obspy.org/t/variable-piercing-point-ppoints-depth/1577)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 1\
**Last updated:** [December 2, 2022, 1:12pm UTC](https://discourse.obspy.org/t/variable-piercing-point-ppoints-depth/1577 "2022-12-02T13:12:11Z")

</div>

Hi (again) Tom!! I was looking at the stacked section in your dissertation, and I was wondering if you instead of using a fixed ppoints you used ppoints that change with depth (considering that it was a subduccion z…

---

## [Batch processing: get the profile](https://discourse.obspy.org/t/batch-processing-get-the-profile/1558)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 1\
**Last updated:** [November 19, 2022, 10:21pm UTC](https://discourse.obspy.org/t/batch-processing-get-the-profile/1558 "2022-11-19T22:21:11Z")

</div>

Thanks you Tom. Just one last question, now I am using the batch option. I manage to stack the data, but to generate the profile I don’t get any folder from it. It loads something but not folder is created. How I shou…

---

## [Time vs. slowness](https://discourse.obspy.org/t/time-vs-slowness/1556)

<div class="topic-metadata">

**Author:** [@Qian](https://discourse.obspy.org/u/Qian)\
**Replies:** 1\
**Last updated:** [November 18, 2022, 1:17pm UTC](https://discourse.obspy.org/t/time-vs-slowness/1556 "2022-11-18T13:17:23Z")

</div>

Hi Mr. Tom Richter. I am working with your scripts to generate stacks of RF. I want to generate a stacked section, but instead of time vs. distance, I want to plot time vs. slowness (as in the figure). Do you thing …

---

## [Piercing points depth](https://discourse.obspy.org/t/piercing-points-depth/1545)

<div class="topic-metadata">

**Author:** [@rafa0110](https://discourse.obspy.org/u/rafa0110)\
**Replies:** 1\
**Last updated:** [November 5, 2022, 9:17pm UTC](https://discourse.obspy.org/t/piercing-points-depth/1545 "2022-11-05T21:17:50Z")

</div>

Hi Tom! I have two questions: I am wondering which criteria can I use for choosing for depth for stacking piercing points (in ppoints = stream.ppoints(70)). You choose 70 km depth for your example, but I guess this s…

---

## [Remove\_response of the instrument](https://discourse.obspy.org/t/remove-response-of-the-instrument/1489)

<div class="topic-metadata">

**Author:** [@rafa0110](https://discourse.obspy.org/u/rafa0110)\
**Replies:** 1\
**Last updated:** [September 6, 2022, 11:08am UTC](https://discourse.obspy.org/t/remove-response-of-the-instrument/1489 "2022-09-06T11:08:06Z")

</div>

Hi Tom! I am using your scripts for working on my master thesis. At the moment, I am using script number 2 to obtain the stacked PRF. I have a question: How to include the ‘remove\_response’ function on the scripts? Th…

---

## [Noise before P-onset in iterative deconvolution](https://discourse.obspy.org/t/noise-before-p-onset-in-iterative-deconvolution/1385)

<div class="topic-metadata">

**Author:** [@Thea92](https://discourse.obspy.org/u/Thea92)\
**Replies:** 2\
**Last updated:** [March 29, 2022, 8:19am UTC](https://discourse.obspy.org/t/noise-before-p-onset-in-iterative-deconvolution/1385 "2022-03-29T08:19:16Z")

</div>

Hi all, for some receiver function techniques (e.g waterlevel), the noise ahead the P-onset can be used to estimate the reliability of the peaks at 0+ seconds. Do you know, how the iterative deconvolution is working an…

---

## [RF calculation using PREM model](https://discourse.obspy.org/t/rf-calculation-using-prem-model/1357)

<div class="topic-metadata">

**Author:** [@Thea92](https://discourse.obspy.org/u/Thea92)\
**Replies:** 3\
**Last updated:** [January 27, 2022, 8:20am UTC](https://discourse.obspy.org/t/rf-calculation-using-prem-model/1357 "2022-01-27T08:20:05Z")

</div>

Hi all, at the moment I’m trying to calculate iterative RFs by using the PREM velocity model. For some distances I get ‘empty’ receiver functions, it is happening when doing: rfstats(stream, phase=phase, dist\_range=dis…

---

## [Peak at zero in Q component receiver function](https://discourse.obspy.org/t/peak-at-zero-in-q-component-receiver-function/1308)

<div class="topic-metadata">

**Author:** [@sahilansari3](https://discourse.obspy.org/u/sahilansari3)\
**Replies:** 3\
**Last updated:** [November 9, 2021, 9:47am UTC](https://discourse.obspy.org/t/peak-at-zero-in-q-component-receiver-function/1308 "2021-11-09T09:47:52Z")

</div>

Hello everyone. I have made my receiver function using the rf package. According to my knowledge in the Q receiver function, there should not be any peak at zero but I got a peak at zero. What are the possible reasons fo…

---

## [Gauss Parameter on Water Level Deconvolution](https://discourse.obspy.org/t/gauss-parameter-on-water-level-deconvolution/1261)

<div class="topic-metadata">

**Author:** [@Indra\_Siregar](https://discourse.obspy.org/u/Indra_Siregar)\
**Replies:** 6\
**Last updated:** [June 21, 2021, 11:26am UTC](https://discourse.obspy.org/t/gauss-parameter-on-water-level-deconvolution/1261 "2021-06-21T11:26:11Z")

</div>

Dear all, Let me ask about the Gaussian parameter on the water-level method. How is the relationship between the gauss parameter and the alpha value? I still don’t understand, because generally, this method uses alpha v…

---

## [How to use trim and multitaper deconvolution for SRF calculation](https://discourse.obspy.org/t/how-to-use-trim-and-multitaper-deconvolution-for-srf-calculation/1244)

<div class="topic-metadata">

**Author:** [@BSPRao](https://discourse.obspy.org/u/BSPRao)\
**Replies:** 6\
**Last updated:** [June 8, 2021, 6:24am UTC](https://discourse.obspy.org/t/how-to-use-trim-and-multitaper-deconvolution-for-srf-calculation/1244 "2021-06-08T06:24:30Z")

</div>

Hello, I tried to use your rf python code to compute the S receiver functions to decipher the Moho and LAB beneath my study region. I downloaded your code and installed all the necessary packages and it is working well.…

---

## [How to import my own SAC data and calculate SRF](https://discourse.obspy.org/t/how-to-import-my-own-sac-data-and-calculate-srf/1229)

<div class="topic-metadata">

**Author:** [@geoAP90](https://discourse.obspy.org/u/geoAP90)\
**Replies:** 11\
**Last updated:** [May 20, 2021, 5:45pm UTC](https://discourse.obspy.org/t/how-to-import-my-own-sac-data-and-calculate-srf/1229 "2021-05-20T17:45:55Z")

</div>

I need to understand how should I import my own data and make the figures I have a set of calculated and selected P-RF. Can I import the .eqr files and just plot the figures using RF software?

---

## [Receiver functions for PKP-arrivals](https://discourse.obspy.org/t/receiver-functions-for-pkp-arrivals/1233)

<div class="topic-metadata">

**Author:** [@Thea92](https://discourse.obspy.org/u/Thea92)\
**Replies:** 4\
**Last updated:** [May 19, 2021, 7:38am UTC](https://discourse.obspy.org/t/receiver-functions-for-pkp-arrivals/1233 "2021-05-19T07:38:02Z")

</div>

Hello everyone, is it possible to calculate Receiver Functions with the PKP-phases using the rf-module? When I try to create the RFstream with stream = RFStream(zrt) and rfstats(stream), the traces are not added to the…

---

## [Trying Notebook2](https://discourse.obspy.org/t/trying-notebook2/1227)

<div class="topic-metadata">

**Author:** [@geoAP90](https://discourse.obspy.org/u/geoAP90)\
**Replies:** 5\
**Last updated:** [May 4, 2021, 1:36pm UTC](https://discourse.obspy.org/t/trying-notebook2/1227 "2021-05-04T13:36:19Z")

</div>

I have just started with the RF package and trying to follow the steps from Rf documentation. I tried the examples given in notebook1 and could make the figures easily but I am not getting any figures while running the n…

---

## [RF iterative deconvolution](https://discourse.obspy.org/t/rf-iterative-deconvolution/1221)

<div class="topic-metadata">

**Author:** [@Thea92](https://discourse.obspy.org/u/Thea92)\
**Replies:** 15\
**Last updated:** [April 15, 2021, 2:13pm UTC](https://discourse.obspy.org/t/rf-iterative-deconvolution/1221 "2021-04-15T14:13:07Z")

</div>

Hi there, I’m using this nice rf tool with obspy. I’m trying to deconvolve the data using the iterative deconvolution. I followed the documentation: RF=rf.deconvolve.deconvolve(RF, method='iterative', gauss=0.25, itmax…

---

## [RF using channels BHZ, BH1, and BH2](https://discourse.obspy.org/t/rf-using-channels-bhz-bh1-and-bh2/1140)

<div class="topic-metadata">

**Author:** [@jritsema](https://discourse.obspy.org/u/jritsema)\
**Replies:** 1\
**Last updated:** [November 26, 2020, 8:07am UTC](https://discourse.obspy.org/t/rf-using-channels-bhz-bh1-and-bh2/1140 "2020-11-26T08:07:40Z")

</div>

Is it possible to use stream.rf() when the horizontal components are called BH1 and BH2 and when CMPAZ is defined in the SAC headers?

---

## [Trim2 of traces with respect to a theoretical arrival time (from taup)](https://discourse.obspy.org/t/trim2-of-traces-with-respect-to-a-theoretical-arrival-time-from-taup/1139)

<div class="topic-metadata">

**Author:** [@jritsema](https://discourse.obspy.org/u/jritsema)\
**Replies:** 2\
**Last updated:** [November 25, 2020, 1:53pm UTC](https://discourse.obspy.org/t/trim2-of-traces-with-respect-to-a-theoretical-arrival-time-from-taup/1139 "2020-11-25T13:53:05Z")

</div>

I have a stream with three component traces. The field stream.stats.onset is filled with a time close to the theoretical P-time (which is filled in stats.sac.t0). I would like to trim the traces with respect to the theor…

[Next page](https://discourse.obspy.org/c/obspy-related-projects/rf/14.md?page=1)