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()