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

Question

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

由于函数形式需要放入Ricker小波，所以我一直无法在此示例基础上使用pywavelets创建自定义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