Python sacp2 Examples

Example #1

File: egplot.py Project: pysmo/aimbat

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

Example #2

File: egalign2.py Project: XStargate/aimbat

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

Example #3

#!/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()

Example #4

File: egalign2.py Project: unmover/aimbat

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

Example #5

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