code: def stChromaFeaturesInit(nfft, fs):
      
      freqs = numpy.array([((f + 1) * fs) / (2 * nfft) for f in range(nfft)])    
      Cp = 27.50    
      nChroma = numpy.round(12.0 * numpy.log2(freqs / Cp)).astype(int)
  
      nFreqsPerChroma = numpy.zeros((nChroma.shape[0], ))
  
      uChroma = numpy.unique(nChroma)
      for u in uChroma:
          idx = numpy.nonzero(nChroma == u)
          nFreqsPerChroma[idx] = idx[0].shape
      
      return nChroma, nFreqsPerChroma
docs: This function initializes the chroma matrices used in the calculation of the chroma features