PitchInFrames

# https://stackoverflow.com/questions/43877971/librosa-pitch-tracking-stft
import librosa
import peakutils
import numpy as np
import matplotlib.pyplot as plt

def detect_pitch(y, sr, t):
  index = magnitudes[:, t].argmax()
  pitch = pitches[index, t]

  return pitch



filename = '/home/nexttest_20200729.wav'
y, sr = librosa.load(filename)

y = y[:]


n_fft = 2048

#pitches, magnitudes = librosa.core.piptrack(y=y, sr=sr, fmin=75, fmax=1600)
pitches, magnitudes = librosa.core.piptrack(y=y, sr=sr, n_fft=n_fft)



pitchy = detect_pitch(y,sr,2)
print("Der Pitch:", pitchy)

print(len(magnitudes), magnitudes.shape, len(pitches), pitches.shape, len(y), y.shape)

for i in range(0,pitches.shape[1]):
    print(i,detect_pitch(y,sr,i))

for i in range(0,pitches.shape[0]):
    print("Nummer:", i, " Magnitude (2):", magnitudes[i,2], " Pitch (2):", pitches[i,2])



indexes = peakutils.indexes(y, thres=0.05/max(y), min_dist=10240)
x = np.linspace(0, len(y)-1, len(y))
print("Shape x:", x.shape)
print("Shape y:", y.shape)
peaks_x = peakutils.interpolate(x, y, ind=indexes)
peaks_x_norm = []
for peak in peaks_x:
    print("Peak:", peak)
    t = round(peak/n_fft)
    print("T:", t)
    print("Pitch:", detect_pitch(y, sr, t))
    peaks_x_norm.append(round(peak))

plt.plot(y)
plt.plot(peaks_x_norm, y[peaks_x_norm], "x")
plt.plot(indexes, y[indexes], "o")
plt.plot(np.zeros_like(y), "--", color="gray")
plt.show()

Output

Peak: 154.17451740228142
T: 0
Pitch: 439.82
Peak: 10973.00106524168
T: 5
Pitch: 512.3134
Peak: 22527.197688260196
T: 11
Pitch: 1430.7959
Peak: 33025.386784057075
T: 16
Pitch: 1579.7072
Peak: 43491.277796775175
T: 21
Pitch: 1424.7646
Peak: 54783.898177332376
T: 27
Pitch: 1471.8829
Peak: 65856.94867421205
T: 32
Pitch: 1568.4641
Peak: 77054.0199016538
T: 38
Pitch: 1467.3224
Peak: 87664.68073674176
T: 43
Pitch: 1562.4485
Peak: 98426.02078009519
T: 48
Pitch: 1562.9667
Peak: 109674.60057979365
T: 54
Pitch: 1661.3
Peak: 120802.04548626121
T: 59
Pitch: 2495.5845
Peak: 131827.88590985513
T: 64
Pitch: 776.8133
Peak: 142764.4072741114
T: 70
Pitch: 822.5438
Peak: 154569.187918604
T: 75
Pitch: 1486.1099
Peak: 164973.65818289487
T: 81
Pitch: 1506.1406
Peak: 175918.7137020767
T: 86
Pitch: 496.76382
Peak: 186963.67849766312
T: 91
Pitch: 1658.482
Peak: 199027.74761241846
T: 97
Pitch: 1669.4775
Peak: 209480.82419132258
T: 102
Pitch: 1478.1989
Peak: 220598.00529181876
T: 108
Pitch: 691.4572
Peak: 236615.01312228944
T: 116
Pitch: 1475.9712
Peak: 252372.90489012047
T: 123
Pitch: 467.1594
Peak: 263993.4979628992
T: 129
Pitch: 461.03122
Peak: 275336.3161003483
T: 134
Pitch: 1476.0253
Peak: 286737.5135514725
T: 140
Pitch: 1663.2999
Peak: 297230.0513971054
T: 145
Pitch: 832.98004
Peak: 308239.6564216623
T: 151
Pitch: 498.58484
Peak: 319255.09065255866
T: 156
Pitch: 1634.6113
Peak: 330687.2112069934
T: 161
Pitch: 1658.923
Peak: 341351.4292968454
T: 167
Pitch: 1662.9713
Peak: 351593.99030480173
T: 172
Pitch: 499.98648
Peak: 363856.46433673584
T: 178
Pitch: 1655.4417
Peak: 374647.0
T: 183
Pitch: 622.9731
Peak: 385874.4950590556
T: 188
Pitch: 829.49506
Peak: 396394.337723163
T: 194
Pitch: 767.05365
Peak: 407165.31716953753
T: 199
Pitch: 1648.6659
Peak: 418683.6821068189
T: 204
Pitch: 1662.068
Peak: 429294.7974067266
T: 210
Pitch: 1656.3007
Peak: 439820.8522538731
T: 215
Pitch: 351.36328
Peak: 451916.9050249872
T: 221
Pitch: 1557.4805
Peak: 463082.42422536603
T: 226
Pitch: 830.7519
Peak: 473514.33764882694
T: 231
Pitch: 2497.5059
Peak: 483873.77526077995
T: 236
Pitch: 772.7677
Peak: 500969.79485515563
T: 245
Pitch: 1657.3533
Peak: 515182.1932971365
T: 252
Pitch: 1247.202
Peak: 528960.5438554683
T: 258
Pitch: 350.50412
Peak: 540369.991649482
T: 264
Pitch: 2487.0857

Dieses Blog durchsuchen

Es gibt ...

PitchInFrames

Code

Output

Kommentare

Kommentar veröffentlichen

Beliebte Posts aus diesem Blog

·

Es brennt.

Bye, bye Nord Stream 2!