pyACA.FeatureSpectralTonalPowerRatio Namespace Reference

Functions
	FeatureSpectralTonalPowerRatio (X, f_s, G_T=5e-4)
	computes the tonal power ratio from the magnitude spectrum

Function Documentation

FeatureSpectralTonalPowerRatio()

FeatureSpectralTonalPowerRatio	(		X,
			f_s,
			G_T = 5e-4 )

computes the tonal power ratio from the magnitude spectrum

Parameters

X	spectrogram (dimension FFTLength X Observations)
f_s	sample rate of audio data
G_T	energy threshold

Returns

vtpr: tonal power ratio

Definition at line 14 of file FeatureSpectralTonalPowerRatio.py.

def FeatureSpectralTonalPowerRatio(X, f_s, G_T=5e-4):
 
    isSpectrum = X.ndim == 1
    if isSpectrum:
        X = np.expand_dims(X, axis=1)
 
    X = X**2
 
    fSum = X.sum(axis=0)
    vtpr = np.zeros(fSum.shape)
 
    for n in range(0, X.shape[1]):
        if fSum[n] < G_T:
            continue
 
        # find local maxima above the threshold
        afPeaks = find_peaks(X[:, n], height=G_T)
 
        if not afPeaks[0].size:
            continue
 
        # calculate ratio
        vtpr[n] = X[afPeaks[0], n].sum() / fSum[n]
 
    return np.squeeze(vtpr) if isSpectrum else vtpr