Documentation
RfPy
is a package containing Python tools for calculating teleseismic
receiver functions. Methods are available to plot and post-process
receiver function data for use in crustal and upper mantle studies.
Current post processing methods include: H-k stacking,
Harmonic decomposition and CCP imaging. The code uses
the StDb
package for querying and building a station database
used in command-line scripts.
Note
RfPy
was written independently of the Python software rf,
although it is possible that some classes and methods defined here might be applied
to rf
objects, since both are heavily based on obspy. The main
differences between RfPy
and rf
are as follows:
RfPy
employs either a Wiener or multitaper spectral deconvolution technique, whereasrf
uses either a water-level spectral deconvolution or a time-domain deconvolution.RfData
objects are used to calculate single-station and single-event receiver functions, whereasrf
can handle multiple stations at once.
Warning
RfPy
was recently updated to fix a problem when running the scripts under Windows OS. The consequence is that version 0.1.1
will throw an error if the extension .py is specified when calling the scripts. The accompanying documentation also uses version 0.2.1
of StDb
in the Tutorials section.