[英]Create a band-pass filter via Scipy in Python?
有沒有一種方法可以在Python 3.6中通過scipy
或librosa
為16KHz wav文件創建快速帶通濾波器,以濾除300-3400Hz的人聲頻帶之外的噪聲? 這是一個樣本wav文件 ,具有低頻背景噪聲。
更新:是的,我已經看到/嘗試過如何使用Scipy.signal.butter實現帶通Butterworth濾波器 。 不幸的是,過濾后的聲音嚴重變形。 本質上,整個代碼執行此操作:
lo,hi=300,3400
sr,y=wavfile.read(wav_file)
b,a=butter(N=6, Wn=[2*lo/sr, 2*hi/sr], btype='band')
x = lfilter(b,a,y)
sounddevice.play(x, sr) # playback
我在做什么錯,或者如何改善此問題,以便正確濾除背景噪音。
這是使用上面的鏈接的原始文件和過濾文件的可視化。 可視化看起來很合理,但是聽起來很可怕:(如何解決?
顯然,寫入未標准化的64位浮點數據時會發生此問題。 通過將x
轉換為16位或32位整數,或將x
標准化為[-1,1]范圍並轉換為32個浮點數,我得到了一個聽起來合理的輸出文件。
我不使用sounddevice
; 相反,我將過濾后的數據保存到新的WAV文件中並進行播放。 以下是對我有用的變體:
# Convert to 16 integers
wavfile.write('off_plus_noise_filtered.wav', sr, x.astype(np.int16))
要么...
# Convert to 32 bit integers
wavfile.write('off_plus_noise_filtered.wav', sr, x.astype(np.int32))
要么...
# Convert to normalized 32 bit floating point
normalized_x = x / np.abs(x).max()
wavfile.write('off_plus_noise_filtered.wav', sr, normalized_x.astype(np.float32))
輸出整數時,可以按比例放大這些值,以最大程度地減少由於截斷浮點值而導致的精度損失:
x16 = (normalized_x * (2**15-1)).astype(np.int16)
wavfile.write('off_plus_noise_filtered.wav', sr, x16)
以下代碼用於從此處生成帶通濾波器: https : //scipy.github.io/old-wiki/pages/Cookbook/ButterworthBandpass
from scipy.signal import butter, lfilter
def butter_bandpass(lowcut, highcut, fs, order=5):
nyq = 0.5 * fs
low = lowcut / nyq
high = highcut / nyq
b, a = butter(order, [low, high], btype='band')
return b, a
def butter_bandpass_filter(data, lowcut, highcut, fs, order=5):
b, a = butter_bandpass(lowcut, highcut, fs, order=order)
y = lfilter(b, a, data)
return y
if __name__ == "__main__":
import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import freqz
# Sample rate and desired cutoff frequencies (in Hz).
fs = 5000.0
lowcut = 500.0
highcut = 1250.0
# Plot the frequency response for a few different orders.
plt.figure(1)
plt.clf()
for order in [3, 6, 9]:
b, a = butter_bandpass(lowcut, highcut, fs, order=order)
w, h = freqz(b, a, worN=2000)
plt.plot((fs * 0.5 / np.pi) * w, abs(h), label="order = %d" % order)
plt.plot([0, 0.5 * fs], [np.sqrt(0.5), np.sqrt(0.5)],
'--', label='sqrt(0.5)')
plt.xlabel('Frequency (Hz)')
plt.ylabel('Gain')
plt.grid(True)
plt.legend(loc='best')
# Filter a noisy signal.
T = 0.05
nsamples = T * fs
t = np.linspace(0, T, nsamples, endpoint=False)
a = 0.02
f0 = 600.0
x = 0.1 * np.sin(2 * np.pi * 1.2 * np.sqrt(t))
x += 0.01 * np.cos(2 * np.pi * 312 * t + 0.1)
x += a * np.cos(2 * np.pi * f0 * t + .11)
x += 0.03 * np.cos(2 * np.pi * 2000 * t)
plt.figure(2)
plt.clf()
plt.plot(t, x, label='Noisy signal')
y = butter_bandpass_filter(x, lowcut, highcut, fs, order=6)
plt.plot(t, y, label='Filtered signal (%g Hz)' % f0)
plt.xlabel('time (seconds)')
plt.hlines([-a, a], 0, T, linestyles='--')
plt.grid(True)
plt.axis('tight')
plt.legend(loc='upper left')
plt.show()
看看這是否有助於您的事業。 您可以在此處指定所需的頻率:
# Sample rate and desired cutoff frequencies (in Hz).
fs = 5000.0
lowcut = 500.0
highcut = 1250.0
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.