axp2 = fig.add_axes(rectp2) axp1 = fig.add_axes(rectp1) axp0 = fig.add_axes(rectp0) # read data and plot gsac, opts = getDataOpts() # prs opts.ynorm = .95 saclist = gsac.saclist prs = sacprs(saclist, opts, axp0) # p1 opts.ynorm = 1.7 p1 = sacp1(saclist, opts, axp1) # p2 opts.reltime = 0 p2 = sacp2(saclist, opts, axp2) # set x limits axp0.set_xlim(625, 762) axp1.set_xlim(625, 762) axp2.set_xlim(-45, 65) # numbering axs = [axp0, axp1, axp2] labs = 'abc' for ax, lab in zip(axs, labs): tt = '(' + lab + ')' trans = transforms.blended_transform_factory(ax.transAxes, ax.transAxes) ax.text(-.05, 1, tt, transform=trans, va='center', ha='right', size=16) fig.savefig('egplot.pdf', format='pdf') plt.show()
if __name__ == "__main__": gsac, opts = load() saclist = gsac.saclist xxlim = -20, 20 reltimes = [0, 1, 2, 3] tlabs = "abcd" npick = len(reltimes) axs = axes2(npick) for i in range(npick): opts.reltime = reltimes[i] tpick = "t" + str(opts.reltime) ax = axs[i] sacp2(saclist, opts, ax) # ax.text(0.03,0.93, 'Aligned by '+ tpick.upper() , fontsize=12, # horizontalalignment='left', # verticalalignment='top', # transform=ax.transAxes) ax.set_xlim(xxlim) tt = "(" + tlabs[i] + ")" trans = transforms.blended_transform_factory(ax.transAxes, ax.transAxes) ax.text(-0.05, 1, tt, transform=trans, va="center", ha="right", size=16) savefig("egalignp2.pdf", format="pdf") show()
#!/usr/bin/env python #------------------------------------------------ # Filename: sacp2.py # Author: Xiaoting Lou # Email: [email protected] # # Copyright (c) 2009-2012 Xiaoting Lou #------------------------------------------------ """ Python script fo SAC p2 style of plotting. :copyright: Xiaoting Lou :license: GNU General Public License, Version 3 (GPLv3) http://www.gnu.org/licenses/gpl.html """ from matplotlib.pyplot import show from pysmo.aimbat.plotphase import getDataOpts, getAxes, sacp2 gsac, opts = getDataOpts() axss = getAxes(opts) ssg = sacp2(gsac.saclist, opts, axss) show()
gsac.stkdh = SacDataHdrs(opts.fstack, opts.delta) opts.pppara = pppara opts.ccpara = ccpara return gsac, opts if __name__ == '__main__': gsac, opts = load() saclist = gsac.saclist xxlim = -20, 20 reltimes = [0, 1, 2, 3] tlabs = 'abcd' npick = len(reltimes) axs = axes2(npick) for i in range(npick): opts.reltime = reltimes[i] tpick = 't' + str(opts.reltime) ax = axs[i] sacp2(saclist, opts, ax) # ax.text(0.03,0.93, 'Aligned by '+ tpick.upper() , fontsize=12, # horizontalalignment='left', # verticalalignment='top', # transform=ax.transAxes) ax.set_xlim(xxlim) tt = '(' + tlabs[i] + ')' trans = transforms.blended_transform_factory(ax.transAxes, ax.transAxes) ax.text(-.05, 1, tt, transform=trans, va='center', ha='right', size=16) plt.savefig('egalignp2.pdf', format='pdf') plt.show()
rectp0 = [.09, .620, .8, .36] axp2 = fig.add_axes(rectp2) axp1 = fig.add_axes(rectp1) axp0 = fig.add_axes(rectp0) # read data and plot gsac, opts = getDataOpts() # prs opts.ynorm = .95 saclist = gsac.saclist prs = sacprs(saclist, opts, axp0) # p1 opts.ynorm = 1.7 p1 = sacp1(saclist, opts, axp1) # p2 opts.reltime = 0 p2 = sacp2(saclist, opts, axp2) # set x limits axp0.set_xlim(625, 762) axp1.set_xlim(625, 762) axp2.set_xlim(-45, 65) # numbering axs = [axp0, axp1, axp2] labs = 'abc' for ax, lab in zip(axs, labs): tt = '(' + lab + ')' trans = transforms.blended_transform_factory(ax.transAxes, ax.transAxes) ax.text(-.05, 1, tt, transform=trans, va='center', ha='right', size=16) fig.savefig('egplot.pdf', format='pdf') plt.show()