вычисление FFT в кадрах и запись в файл - PullRequest
Новая
       107

вычисление FFT в кадрах и запись в файл

0 голосов
/ 19 апреля 2020

Я новичок в python, я пытаюсь получить значение FFT загруженного файла WAV и вернуть FFT каждого кадра в каждой строке текстового файла (используя GCP)

, используя scipy или librosa

Требуемая частота кадров составляет 30 кадров в секунду

. Волновой файл будет иметь частоту дискретизации 48 КБ

, поэтому мои вопросы

  1. как разделить выборки для всего файла WAV на выборки каждого кадра
  2. Как добавить пустые выборки, чтобы длина выборок кадра была равна 2 (как 48000/30 = 1600 добавить 448 пустых образцов, чтобы сделать его 2048)
  3. как мне нормализовать результирующий массив FFT на [-1,1]?

1 Ответ

1 голос
/ 20 апреля 2020

Вы можете использовать pyaudio с обратным вызовом, чтобы получить то, что вы делаете. 

import pyaudio
import wave
import time
import struct
import sys
import numpy as np

if len(sys.argv) < 2:
   print("Plays a wave file.\n\nUsage: %s filename.wav" % sys.argv[0])
   sys.exit(-1)

wf = wave.open(sys.argv[1], 'rb')

# instantiate PyAudio (1)
p = pyaudio.PyAudio()

def callback_test(data, frame_count, time_info, status):
    frame_count =1024
    elm = wf.readframes(frame_count) # read n frames
    da_i = np.frombuffer(elm, dtype='<i2') # convert to little endian int pairs
    da_fft = np.fft.rfft(da_i) # fast fourier transform for real values

    da_ifft = np.fft.irfft(da_fft)  # inverse fast fourier transform for real values 
    da_i = da_ifft.astype('<i2') # convert to little endian int pairs
    da_m = da_i.tobytes() # convert to bytes 
    return (da_m, pyaudio.paContinue)

# open stream using callback (3)
stream = p.open(format=p.get_format_from_width(wf.getsampwidth()),
                channels=wf.getnchannels(),
                rate=wf.getframerate(),# sampling frequency
                output=True,
                stream_callback=callback_test)

# # start the stream (4)
stream.start_stream()

# # wait for stream to finish (5)
while stream.is_active():
    time.sleep(0.1)

# # stop stream (6)
stream.stop_stream()
stream.close()
wf.close()

# close PyAudio (7)
p.terminate()

Пожалуйста, обратитесь к этим ссылкам для дальнейшего изучения:

https://people.csail.mit.edu/hubert/pyaudio/docs/#example -callback-mode- audio-io и Python изменение высоты тона wav-файла

...