def main():
args = parseargs()
if args.l:
# Load layout
font = load_attribute(args.l + '.font')
matplotlib.rc('font', **font)
graph = None
if args.show or args.s:
fig = plt.figure()
ax_main, ax_res = Grid(fig, rect=111, nrows_ncols=(2, 1),
axes_pad=0.2, label_mode='L')
# ax_main.get_xaxis().set_visible(False)
# div = make_axes_locatable(ax_main)
# ax_res = div.append_axes('bottom', pad=0.05, label='Residuals',
# size='50%')
# ax_res = fig.add_subplot(212, label='Residuals')
graph = (fig, (ax_main, ax_res))
f_class = args.f()
control = (None, None, None)
colors = iter(cm.brg(np.linspace(0, 1, len(args.i) + bool(args.c))))
if args.c:
_, dat = unpack_data(args.c, args.delimiter).__next__()
data = [('control', dat)]
control = analysis_update(data, f_class, colors.__next__(), graph,
report=False, control=control, notify=False)
for path in args.i:
color = colors.__next__()
try:
data = unpack_data(path, args.delimiter, filtr=args.filtr,
split_column=args.split)
analysis_update(data, f_class, color, graph, args.r,
control=control, notify=True)
except Exception as e:
raise e
print ('Could not process', path, os.linesep)
return None
if args.s:
plt.tight_layout()
bname = os.path.basename(path)
name = os.path.splitext(bname)[0]
fig.savefig(os.path.join(args.s, name + '.png'), dpi=160,
bbox_inches='tight', pad_inches=0)
if args.show:
plt.tight_layout()
plt.show()
#REDVOL = ['90']
#REDVOL = ['85']
#REDVOL = ['80']
#REDVOL = ['75']
#REDVOL = ['70']
#REDVOL = ['65']
#fs = 24
#plt.figure(figsize=(14,12))
fs = 10
plt.figure(figsize=(6,4))
plt.rc('text',usetex=True)
plt.rc('font', family='serif')
#for redvol in REDVOL :
for redvol in REDVOL :
color = iter(cm.brg(np.linspace(0.0,0.8, len(CONFIN))))
for confin in CONFIN :
c = next(color)
#(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA) = readMetrics(confin)
(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA, CURV) = readMetrics(redvol,confin)
v = round(REDVOL[0],2)
lam = LAMBDA[0]
# calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0
lamc = 1.0
if (abs(v - 1.00) < 1e-12) :
lamc = 1.000000000000000000000000
if (abs(v - 0.99) < 1e-12) :
lamc = 1.090275107548953613318768
if (abs(v - 0.98) < 1e-12) :
def main():
args = parseargs()
if args.l:
# Load layout
font = load_attribute(args.l + '.font')
matplotlib.rc('font', **font)
graph = None
if args.show or args.s:
fig = plt.figure()
ax_main, ax_res = Grid(fig,
rect=111,
nrows_ncols=(2, 1),
axes_pad=0.2,
label_mode='L')
# ax_main.get_xaxis().set_visible(False)
# div = make_axes_locatable(ax_main)
# ax_res = div.append_axes('bottom', pad=0.05, label='Residuals',
# size='50%')
# ax_res = fig.add_subplot(212, label='Residuals')
graph = (fig, (ax_main, ax_res))
f_class = args.f()
control = (None, None, None)
colors = iter(cm.brg(np.linspace(0, 1, len(args.i) + bool(args.c))))
if args.c:
_, dat = unpack_data(args.c, args.delimiter).__next__()
data = [('control', dat)]
control = analysis_update(data,
f_class,
colors.__next__(),
graph,
report=False,
control=control,
notify=False)
for path in args.i:
color = colors.__next__()
try:
data = unpack_data(path,
args.delimiter,
filtr=args.filtr,
split_column=args.split)
analysis_update(data,
f_class,
color,
graph,
args.r,
control=control,
notify=True)
except Exception as e:
raise e
print('Could not process', path, os.linesep)
return None
if args.s:
plt.tight_layout()
bname = os.path.basename(path)
name = os.path.splitext(bname)[0]
fig.savefig(os.path.join(args.s, name + '.png'),
dpi=160,
bbox_inches='tight',
pad_inches=0)
if args.show:
plt.tight_layout()
plt.show()
x_real_validate = np.zeros(len(A1[0]))
x_real_validate[get_index(1560+float(signals[i])/2.0, wavelengths)] = 0.8
x_real_validate[get_index(1560-float(signals[i])/2.0, wavelengths)] = 1.0
x_real_vals.append(x_real_validate)
""" ---------------------------------------------------------------------------
PLOT THE RESULTS FROM x_vals & x_real_vals BELOW
--------------------------------------------------------------------------- """
font = {'size' : 16}
matplotlib.rc('font', **font)
f, ax = plt.subplots(ncols=1,nrows=len(x_vals), figsize=(6,7))
for i in xrange(len(x_vals)):
dl = float(signals[i])
plt.subplot(len(x_vals), 1, i+1)
color=cm.brg(np.linspace(0.0,0.5,len(x_vals)))
plt.plot(wavelengths, x_vals[i]/max(x_vals[i]), color=color[i], label="spacing = "+str(dl))
# plt.plot(wavelengths, x_real_vals[i]/max(x_real_vals[i]), 'ko')
if i==len(x_vals)//2:
plt.ylabel("Intensity [a.u.]")
plt.xlim([1550,1570])
plt.ylim([-0.05, 1.05])
frame = plt.gca()
if i!=len(x_vals)-1:
frame.axes.get_xaxis().set_ticklabels([])
frame.axes.get_yaxis().set_ticklabels([])
plt.yticks(np.linspace(0,1.,6))
hl, hw = 0.7, 0.07
arrow_height = 0.5
# if i==2:
# plt.text(1555.0, arrow_height, '400 pm', horizontalalignment='center', verticalalignment='center')
REDVOL = ['65','70','75','80','85','90','95']
#REDVOL = ['65','66','67','68','69','70','71','72','73','74','75','76','77','78','79','80','81','82','83','84','85','86','87','88','89','90','91','92','93','94','95','96','97','98','99','100']
CAPNUM = ['100']
CAPNUM = ['1']
#fs = 24
#plt.figure(figsize=(14,12))
fs = 10
plt.figure(figsize=(6,4))
plt.rc('text',usetex=True)
plt.rc('font', family='serif')
#for redvol in REDVOL :
for capnum in CAPNUM :
color = iter(cm.brg(np.linspace(0.0,0.8, len(REDVOL))))
for redvol in REDVOL :
c = next(color)
#(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA) = readMetrics(confin)
(STR, EPS, VOLRAT, LAMBDA, AREA, VLME, DP) = readMetrics(redvol,capnum)
v = round(float(redvol)/100.0,2)
# calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0
lamc = 1.0
if (abs(v - 0.99) < 1e-12) :
lamc = 1.090275107548953613318768
if (abs(v - 0.98) < 1e-12) :
lamc = 1.133537870150389282270677
if (abs(v - 0.97) < 1e-12) :
lamc = 1.169545928031930431602916
from readMetrics import readMetrics import numpy as np from matplotlib import pyplot as plt from matplotlib.pyplot import cm from scipy.interpolate import splrep, splev, interp1d from math import log10, floor MAXRAD = ['90', '92', '94', '96'] MAXRAD = ['90', '99'] #MAXRAD = ['80','82','84','86','88','90', '92', '94', '96', '98'] MAXRAD = ['90','91','92','93','94','95','96','97','98','99'] #MAXRAD = ['99'] color = iter(cm.brg(np.linspace(0.0,0.8, len(MAXRAD)))) fs = 24 plt.figure(figsize=(14,12)) for maxrad in MAXRAD : c = next(color) (STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, LENGTH, DP, ETA) = readMetrics(maxrad) v = round(REDVOL[0],2) lam = LAMBDA[0] # calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0 if (abs(v - 0.95) < 1e-12) : lamc = 1.2324 if (abs(v - 0.90) < 1e-12) : lamc = 1.3712 if (abs(v - 0.85) < 1e-12) : lamc = 1.5050 if (abs(v - 0.80) < 1e-12) :
vois = {}
for filename in sys.argv[1:]:
# Opens file and extract MNI coordinates
# ----------------------------------------------------------
f = open(filename, "r")
tokens = [line.split() for line in f.readlines()[1:]]
coords = [x[2:5] for x in tokens]
mni = [[float(y) for y in x] for x in coords]
name = tokens[0][1]
vois[name] = mni
#cols = cm.tab10(linspace(0,1,len(vois.keys())))
#cols = cm.rainbow(linspace(0,1,len(vois.keys())))
cols = cm.brg(linspace(0, 1, len(vois.keys())))
# Adds transparency
# ----------------------------------------------------------
for c in cols:
c[3] = 0.7 # Sets alpha
# Plots and saves
# ----------------------------------------------------------
# Calculates the proper size of the marker
L = [len(vois[x]) for x in vois.keys()]
N = numpy.max(L)
msize = 100
