使用python_speech_features提取音频文件特征

1. 读取wav文件

使用scipy.io.wavfile

import scipy.io.wavfile as wav
fs, signal = wav.read(filename)

fs是wav文件的采样率，signal是wav文件的内容，filename是要读取的音频文件的路径。我们将signal绘制出来就是下图这个样子。

image.png

2. 使用python_speech_features提取特征

①MFCC：

默认提取的特征维度是13，通常的做法是将该特征进行一阶差分和二阶差分，并将结果进行合并。

from python_speech_features import *
import numpy as np
def get_mfcc(data, fs):
    wav_feature =  mfcc(data, fs)
    d_mfcc_feat = delta(wav_feature, 1)
    d_mfcc_feat2 = delta(wav_feature, 2)
    feature = np.hstack((wav_feature, d_mfcc_feat, d_mfcc_feat2))
    return feature

参数介绍：
内容来源于 金泽夕
https://www.cnblogs.com/zhuimengzhe/p/10223510.html

mfcc：

python_speech_features.base.fbank(signal, samplerate=16000, winlen=0.025, winstep=0.01, nfilt=26, nfft=512, lowfreq=0, highfreq=None, preemph=0.97, winfunc=<function >)

• signal - 需要用来计算特征的音频信号，应该是一个N*1的数组
• samplerate - 我们用来工作的信号的采样率
• winlen - 分析窗口的长度，按秒计，默认0.025s(25ms)
• winstep - 连续窗口之间的步长，按秒计，默认0.01s（10ms）
• numcep - 倒频谱返回的数量，默认13
• nfilt - 滤波器组的滤波器数量，默认26
• nfft - FFT的大小，默认512
• lowfreq - 梅尔滤波器的最低边缘，单位赫兹，默认为0
• highfreq - 梅尔滤波器的最高边缘，单位赫兹，默认为采样率/2
• preemph - 应用预加重过滤器和预加重过滤器的系数，0表示没有过滤器，默认0.97
• ceplifter - 将升降器应用于最终的倒谱系数。 0没有升降机。默认值为22。
• appendEnergy - 如果是true，则将第0个倒谱系数替换为总帧能量的对数。
• winfunc - 分析窗口应用于每个框架。 默认情况下不应用任何窗口。 你可以在这里使用numpy窗口函数 例如：winfunc=numpy.hamming
  返回： 一个大小为numcep的numpy数组，包含着特征，每一行都包含一个特征向量。

delta：

python_speech_features.base.delta(feat, N)

• feat - 一个大小为特征数量的numpy数组，每一行都有一个特征向量
• N - 对于每一帧，计算delta特征根据前后N帧
• 返回：一个大小为特征数量的numpy数组，包含有delta特征，每一行都有一个delta向量

②logfbank

def get_fbank(data, fs):
    wav_feature = logfbank(data, fs)
    return wav_feature

参数介绍：

python_speech_features.base.logfbank(signal, samplerate=16000, winlen=0.025, winstep=0.01, nfilt=26, nfft=512, lowfreq=0, highfreq=None, preemph=0.97)

• signal - 需要用来计算特征的音频信号，应该是一个N*1的数组
• samplerate - 我们用来工作的信号的采样率
• winlen - 分析窗口的长度，按秒计，默认0.025s(25ms)
• winstep - 连续窗口之间的步长，按秒计，默认0.01s（10ms）
• nfilt - 滤波器组的滤波器数量，默认26
• nfft - FFT的大小，默认512
• lowfreq - 梅尔滤波器的最低边缘，单位赫兹，默认为0
• highfreq - 梅尔滤波器的最高边缘，单位赫兹，默认为采样率/2
• preemph - 应用预加重过滤器和预加重过滤器的系数，0表示没有过滤器，默认0.97
• 返回： 一个包含特征的大小为nfilt的numpy数组，每一行都有一个特征向量