-
Notifications
You must be signed in to change notification settings - Fork 2
/
testDMV.py
114 lines (98 loc) · 3.01 KB
/
testDMV.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
# %% Imports
from readDMV import readDMV
import matplotlib.pyplot as plt
# %% RNC
c1 = readDMV('/Users/vonw/data/paeri/raw/AE160201/160201C1.RNC')
plt.figure()
c1.mean_rad[0].plot()
c1.mean_rad[-1].plot()
plt.show()
c2 = readDMV('/Users/vonw/data/paeri/raw/AE160201/160201C2.RNC')
plt.figure()
c2.mean_rad[0].plot()
c2.mean_rad[-1].plot()
plt.show()
# %% RLC
b1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602B1.RLC')
plt.figure()
b1.atmosphericRadiance[0].plot()
b1.atmosphericRadiance[-1].plot()
plt.show()
f1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602F1.RLC')
plt.figure()
f1.atmosphericRadiance[0].plot()
f1.atmosphericRadiance[-1].plot()
plt.show()
b2 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602B2.RLC')
plt.figure()
b2.atmosphericRadiance[0].plot()
b2.atmosphericRadiance[-1].plot()
plt.show()
f2 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602F2.RLC')
plt.figure()
f2.atmosphericRadiance[0].plot()
f2.atmosphericRadiance[-1].plot()
plt.show()
c1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602C1.RLC')
plt.figure()
c1.averageRadiance[0].plot()
c1.averageRadiance[-1].plot()
plt.show()
c2 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602C2.RLC')
plt.figure()
c2.averageRadiance[0].plot()
c2.averageRadiance[-1].plot()
plt.show()
# %% CXS
b1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602B1.CXS')
plt.figure()
b1.Ch1BackwardScanRealPartCounts[0].plot()
b1.Ch1BackwardScanImagPartCounts[0].plot()
b1.Ch1BackwardScanRealPartCounts[-1].plot()
b1.Ch1BackwardScanImagPartCounts[-1].plot()
plt.show()
f1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602F1.CXS')
plt.figure()
f1.Ch1ForwardScanRealPartCounts[0].plot()
f1.Ch1ForwardScanImagPartCounts[0].plot()
f1.Ch1ForwardScanRealPartCounts[-1].plot()
f1.Ch1ForwardScanImagPartCounts[-1].plot()
plt.show()
b2 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602B2.CXS')
plt.figure()
b2.Ch2BackwardScanRealPartCounts[0].plot()
b2.Ch2BackwardScanImagPartCounts[0].plot()
b2.Ch2BackwardScanRealPartCounts[-1].plot()
b2.Ch2BackwardScanImagPartCounts[-1].plot()
plt.show()
f2 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602F2.CXS')
plt.figure()
f2.Ch2ForwardScanRealPartCounts[0].plot()
f2.Ch2ForwardScanImagPartCounts[0].plot()
f2.Ch2ForwardScanRealPartCounts[-1].plot()
f2.Ch2ForwardScanImagPartCounts[-1].plot()
plt.show()
# %% CXV
b1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602B1.CXV')
plt.figure()
b1.RealPartCounts[0].plot()
b1.ImagPartCounts[0].plot()
b1.RealPartCounts[-1].plot()
b1.ImagPartCounts[-1].plot()
plt.show()
f1 = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602F1.CXV')
plt.figure()
f1.RealPartCounts[0].plot()
f1.ImagPartCounts[0].plot()
f1.RealPartCounts[-1].plot()
f1.ImagPartCounts[-1].plot()
plt.show()
# %% SUM
sm = readDMV('/Users/vonw/data/paeri/reprocessed/AE110602/110602.SUM')
plt.figure()
sm.ABBapexTemp.plot()
plt.show()
plt.figure()
sm.HBB2minNENestimateNo2ch1[0].plot()
sm.HBB2minNENestimateNo2ch1[-1].plot()
plt.show()