S, sig95, gs, signif_g, time_base, scalemin=.1, scalemax=128, ylabel='Echo Intens.', plot_percentile=True) plt.savefig(fig_name_base + '_wave3' + str(depth[level]).replace('.0', 'm') + '.png', bbox_inches='tight', dpi=(100)) plt.close() """----------------------- plotting power spectrum FFT --------------------------------""" (plt, fig) = wavelet_analy_plot.fft_power_spec(x, time_base, Fs=24) plt.savefig(fig_name_base + '_FFTspec' + str(depth[level]).replace('.0', 'm') + '.png', bbox_inches='tight', dpi=(100)) plt.close() """ # Do FFT analysis of array sp = np.fft.fft(x) # Getting the related frequencies freq = np.fft.fftfreq(t.shape[-1], d=.25) pyy = sp*np.conj(sp) """ """----------------- plot scale averaged power timeseries -----------------------------""" #Need to know what scales to plot
plt, fig = wavelet_analy_plot.plot_wavetransf_time(x, wa, T, S, sig95, gs, signif_g, time_base, ylabel='Echo Intens.', plot_percentile=True) plt.savefig(fig_name_base + '_wave2' + str(depth[level]).replace('.0','m') + '.png', bbox_inches='tight', dpi = (100)) plt.close() """----------------- plotting contours w/global and timeseries ----------""" """----------------- zoom in to specified scales ----------""" plt, fig = wavelet_analy_plot.plot_wavetransf_time_zoom(x, wa, T, S, sig95, gs, signif_g, time_base, scalemin=.1, scalemax=128, ylabel='Echo Intens.', plot_percentile=True) plt.savefig(fig_name_base + '_wave3' + str(depth[level]).replace('.0','m') + '.png', bbox_inches='tight', dpi = (100)) plt.close() """----------------------- plotting power spectrum FFT --------------------------------""" (plt, fig) = wavelet_analy_plot.fft_power_spec(x, time_base, Fs=24) plt.savefig(fig_name_base + '_FFTspec' + str(depth[level]).replace('.0','m') + '.png', bbox_inches='tight', dpi = (100)) plt.close() """ # Do FFT analysis of array sp = np.fft.fft(x) # Getting the related frequencies freq = np.fft.fftfreq(t.shape[-1], d=.25) pyy = sp*np.conj(sp) """ """----------------- plot scale averaged power timeseries -----------------------------""" #Need to know what scales to plot scales_bin = [2, 8]
"""----------------------------- plot setup ------------------------------------------""" T, S = np.meshgrid(t, scales) """----------- plotting WaveTransform Power with confidence interval contour ----------""" """----------------- plotting contours w/global and timeseries ----------""" """----------------- zoom in to specified scales ----------""" plt, fig = wavelet_analy_plot.plot_wavetransf_time_zoom(data_1['anom'], wa, T, S, sig95, gs, signif_g, data_1['time_base'], scalemin=.1, scalemax=64, ylabel=par_1[1], plot_percentile=True) plt.savefig(fig_name_base + '_wave3' + depth + '.png', bbox_inches='tight', dpi = (100)) plt.close() """----------------------- plotting power spectrum FFT --------------------------------""" (plt, fig) = wavelet_analy_plot.fft_power_spec(data_1['anom'], data_1['time_base'], Fs=24) plt.savefig(fig_name_base + '_FFTspec' + depth + '.png', bbox_inches='tight', dpi = (100)) plt.close()