playground

2 files changed, 209 insertions, 0 deletions

diff --git a/labjackStreamTest.py b/labjackStreamTest.py
new file mode 100644
index 0000000..d4c6a9e
--- /dev/null
+++ b/labjackStreamTest.py
@@ -0,0 +1,148 @@
+import sys
+import traceback
+from datetime import datetime
+
+import u3
+import u6
+import ue9
+
+
+# MAX_REQUESTS is the number of packets to be read.
+MAX_REQUESTS = 75
+# SCAN_FREQUENCY is the scan frequency of stream mode in Hz
+SCAN_FREQUENCY = 5000
+
+d = None
+
+###############################################################################
+# U3
+# Uncomment these lines to stream from a U3
+###############################################################################
+'''
+# At high frequencies ( >5 kHz), the number of samples will be MAX_REQUESTS
+# times 48 (packets per request) times 25 (samples per packet).
+d = u3.U3()
+
+# To learn the if the U3 is an HV
+d.configU3()
+
+# For applying the proper calibration to readings.
+d.getCalibrationData()
+
+# Set the FIO0 and FIO1 to Analog (d3 = b00000011)
+d.configIO(FIOAnalog=3)
+
+print("Configuring U3 stream")
+d.streamConfig(NumChannels=2, PChannels=[0, 1], NChannels=[31, 31], Resolution=3, ScanFrequency=SCAN_FREQUENCY)
+'''
+
+###############################################################################
+# U6
+# Uncomment these lines to stream from a U6
+###############################################################################
+'''
+# At high frequencies ( >5 kHz), the number of samples will be MAX_REQUESTS
+# times 48 (packets per request) times 25 (samples per packet).
+d = u6.U6()
+
+# For applying the proper calibration to readings.
+d.getCalibrationData()
+
+print("Configuring U6 stream")
+
+d.streamConfig(NumChannels=2, ChannelNumbers=[0, 1], ChannelOptions=[0, 0], SettlingFactor=1, ResolutionIndex=1, ScanFrequency=SCAN_FREQUENCY)
+'''
+
+###############################################################################
+# UE9
+# Uncomment these lines to stream from a UE9
+###############################################################################
+# At 96 Hz or higher frequencies, the number of samples will be MAX_REQUESTS
+# times 8 (packets per request) times 16 (samples per packet).
+# Currently over ethernet packets per request is 1.
+d = ue9.UE9()
+#d = ue9.UE9(ethernet=True, ipAddress="192.168.1.209")  # Over TCP/ethernet connect to UE9 with IP address 192.168.1.209
+
+# For applying the proper calibration to readings.
+c=d.getCalibrationData()
+
+# by evmik
+# fixing missing slope for gain '0'
+c['AINSlopes']['0']= 0.0000775030
+
+
+print("Configuring UE9 stream")
+
+d.streamConfig(NumChannels=2, ChannelNumbers=[0, 1], ChannelOptions=[0, 0], SettlingTime=0, Resolution=12, ScanFrequency=SCAN_FREQUENCY)
+'''
+'''
+
+if d is None:
+    print("""Configure a device first.
+Please open streamTest.py in a text editor and uncomment the lines for your device.
+
+Exiting...""")
+    sys.exit(0)
+
+try:
+    print("Start stream")
+    d.streamStart()
+    start = datetime.now()
+    print("Start time is %s" % start)
+
+    missed = 0
+    dataCount = 0
+    packetCount = 0
+
+    for r in d.streamData():
+        if r is not None:
+            # Our stop condition
+            if dataCount >= MAX_REQUESTS:
+                break
+
+            if r["errors"] != 0:
+                print("Errors counted: %s ; %s" % (r["errors"], datetime.now()))
+
+            if r["numPackets"] != d.packetsPerRequest:
+                print("----- UNDERFLOW : %s ; %s" %
+                      (r["numPackets"], datetime.now()))
+
+            if r["missed"] != 0:
+                missed += r['missed']
+                print("+++ Missed %s" % r["missed"])
+
+            # Comment out these prints and do something with r
+            print("Average of %s AIN0, %s AIN1 readings: %s, %s" %
+                  (len(r["AIN0"]), len(r["AIN1"]), sum(r["AIN0"])/len(r["AIN0"]), sum(r["AIN1"])/len(r["AIN1"])))
+
+            dataCount += 1
+            packetCount += r['numPackets']
+        else:
+            # Got no data back from our read.
+            # This only happens if your stream isn't faster than the USB read
+            # timeout, ~1 sec.
+            print("No data ; %s" % datetime.now())
+except:
+    print("".join(i for i in traceback.format_exc()))
+finally:
+    stop = datetime.now()
+    d.streamStop()
+    print("Stream stopped.\n")
+    d.close()
+
+    sampleTotal = packetCount * d.streamSamplesPerPacket
+
+    scanTotal = sampleTotal / 2  # sampleTotal / NumChannels
+    print("%s requests with %s packets per request with %s samples per packet = %s samples total." %
+          (dataCount, (float(packetCount)/dataCount), d.streamSamplesPerPacket, sampleTotal))
+    print("%s samples were lost due to errors." % missed)
+    sampleTotal -= missed
+    print("Adjusted number of samples = %s" % sampleTotal)
+
+    runTime = (stop-start).seconds + float((stop-start).microseconds)/1000000
+    print("The experiment took %s seconds." % runTime)
+    print("Actual Scan Rate = %s Hz" % SCAN_FREQUENCY)
+    print("Timed Scan Rate = %s scans / %s seconds = %s Hz" %
+          (scanTotal, runTime, float(scanTotal)/runTime))
+    print("Timed Sample Rate = %s samples / %s seconds = %s Hz" %
+          (sampleTotal, runTime, float(sampleTotal)/runTime))
diff --git a/panel.py b/panel.py
new file mode 100644
index 0000000..a9fab18
--- /dev/null
+++ b/panel.py
@@ -0,0 +1,61 @@
+import tkinter as tk
+from tkinter import *
+from tkinter import ttk
+
+import matplotlib.pyplot as plt
+from matplotlib.figure import Figure
+from matplotlib.backends.backend_tkagg import (FigureCanvasTkAgg, NavigationToolbar2Tk)
+
+cnt=1
+
+class Experiment:
+
+    def __init__(self, root):
+        self.gui(root)
+
+    def gui(self, root): 
+
+        self.bStart=Button(root,text="START",command=self.hello,font=('Arial','24'))
+        self.bStart.pack()
+
+        self.bStop=Button(root,text="STOP",command=exit,font=('Arial','24'))
+        self.bStop.pack()
+
+        self.fig=plt.figure()
+
+        self.ax = self.fig.add_subplot(1,1,1)
+        self.ax.set_xlim([0,20])
+        self.ax.set_ylim([0,20])
+        self.ax.plot([i for i in range(10)],[i for i in range(10)])
+
+        self.canvas = FigureCanvasTkAgg(self.fig, master = root)
+        self.canvas.draw()
+
+        # placing the canvas on the Tkinter window
+        self.canvas.get_tk_widget().pack()
+               
+        # creating the Matplotlib toolbar
+        self.toolbar = NavigationToolbar2Tk(self.canvas, root)
+        self.toolbar.update()
+                             
+        # placing the toolbar on the Tkinter window
+        self.canvas.get_tk_widget().pack()
+
+    def hello(self):
+        global cnt
+        cnt = cnt + 1
+        self.ax.plot([i for i in range(20)],[i/cnt for i in range(20)])
+        self.canvas.draw()
+        print("Hello ")
+        print(cnt)
+
+ 
+if __name__ == '__main__': 
+    cnt=1
+
+    root=Tk()
+
+    Experiment(root) 
+
+    root.mainloop()
+