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)
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)
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)
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() 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]
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() #plt.axis('off') xmin = -0.1 xmax = 1.1 ymin = -0.1 ymax = 1