0

I am using SpectrumLab for logging meteors with an SDR. SpectrumLab records waterfall screenshots and a wav file of an event. I am trying to reproduce the waterfall screenshot of SpectrumLab from the wave file but the pattern looks different:

import librosa
import matplotlib.pyplot as plt
import librosa.display
audio_data = 'event20230131_101027_26.wav'
x , sr = librosa.load(audio_data)
plt.figure(figsize=(14, 5))
librosa.display.waveshow(x, sr=sr)

X = librosa.stft(x)
Xdb = librosa.amplitude_to_db(abs(X))
plt.figure(figsize=(14, 5))
librosa.display.specshow(Xdb, sr=sr, x_axis='time', y_axis='log')
plt.colorbar()
plt.show()

Screenshot of SpectrumLab: enter image description here

Screenshot generated with librosa: enter image description here

Wave file which was recorded together with the screenshot in SpectrumLab:

Wavefile

honeymoon
  • 2,400
  • 5
  • 34
  • 43
  • 1
    What exactly does the SpectrumLab screenshot show? I assume the y-axis is radio frequency, which is far different from the audio frequencies that you have in the wav file. It's not obvious how the SpectrumLab input is converted to a wav file, so it's worth looking into how that is performed before expecting to see the same output. – Matt Pitkin Feb 03 '23 at 10:40
  • Yes, y-axis is radio frequency. – honeymoon Feb 03 '23 at 11:38

0 Answers0