Seispy is a Python module for processing seismological data and calculating Receiver Functions. The advanced functions are available to improve the Obspy.
- Python >= 3.8
- ObsPy >= 1.2.0
- NumPy >= 1.16
- SciPy >= 1.2.0
- matplotlib >= 3.6.0
- PySide6 >= 6.3.0
- scikits.bootstrap >= 1.0.0
See Seispy documentation in detail.
seispy.distaz
: Calculate distance and azimuth credited by the lithospheric seismology program at USC, butnumpy.ndarray
operations are supported.seispy.geo
: Tiny codes of geophysics.seispy.decon
: Functions of deconvolution transferred from iwbailey/processRFmatlab including- Iterative time domain deconvolution method (Ligorría and Ammon 1999 BSSA).
- Water level frequency domain deconvolution method (CJ. Ammon 1991 BSSA)
seispy.rf
: Procedure for RF calculation. The functions ofmatch_eq
,search_eq
invokedobspy.core.UTCDateTime
andobspy.clients
from the Obspy.seispy.eq
: RF processing for each event, which invokedobspy.io.sac
,obspy.signal
,obspy.taup
andobspy.core.Stream
from the Obspy.seispy.hk
: H-k stacking for single station (Zhu and Kanamori 2000 JGR).seispy.rfani
: A joint method for crustal anisotropic calculation (Liu and Niu 2011 GJI).seispy.slantstack
: Slant stacking for single station (Tauzin et al., 2008)seispy.rfcorrect
: Subsequent process of RFs including moveout correction and time to depth conversion (1D and 3D) (see Xu et al., 2018 EPSL)seispy.ccpprofile
: CCP stacking along a profile.seispy.ccp3d
: 3-D CCP stacking with extracting depth D410 and D660.
prf
: Calculate PRFs for a station.pickrf
: Pick PRFs with virtual quality control after the calculation.plotrt
: Plot PRFs with R and T components order by back-azimuth.plotr
: Plot PRFs with R component order by back-azimuth.hk
: H-Kappa stacking for estimating Moho depth and crustal Vp/Vs.rf2depth
: Convert PRFs to depth axis.ccp_profile
: Stack PRFs along a profile with a CCP stacking method.ccp3d
: Stack PRFs with spaced bins.rfani
: Estimating crustal anisotropy with a joint method.rfharmo
: Harmonic decomposition to extract constant component of RF and plot dip/anisotropic components.
veltxt2mod
: Create 3D velocity model withnumpy.lib.npyio.NpzFile
format from a ASCII table file.ndk2dat
: Convert the GCMT catalog file ("ndk" format) to the list file for theprf
command.updatecatalog
: Automatically update the GCMT catalog.setpar
: Set up the values in configure files.