Python librosa implementation of https://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=8865624

cmc.py

import librosa
import scipy
import numpy as np


def cmc(y, sr, P=12, n_level=7, eps=1e-12, **kwargs):
    log_abs_power_cqt = np.log(
        np.abs(librosa.core.cqt(y=y, sr=sr, **kwargs))**2 + eps)
    cmc_levels = []
    for _ in range(n_level):
        octave = scipy.fftpack.dct(log_abs_power_cqt,
                                   norm='ortho',
                                   axis=0,
                                   type=2)
        cmc = octave[:P]
        log_abs_power_cqt = octave[P:]  # residual to next dct
        cmc_levels.append(cmc)
    return np.concatenate(cmc_levels, 0).T