#!/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()