- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我目前有几千个音频剪辑需要通过机器学习进行分类。
经过一番挖掘,我发现如果对音频进行短时傅立叶变换,它就会变成二维图像,因此我可以对这些图像而不是音频文件本身使用各种图像分类算法。
为此我找到了一个python package这就是 STFT,我所需要的就是绘制它,这样我就可以获得图像。为了绘图,我发现 this github repo非常有用。
最后我的代码如下:
import stft
import scipy
import scipy.io.wavfile as wav
import matplotlib.pylab as pylab
def save_stft_image(source_filename, destination_filename):
fs, audio = wav.read(source_filename)
X = stft.spectrogram(audio)
print X.shape
fig = pylab.figure()
ax = pylab.Axes(fig, [0,0,1,1])
ax.set_axis_off()
fig.add_axes(ax)
pylab.imshow(scipy.absolute(X[:][:][0].T), origin='lower', aspect='auto', interpolation='nearest')
pylab.savefig(destination_filename)
save_stft_image("Example.wav","Example.png")
代码有效,但是我观察到,当 print X.shape
行执行时,我得到 (513L, 943L, 2L)
。所以结果是 3 维的。因此,当我只写 X[:][:][0]
或 X[:][:][1]
时,我会得到一个图像。
我一直在阅读 STFT 的“冗余”,您可以删除一半,因为您不需要它。第三维度是多余的还是我在这里做错了什么?如果是这样,我该如何正确绘制它?
谢谢。
编辑:所以新的代码和输出是:
import stft
import os
import scipy
import scipy.io.wavfile as wav
import matplotlib.pylab as pylab
def save_stft_image(source_filename, destination_filename):
fs, audio = wav.read(source_filename)
audio = scipy.mean(audio, axis = 1)
X = stft.spectrogram(audio)
print X.shape
fig = pylab.figure()
ax = pylab.Axes(fig, [0,0,1,1])
ax.set_axis_off()
fig.add_axes(ax)
pylab.imshow(scipy.absolute(X.T), origin='lower', aspect='auto', interpolation='nearest')
pylab.savefig(destination_filename)
save_stft_image("Example.wav","Example.png")
在左边我看到一列几乎看不见的颜色。我正在研究的声音是呼吸声,因此它们的频率非常低。也许这就是为什么可视化是一个非常薄的颜色列。
最佳答案
您可能有立体声音频文件?因此 X[:][:][0]
和 X[:][:][1]
对应于每个 channel 。
您可以通过scipy.mean(audio, axis=1)将多声道转换为单声道
。
关于python - 在Python中绘制音频频谱图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36893767/
我对复值频谱执行 iFFT 并通过让第一个样本归零来更改相应的时域信号。最后,我通过 FFT 将其转换回频域。 我想知道这里使用两侧(对称)频谱或仅使用一侧频谱(仅正频率)之间的(物理)差异在哪里,因
我正在尝试提高使用 cython 计算 Jonswap 频谱的性能。但它似乎比原始代码慢得多。我该如何改进? 赛通代码: from libc.math cimport exp from libc.st
我目前可以使用 JavaScript Web Audio API 播放音轨。在播放此轨道时,我可以提取 FFT 频谱数据、峰值、RMS 值等。 但是,出于我的应用程序的目的,我需要能够在开始播放轨道之
我正在尝试将频谱保存在我的 FMOD_DSP_PARAMETER_FFT 中,但我只收到充满零的频谱,如果你能看到我的错误我会同意,我认为我没有很好地将 DSP 连接到 channel 或类似的东西,
我有一个复杂的信号,我想对其进行 FFT。使用 Mathematica 我得到以下结果: (* Some acquisition params *) fS = 100. 10^6; time = 10
我是一名优秀的程序员,十分优秀!