##// END OF EJS Templates
DBS and SA scripts corrections
Julio Valdez -
r763:df4875709e0e
parent child
Show More
@@ -1,169 +1,169
1 1 # DIAS 19 Y 20 FEB 2014
2 2 # Comprobacion de Resultados DBS con SA
3 3
4 4 import os, sys
5 5
6 6 path = os.path.split(os.getcwd())[0]
7 sys.path.append(path)
7 path = os.path.split(path)[0]
8 8
9 from controller import *
9 sys.path.insert(0, path)
10
11 from schainpy.controller import Project
10 12
11 13 desc = "DBS Experiment Test"
12 14 filename = "DBStest.xml"
13 15
14 16 controllerObj = Project()
15 17
16 18 controllerObj.setup(id = '191', name='test01', description=desc)
17 19
18 20 #Experimentos
19 21
20 22 #2014050 19 Feb 2014
21 23 # path = '/home/soporte/Documents/MST_Data/DBS/d2014050'
22 24 # pathFigure = '/home/soporte/workspace/Graficos/DBS/d2014050p/'
23 25 # xmin = '15.5'
24 26 # xmax = '23.99999999'
25 27 # startTime = '17:25:00'
26 28 # filehdf5 = "DBS_2014050.hdf5"
27 29
28 30 #2014051 20 Feb 2014
29 path = '/home/soporte/Data/MST/DBS/d2014051'
30 pathFigure = '/home/soporte/workspace/Graficos/DBS/prueba1/'
31 path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/MST/DBS/d2014051'
32 # path = '/media/joscanoa/disco4/Data/2014/DBS_SA JAN 2014/DBS_SA/250/d2014050'
33 pathFigure = '/home/joscanoa/Pictures/testHDF5/figure'
34 xmax = '4'
31 35 xmin = '0'
32 xmax = '7.5'
33 36 startTime = '00:00:00'
34 37 filehdf5 = "DBS_2014051.hdf5"
35 38
36 39
37 40
38 41 #------------------------------------------------------------------------------------------------
39 42 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
40 43 path=path,
41 44 startDate='2014/01/31',
42 45 endDate='2014/03/31',
43 46 startTime=startTime,
44 47 endTime='23:59:59',
45 48 online=0,
46 49 delay=5,
47 50 walk=0)
48 51
49 52 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
50 53
51 54
52 55 #--------------------------------------------------------------------------------------------------
53 56
54 57 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
55 58
56 59 opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
57 60
58 61 opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other')
59 62 opObj11.addParameter(name='n', value='256', format='int')
60 63 # opObj11.addParameter(name='n', value='16', format='int')
61 64
62 65 opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex')
63 66 opObj11.addParameter(name='minIndex', value='10', format='float')
64 67 opObj11.addParameter(name='maxIndex', value='60', format='float')
65 68
66 69 #---------------------------------------------------------------------------------------------------
67 70
68 71 procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId())
69 72 procUnitConfObj1.addParameter(name='nFFTPoints', value='64', format='int')
70 73 procUnitConfObj1.addParameter(name='nProfiles', value='64', format='int')
71 74 # procUnitConfObj1.addParameter(name='ippFactor', value='2', format='int')
72 75 procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(0,1),(2,1)', format='pairsList')
73 76
74 77 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
75 78 opObj11.addParameter(name='n', value='5', format='int')
76 79
80 # opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
81 # opObj14.addParameter(name='id', value='1', format='int')
82 # opObj14.addParameter(name='wintitle', value='Con interf', format='str')
83 # opObj14.addParameter(name='save', value='0', format='bool')
84 # opObj14.addParameter(name='figpath', value=pathFigure, format='str')
85 # opObj14.addParameter(name='zmin', value='5', format='int')
86 # opObj14.addParameter(name='zmax', value='90', format='int')
87
88 opObj12 = procUnitConfObj1.addOperation(name='removeInterference')
89 opObj13 = procUnitConfObj1.addOperation(name='removeDC')
90 opObj13.addParameter(name='mode', value='1', format='int')
91 #
77 92 opObj14 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other')
78 93 opObj14.addParameter(name='id', value='1', format='int')
79 94 opObj14.addParameter(name='wintitle', value='Con interf', format='str')
80 opObj14.addParameter(name='save', value='1', format='bool')
95 opObj14.addParameter(name='save', value='0', format='bool')
81 96 opObj14.addParameter(name='figpath', value=pathFigure, format='str')
82 97 opObj14.addParameter(name='zmin', value='5', format='int')
83 98 opObj14.addParameter(name='zmax', value='90', format='int')
84 99
85 opObj12 = procUnitConfObj1.addOperation(name='removeInterference')
86 opObj13 = procUnitConfObj1.addOperation(name='removeDC')
87 opObj13.addParameter(name='mode', value='1', format='int')
88 100
89 101 opObj12 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other')
90 102 opObj12.addParameter(name='id', value='2', format='int')
91 103 opObj12.addParameter(name='wintitle', value='RTI Plot', format='str')
92 104 opObj12.addParameter(name='save', value='1', format='bool')
93 105 opObj12.addParameter(name='figpath', value = pathFigure, format='str')
94 106 opObj12.addParameter(name='xmin', value=xmin, format='float')
95 107 opObj12.addParameter(name='xmax', value=xmax, format='float')
96 108 opObj12.addParameter(name='zmin', value='5', format='int')
97 109 opObj12.addParameter(name='zmax', value='90', format='int')
98 110
99 111 #--------------------------------------------------------------------------------------------------
100 112
101 113 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
102 114 opObj20 = procUnitConfObj2.addOperation(name='GetMoments')
103 115
104 opObj21 = procUnitConfObj2.addOperation(name='MomentsPlot', optype='other')
105 opObj21.addParameter(name='id', value='3', format='int')
106 opObj21.addParameter(name='wintitle', value='Moments Plot', format='str')
107 opObj21.addParameter(name='save', value='1', format='bool')
108 opObj21.addParameter(name='figpath', value=pathFigure, format='str')
109 opObj21.addParameter(name='zmin', value='5', format='int')
110 opObj21.addParameter(name='zmax', value='90', format='int')
111
112 opObj21 = procUnitConfObj2.addOperation(name='ParametersPlot', optype='other')
113 opObj21.addParameter(name='id', value='5', format='int')
114 opObj21.addParameter(name='wintitle', value='Radial Velocity Plot', format='str')
115 opObj21.addParameter(name='save', value='1', format='bool')
116 opObj21.addParameter(name='figpath', value=pathFigure, format='str')
117 opObj21.addParameter(name='SNRmin', value='-10', format='int')
118 opObj21.addParameter(name='SNRmax', value='60', format='int')
119 opObj21.addParameter(name='SNRthresh', value='0', format='float')
120 opObj21.addParameter(name='xmin', value=xmin, format='float')
121 opObj21.addParameter(name='xmax', value=xmax, format='float')
122
123 opObj21 = procUnitConfObj2.addOperation(name='ParametersPlot', optype='other')
124 opObj21.addParameter(name='id', value='6', format='int')
125 opObj21.addParameter(name='wintitle', value='Spectral width Plot', format='str')
126 opObj21.addParameter(name='save', value='1', format='bool')
127 opObj21.addParameter(name='figpath', value=pathFigure, format='str')
128 opObj21.addParameter(name='SNRmin', value='-10', format='int')
129 opObj21.addParameter(name='SNRmax', value='60', format='int')
130 opObj21.addParameter(name='SNRthresh', value='0', format='float')
131 opObj21.addParameter(name='xmin', value=xmin, format='float')
132 opObj21.addParameter(name='xmax', value=xmax, format='float')
133 opObj21.addParameter(name='zmin', value=0, format='float')
134 opObj21.addParameter(name='paramIndex', value=2, format='int')
135 opObj21.addParameter(name='onlyPositive', value=1, format='bool')
116 # opObj21 = procUnitConfObj2.addOperation(name='MomentsPlot', optype='other')
117 # opObj21.addParameter(name='id', value='3', format='int')
118 # opObj21.addParameter(name='wintitle', value='Moments Plot', format='str')
119 # opObj21.addParameter(name='save', value='1', format='bool')
120 # opObj21.addParameter(name='figpath', value=pathFigure, format='str')
121 # opObj21.addParameter(name='zmin', value='5', format='int')
122 # opObj21.addParameter(name='zmax', value='90', format='int')
123
124 # opObj21 = procUnitConfObj2.addOperation(name='ParametersPlot', optype='other')
125 # opObj21.addParameter(name='id', value='5', format='int')
126 # opObj21.addParameter(name='wintitle', value='Radial Velocity Plot', format='str')
127 # opObj21.addParameter(name='save', value='0', format='bool')
128 # opObj21.addParameter(name='figpath', value=pathFigure, format='str')
129 # opObj21.addParameter(name='SNRmin', value='-10', format='int')
130 # opObj21.addParameter(name='SNRmax', value='60', format='int')
131 # opObj21.addParameter(name='channelList', value='0,2', format='intlist')
132 # opObj21.addParameter(name='SNR', value='1', format='bool')
133 # opObj21.addParameter(name='SNRthresh', value='0', format='float')
134 # opObj21.addParameter(name='xmin', value=xmin, format='float')
135 # opObj21.addParameter(name='xmax', value=xmax, format='float')
136 136
137 137 opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
138 138 opObj22.addParameter(name='technique', value='DBS', format='str')
139 139 opObj22.addParameter(name='correctAzimuth', value='51.06', format='float')
140 140 opObj22.addParameter(name='correctFactor', value='-1', format='float')
141 141 opObj22.addParameter(name='dirCosx', value='0.041016, 0, -0.054688', format='floatlist')
142 142 opObj22.addParameter(name='dirCosy', value='-0.041016, 0.025391, -0.023438', format='floatlist')
143 # opObj22.addParameter(name='horizontalOnly', value='1', format='bool')
144 # opObj22.addParameter(name='channelList', value='1,2', format='intlist')
145
143 # opObj22.addParameter(name='horizontalOnly', value='0', format='bool')
144 # opObj22.addParameter(name='channelList', value='1,2,3', format='intlist')
145
146 146 opObj23 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
147 opObj23.addParameter(name='id', value='4', format='int')
147 opObj23.addParameter(name='id', value='3', format='int')
148 148 opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
149 149 opObj23.addParameter(name='save', value='1', format='bool')
150 150 opObj23.addParameter(name='figpath', value = pathFigure, format='str')
151 151 opObj23.addParameter(name='zmin', value='-10', format='int')
152 152 opObj23.addParameter(name='zmax', value='10', format='int')
153 153 opObj23.addParameter(name='zmin_ver', value='-80', format='float')
154 154 opObj23.addParameter(name='zmax_ver', value='80', format='float')
155 155 opObj23.addParameter(name='SNRmin', value='-10', format='int')
156 156 opObj23.addParameter(name='SNRmax', value='60', format='int')
157 157 opObj23.addParameter(name='SNRthresh', value='0', format='float')
158 158 opObj23.addParameter(name='xmin', value=xmin, format='float')
159 159 opObj23.addParameter(name='xmax', value=xmax, format='float')
160 160
161 161 #--------------------------------------------------------------------------------------------------
162 162 print "Escribiendo el archivo XML"
163 163 controllerObj.writeXml(filename)
164 164 print "Leyendo el archivo XML"
165 165 controllerObj.readXml(filename)
166 166
167 167 controllerObj.createObjects()
168 168 controllerObj.connectObjects()
169 169 controllerObj.run() No newline at end of file
@@ -1,139 +1,141
1 1 # DIAS 19 Y 20 FEB 2014
2 2 # Comprobacion de Resultados DBS con SA
3 3
4 4 import os, sys
5 5
6 6 path = os.path.split(os.getcwd())[0]
7 sys.path.append(path)
7 path = os.path.split(path)[0]
8 8
9 from controller import *
9 sys.path.insert(0, path)
10
11 from schainpy.controller import Project
10 12
11 13 desc = "SA Experiment Test"
12 14 filename = "SA2014050.xml"
13 15
14 16 controllerObj = Project()
15 17
16 18 controllerObj.setup(id = '191', name='test01', description=desc)
17 19
18 20
19 21 #Experimentos
20 22
21 23 #2014050 19 Feb 2014
22 path = '/home/soporte/Data/MST/SA/d2014050'
23 pathFigure = '/home/soporte/workspace/Graficos/SA/prueba1/'
24 path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/MST/SA/d2014050'
25 pathFigure = '/media/joscanoa/84A65E64A65E5730/soporte/workspace/Graficos/SA/prueba1/'
24 26 xmin = '15.5'
25 27 xmax = '24'
26 28 startTime = '15:30:00'
27 29 filehdf5 = "SA_2014050.hdf5"
28 30
29 31 #2014051 20 Feb 2014
30 32 # path = '/home/soporte/Data/MST/SA/d2014051'
31 33 # pathFigure = '/home/soporte/workspace/Graficos/SA/new/'
32 34 # xmin = '0.0'
33 35 # xmax = '8.0'
34 36 # startTime = '00:00:00'
35 37 # filehdf5 = "SA_2014051.hdf5"
36 38
37 39 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
38 40 path=path,
39 41 startDate='2014/01/01',
40 42 endDate='2014/03/31',
41 43 startTime=startTime,
42 44 endTime='23:59:59',
43 45 online=0,
44 46 delay=5,
45 47 walk=0)
46 48
47 49 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
48 50
49 51
50 52 #--------------------------------------------------------------------------------------------------
51 53
52 54 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
53 55
54 56 opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
55 57
56 58 opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other')
57 59 opObj11.addParameter(name='n', value='600', format='int')
58 60 # opObj11.addParameter(name='n', value='10', format='int')
59 61
60 62 opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex')
61 63 opObj11.addParameter(name='minIndex', value='10', format='float')
62 64 opObj11.addParameter(name='maxIndex', value='60', format='float')
63 65 #---------------------------------------------------------------------------------------------------
64 66 procUnitConfObj1 = controllerObj.addProcUnit(datatype='CorrelationProc', inputId=procUnitConfObj0.getId())
65 67 # procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(1,0),(2,3)', format='pairsList')
66 68 procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList')
67 69 procUnitConfObj1.addParameter(name='fullT', value='1', format='bool')
68 70 procUnitConfObj1.addParameter(name='removeDC', value='1', format='bool')
69 71 #procUnitConfObj1.addParameter(name='lagT', value='0,1,2,3', format='intlist')
70 72
71 73 opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
72 74 opObj12.addParameter(name='id', value='1', format='int')
73 75 opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
74 76 opObj12.addParameter(name='save', value='1', format='bool')
75 77 opObj12.addParameter(name='zmin', value='0', format='int')
76 78 opObj12.addParameter(name='zmax', value='1', format='int')
77 79 opObj12.addParameter(name='figpath', value = pathFigure, format='str')
78 80
79 81 opObj12 = procUnitConfObj1.addOperation(name='removeNoise')
80 82 opObj12.addParameter(name='mode', value='2', format='int')
81 83 opObj12 = procUnitConfObj1.addOperation(name='calculateNormFactor')
82 84
83 85 opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other')
84 86 opObj12.addParameter(name='id', value='2', format='int')
85 87 opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str')
86 88 opObj12.addParameter(name='save', value='1', format='bool')
87 89 opObj12.addParameter(name='zmin', value='0', format='int')
88 90 opObj12.addParameter(name='zmax', value='1', format='int')
89 91 opObj12.addParameter(name='figpath', value = pathFigure, format='str')
90 92
91 93 #---------------------------------------------------------------------------------------------------
92 94 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
93 95 opObj20 = procUnitConfObj2.addOperation(name='GetLags')
94 96
95 97 opObj21 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
96 98 opObj21.addParameter(name='technique', value='SA', format='str')
97 99 # opObj21.addParameter(name='correctFactor', value='-1', format='float')
98 100 opObj21.addParameter(name='positionX', value='36,0,36,0', format='floatlist')
99 101 opObj21.addParameter(name='positionY', value='36,0,0,36', format='floatlist')
100 102 opObj21.addParameter(name='azimuth', value='51.06', format='float')
101 103 opObj21.addParameter(name='crosspairsList', value='(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList')#COrregir
102 104 #
103 105 opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
104 106 opObj22.addParameter(name='id', value='4', format='int')
105 107 opObj22.addParameter(name='wintitle', value='Wind Profiler', format='str')
106 108 opObj22.addParameter(name='save', value='1', format='bool')
107 109 opObj22.addParameter(name='figpath', value = pathFigure, format='str')
108 110 opObj22.addParameter(name='zmin', value='-15', format='int')
109 111 opObj22.addParameter(name='zmax', value='15', format='int')
110 112 opObj22.addParameter(name='zmin_ver', value='-80', format='float')
111 113 opObj22.addParameter(name='zmax_ver', value='80', format='float')
112 114 opObj22.addParameter(name='SNRmin', value='-20', format='int')
113 115 opObj22.addParameter(name='SNRmax', value='40', format='int')
114 116 opObj22.addParameter(name='SNRthresh', value='-3.5', format='float')
115 117 opObj22.addParameter(name='xmin', value=xmin, format='float')
116 118 opObj22.addParameter(name='xmax', value=xmax, format='float')
117 119 # #-----------------------------------------------------------------------------------
118 120 #
119 121 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId())
120 122 # procUnitConfObj2.addParameter(name='nFFTPoints', value='128', format='int')
121 123 # procUnitConfObj2.addParameter(name='nProfiles', value='128', format='int')
122 124 # procUnitConfObj2.addParameter(name='pairsList', value='(0,0),(0,1),(2,1)', format='pairsList')
123 125 #
124 126 # opObj22 = procUnitConfObj2.addOperation(name='SpectraPlot', optype='other')
125 127 # opObj22.addParameter(name='id', value='5', format='int')
126 128 # opObj22.addParameter(name='wintitle', value='Spectra Plot', format='str')
127 129 # opObj22.addParameter(name='save', value='1', format='bool')
128 130 # opObj22.addParameter(name='figpath', value = pathFigure, format='str')
129 131
130 132 #-----------------------------------------------------------------------------------
131 133
132 134 print "Escribiendo el archivo XML"
133 135 controllerObj.writeXml(filename)
134 136 print "Leyendo el archivo XML"
135 137 controllerObj.readXml(filename)
136 138
137 139 controllerObj.createObjects()
138 140 controllerObj.connectObjects()
139 141 controllerObj.run() No newline at end of file
General Comments 0
You need to be logged in to leave comments. Login now