def Envelope(syn, obs, nt, dt, eps=0.05):
# envelope difference
# (Yuan et al 2015, eq 16)
esyn = abs(_analytic(syn))
eobs = abs(_analytic(obs))
if esyn.all() == 0.0:
etmp = 0.0
else:
etmp = (esyn - eobs)/(esyn + eps*esyn.max())
wadj = etmp*syn - _np.imag(_analytic(etmp*_np.imag(_analytic(syn))))
return wadj
评论列表
文章目录