#!/usr/bin/python2
# -*- coding: utf-8 -*-
import matplotlib
matplotlib.rcParams['backend'] = 'Agg'
matplotlib.rc('font', size='10')
matplotlib.rc('text', usetex=True) # FIXME !! nicer but slower !!!
#import pygolem_lite
from numpy import *
#from pygolem_lite.config import *
from pygolem_lite.modules import *
#from matplotlib.pyplot import *
import time
from shutil import copy
import os, sys
from scipy.interpolate import interp1d
def getdata():
for i in range(32):
try:
tvec, data = load_adv("PapouchSt_%02i" % i )
if 'tvec_0' in locals():
data = interp1d(tvec, squeeze( data), copy=False, bounds_error=False, fill_value=0)(tvec_0) # interplate to the same time resolution !!!
data_0 = vstack([data_0.T, data]).T
else:
data_0 = data
tvec_0 = tvec
except IOError:
pass
if 'tvec_0' in locals():
save_adv('Papouch_All', tvec_0, data_0)
else:
print "Missing data Papouch"
def graphs(file_type):
name = "Papouch Stockel"
# graph all
channels = [ 'ch '+str(i) for i in range(1,13)] # channels from 1 to 12 !!
# integrated
out = get_data('Papouch_All', channels, "",xlim=[0,40] ,integrated=True )
multiplot(out, name + ' Integrated' , 'graph_int', (6,len(out)*2), 100, 'vertical', file_type)
os.system("convert -resize 150x120\! graph_int.png icon.png")
# graph all
out = get_data('Papouch_All', channels , "",xlim=[0,40] ,integrated=False, reduction = False)
multiplot(out, name , 'graph', (6,len(out)*2), 100, 'vertical', file_type)
os.system("convert -resize 150x120\! graph.png icon.png")
# icon
#out = get_data('Papouch_All', "", "" , xlabel = "",integrated=False )
#paralel_multiplot(out, "" , 'icon', (4,3), 40)
def main():
if sys.argv[1] == "acquisition":
getdata()
if sys.argv[1] == "plots":
graphs('png')
#graphs('svg')
saveconst('status', 0)
if __name__ == "__main__":
main()
