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