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