Mel Powers

mel_powers.py

from essentia.standard import MelBands

# Compute the mel band powers
melbands = MelBands(lowFrequencyBound=f_low,
                    highFrequencyBound=f_high,
                    inputSize=1024,
                    numberBands=n_bands,
                    type='magnitude',  # We already computed the power.
                    sampleRate=44100)
mels = melbands(spec)

# Plot the mel band powers
plt.bar(np.arange(len(mels)), mels, align='center')
plt.title('Mel Band Powers')
plt.xlabel('Mel Bands')
plt.ylabel('Power')