Changeset 636 in ETALON for CLIO

Timestamp:

Feb 6, 2017, 6:31:47 PM (7 years ago)

Author:

delerue

Message:

CLIO analysis script updated

File:

• CLIO/control/read_data.py (modified) (1 diff)

CLIO/control/read_data.py

@@ -6 +6 @@
 """
 
+
 import numpy as np
-#regexp = r"([0-9A-F]+)(...)"  # match [digits, whitespace, anything]
-#regexp = r"([0-9A-F]+)"  # match [digits, whitespace, anything]
-#output = np.fromregex('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-04-38.txt', regexp, [('num', np.int64)])
-#output = np.fromregex('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-04-38_head5.txt', regexp, [('num', np.uint32)])
-#output = np.fromregex('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-04-38_head5.txt', regexp, [('num', 'S8')])
-#print(output)
-#print(output['num'])
+import matplotlib.pyplot as plt
 
+def read_data(filename):
+    alldata=np.loadtxt(filename, dtype='int16', comments='#', delimiter='  ', usecols=range(1,33), converters={ _:lambda s:  int(s , 16)  for _ in range(0,33) })
+    #alldata=np.loadtxt('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-40-25_tail20.txt', dtype='int16', delimiter='  ', usecols=range(1,33), converters={ _:lambda s:  int('0'+s , 16)  for _ in range(0,33) })
+    #print(alldata)
+    data_length=len(alldata[:,1])
+    print(alldata[:,1])
+    print(range(1,data_length+1))
 
-#alldata=np.loadtxt('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-04-38_head5.txt', dtype='int16', comments='#', delimiter='  ', usecols=range(1,33), converters={ _:lambda s:  int(s , 16)  for _ in range(0,33) })
-
-alldata=np.loadtxt('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-40-25_tail20.txt', dtype='int16', delimiter='  ', usecols=range(1,33), converters={ _:lambda s:  int('0'+s , 16)  for _ in range(0,33) })
-print(alldata)
+    fig=plt.figure(figsize=(6, 4))
+    for idata in range(0,32):
+        txtline='r--'
+        if idata == 1:
+            txtline='bs'
+        if idata == 2:
+            txtline='g^'
+        
+        plt.plot(range(1,data_length+1),alldata[:,idata], txtline)
+    plt.show()
+    imagename=filename.replace('.txt','.png')
+    print(imagename)
+    fig.savefig(imagename)
+    plt.close(fig)
+    return;
+read_data('/Users/delerue/Downloads/data_CLIO_2017-01-23_22-04-38_head5.txt')