RF package#
seispy.eq module#
- class seispy.eq.EQ(pathname, datestr, suffix='SAC')[source]#
Bases:
object
Methods
arr_correct
([write_to_sac])offset = sac.b - real o
channel_correct
([switchEN, reverseE, reverseN])_summary_
deconvolute
(shift, time_after[, f0, method, ...])Deconvolution
trim
(time_before, time_after[, isreturn])offset = sac.b - real o
cleanstream
decon_p
decon_s
detrend
filter
fix_channel_name
from_stream
get_arrival
get_time_offset
judge_rf
phase_trigger
readstream
rotate
s_condition
saverf
search_baz
search_inc
set_comp
snr
write
- channel_correct(switchEN=False, reverseE=False, reverseN=False)[source]#
_summary_
- Parameters:
switchEN (bool, optional) – _description_, defaults to False
reverseE (bool, optional) – _description_, defaults to False
reverseN (bool, optional) – _description_, defaults to False
- deconvolute(shift, time_after, f0=2.0, method='iter', only_r=False, itmax=400, minderr=0.001, wlevel=0.05, target_dt=None)[source]#
Deconvolution
- Parameters:
- shiftfloat
Time shift before P arrival
- time_afterfloat
Time length after P arrival
- f0float or list, optional
Gaussian factors, by default 2.0
- methodstr, optional
method for deconvolution in
iter
orwater
, by defaultiter
- only_rbool, optional
Whether only calculate RF in prime component, by default False
- itmaxint, optional
Maximum iterative number, valid for method of
iter
, by default 400- minderrfloat, optional
Minium residual error, valid for method of
iter
, by default 0.001- wlevelfloat, optional
Water level, valid for method of
water
, by default 0.05- target_dtNone or float, optional
Time delta for resampling, by default None
seispy.rf module#
- class seispy.rf.RF(cfg_file=None, log=None)[source]#
Bases:
object
- Attributes:
- date_begin
- date_end
Methods
baz_correct
cal_phase
channel_correct
deconv
detrend
drop_eq_snr
filter
load_stainfo
loadpjt
match_eq
pick
rotate
save_raw_data
savepjt
saverf
search_eq
trim
- property date_begin#
- property date_end#
- seispy.rf.match_eq(eq_lst, pathname, stla, stlo, logger, ref_comp='Z', suffix='SAC', offset=None, tolerance=210, dateformat='%Y.%j.%H.%M.%S')[source]#
- seispy.rf.read_catalog(logpath: str, b_time, e_time, stla: float, stlo: float, magmin=5.5, magmax=10.0, dismin=30.0, dismax=90.0, depthmin=0, depthmax=800)[source]#
Read local catalog with seispy or QUAKEML format
- Parameters:
logpath (str) – Path to catalogs
b_time (obspy.UTCDateTime) – Start time
e_time (obspy.UTCDateTime) – End time
stla (float) – Station latitude
stlo (float) – Station longitude
magmin (float, optional) – Minimum magnitude, defaults to 5.5
magmax (float, optional) – Maximum magnitude, defaults to 10
dismin (float, optional) – Minimum distance, defaults to 30
dismax (float, optional) – Maximum distance, defaults to 90
- Returns:
list of earthquakes
- Return type:
pandas.DataFrame
seispy.para module#
- class seispy.para.RFPara[source]#
Bases:
object
- Attributes:
- catalogpath
- criterion
- datapath
- decon_method
- rfpath
- server
Methods
get_para
read_para
- property catalogpath#
- property criterion#
- property datapath#
- property decon_method#
- property rfpath#
- property server#
- class seispy.para.StaInfo[source]#
Bases:
object
Methods
get_station_from_ws
([server])Get station information from web-service with given network and station or other optional condition.
get_stainfo
load_stainfo