Post-processing package

seispy.hk module

seispy.hk.ci(allstack, h, kappa, ev_num)

Search best H and kappa from stacked matrix. Calculate error for H and kappa :param allstack: stacked HK matrix :param h: 1-D array of H :param kappa: 1-D array of kappa :param ev_num: event number :return:

seispy.hk.hk()

seispy.hk.hk_sta_test()

seispy.hk.hksta(hpara, isplot=False, isdisplay=False)

seispy.hk.hkstack(seis, t0, dt, p, h, kappa, vp=6.3, weight=(0.7, 0.2, 0.1))

seispy.hk.hktest()

seispy.hk.plot(stack, allstack, h, kappa, besth, bestk, cvalue, cmap=<matplotlib.colors.ListedColormap object>, title=None, path=None)

seispy.hk.print_result(besth, bestk, maxhsig, maxksig, print_comment=True)

seispy.hk.time2idx(times, ti0, dt)

seispy.hk.tppps(depth, eta_p, eta_s)

seispy.hk.tps(depth, eta_p, eta_s)

seispy.hk.tpsps(depth, eta_s)

seispy.hk.transarray(array, axis=0)

seispy.hk.vslow(v, rayp)

class seispy.hkpara.HKPara

Bases: object

Attributes:

hrange

krange

property hrange

property krange

seispy.hkpara.hkpara(cfg_file)

seispy.slantstack module

class seispy.slantstack.SlantStack(seis, timeaxis, dis)

Bases: object

Methods

plot([cmap, xlim, vmin, vmax, figpath, colorbar])

Imaging for slant stacking

stack
syn_tps

plot(cmap='jet', xlim=None, vmin=None, vmax=None, figpath=None, colorbar=True)

Imaging for slant stacking

Parameters:

cmapstr, optional

The colormap to use, by default ‘jet’

xlim_type_, optional

Set the x limits of the current axes., by default None

vmin_type_, optional

Normalize the minium value data to the vmin, by default None

vmax_type_, optional

Normalize the maximum value data to the vmax, by default None

figpath_type_, optional

Output path to the image, by default None

colorbarbool, optional

Wether plot the colorbar, by default True

stack(ref_dis=None, rayp_range=None, tau_range=None)

syn_tps(phase_list, velmodel='iasp91', focal_dep=10)

seispy.rfani module

class seispy.rfani.RFAni(sacdatar, tb, te, tlen=3, val=10, rayp=0.06, model='iasp91')

Bases: object

Methods

plot_polar([cmap, show, outpath])

Polar map of crustal anisotropy.

baz_stack
cut_energy_waveform
init_ani_para
joint_ani
plot_correct
plot_stack_baz
radial_energy_max
rotate_to_fast_slow
search_peak
search_peak_amp
search_peak_list
xyz2grd

baz_stack(val=10)

cut_energy_waveform(idx, nb, ne)

init_ani_para()

joint_ani(weight=[0.5, 0.3, 0.2])

plot_correct(fvd=0, dt=0.44, enf=80, outpath=None)

plot_polar(cmap=<matplotlib.colors.ListedColormap object>, show=False, outpath='./')

Polar map of crustal anisotropy. See Liu and Niu (2012, doi: 10.1111/j.1365-246X.2011.05249.x) in detail.

Parameters:

• cmap (optional) – Colormap of matplotlib, defaults to ‘rainbow’

• show (bool, optional) – If show the polar map in the Matplotlib window, defaults to True

• outpath (str, optional) – Output path to saving the figure. If show the figure in the Matplotlib window, this option will be invalid, defaults to current directory.

plot_stack_baz(enf=60, outpath='./')

radial_energy_max()

rotate_to_fast_slow()

search_peak(energy, opt='max')

search_peak_amp()

search_peak_list(energy, opt='max')

xyz2grd(energy)

seispy.rfani.average_delay(fd, td)

seispy.rfani.joint_stack(energy_r, energy_cc, energy_tc, weight=[0.4, 0.3, 0.3])

seispy.harmonics module

class seispy.harmonics.Harmonics(rfsta, tmin=-5, tmax=10)

Bases: object

Methods

cut_trace()	Trim RFs from tmin to tmax
plot([outpath, enf])	Plot harmonic and unmodeled components
write_constant([out_sac_path])	Write constant component to SAC file.

harmo_trans

cut_trace()

Trim RFs from tmin to tmax

harmo_trans()

plot(outpath='./', enf=2.0)

Plot harmonic and unmodeled components

Parameters:

• outpath (str, optional) – Output path, defaults to ‘./’

• enf (float, optional) – Amplification factor, defaults to 2.0

write_constant(out_sac_path='./')

Write constant component to SAC file.

Parameters:

out_sac_path (str, optional) – Output path, defaults to ‘./’

previous

RF package

next

CCP package