#
# hello.load.45.py
#
# receive and display loading step response
# hello.step.45.py serial_port
#
# Neil Gershenfeld
# CBA MIT 10/29/10
#
# (c) Massachusetts Institute of Technology 2010
# Permission granted for experimental and personal use;
# license for commercial sale available from MIT
#

from Tkinter import *
from time import gmtime, strftime, time
import serial

WINDOW = 600 # window size
eps = 0.5 # filter time constant
filter1 = 0.0 # filtered value
filter2 = 0.0 # filtered value
filter3 = 0.0 # filtered value

# Define parameters for file read and write
RECSAMPLE = 2 # record every 3 minutes
sample1 = time() # time of last sample
sample2 = time() # time of new sample
samplestarttime = strftime("%a, %d %b %Y %H:%M:%S +0000", gmtime()) #log the start time of the sample
samplecurrenttime = strftime("%a, %d %b %Y %H:%M:%S +0000", gmtime()) #log the current sample time of the sample
f1sum = 0
f2sum = 0
f3sum = 0
index = 0
f1mean = 0
f2mean = 0
f3mean = 0

def idle(parent,canvas):
   global filter1, filter2, filter3, eps
   global sample1, sample2, samplestarttime, samplecurrenttime
   global f1sum, f2sum, f3sum, index, f1mean, f2mean, f3mean
   #
   # idle routine
   #
   byte2 = 0
   byte3 = 0
   byte4 = 0
   ser.flush()
   #
   # find framing 
   #
   while 1:
      byte1 = byte2
      byte2 = byte3
      byte3 = byte4
      byte4 = ord(ser.read())
      if ((byte1 == 1) & (byte2 == 2) & (byte3 == 3) & (byte4 == 4)):
         break
   #
   # read and plot
   #
   # read three sets of low and high bytes corresponding to 1ms, 10ms, and 100ms samples
   up_low1 = ord(ser.read())
   up_high1 = ord(ser.read())
   down_low1 = ord(ser.read())
   down_high1 = ord(ser.read())
   up_low2 = ord(ser.read())
   up_high2 = ord(ser.read())
   down_low2 = ord(ser.read())
   down_high2 = ord(ser.read())
   up_low3 = ord(ser.read())
   up_high3 = ord(ser.read())
   down_low3 = ord(ser.read())
   down_high3 = ord(ser.read())
   # format bytes read into first 16 bits and then 32 bits
   up_value1 = 256*up_high1 + up_low1
   down_value1 = 256*down_high1 + down_low1
   up_value2 = 256*up_high2 + up_low2
   down_value2 = 256*down_high2 + down_low2
   up_value3 = 256*up_high3 + up_low3
   down_value3 = 256*down_high3 + down_low3
   value1 = (up_value1 + (1023 - down_value1))/2.0
   value2 = (up_value2 + (1023 - down_value2))/2.0
   value3 = (up_value3 + (1023 - down_value3))/2.0
   filter1 = (1-eps)*filter1 + eps*value1
   filter2 = (1-eps)*filter2 + eps*value2
   filter3 = (1-eps)*filter3 + eps*value3
   # converting three rc filter values into integers
   x1 = int(.2*WINDOW + (.9-.2)*WINDOW*filter1/1023.0)
   x2 = int(.2*WINDOW + (.9-.2)*WINDOW*filter2/1023.0)
   x3 = int(.2*WINDOW + (.9-.2)*WINDOW*filter3/1023.0)
   canvas.itemconfigure("text1",text="%.1f"%filter1)
   canvas.itemconfigure("text2",text="%.1f"%filter2)
   canvas.itemconfigure("text3",text="%.1f"%filter3)
   canvas.coords('rect11',.2*WINDOW,.05*WINDOW,x1,.2*WINDOW)
   canvas.coords('rect12',x1,.05*WINDOW,.9*WINDOW,.2*WINDOW)
   canvas.coords('rect21',.2*WINDOW,.3*WINDOW,x2,.45*WINDOW)
   canvas.coords('rect22',x2,.3*WINDOW,.9*WINDOW,.45*WINDOW)
   canvas.coords('rect31',.2*WINDOW,.55*WINDOW,x3,.7*WINDOW)
   canvas.coords('rect32',x3,.55*WINDOW,.9*WINDOW,.7*WINDOW)
   canvas.update()

   # record samples in a file
   sample2 = time(); 
   samplelapse = sample2 - sample1
   if samplelapse > RECSAMPLE:
      index = index + 1
      f = open('logsamples.txt', 'a')
      f1 = int(filter1)
      f2 = int(filter2)
      f3 = int(filter3)
      f1sum = f1sum + f1
      f2sum = f2sum + f2
      f3sum = f3sum + f3

      f1mean = f1sum/index 
      f2mean = f2sum/index 
      f3mean = f3sum/index 
      value = (f1, f2, f3,strftime("%a, %d %b %Y %H:%M:%S +0000", gmtime()))
      samplecurrenttime = strftime("%a, %d %b %Y %H:%M:%S +0000", gmtime()) #log the current sample time of the sample
      s = str(value) + '\n' 
      f.write(s)
      f.close()
      sample1 = sample2

   canvas.itemconfigure("text5",text="%s"%samplestarttime)
   canvas.itemconfigure("text15",text="%s"%samplecurrenttime)
   canvas.itemconfigure("text7",text="%d"%RECSAMPLE)
   canvas.itemconfigure("text9",text="%d"%f1mean)
   canvas.itemconfigure("text11",text="%d"%f2mean)
   canvas.itemconfigure("text13",text="%d"%f3mean)

   parent.after_idle(idle,parent,canvas)

#
#  check command line arguments
#
if (len(sys.argv) != 2):
   print "command line: python history.load.45.py serial_port"
   sys.exit()
port = sys.argv[1]
#
# open serial port
#
ser = serial.Serial(port,9600)
ser.setDTR()
#
# set up GUI
#
root = Tk()
root.title('history.load.45.py (q to exit)')
root.bind('q','exit')
canvas = Canvas(root, width=WINDOW, height=1.5*WINDOW, background='white')
#
canvas.create_text(.1*WINDOW,.125*WINDOW,text="1",font=("Helvetica", 24),tags="text1",fill="#0000b0")
canvas.create_rectangle(.2*WINDOW,.05*WINDOW,.3*WINDOW,.2*WINDOW, tags='rect11', fill='#b00000')
canvas.create_rectangle(.3*WINDOW,.05*WINDOW,.9*WINDOW,.2*WINDOW, tags='rect12', fill='#0000b0')
#
canvas.create_text(.1*WINDOW,.375*WINDOW,text="2",font=("Helvetica", 24),tags="text2",fill="#0000b0")
canvas.create_rectangle(.2*WINDOW,.3*WINDOW,.3*WINDOW,.45*WINDOW, tags='rect21', fill='#b00000')
canvas.create_rectangle(.3*WINDOW,.3*WINDOW,.9*WINDOW,.45*WINDOW, tags='rect22', fill='#0000b0')
#
canvas.create_text(.1*WINDOW,.625*WINDOW,text="3",font=("Helvetica", 24),tags="text3",fill="#0000b0")
canvas.create_rectangle(.2*WINDOW,.55*WINDOW,.3*WINDOW,.7*WINDOW, tags='rect31', fill='#b00000')
canvas.create_rectangle(.3*WINDOW,.55*WINDOW,.9*WINDOW,.7*WINDOW, tags='rect32', fill='#0000b0')

# Display Data Statistics
canvas.create_text(.3*WINDOW,.8*WINDOW,text="Data Collection Start Time:",font=("Helvetica", 18),tags="text4",fill="#b00000")
canvas.create_text(.75*WINDOW,.8*WINDOW,text="1",font=("Helvetica", 18),tags="text5",fill="#b00000")

canvas.create_text(.3*WINDOW,.875*WINDOW,text="Data Collection End Time:",font=("Helvetica", 18),tags="text14",fill="#b00000")
canvas.create_text(.75*WINDOW,.875*WINDOW,text="2",font=("Helvetica", 18),tags="text15",fill="#b00000")

canvas.create_text(.3*WINDOW,.95*WINDOW,text="Time Between Samples (secs):",font=("Helvetica", 18),tags="text6",fill="#b00000")
canvas.create_text(.75*WINDOW,.95*WINDOW,text="2",font=("Helvetica", 18),tags="text7",fill="#b00000")

canvas.create_text(.3*WINDOW,1.025*WINDOW,text="Sample1 (mean):",font=("Helvetica", 24),tags="text8",fill="#b00000")
canvas.create_text(.65*WINDOW,1.025*WINDOW,text="2",font=("Helvetica", 24),tags="text9",fill="#b00000")

canvas.create_text(.3*WINDOW,1.1*WINDOW,text="Sample2 (mean):",font=("Helvetica", 24),tags="text10",fill="#b00000")
canvas.create_text(.65*WINDOW,1.1*WINDOW,text="2",font=("Helvetica", 24),tags="text11",fill="#b00000")

canvas.create_text(.3*WINDOW,1.175*WINDOW,text="Sample3 (mean):",font=("Helvetica", 24),tags="text12",fill="#b00000")
canvas.create_text(.65*WINDOW,1.175*WINDOW,text="2",font=("Helvetica", 24),tags="text13",fill="#b00000")

canvas.pack()
#
# start idle loop
#
root.after(100,idle,root,canvas)
root.mainloop()