hezar.utils.audio_utils module

Common audio utils taken from transformers.audio_utils

hezar.utils.audio_utils.amplitude_to_db(spectrogram: ndarray, reference: float = 1.0, min_value: float = 1e-05, db_range: float | None = None) → ndarray

Converts an amplitude spectrogram to the decibel scale. This computes 20 * log10(spectrogram / reference), using basic logarithm properties for numerical stability.

The motivation behind applying the log function on the (mel) spectrogram is that humans do not hear loudness on a linear scale. Generally to double the perceived volume of a sound we need to put 8 times as much energy into it. This means that large variations in energy may not sound all that different if the sound is loud to begin with. This compression operation makes the (mel) spectrogram features match more closely what humans actually hear.

Parameters:

• spectrogram (np.ndarray) – The input amplitude (mel) spectrogram.

• reference (float, optional, defaults to 1.0) – Sets the input spectrogram value that corresponds to 0 dB. For example, use np.max(spectrogram) to set the loudest part to 0 dB. Must be greater than zero.

• min_value (float, optional, defaults to 1e-5) – The spectrogram will be clipped to this minimum value before conversion to decibels, to avoid taking log(0). The default of 1e-5 corresponds to a minimum of -100 dB. Must be greater than zero.

• db_range (float, optional) – Sets the maximum dynamic range in decibels. For example, if db_range = 80, the difference between the peak value and the smallest value will never be more than 80 dB. Must be greater than zero.

Returns:

the spectrogram in decibels

Return type:

np.ndarray

hezar.utils.audio_utils.hertz_to_mel(freq: float | ndarray, mel_scale: str = 'htk') → float | ndarray

Convert frequency from hertz to mels.

Parameters:

• freq (float or np.ndarray) – The frequency, or multiple frequencies, in hertz (Hz).

• mel_scale (str, optional, defaults to “htk”) – The mel frequency scale to use, “htk” or “slaney”.

Returns:

The frequencies on the mel scale.

Return type:

float or np.ndarray

hezar.utils.audio_utils.load_audio_files(paths: str | List[str], sampling_rate: int = 16000)

hezar.utils.audio_utils.mel_filter_bank(num_frequency_bins: int, num_mel_filters: int, min_frequency: float, max_frequency: float, sampling_rate: int, norm: str | None = None, mel_scale: str = 'htk') → ndarray

Creates a frequency bin conversion matrix used to obtain a mel spectrogram. This is called a mel filter bank, and various implementation exist, which differ in the number of filters, the shape of the filters, the way the filters are spaced, the bandwidth of the filters, and the manner in which the spectrum is warped. The goal of these features is to approximate the non-linear human perception of the variation in pitch with respect to the frequency.

Different banks of mel filters were introduced in the literature. The following variations are supported:

• MFCC FB-20: introduced in 1980 by Davis and Mermelstein, it assumes a sampling frequency of 10 kHz and a speech bandwidth of [0, 4600] Hz.

• MFCC FB-24 HTK: from the Cambridge HMM Toolkit (HTK) (1995) uses a filter bank of 24 filters for a speech bandwidth of [0, 8000] Hz. This assumes sampling rate ≥ 16 kHz.

• MFCC FB-40: from the Auditory Toolbox for MATLAB written by Slaney in 1998, assumes a sampling rate of 16 kHz and speech bandwidth of [133, 6854] Hz. This version also includes area normalization.

• HFCC-E FB-29 (Human Factor Cepstral Coefficients) of Skowronski and Harris (2004), assumes a sampling rate of 12.5 kHz and speech bandwidth of [0, 6250] Hz.

This code is adapted from torchaudio and librosa. Note that the default parameters of torchaudio’s melscale_fbanks implement the “htk” filters while librosa uses the “slaney” implementation.

Parameters:

• num_frequency_bins (int) – Number of frequencies used to compute the spectrogram (should be the same as in stft).

• num_mel_filters (int) – Number of mel filters to generate.

• min_frequency (float) – Lowest frequency of interest in Hz.

• max_frequency (float) – Highest frequency of interest in Hz. This should not exceed sampling_rate / 2.

• sampling_rate (int) – Sample rate of the audio waveform.

• norm (str, optional) – If “slaney”, divide the triangular mel weights by the width of the mel band (area normalization).

• mel_scale (str, optional, defaults to “htk”) – The mel frequency scale to use, “htk” or “slaney”.

Returns:

Triangular filter bank matrix. This is a projection matrix to go from a spectrogram to a mel spectrogram.

Return type:

np.ndarray of shape (num_frequency_bins, num_mel_filters)

hezar.utils.audio_utils.mel_to_hertz(mels: float | ndarray, mel_scale: str = 'htk') → float | ndarray

Convert frequency from mels to hertz.

Parameters:

• mels (float or np.ndarray) – The frequency, or multiple frequencies, in mels.

• mel_scale (str, optional, “htk”) – The mel frequency scale to use, “htk” or “slaney”.

Returns:

The frequencies in hertz.

Return type:

float or np.ndarray

hezar.utils.audio_utils.power_to_db(spectrogram: ndarray, reference: float = 1.0, min_value: float = 1e-10, db_range: float | None = None) → ndarray

Converts a power spectrogram to the decibel scale. This computes 10 * log10(spectrogram / reference), using basic logarithm properties for numerical stability.

The motivation behind applying the log function on the (mel) spectrogram is that humans do not hear loudness on a linear scale. Generally to double the perceived volume of a sound we need to put 8 times as much energy into it. This means that large variations in energy may not sound all that different if the sound is loud to begin with. This compression operation makes the (mel) spectrogram features match more closely what humans actually hear.

Based on the implementation of librosa.power_to_db.

Parameters:

• spectrogram (np.ndarray) – The input power (mel) spectrogram. Note that a power spectrogram has the amplitudes squared!

• reference (float, optional, defaults to 1.0) – Sets the input spectrogram value that corresponds to 0 dB. For example, use np.max(spectrogram) to set the loudest part to 0 dB. Must be greater than zero.

• min_value (float, optional, defaults to 1e-10) – The spectrogram will be clipped to this minimum value before conversion to decibels, to avoid taking log(0). The default of 1e-10 corresponds to a minimum of -100 dB. Must be greater than zero.

• db_range (float, optional) – Sets the maximum dynamic range in decibels. For example, if db_range = 80, the difference between the peak value and the smallest value will never be more than 80 dB. Must be greater than zero.

Returns:

the spectrogram in decibels

Return type:

np.ndarray

hezar.utils.audio_utils.spectrogram(waveform: ~numpy.ndarray, window: ~numpy.ndarray, frame_length: int, hop_length: int, fft_length: int | None = None, power: float | None = 1.0, center: bool = True, pad_mode: str = 'reflect', onesided: bool = True, preemphasis: float | None = None, mel_filters: ~numpy.ndarray | None = None, mel_floor: float = 1e-10, log_mel: str | None = None, reference: float = 1.0, min_value: float = 1e-10, db_range: float | None = None, dtype: ~numpy.dtype = <class 'numpy.float32'>) → ndarray

Calculates a spectrogram over one waveform using the Short-Time Fourier Transform.

This function can create the following kinds of spectrograms:

• amplitude spectrogram (power = 1.0)

• power spectrogram (power = 2.0)

• complex-valued spectrogram (power = None)

• log spectrogram (use log_mel argument)

• mel spectrogram (provide mel_filters)

• log-mel spectrogram (provide mel_filters and log_mel)

How this works:

1. The input waveform is split into frames of size frame_length that are partially overlapping by frame_length - hop_length samples.

2. Each frame is multiplied by the window and placed into a buffer of size fft_length.

3. The DFT is taken of each windowed frame.

4. The results are stacked into a spectrogram.

We make a distinction between the following “blocks” of sample data, each of which may have a different lengths:

• The analysis frame. This is the size of the time slices that the input waveform is split into.

• The window. Each analysis frame is multiplied by the window to avoid spectral leakage.

• The FFT input buffer. The length of this determines how many frequency bins are in the spectrogram.

In this implementation, the window is assumed to be zero-padded to have the same size as the analysis frame. A padded window can be obtained from window_function(). The FFT input buffer may be larger than the analysis frame, typically the next power of two.

Note: This function is not optimized for speed yet. It should be mostly compatible with librosa.stft and torchaudio.functional.transforms.Spectrogram, although it is more flexible due to the different ways spectrograms can be constructed.

Parameters:

• waveform (np.ndarray of shape (length,)) – The input waveform. This must be a single real-valued, mono waveform.

• window (np.ndarray of shape (frame_length,)) – The windowing function to apply, including zero-padding if necessary. The actual window length may be shorter than frame_length, but we’re assuming the array has already been zero-padded.

• frame_length (int) – The length of the analysis frames in samples. With librosa this is always equal to fft_length but we also allow smaller sizes.

• hop_length (int) – The stride between successive analysis frames in samples.

• fft_length (int, optional) – The size of the FFT buffer in samples. This determines how many frequency bins the spectrogram will have. For optimal speed, this should be a power of two. If None, uses frame_length.

• power (float, optional, defaults to 1.0) – If 1.0, returns the amplitude spectrogram. If 2.0, returns the power spectrogram. If None, returns complex numbers.

• center (bool, optional, defaults to True) – Whether to pad the waveform so that frame t is centered around time t * hop_length. If False, frame t will start at time t * hop_length.

• pad_mode (str, optional, defaults to “reflect”) – Padding mode used when center is True. Possible values are: “constant” (pad with zeros), “edge” (pad with edge values), “reflect” (pads with mirrored values).

• onesided (bool, optional, defaults to True) – If True, only computes the positive frequencies and returns a spectrogram containing fft_length // 2 + 1 frequency bins. If False, also computes the negative frequencies and returns fft_length frequency bins.

• preemphasis (float, optional) – Coefficient for a low-pass filter that applies pre-emphasis before the DFT.

• mel_filters (np.ndarray of shape (num_freq_bins, num_mel_filters), optional) – The mel filter bank. If supplied, applies a this filter bank to create a mel spectrogram.

• mel_floor (float, optional, defaults to 1e-10) – Minimum value of mel frequency banks.

• log_mel (str, optional) – How to convert the spectrogram to log scale. Possible options are: None (don’t convert), “log” (take the natural logarithm) “log10” (take the base-10 logarithm), “dB” (convert to decibels). Can only be used when power is not None.

• reference (float, optional, defaults to 1.0) – Sets the input spectrogram value that corresponds to 0 dB. For example, use np.max(spectrogram) to set the loudest part to 0 dB. Must be greater than zero.

• min_value (float, optional, defaults to 1e-10) – The spectrogram will be clipped to this minimum value before conversion to decibels, to avoid taking log(0). For a power spectrogram, the default of 1e-10 corresponds to a minimum of -100 dB. For an amplitude spectrogram, the value 1e-5 corresponds to -100 dB. Must be greater than zero.

• db_range (float, optional) – Sets the maximum dynamic range in decibels. For example, if db_range = 80, the difference between the peak value and the smallest value will never be more than 80 dB. Must be greater than zero.

• dtype (np.dtype, optional, defaults to np.float32) – Data type of the spectrogram tensor. If power is None, this argument is ignored and the dtype will be np.complex64.

Returns:

nd.array containing a spectrogram of shape (num_frequency_bins, length) for a regular spectrogram or shape (num_mel_filters, length) for a mel spectrogram.

hezar.utils.audio_utils.window_function(window_length: int, name: str = 'hann', periodic: bool = True, frame_length: int | None = None, center: bool = True) → ndarray

Returns an array containing the specified window. This window is intended to be used with stft.

The following window types are supported:

• “boxcar”: a rectangular window

• “hamming”: the Hamming window

• “hann”: the Hann window

Parameters:

• window_length (int) – The length of the window in samples.

• name (str, optional, defaults to “hann”) – The name of the window function.

• periodic (bool, optional, defaults to True) – Whether the window is periodic or symmetric.

• frame_length (int, optional) – The length of the analysis frames in samples. Provide a value for frame_length if the window is smaller than the frame length, so that it will be zero-padded.

• center (bool, optional, defaults to True) – Whether to center the window inside the FFT buffer. Only used when frame_length is provided.

Returns:

np.ndarray of shape (window_length,) or (frame_length,) containing the window.