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

Question

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

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

Answer 1

抱歉，我没有关于 pywavelet 的直接答案，但以防万一出现另一个答案需要一些时间，您可能需要考虑我用来创建 Ricker 小波的旧配方（使用名为的软件构建合成地震图地球仪2）：

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