##// END OF EJS Templates
schain
RSS
to be run with jroIO_kamisr.py
Karim Valdez -
r666:4bdc5a5e1251
parent child
Show More
Show file before | Show file after | Hide comments
@@ -0,0 +1,222
1 #! /usr/bin/python
2 #! /usr/bin/env python
3
4 import os, sys
5 import time
6
7 path = os.path.dirname(os.getcwd())
8 path = os.path.dirname(path)
9 sys.path.insert(0, path)
10
11 from schainpy.controller import Project
12
13 desc = "AMISR Experiment"
14
15 filename = "amisr_reader.xml"
16
17 controllerObj = Project()
18
19 controllerObj.setup(id = '191', name='eej_proc', description=desc)
20
21
22 # path = os.path.join(os.environ['HOME'],'amisr')
23 # path = '/media/signalchain/HD-PXU2/AMISR_JULIA_MODE'
24 # path = '/media/soporte/E9F4-F053/AMISR/Data/NoiseTest/EEJ'
25 # path = '/media/soporte/E9F4-F053/AMISR/Data/EEJ'
26 # path = '/mnt/data_amisr'
27 # path = '/media/soporte/E9F4-F053/AMISR/Data/winds'
28 path = '/mnt/data_amisr'
29
30 #path = '/media/soporte/AMISR_104'
31 #figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/esf')
32 #figpath = '/media/soporte/E9F4-F053/AMISR/Data/JULIA/ESF'
33 figpath = '/home/soporte/Data/winds'
34 remotefolder = "/home/wmaster/graficos"
35 #path = '/media/soporte/AMISR_104'
36 #figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/eej')
37 #figpath = '/media/soporte/E9F4-F053/AMISR/Data/winds/plots'
38
39 xmin = '08'
40 xmax = '18'
41 dbmin = '50' #'60'#'55' #'40' #noise esf eej
42 dbmax = '80' #'70' #'55'
43
44 ippFactor = '5'
45
46 code = '1,-1,-1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,1,-1,1,1,-1,1'
47 nCode = '1'
48 nBaud = '28'
49
50
51 today = time.strftime("%Y/%m/%d")
52
53
54 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
55 path=path,
56 startDate='2015/11/12', #'2014/10/07',
57 endDate='2015/11/12', #'2014/10/07',
58 startTime='00:00:00',#'07:00:00',
59 endTime='23:59:59',#'15:00:00',
60 walk=0,
61 # code = code,
62 # nCode = nCode,
63 # nBaud = nBaud,
64 timezone='lt',
65 online=1)
66
67 #AMISR Processing Unit
68
69 #Voltage Processing Unit
70 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
71 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
72 opObj10.addParameter(name='frequency', value='445e6', format='float')
73
74
75
76 opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
77 opObj12.addParameter(name='minHei', value='0', format='float')
78 opObj12.addParameter(name='maxHei', value='10', format='float')
79
80
81
82 #Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
83 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
84 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=64, format='int')
85 procUnitConfObjSpectraBeam0.addParameter(name='ippFactor', value=ippFactor, format='int')
86 #
87 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
88 opObj11.addParameter(name='n', value='64', format='int')
89 #opObj11.addParameter(name='timeInterval', value='30', format='float')
90
91
92 # # #RemoveDc
93 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
94 #
95 #Noise Estimation
96 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
97 opObj11.addParameter(name='minHei', value='5', format='float')
98 opObj11.addParameter(name='maxHei', value='9', format='float')
99
100
101 #SpectraPlot
102 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='other')
103 opObj11.addParameter(name='id', value='1', format='int')
104 opObj11.addParameter(name='wintitle', value='Winds AMISR', format='str')
105 opObj11.addParameter(name='zmin', value=dbmin, format='int')
106 opObj11.addParameter(name='zmax', value=dbmax, format='int')
107 opObj11.addParameter(name='save', value='1', format='bool')
108 opObj11.addParameter(name='figpath', value = figpath, format='str')
109 opObj11.addParameter(name='ftp', value='1', format='int')
110 opObj11.addParameter(name='wr_period', value='2', format='int')
111 opObj11.addParameter(name='exp_code', value='21', format='int')
112 opObj11.addParameter(name='sub_exp_code', value='3', format='int')
113 opObj11.addParameter(name='ftp_wei', value='0', format='int')
114 opObj11.addParameter(name='plot_pos', value='0', format='int')
115
116 # #RTIPlot
117 # #title0 = 'RTI AMISR Beam 0'
118 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='other')
119 opObj11.addParameter(name='id', value='2', format='int')
120 opObj11.addParameter(name='wintitle', value='Winds AMISR', format='str')
121 opObj11.addParameter(name='showprofile', value='0', format='int')
122 opObj11.addParameter(name='xmin', value=xmin, format='float')
123 opObj11.addParameter(name='xmax', value=xmax, format='float')
124 opObj11.addParameter(name='zmin', value=dbmin, format='int')
125 opObj11.addParameter(name='zmax', value=dbmax, format='int')
126 opObj11.addParameter(name='save', value='1', format='bool')
127 opObj11.addParameter(name='figpath', value = figpath, format='str')
128 opObj11.addParameter(name='ftp', value='1', format='int')
129 opObj11.addParameter(name='wr_period', value='2', format='int')
130 opObj11.addParameter(name='exp_code', value='21', format='int')
131 opObj11.addParameter(name='sub_exp_code', value='3', format='int')
132 opObj11.addParameter(name='ftp_wei', value='0', format='int')
133 opObj11.addParameter(name='plot_pos', value='0', format='int')
134 #
135 #
136 # #Noise
137 #title0 = 'RTI AMISR Beam 0'
138 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='Noise', optype='other')
139 opObj11.addParameter(name='id', value='3', format='int')
140 # opObj11.addParameter(name='wintitle', value=title0, format='str')
141 opObj11.addParameter(name='showprofile', value='0', format='int')
142 opObj11.addParameter(name='xmin', value=xmin, format='float')
143 opObj11.addParameter(name='xmax', value=xmax, format='float')
144 opObj11.addParameter(name='ymin', value=dbmin, format='int')
145 opObj11.addParameter(name='ymax', value=dbmax, format='int')
146 opObj11.addParameter(name='save', value='1', format='bool')
147 opObj11.addParameter(name='figpath', value = figpath, format='str')
148 opObj11.addParameter(name='ftp', value='1', format='int')
149 opObj11.addParameter(name='wr_period', value='2', format='int')
150 opObj11.addParameter(name='exp_code', value='21', format='int')
151 opObj11.addParameter(name='sub_exp_code', value='3', format='int')
152 opObj11.addParameter(name='ftp_wei', value='0', format='int')
153 opObj11.addParameter(name='plot_pos', value='0', format='int')
154
155 #For saving Pdata (doesn't work with amisr data yet!)
156 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraWriter', optype='other')
157 # opObj11.addParameter(name='path', value=figpath)
158 # opObj11.addParameter(name='blocksPerFile', value='100', format='int')
159 # opObj11.addParameter(name='datatype', value="4", format="int") #size of data to be saved
160
161 # opObj22.addParameter(name='azimuth', value='0,-90,0,90,180', format='floatlist')
162 # opObj22.addParameter(name='elevation', value='75.6,75.6,90,75.60,75.6', format='floatlist')
163
164 #Parameters Process
165 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObjSpectraBeam0.getId())
166 opObj20 = procUnitConfObj2.addOperation(name='GetMoments')
167
168 opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
169 opObj22.addParameter(name='technique', value='DBS', format='str')
170 opObj22.addParameter(name='correctAzimuth', value='51.06', format='float')
171 opObj22.addParameter(name='correctFactor', value='-1', format='float')
172 opObj22.addParameter(name='azimuth', value='0,-90,0,90,180', format='floatlist')
173 opObj22.addParameter(name='elevation', value='75.6,75.6,90,75.60,75.60', format='floatlist')
174
175 #WindProfilerPlot
176 opObj23 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
177 opObj23.addParameter(name='id', value='4', format='int')
178 opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
179 opObj23.addParameter(name='save', value='1', format='bool')
180 opObj23.addParameter(name='figpath', value = figpath, format='str')
181 opObj23.addParameter(name='zmin', value='-20', format='int')
182 opObj23.addParameter(name='zmax', value='20', format='int')
183 opObj23.addParameter(name='zmin_ver', value='-100', format='float')
184 opObj23.addParameter(name='zmax_ver', value='100', format='float')
185 opObj23.addParameter(name='SNRmin', value='-20', format='int')
186 opObj23.addParameter(name='SNRmax', value='30', format='int')
187 opObj23.addParameter(name='SNRthresh', value='-50', format='float')
188 opObj23.addParameter(name='xmin', value=xmin, format='float')
189 opObj23.addParameter(name='xmax', value=xmax, format='float')
190 opObj23.addParameter(name='ftp', value='1', format='int')
191 opObj23.addParameter(name='wr_period', value='2', format='int')
192 opObj23.addParameter(name='exp_code', value='21', format='int')
193 opObj23.addParameter(name='sub_exp_code', value='3', format='int')
194 opObj23.addParameter(name='ftp_wei', value='0', format='int')
195 opObj23.addParameter(name='plot_pos', value='0', format='int')
196
197 #---------------------------------------------------------------------------------------------
198
199 procUnitConfObj2 = controllerObj.addProcUnit(name='SendToServer')
200 procUnitConfObj2.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
201 procUnitConfObj2.addParameter(name='username', value='wmaster', format='str')
202 procUnitConfObj2.addParameter(name='password', value='mst2010vhf', format='str')
203 procUnitConfObj2.addParameter(name='localfolder', value=figpath, format='str')
204 procUnitConfObj2.addParameter(name='remotefolder', value=remotefolder, format='str')
205 procUnitConfObj2.addParameter(name='ext', value='.png', format='str')
206 procUnitConfObj2.addParameter(name='period', value=60, format='int')
207 procUnitConfObj2.addParameter(name='protocol', value='ftp', format='str')
208 #-----------------------------------------------------------------------------------------------
209
210
211 # print "Escribiendo el archivo XML"
212 # controllerObj.writeXml(path +'/'+filename)
213 # print "Leyendo el archivo XML"
214 # controllerObj.readXml(path +'/'+filename)
215
216 controllerObj.createObjects()
217 controllerObj.connectObjects()
218 controllerObj.run()
219
220 #21 3 pm
221
222
Show file before | Show file after | Hide comments
@@ -0,0 +1,176
1 #! /usr/bin/python
2 #! /usr/bin/env python
3
4 import os, sys
5 import time
6
7 path = os.path.dirname(os.getcwd())
8 path = os.path.dirname(path)
9 sys.path.insert(0, path)
10
11 from schainpy.controller import Project
12
13 desc = "AMISR Experiment"
14
15 filename = "amisr_reader.xml"
16
17 controllerObj = Project()
18
19 controllerObj.setup(id = '191', name='eej_proc', description=desc)
20
21
22 # path = os.path.join(os.environ['HOME'],'amisr')
23 # path = '/media/signalchain/HD-PXU2/AMISR_JULIA_MODE'
24 # path = '/media/soporte/E9F4-F053/AMISR/Data/NoiseTest/EEJ'
25 # path = '/media/soporte/E9F4-F053/AMISR/Data/EEJ'
26 # path = '/mnt/data_amisr'
27 # path = '/media/soporte/E9F4-F053/AMISR/Data/winds'
28 path = '/mnt/data_amisr'
29 #path = '/media/soporte/AMISR_104'
30 #figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/eej')
31 figpath = '/home/soporte/Data/winds/plots'
32 remotefolder = "/home/wmaster/graficos"
33
34 xmin = '08'
35 xmax = '18'
36 dbmin = '50' #'60'#'55' #'40' #noise esf eej
37 dbmax = '80' #'70' #'55'
38
39 #to consider that each real IPP is actually original_IPP * nchannels
40 ippFactor = '5'
41
42 # code = '1,-1,-1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,1,-1,1,1,-1,1'
43 # nCode = '1'
44 # nBaud = '28'
45
46
47 today = time.strftime("%Y/%m/%d")
48
49
50 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
51 path=path,
52 startDate='2016/10/28', #'2014/10/07',
53 endDate='2014/10/28', #'2014/10/07',
54 startTime='00:00:00',#'07:00:00',
55 endTime='23:59:59',#'15:00:00',
56 walk=0,
57 # code = code,
58 # nCode = nCode,
59 # nBaud = nBaud,
60 timezone='lt',
61 online=1)
62
63 #AMISR Processing Unit
64
65 #Voltage Processing Unit
66 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
67 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
68 opObj10.addParameter(name='frequency', value='445e6', format='float')
69
70
71
72 opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
73 opObj12.addParameter(name='minHei', value='0', format='float')
74 opObj12.addParameter(name='maxHei', value='10', format='float')
75
76
77
78 #Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
79 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
80 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=32, format='int')
81 procUnitConfObjSpectraBeam0.addParameter(name='ippFactor', value=ippFactor, format='int')
82 #
83 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
84 opObj11.addParameter(name='n', value='16', format='int')
85 #opObj11.addParameter(name='timeInterval', value='30', format='float')
86
87
88 # # #RemoveDc
89 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
90 #
91 #Noise Estimation
92 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
93 opObj11.addParameter(name='minHei', value='5', format='float')
94 opObj11.addParameter(name='maxHei', value='9', format='float')
95
96
97 #SpectraPlot
98 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='other')
99 opObj11.addParameter(name='id', value='1', format='int')
100 opObj11.addParameter(name='wintitle', value='Winds AMISR', format='str')
101 opObj11.addParameter(name='zmin', value=dbmin, format='int')
102 opObj11.addParameter(name='zmax', value=dbmax, format='int')
103 opObj11.addParameter(name='save', value='1', format='bool')
104 opObj11.addParameter(name='figpath', value = figpath, format='str')
105 opObj11.addParameter(name='ftp', value='1', format='int')
106 opObj11.addParameter(name='wr_period', value='2', format='int')
107 opObj11.addParameter(name='exp_code', value='21', format='int')
108 opObj11.addParameter(name='sub_exp_code', value='3', format='int')
109 opObj11.addParameter(name='ftp_wei', value='0', format='int')
110 opObj11.addParameter(name='plot_pos', value='0', format='int')
111
112 # #RTIPlot
113 # #title0 = 'RTI AMISR Beam 0'
114 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='other')
115 opObj11.addParameter(name='id', value='2', format='int')
116 opObj11.addParameter(name='wintitle', value='Winds AMISR', format='str')
117 opObj11.addParameter(name='showprofile', value='0', format='int')
118 opObj11.addParameter(name='xmin', value=xmin, format='float')
119 opObj11.addParameter(name='xmax', value=xmax, format='float')
120 opObj11.addParameter(name='zmin', value=dbmin, format='int')
121 opObj11.addParameter(name='zmax', value=dbmax, format='int')
122 opObj11.addParameter(name='save', value='1', format='bool')
123 opObj11.addParameter(name='figpath', value = figpath, format='str')
124 opObj11.addParameter(name='ftp', value='1', format='int')
125 opObj11.addParameter(name='wr_period', value='2', format='int')
126 opObj11.addParameter(name='exp_code', value='21', format='int')
127 opObj11.addParameter(name='sub_exp_code', value='3', format='int')
128 opObj11.addParameter(name='ftp_wei', value='0', format='int')
129 opObj11.addParameter(name='plot_pos', value='0', format='int')
130 #
131 #
132 # #Noise
133 #title0 = 'RTI AMISR Beam 0'
134 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='Noise', optype='other')
135 opObj11.addParameter(name='id', value='3', format='int')
136 # opObj11.addParameter(name='wintitle', value=title0, format='str')
137 opObj11.addParameter(name='showprofile', value='0', format='int')
138 opObj11.addParameter(name='xmin', value=xmin, format='float')
139 opObj11.addParameter(name='xmax', value=xmax, format='float')
140 opObj11.addParameter(name='ymin', value=dbmin, format='int')
141 opObj11.addParameter(name='ymax', value=dbmax, format='int')
142 opObj11.addParameter(name='save', value='1', format='bool')
143 opObj11.addParameter(name='figpath', value = figpath, format='str')
144
145
146 #For saving Pdata (doesn't work with amisr data yet!)
147 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraWriter', optype='other')
148 # opObj11.addParameter(name='path', value=figpath)
149 # opObj11.addParameter(name='blocksPerFile', value='100', format='int')
150 # opObj11.addParameter(name='datatype', value="4", format="int") #size of data to be saved
151
152
153 # procUnitConfObj2 = controllerObj.addProcUnit(name='SendToServer')
154 # procUnitConfObj2.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
155 # procUnitConfObj2.addParameter(name='username', value='wmaster', format='str')
156 # procUnitConfObj2.addParameter(name='password', value='mst2010vhf', format='str')
157 # procUnitConfObj2.addParameter(name='localfolder', value=pathFigure, format='str')
158 # procUnitConfObj2.addParameter(name='remotefolder', value=remotefolder, format='str')
159 # procUnitConfObj2.addParameter(name='ext', value='.png', format='str')
160 # procUnitConfObj2.addParameter(name='period', value=5, format='int')
161 # procUnitConfObj2.addParameter(name='protocol', value='ftp', format='str')
162 #-----------------------------------------------------------------------------------------------
163
164
165 # print "Escribiendo el archivo XML"
166 # controllerObj.writeXml(path +'/'+filename)
167 # print "Leyendo el archivo XML"
168 # controllerObj.readXml(path +'/'+filename)
169
170 controllerObj.createObjects()
171 controllerObj.connectObjects()
172 controllerObj.run()
173
174 #21 3 pm
175
176
General Comments 0
You need to be logged in to leave comments. Login now