如何使用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