def psd_plot(station): ppsd = PPSD.load("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h.pkl.bz2" % station) ppsd.plot(show=False, show_coverage=False, show_percentiles=True, period_lim=(0.05, 1000), max_percentage=10) fig = mpl.pyplot.figure(1) ax = fig.axes[0] ax.set_ylim((-200, -50)) ax.set_xticklabels(['%s' % t for t in ax.get_xticks()]) ax.set_xlabel('period (s)') ax.set_ylabel('amplitude (dB)') col = ax.collections[0] bbox = col.colorbar.ax.get_position() fig.delaxes(col.colorbar.ax) if COLORBAR: fig.colorbar(col, cax=fig.add_axes( [bbox.x0 + 0.1, bbox.y0, bbox.width, bbox.height]), extend='max') #col.colorbar.set_ticks(range(8)) col.colorbar.set_label('probability (%)') bbox = ax.get_position() ax.set_position( [bbox.x0, bbox.y0, bbox.width + 0.2 - 0.1 * COLORBAR, bbox.height]) fig.savefig("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h_cb.png" % station, dpi=150, bbox_inches='tight') mpl.pyplot.close()
def correct_sensitivity(station): ppsd = PPSD.load("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h_wrong_gain.pkl" % station) dB = 10 * np.log10(2516580000.0 / 629145000.0) if station == 'PB14': dB = 10 * np.log10(588000000.0 / 629145000.0) ppsd.spec_bins = ppsd.spec_bins - dB ppsd.yedges = ppsd.yedges - dB ppsd.save("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h.pkl" % station)
def correct_sensitivity(station): ppsd = PPSD.load( "/home/richter/Results/IPOC/PPSD/ppsd_%s_6h_wrong_gain.pkl" % station) dB = 10 * np.log10(2516580000.0 / 629145000.0) if station == 'PB14': dB = 10 * np.log10(588000000.0 / 629145000.0) ppsd.spec_bins = ppsd.spec_bins - dB ppsd.yedges = ppsd.yedges - dB ppsd.save("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h.pkl" % station)
def psd_plot(station): ppsd = PPSD.load("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h.pkl.bz2" % station) ppsd.plot(show=False, show_coverage=False, show_percentiles=True, period_lim=(0.05, 1000), max_percentage=10) fig = mpl.pyplot.figure(1) ax = fig.axes[0] ax.set_ylim((-200, -50)) ax.set_xticklabels(['%s' % t for t in ax.get_xticks()]) ax.set_xlabel('period (s)') ax.set_ylabel('amplitude (dB)') col = ax.collections[0] bbox = col.colorbar.ax.get_position() fig.delaxes(col.colorbar.ax) if COLORBAR: fig.colorbar(col, cax=fig.add_axes([bbox.x0 + 0.1, bbox.y0, bbox.width, bbox.height]), extend='max') #col.colorbar.set_ticks(range(8)) col.colorbar.set_label('probability (%)') bbox = ax.get_position() ax.set_position([bbox.x0, bbox.y0, bbox.width + 0.2 - 0.1 * COLORBAR, bbox.height]) fig.savefig("/home/richter/Results/IPOC/PPSD/ppsd_%s_6h_cb.png" % station, dpi=150, bbox_inches='tight') mpl.pyplot.close()