Exemplos de terntransform em Python

Exemplo n.º 1

Arquivo: histplots.py Projeto: bpostlethwaite/Masters-Thesis

def plottern(plt, alpha, vp):
    gran = np.array([6.208, poisson(0.25, True)])
    mafc = np.array([6.942, poisson(0.283, True)])
    gray = np.array([6.302, poisson(0.257, True)])
    vp = np.array([vp, alpha])

    a = gran
    b = mafc
    c = gray
    d = vp

    endmembers = ("Granite Gneiss", "Mafic Granulite", "Gray Gneiss")
    datalabels = ["Vp", r"$\sigma$"]
    colors = ["b", "g", "c", "m", "r"]
    lines = terntransform(a, b, c, d)

    step = 0.2

    for line, color, label in zip(lines, colors, datalabels):
        tern.plot(line, color=color, linewidth=2.0, label=label)

    ip = baryIntersect(lines)
    tern.plotIntersect(ip)

    tern.gridlines(step, "0.7")
    tern.draw_boundary()
    tern.addlabels(endmembers)
    plt.legend()

    # pyplot.box(on='off')
    h = plt.gca()
    # plt.axis('off')
    xmin = -0.1
    xmax = 1.1
    ymin = -0.1
    ymax = 1
    v = [xmin, xmax, ymin, ymax]
    plt.axis(v)
    h.axes.get_xaxis().set_visible(False)
    h.axes.get_yaxis().set_visible(False)

Exemplo n.º 2

Arquivo: histplots.py Projeto: padhys/Masters-Thesis

def plottern(plt, alpha, vp):
    gran = np.array([6.208, 0.25])
    mafc = np.array([6.942, 0.283])
    gray = np.array([6.302, 0.257])
    vp = np.array([vp, alpha])

    a = gran
    b = mafc
    c = gray
    d = vp

    endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
    datalabels = ['Vp',r'$\sigma$']
    colors = ['b','g','c','m','r']
    lines = terntransform(a, b, c, d)

    step = 0.2

    for line,color,label in zip(lines,colors,datalabels):
        tern.plot(line, color = color, linewidth=2.0, label = label)

    ip = baryIntersect(lines)
    tern.plotIntersect(ip)

    tern.gridlines(step, '0.7')
    tern.draw_boundary()
    tern.addlabels(endmembers)
    plt.legend()

        #pyplot.box(on='off')
    h = plt.gca()
    #plt.axis('off')
    xmin = -0.1
    xmax = 1.1
    ymin = -0.1
    ymax = 1
    v = [xmin, xmax, ymin, ymax]
    plt.axis(v)
    h.axes.get_xaxis().set_visible(False)
    h.axes.get_yaxis().set_visible(False)

Exemplo n.º 3

def plottern(plt, alpha, vp):
    gran = np.array([6.208, poisson(0.25, True)])
    mafc = np.array([6.942, poisson(0.283, True)])
    gray = np.array([6.302, poisson(0.257, True)])
    vp = np.array([vp, alpha])

    a = gran
    b = mafc
    c = gray
    d = vp

    endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
    datalabels = ['Vp',r'$\sigma$']
    colors = ['b','g','c','m','r']
    lines = terntransform(a, b, c, d)

    step = 0.2

    for line,color,label in zip(lines,colors,datalabels):
        tern.plot(line, color = color, linewidth=2.0, label = label)

    ip = baryIntersect(lines)
    tern.plotIntersect(ip)

    tern.gridlines(step, '0.7')
    tern.draw_boundary()
    tern.addlabels(endmembers)
    plt.legend()

        #pyplot.box(on='off')
    h = plt.gca()
    #plt.axis('off')
    xmin = -0.1
    xmax = 1.1
    ymin = -0.1
    ymax = 1
    v = [xmin, xmax, ymin, ymax]
    plt.axis(v)
    h.axes.get_xaxis().set_visible(False)
    h.axes.get_yaxis().set_visible(False)

Exemplo n.º 4

Arquivo: terntest.py Projeto: padhys/Masters-Thesis

gran = np.array([6.208, 0.25])
mafc = np.array([6.942, 0.283])
gray = np.array([6.302, 0.257])
vp = np.array([6.316, 0.2585])

a = gran
b = mafc
c = gray
d = vp

endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
datalabels = ['Vp',r'$\sigma$']
colors = ['b','g','c','m','r']
lines = terntransform(a, b, c, d)

step = 0.2



plt.figure()

for line,color,label in zip(lines,colors,datalabels):
    tern.plot(line, color = color, linewidth=2.0, label = label)

ip = baryIntersect(lines)
tern.plotIntersect(ip)

tern.gridlines(step, '0.7')
tern.draw_boundary()

Exemplo n.º 5

Arquivo: terntest.py Projeto: bpostlethwaite/Masters-Thesis

# Select Assemblage
data = lithdata[3:6]
endmembers = endmembers[3:6]

# Data

plt.figure()
for ind, region in enumerate(regions):

    data = lithdata[3:6]
    data.append(region[1])

    datalabels = ['Vp',r'$\sigma$']
    colors = ['b','g','c','m','r','y']

    lines = terntransform( *data )

    for line, color, label in zip(lines, colors, datalabels):
        tern.plot(line, color = color, linewidth=2.0, label = label if (ind == 0) else None)

    ip = baryIntersect(lines)
    tern.plotIntersect(ip, marker = 'o'+colors[ind], label = region[0])

step = 0.2
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()

    #pyplot.box(on='off')
h = plt.gca()

Exemplo n.º 6

Arquivo: terntest.py Projeto: bpostlethwaite/Masters-Thesis

# Select Assemblage
data = lithdata[3:6]
endmembers = endmembers[3:6]

# Data

plt.figure()
for ind, region in enumerate(regions):

    data = lithdata[3:6]
    data.append(region[1])

    datalabels = ['Vp', r'$\sigma$']
    colors = ['b', 'g', 'c', 'm', 'r', 'y']

    lines = terntransform(*data)

    for line, color, label in zip(lines, colors, datalabels):
        tern.plot(line,
                  color=color,
                  linewidth=2.0,
                  label=label if (ind == 0) else None)

    ip = baryIntersect(lines)
    tern.plotIntersect(ip, marker='o' + colors[ind], label=region[0])

step = 0.2
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()