Spectral Analysis, Seismology, Time Frequency Analysis, Signal Processing, Data Analysis.
fmax = 25. freq_res = (nw/twin)*2 dt = tr.stats.delta t_vals,freq_vals,quad_vals,thomp_vals = spec.spectrogram(noisy_signal,dt,twin = twin, olap = olap, nw = nw, fmax = fmax) max_qi_vals = np.max(quad_vals, axis=1) log_psd_db = 10*np.log10(max_qi_vals) fig, ax = plt.subplots(1, 1, figsize=(10, 6),