如何使用pywavelets定义自定义Ricker小波？

Question

我一直在努力使用pywavelets定义自定义小波。我想使用ricker wavelet执行DWT。然而，ricker小波似乎不是pywavelets提供的内置小波族的一部分。官方的pywavelet文档给出了如何定义自定义Haar小波的示例。 http://www.pybytes.com/pywavelets/ref/wavelets.html

由于需要为Ricker小波设置函数形式，我无法在此示例的基础上使用pywavelet创建自定义Ricker小波。

Answer 1

很抱歉，我没有关于pywavelet的直接答案，但是为了防止出现另一个答案需要花费一些时间你可能想要考虑我用来制作Ricker小波的旧配方（用于构建合成地震图一个名为GEOMS2）的软件：

def ricker_wavelet(f, size, dt=1):
    t = np.int_(np.linspace(-size, size, (2*size+1)/dt))
    y = (1.0 - 2.0*(np.pi**2)*(f**2)*(t**2)) * np.exp(-(np.pi**2)*(f**2)*(t**2))
    data = np.hstack((t[:,np.newaxis],y[:,np.newaxis]))
    return data