##// END OF EJS Templates
schain
RSS
f_
joabAM -
r1374:7243e1f8c2ee
parent child
Show More
Side by Side Unified
Context file:
r1374:7243e1f8c2ee Loading diff...:
Show file before | Show file after | Hide comments
@@ -132,13 +132,14 class AMISRReader(ProcessingUnit):
132
132
133 def readAMISRHeader(self,fp):
133 def readAMISRHeader(self,fp):
134
134
135 if self.isConfig :
135 if self.isConfig and (not self.flagNoMoreFiles):
136 newShape = fp.get('Raw11/Data/Samples/Data').shape
136 newShape = fp.get('Raw11/Data/Samples/Data').shape[1:]
137 if self.dataShape != newShape:
137 if self.dataShape != newShape and newShape != None:
138 print("\nERROR, NEW FILE IS A DIFFERENT EXPERIMENT \n")
138 print("\nNEW FILE HAS A DIFFERENT SHAPE")
139 print(self.dataShape,newShape,"\n")
139 return 0
140 return 0
140 else:
141 else:
141 self.dataShape = fp.get('Raw11/Data/Samples/Data').shape
142 self.dataShape = fp.get('Raw11/Data/Samples/Data').shape[1:]
142
143
143
144
144 header = 'Raw11/Data/RadacHeader'
145 header = 'Raw11/Data/RadacHeader'
Show file before | Show file after | Hide comments
@@ -15,24 +15,25 def main():
15 desc = "AMISR EEJ Experiment"
15 desc = "AMISR EEJ Experiment"
16 filename = "amisr_reader.xml"
16 filename = "amisr_reader.xml"
17 xmin = '07'
17 xmin = '07'
18 xmax = '18' #-> el plot genera +1 en la hora, es decir aparece 18 como máximo
18 xmax = '18'
19 xrange = int(xmax)-int(xmin)
19 ymin = '0'
20 ymin = '0'
20 ymax = '300'
21 ymax = '300'
21 dbmin = '45' #'60'#'55' #'40' #noise esf eej
22 dbmin = '45' #'60'#'55' #'40' #noise esf eej
22 dbmax = '65' #'70' #'55'
23 dbmax = '65' #'70' #'55'
23 showSPC = '1' #view plot Spectra
24 showSPC = '0' #view plot Spectra
24 showRTI = '1' #view plot RTI
25 showRTI = '0' #view plot RTI
25 showNOISE = '1' #view plot NOISE
26 showNOISE = '0' #view plot NOISE
26 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
27 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
27 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'
28 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'
28 nCode = '1'
29 nCode = '1'
29 nBaud = '28'
30 nBaud = '28'
30 nosamp = '1' # oversample for EEJ
31 nosamp = '1' # oversample for EEJ
31 today = time.strftime("%Y/%m/%d")
32 today = time.strftime("%Y/%m/%d")
32 #startDate=today
33 startDate=today
33 #endDate=today
34 endDate=today
34 startDate='2021/07/11'
35 #startDate='2021/07/11'
35 endDate='2021/07/11'
36 #endDate='2021/07/11'
36 #inPath= '/home/soporte/dataAMISR_test/'
37 #inPath= '/home/soporte/dataAMISR_test/'
37 inPath= '/home/soporte/dataAMISR/'
38 inPath= '/home/soporte/dataAMISR/'
38 inPath= '/media/soporte/UARS_4T_D02/AMISR_DATA/2021/'
39 inPath= '/media/soporte/UARS_4T_D02/AMISR_DATA/2021/'
@@ -63,8 +64,8 def main():
63 path=inPath,
64 path=inPath,
64 startDate=startDate,#startDate, #'2014/10/07',
65 startDate=startDate,#startDate, #'2014/10/07',
65 endDate=endDate, #endDate '2014/10/07',
66 endDate=endDate, #endDate '2014/10/07',
66 startTime='17:55:30',#'07:00:00',
67 startTime='07:00:00',#'07:00:00',
67 endTime='19:00:00',#'15:00:00',
68 endTime='18:00:00',#'15:00:00',
68 walk=1,
69 walk=1,
69 code = code,
70 code = code,
70 nCode = nCode,
71 nCode = nCode,
@@ -143,10 +144,10 def main():
143 opObj14.addParameter(name='wintitle', value='title0', format='str')
144 opObj14.addParameter(name='wintitle', value='title0', format='str')
144 opObj14.addParameter(name='showprofile', value='0', format='int')
145 opObj14.addParameter(name='showprofile', value='0', format='int')
145 opObj14.addParameter(name='tmin', value=xmin, format='int')
146 opObj14.addParameter(name='tmin', value=xmin, format='int')
146 opObj14.addParameter(name='tmax', value=xmax, format='int')
147 opObj14.addParameter(name='xrange', value=xrange, format='int')
147 opObj14.addParameter(name='ymin', value=dbmin, format='int')
148 opObj14.addParameter(name='ymin', value=dbmin, format='int')
148 opObj14.addParameter(name='ymax', value=dbmax, format='int')
149 opObj14.addParameter(name='ymax', value=dbmax, format='int')
149 opObj14.addParameter(name='save', value=outPath, format='str')
150 opObj14.addParameter(name='save', value=outPath+'/plots', format='str')
150 opObj14.addParameter(name='localtime', value=localtime,format='int')
151 opObj14.addParameter(name='localtime', value=localtime,format='int')
151 opObj14.addParameter(name='show', value = showNOISE, format='int')
152 opObj14.addParameter(name='show', value = showNOISE, format='int')
152
153
@@ -156,7 +157,7 def main():
156 opObj15.addParameter(name='localtime', value=localtime,format='int')
157 opObj15.addParameter(name='localtime', value=localtime,format='int')
157 opObj15.addParameter(name='wintitle', value='RTI', format='str')
158 opObj15.addParameter(name='wintitle', value='RTI', format='str')
158 opObj15.addParameter(name='tmin', value=xmin, format='int')
159 opObj15.addParameter(name='tmin', value=xmin, format='int')
159 opObj15.addParameter(name='tmax', value=xmax, format='int') #max value =23
160 opObj15.addParameter(name='xrange', value=xrange, format='int')
160 opObj15.addParameter(name='ymin', value=ymin, format='int')
161 opObj15.addParameter(name='ymin', value=ymin, format='int')
161 opObj15.addParameter(name='zmin', value=dbmin, format='int')
162 opObj15.addParameter(name='zmin', value=dbmin, format='int')
162 opObj15.addParameter(name='zmax', value=dbmax, format='int')
163 opObj15.addParameter(name='zmax', value=dbmax, format='int')
@@ -169,19 +170,19 def main():
169
170
170 ##.......................................................................................
171 ##.......................................................................................
171 ##.......................................................................................
172 ##.......................................................................................
172
173 #
173 procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
174 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
174 opObj16 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
175 # opObj16 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
175
176
176
177
177 #Using HDFWriter::::
178 #Using HDFWriter::::
178 #.......................................................................................
179 #.......................................................................................
179 #.......................................................................................
180 #.......................................................................................
180 opObj17 = procUnitConfObj2.addOperation(name='HDFWriter', optype='external')
181 # opObj17 = procUnitConfObj2.addOperation(name='HDFWriter', optype='external')
181 opObj17.addParameter(name='path', value=outPath)
182 # opObj17.addParameter(name='path', value=outPath)
182 opObj17.addParameter(name='blocksPerFile', value='10', format='int')
183 # opObj17.addParameter(name='blocksPerFile', value='10', format='int')
183 opObj17.addParameter(name='metadataList',value='type,inputUnit,heightList',format='list')
184 # opObj17.addParameter(name='metadataList',value='type,inputUnit,heightList',format='list')
184 opObj17.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
185 # opObj17.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
185
186
186
187
187
188
Show file before | Show file after | Hide comments
@@ -15,15 +15,15 def main():
15 desc = "AMISR EEJ Experiment"
15 desc = "AMISR EEJ Experiment"
16 filename = "amisr_reader.xml"
16 filename = "amisr_reader.xml"
17 xmin = '07'
17 xmin = '07'
18 xmax = '18' #-> el plot genera +1 en la hora, es decir aparece 18 como máximo
18 xmax = '18' #
19 ymin = '0'
19 xrange = int(xmax)-int(xmin)
20 ymax = '300'
20
21 dbmin = '45' #'60'#'55' #'40' #noise esf eej
21 dbmin = '45' #'60'#'55' #'40' #noise esf eej
22 dbmax = '65' #'70' #'55'
22 dbmax = '65' #'70' #'55'
23 showSPC = '1' #view plot Spectra
23 showSPC = '1' #view plot Spectra
24 showRTI = '1' #view plot RTI
24 showRTI = '1' #view plot RTI
25 showNOISE = '1' #view plot NOISE
25 showNOISE = '1' #view plot NOISE
26 localtime='0' #para ajustar el horario en las gráficas '0' para dejar en utc
26 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
27 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'
27 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'
28 nCode = '1'
28 nCode = '1'
29 nBaud = '28'
29 nBaud = '28'
@@ -62,20 +62,20 def main():
62 startDate=startDate,
62 startDate=startDate,
63 endDate=endDate,
63 endDate=endDate,
64 startTime='07:01:30',#'07:00:00',
64 startTime='07:01:30',#'07:00:00',
65 endTime='18:55:00',#'15:00:00',
65 endTime='17:55:00',#'15:00:00',
66 walk=0,
66 walk=0,
67 code = code,
67 code = code,
68 nCode = nCode,
68 nCode = nCode,
69 nBaud = nBaud,
69 nBaud = nBaud,
70 timezone='ut',
70 timezone='lt',
71 online=1)
71 online=1)
72
72
73 #AMISR Processing Unit
73 #AMISR Processing Unit
74 ##.......................................................................................
74 ##.......................................................................................
75 ##.......................................................................................
75 ##.......................................................................................
76 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
76 procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
77 opObj10 = procUnitConfObj0.addOperation(name='setRadarFrequency')
77 opObj10 = procUnitConfObj0.addOperation(name='setAttribute')
78 opObj10.addParameter(name='frequency', value='445e6', format='float')
78 opObj10.addParameter(name='frequency', value='445.0e6')
79
79
80
80
81 opObj01 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
81 opObj01 = procUnitConfObj0.addOperation(name='Decoder', optype='other')
@@ -85,8 +85,6 def main():
85 opObj01.addParameter(name='osamp', value=nosamp, format='int')
85 opObj01.addParameter(name='osamp', value=nosamp, format='int')
86
86
87
87
88 # opObj02 = procUnitConfObj0.addOperation(name='CohInt', optype='other')
89 # opObj02.addParameter(name='n', value='5', format='int')
90
88
91
89
92
90
@@ -100,14 +98,9 def main():
100 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
98 opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
101 opObj11.addParameter(name='n', value='150', format='int') #300 normal value
99 opObj11.addParameter(name='n', value='150', format='int') #300 normal value
102
100
103 # ## Remove DC signal
104 # opObj11 = procUnitConfObj1.addOperation(name='removeDC')
105 ##.......................................................................................
101 ##.......................................................................................
106 ##.......................................................................................
102 ##.......................................................................................
107
103
108 opObj13 = procUnitConfObj1.addOperation(name='getNoise' , optype ='self')
109 opObj13.addParameter(name='minHei', value='100', format='float')
110 opObj13.addParameter(name='maxHei', value='280', format='float')
111
104
112 #
105 #
113 opObj12 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='external')
106 opObj12 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='external')
@@ -124,38 +117,25 def main():
124 opObj12.addParameter(name='show', value = showSPC, format='int')
117 opObj12.addParameter(name='show', value = showSPC, format='int')
125
118
126
119
127
120 opObj14 = procUnitConfObj1.addOperation(name='NoisePlot', optype='external')
128
121 opObj14.addParameter(name='id', value='3', format='int')
129
122 opObj14.addParameter(name='wintitle', value='title0', format='str')
130
123 opObj14.addParameter(name='showprofile', value='0', format='int')
131 ##Generate *.pdata from AMISR data
124 opObj14.addParameter(name='tmin', value=xmin, format='int')
132 ##.......................................................................................
125 opObj14.addParameter(name='xrange', value=xrange, format='int')
133 ##.......................................................................................
126 opObj14.addParameter(name='ymin', value=dbmin, format='int')
134 # opObj13 = procUnitConfObj1.addOperation(name='SpectraWriter', optype='external')
127 opObj14.addParameter(name='ymax', value=dbmax, format='int')
135 # opObj13.addParameter(name='path', value=outPath)
128 opObj14.addParameter(name='save', value=outPath+'/plots', format='str')
136 # opObj13.addParameter(name='blocksPerFile', value='10', format='int')
129 opObj14.addParameter(name='localtime', value=localtime,format='int')
137
130 opObj14.addParameter(name='show', value = showNOISE, format='int')
138
139 # opObj14 = procUnitConfObj1.addOperation(name='NoisePlot', optype='external')
140 # opObj14.addParameter(name='id', value='3', format='int')
141 # opObj14.addParameter(name='wintitle', value='title0', format='str')
142 # opObj14.addParameter(name='showprofile', value='0', format='int')
143 # opObj14.addParameter(name='xmin', value=xmin, format='int')
144 # opObj14.addParameter(name='xmax', value=xmax, format='int')
145 # opObj14.addParameter(name='ymin', value=dbmin, format='int')
146 # opObj14.addParameter(name='ymax', value=dbmax, format='int')
147 # opObj14.addParameter(name='save', value=outPath+'/plots', format='str')
148 # opObj14.addParameter(name='localtime', value=localtime,format='int')
149 # opObj14.addParameter(name='show', value = showNOISE, format='int')
150
131
151
132
152 opObj15 = procUnitConfObj1.addOperation(name='RTIPlot', optype='external')
133 opObj15 = procUnitConfObj1.addOperation(name='RTIPlot', optype='external')
153 opObj15.addParameter(name='id', value='2', format='int')
134 opObj15.addParameter(name='id', value='2', format='int')
154 opObj15.addParameter(name='localtime', value=localtime,format='int')
135 opObj15.addParameter(name='localtime', value=localtime,format='int')
155 opObj15.addParameter(name='wintitle', value='RTI', format='str')
136 opObj15.addParameter(name='wintitle', value='RTI', format='str')
156 opObj15.addParameter(name='xmin', value=xmin, format='int')
137 opObj15.addParameter(name='tmin', value=xmin, format='int')
157 opObj15.addParameter(name='xmax', value=xmax, format='int') #max value =23
138 opObj15.addParameter(name='xrange', value=xrange, format='int')
158 opObj15.addParameter(name='ymin', value=ymin, format='int')
159 opObj15.addParameter(name='zmin', value=dbmin, format='int')
139 opObj15.addParameter(name='zmin', value=dbmin, format='int')
160 opObj15.addParameter(name='zmax', value=dbmax, format='int')
140 opObj15.addParameter(name='zmax', value=dbmax, format='int')
161 opObj15.addParameter(name='showprofile', value='1', format='int')
141 opObj15.addParameter(name='showprofile', value='1', format='int')
@@ -168,30 +148,9 def main():
168 ##.......................................................................................
148 ##.......................................................................................
169 ##.......................................................................................
149 ##.......................................................................................
170 #
150 #
171 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId())
172 # opObj16 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
173 #
174 #
175 # #Using ParamWriter::::
176 # ##.......................................................................................
177 # ##.......................................................................................
178 # opObj17 = procUnitConfObj2.addOperation(name='ParamWriter', optype='external')
179 # opObj17.addParameter(name='path', value=outPath)
180 # opObj17.addParameter(name='blocksPerFile', value='10', format='int')
181 # opObj17.addParameter(name='metadataList',value='type,inputUnit,heightList',format='list')
182 # opObj17.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
183 # opObj17.addParameter(name='mode',value='1',format='int') #'0' channels, '1' parameters, '3' table (for meteors)
184 # ##opObj17.addParameter(name='setType', value ='anything', format='str')
185 #
186
187 ##.......................................................................................
188 ##.......................................................................................
189 #print("Escribiendo el archivo XML",controllerObj.writeXml(path +'/'+filename))
190
151
191 controllerObj.start()
152 controllerObj.start()
192
153
193 #print("Leyendo el archivo XML",controllerObj.readXml(path +'/'+filename))
194
195 ##.......................................................................................
154 ##.......................................................................................
196 ##.......................................................................................
155 ##.......................................................................................
197
156
Show file before | Show file after | Hide comments
@@ -23,24 +23,25 def main():
23 #path = os.path.join(os.environ['HOME'],'amisr')
23 #path = os.path.join(os.environ['HOME'],'amisr')
24 inPath = '/mnt/data_amisr'
24 inPath = '/mnt/data_amisr'
25 inPath= '/home/soporte/dataAMISR/'
25 inPath= '/home/soporte/dataAMISR/'
26 #path= '/home/soporte/dataAMISR_test/'
26 inPath= '/media/soporte/DATA/dataAMISR/'
27
27 outPath = '/home/soporte/Data/ESF'
28 outPath = '/home/soporte/Data/ESF'
28 remotefolder = "/home/wmaster/graficos"
29 remotefolder = "/home/wmaster/graficos"
29 xmin = '19'
30 xmin = '18'
30 xmax = '07'
31 xmax = '31'
31 xrange = (int(xmax)+24)-int(xmin)
32 xrange = int(xmax)-int(xmin)
32 dbmin = '60' #'60'#'55' #'40' #noise esf eej
33 dbmin = '60' #'60'#'55' #'40' #noise esf eej
33 dbmax = '75' #'70' #'55'
34 dbmax = '75' #'70' #'55'
34 showSPC = '0' #view plot Spectra
35 showSPC = '1' #view plot Spectra
35 showRTI = '0' #view plot RTI
36 showRTI = '0' #view plot RTI
36 showNOISE = '0' #view plot NOISE
37 showNOISE = '1' #view plot NOISE
37 # # 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'
38 # # 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'
38 # # code = '1,1,0,1,1,0,1,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,0,0,0'
39 # # code = '1,1,0,1,1,0,1,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,0,0,0'
39 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'
40 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'
40 nCode = '1'
41 nCode = '1'
41 nBaud = '28'
42 nBaud = '28'
42 nosamp = '2' # oversample
43 nosamp = '2' # oversample
43 localtime='0' #para ajustar el horario en las gráficas '0' para dejar en utc
44 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
44
45
45 str0 = datetime.date.today()
46 str0 = datetime.date.today()
46 str1 = str0 + datetime.timedelta(days=1)
47 str1 = str0 + datetime.timedelta(days=1)
@@ -52,6 +53,8 def main():
52 startDate=yesterday
53 startDate=yesterday
53 endDate=today
54 endDate=today
54
55
56 startDate='2019/12/17'
57 endDate='2019/12/18'
55 ##.......................................................................................
58 ##.......................................................................................
56 ##.......................................................................................
59 ##.......................................................................................
57 l = startDate.split('/') #adding day of the year to outPath
60 l = startDate.split('/') #adding day of the year to outPath
@@ -67,21 +70,21 def main():
67 ##.......................................................................................
70 ##.......................................................................................
68
71
69 controllerObj = Project()
72 controllerObj = Project()
70 controllerObj.setup(id = '191', name='esf_proc', description=desc)
73 controllerObj.setup(id = '20', name='esf_proc', description=desc)
71 ##.......................................................................................
74 ##.......................................................................................
72 ##.......................................................................................
75 ##.......................................................................................
73
76
74 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
77 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
75 path=inPath,
78 path=inPath,
76 startDate='2019/12/17',#startDate,#'2016/07/12',
79 startDate=startDate,#startDate,#'2016/07/12',
77 endDate='2019/12/18',#endDate,#'2016/07/13',
80 endDate=endDate,#endDate,#'2016/07/13',
78 startTime='19:01:30',#'07:00:00',
81 startTime='23:40:30',#'07:00:00',
79 endTime='07:00:00',#'15:00:00',
82 endTime='07:00:00',#'15:00:00',
80 walk=0,
83 walk=0,
81 code = code,
84 code = code,
82 nCode = nCode,
85 nCode = nCode,
83 nBaud = nBaud,
86 nBaud = nBaud,
84 timezone='ut',
87 timezone='lt',
85 online=0)
88 online=0)
86
89
87 #AMISR Processing Unit
90 #AMISR Processing Unit
@@ -90,10 +93,10 def main():
90
93
91 #Voltage Processing Unit
94 #Voltage Processing Unit
92 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
95 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
93 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
96 opObj10 = procUnitConfObjBeam0.addOperation(name='setAttribute')
94 #opObj10.addParameter(name='frequency', value='445e6', format='int') #changed on Dec 3, 15:40h
97 opObj10.addParameter(name='frequency', value='445.09e6')
95
96
98
99 '''
97 #Noise--> no code
100 #Noise--> no code
98
101
99 opObj11 = procUnitConfObjBeam0.addOperation(name='Decoder', optype='other')
102 opObj11 = procUnitConfObjBeam0.addOperation(name='Decoder', optype='other')
@@ -118,15 +121,15 def main():
118 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
121 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
119 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=16, format='int')
122 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=16, format='int')
120 #
123 #
121 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
124 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
122 # opObj11.addParameter(name='n', value='60', format='int')
125 opObj11.addParameter(name='n', value='260', format='int')
123 # opObj11.addParameter(name='timeInterval', value='30', format='int')
126 #opObj11.addParameter(name='timeInterval', value='30', format='int')
124 #
127 #
125 #
128 #
126 # # # #RemoveDc
129 # # # #RemoveDc
127 # #opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
130 #opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
128 #
131 #
129 # ##Noise Estimation
132 # ##Noise Estimation #so reconoce nueva version
130 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise', optype='self')
133 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise', optype='self')
131 # opObj11.addParameter(name='minHei', value='100', format='int')
134 # opObj11.addParameter(name='minHei', value='100', format='int')
132 # opObj11.addParameter(name='maxHei', value='280', format='int')
135 # opObj11.addParameter(name='maxHei', value='280', format='int')
@@ -145,7 +148,6 def main():
145 opObj11.addParameter(name='id', value='1', format='int')
148 opObj11.addParameter(name='id', value='1', format='int')
146 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
149 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
147 opObj11.addParameter(name='showprofile', value='1', format='int')
150 opObj11.addParameter(name='showprofile', value='1', format='int')
148 #opObj11.addParameter(name='zmin', value='38', format='int')
149 opObj11.addParameter(name='zmin', value=dbmin, format='int')
151 opObj11.addParameter(name='zmin', value=dbmin, format='int')
150 opObj11.addParameter(name='zmax', value=dbmax, format='int')
152 opObj11.addParameter(name='zmax', value=dbmax, format='int')
151 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
153 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
@@ -164,7 +166,7 def main():
164 opObj11.addParameter(name='id', value='2', format='int')
166 opObj11.addParameter(name='id', value='2', format='int')
165 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
167 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
166 opObj11.addParameter(name='showprofile', value='0', format='int')
168 opObj11.addParameter(name='showprofile', value='0', format='int')
167 opObj11.addParameter(name='xmin', value=xmin, format='int')
169 opObj11.addParameter(name='tmin', value=xmin, format='int')
168 opObj11.addParameter(name='xrange', value=xrange, format='int')
170 opObj11.addParameter(name='xrange', value=xrange, format='int')
169 opObj11.addParameter(name='zmin', value=dbmin, format='int')
171 opObj11.addParameter(name='zmin', value=dbmin, format='int')
170 opObj11.addParameter(name='zmax', value=dbmax, format='int')
172 opObj11.addParameter(name='zmax', value=dbmax, format='int')
@@ -186,7 +188,7 def main():
186 opObj11.addParameter(name='id', value='3', format='int')
188 opObj11.addParameter(name='id', value='3', format='int')
187 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
189 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
188 opObj11.addParameter(name='showprofile', value='0', format='int')
190 opObj11.addParameter(name='showprofile', value='0', format='int')
189 opObj11.addParameter(name='xmin', value=xmin, format='int')
191 opObj11.addParameter(name='tmin', value=xmin, format='int')
190 opObj11.addParameter(name='xrange', value=xrange,format='int')
192 opObj11.addParameter(name='xrange', value=xrange, format='int')
191 #opObj11.addParameter(name='ymin', value=dbmin, format='int')
193 #opObj11.addParameter(name='ymin', value=dbmin, format='int')
192 #opObj11.addParameter(name='ymax', value=dbmax, format='int')
194 #opObj11.addParameter(name='ymax', value=dbmax, format='int')
@@ -198,11 +200,11 def main():
198 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraWriter', optype='external')
200 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraWriter', optype='external')
199 opObj11.addParameter(name='path', value=outPath)
201 opObj11.addParameter(name='path', value=outPath)
200 opObj11.addParameter(name='blocksPerFile', value='10', format='int')
202 opObj11.addParameter(name='blocksPerFile', value='10', format='int')
201
203 '''
202 ##.......................................................................................
204 ##.......................................................................................
203 ##.......................................................................................
205 ##.......................................................................................
204
206
205 #generate moments
207 # #generate moments
206 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObjSpectraBeam0.getId())
208 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObjSpectraBeam0.getId())
207 # opObj20 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
209 # opObj20 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
208 #
210 #
@@ -213,15 +215,15 def main():
213 # opObj12.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
215 # opObj12.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
214 # opObj12.addParameter(name='mode',value='1',format='int')
216 # opObj12.addParameter(name='mode',value='1',format='int')
215 #
217 #
216
218 #
217 #----
219 # ----
218
220
219
221
220 # print "Escribiendo el archivo XML"
222 # print "Escribiendo el archivo XML"
221 # controllerObj.writeXml(filename)
223 # controllerObj.writeXml(outPath+"/"filename)
222 controllerObj.start()
224 controllerObj.start()
223 # print "Leyendo el archivo XML"
225 # print "Leyendo el archivo XML"
224 # controllerObj.readXml(filename)
226 # controllerObj.readXml(outPath+"/"filename)
225
227
226 if __name__ == '__main__':
228 if __name__ == '__main__':
227 import time
229 import time
Show file before | Show file after | Hide comments
@@ -26,14 +26,13 def main():
26 #inPath= '/home/soporte/dataAMISR_test/'
26 #inPath= '/home/soporte/dataAMISR_test/'
27 outPath = '/home/soporte/Data/ESF'
27 outPath = '/home/soporte/Data/ESF'
28 remotefolder = "/home/wmaster/graficos"
28 remotefolder = "/home/wmaster/graficos"
29 xmin = '16'
29 xmin = '18'
30 xmax = '07'
30 xmax = '07'
31 xrange = (int(xmax)+24)-int(xmin)
32 dbmin = '60' #'60'#'55' #'40' #noise esf eej
31 dbmin = '60' #'60'#'55' #'40' #noise esf eej
33 dbmax = '75' #'70' #'55'
32 dbmax = '75' #'70' #'55'
34 showSPC = '1' #view plot Spectra
33 showSPC = '1' #view plot Spectra
35 showRTI = '1' #view plot RTI
34 showRTI = '1' #view plot RTI
36 showNOISE = '0' #view plot NOISE
35 showNOISE = '1' #view plot NOISE
37 # # 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'
36 # # 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'
38 # # code = '1,1,0,1,1,0,1,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,0,0,0'
37 # # code = '1,1,0,1,1,0,1,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,0,0,0'
39 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'
38 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'
@@ -67,7 +66,7 def main():
67 ##.......................................................................................
66 ##.......................................................................................
68
67
69 controllerObj = Project()
68 controllerObj = Project()
70 controllerObj.setup(id = '191', name='esf_proc', description=desc)
69 controllerObj.setup(id = '21', name='esf_proc', description=desc)
71 ##.......................................................................................
70 ##.......................................................................................
72 ##.......................................................................................
71 ##.......................................................................................
73
72
@@ -126,16 +125,10 def main():
126 # # #RemoveDc
125 # # #RemoveDc
127 #opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
126 #opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
128
127
129 ##Noise Estimation
130 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise', optype='self')
131 opObj11.addParameter(name='minHei', value='100', format='int')
132 opObj11.addParameter(name='maxHei', value='280', format='int')
133 #opObj11.addParameter(name='minHei', value='15', format='int')
134 #opObj11.addParameter(name='maxHei', value='20', format='int')
135
128
136 #---
129 #---
137 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='selectChannels')
130 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='selectChannels')
138 # opObj11.addParameter(name='channelList', value='7,8,9', format='intlist')
131 opObj11.addParameter(name='channelList', value='0,2,4,6,8', format='intlist')
139
132
140 ##.......................................................................................
133 ##.......................................................................................
141 ##.......................................................................................
134 ##.......................................................................................
@@ -145,7 +138,6 def main():
145 opObj11.addParameter(name='id', value='1', format='int')
138 opObj11.addParameter(name='id', value='1', format='int')
146 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
139 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
147 opObj11.addParameter(name='showprofile', value='1', format='int')
140 opObj11.addParameter(name='showprofile', value='1', format='int')
148 #opObj11.addParameter(name='zmin', value='38', format='int')
149 opObj11.addParameter(name='zmin', value=dbmin, format='int')
141 opObj11.addParameter(name='zmin', value=dbmin, format='int')
150 opObj11.addParameter(name='zmax', value=dbmax, format='int')
142 opObj11.addParameter(name='zmax', value=dbmax, format='int')
151 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
143 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
@@ -166,8 +158,8 def main():
166 opObj11.addParameter(name='id', value='2', format='int')
158 opObj11.addParameter(name='id', value='2', format='int')
167 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
159 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
168 opObj11.addParameter(name='showprofile', value='0', format='int')
160 opObj11.addParameter(name='showprofile', value='0', format='int')
169 opObj11.addParameter(name='xmin', value=xmin, format='int')
161 opObj11.addParameter(name='tmin', value=xmin, format='int')
170 opObj11.addParameter(name='xrange', value=xrange,format='int')
162 opObj11.addParameter(name='tmax', value=xmax,format='int')
171 opObj11.addParameter(name='zmin', value=dbmin, format='int')
163 opObj11.addParameter(name='zmin', value=dbmin, format='int')
172 opObj11.addParameter(name='zmax', value=dbmax, format='int')
164 opObj11.addParameter(name='zmax', value=dbmax, format='int')
173 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
165 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
@@ -184,40 +176,18 def main():
184 # # #
176 # # #
185 # #Noise
177 # #Noise
186 #title0 = 'RTI AMISR Beam 0'
178 #title0 = 'RTI AMISR Beam 0'
187 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='NoisePlot', optype='external')
179 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='NoisePlot', optype='external')
188 # opObj11.addParameter(name='id', value='3', format='int')
180 opObj11.addParameter(name='id', value='3', format='int')
189 # opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
181 opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
190 # opObj11.addParameter(name='showprofile', value='0', format='int')
182 opObj11.addParameter(name='showprofile', value='0', format='int')
191 # opObj11.addParameter(name='xmin', value=xmin, format='int')
183 opObj11.addParameter(name='tmin', value=xmin, format='int')
192 # opObj11.addParameter(name='xmin', value=xmin, format='int')
184 opObj11.addParameter(name='tmax', value=xmax, format='int')
193 # opObj11.addParameter(name='xrange', value='16',format='int')
194 # #opObj11.addParameter(name='ymin', value=dbmin, format='int')
185 ##opObj11.addParameter(name='ymin', value=dbmin, format='int')
195 # #opObj11.addParameter(name='ymax', value=dbmax, format='int')
186 ##opObj11.addParameter(name='ymax', value=dbmax, format='int')
196 # opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
187 opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
197 # opObj11.addParameter(name='show', value = showNOISE, format='int')
188 opObj11.addParameter(name='show', value = showNOISE, format='int')
198 # opObj11.addParameter(name='localtime', value=localtime,format='int')
189 opObj11.addParameter(name='localtime', value=localtime,format='int')
199 #
200 # #Generate *.pdata from AMISR data
201 # opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraWriter', optype='external')
202 # opObj11.addParameter(name='path', value=outPath, format='str')
203 # opObj11.addParameter(name='blocksPerFile', value='10', format='int')
204 #
205 # ##.......................................................................................
206 # ##.......................................................................................
207 #
208 # #generate moments
209 # procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObjSpectraBeam0.getId())
210 # opObj20 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
211 #
212 # opObj12 = procUnitConfObj2.addOperation(name='ParamWriter', optype='external')
213 # opObj12.addParameter(name='path', value=outPath)
214 # opObj12.addParameter(name='blocksPerFile', value='10', format='int')
215 # opObj12.addParameter(name='metadataList',value='type,inputUnit,heightList',format='list')
216 # opObj12.addParameter(name='dataList',value='moments,data_SNR,utctime',format='list')
217 # opObj12.addParameter(name='mode',value='1',format='int')
218 #
219
190
220 #----
221
191
222
192
223 # print "Escribiendo el archivo XML"
193 # print "Escribiendo el archivo XML"
Show file before | Show file after | Hide comments
@@ -16,7 +16,7 def main():
16
16
17 controllerObj = Project()
17 controllerObj = Project()
18
18
19 controllerObj.setup(id = '10', name='eej_proc', description=desc)
19 controllerObj.setup(id = '100', name='eej_proc', description=desc)
20
20
21
21
22 #path = os.path.join(os.environ['HOME'],'amisr')
22 #path = os.path.join(os.environ['HOME'],'amisr')
@@ -25,16 +25,16 def main():
25 # path = '/media/soporte/E9F4-F053/AMISR/Data/ESF'
25 # path = '/media/soporte/E9F4-F053/AMISR/Data/ESF'
26 #path = '/mnt/data_amisr'
26 #path = '/mnt/data_amisr'
27
27
28 path = '/home/soporte/Data/EEJ/EEJ2021192'
28 path = '/home/soporte/Data/EEJ/EEJ2021192/'
29 path = '/home/soporte/Data/ESF/ESF2019351/'
29 #figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/esf')
30 #figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/esf')
30 #figpath = '/media/soporte/E9F4-F053/AMISR/Data/JULIA/ESF'
31 #figpath = '/media/soporte/E9F4-F053/AMISR/Data/JULIA/ESF'
31 figpath = '/home/soporte/Data/OutTest/EEJ'
32 figpath = '/home/soporte/Data/OutTest/EEJ'
33 figpath = '/home/soporte/Data/OutTest/ESF'
32
34
35 xmin = '18'
36 xmax = '07' #-> el plot genera +1 en la hora, es decir aparece 18 como máximo
33
37
34 xmin = '07'
35 xmax = '18' #-> el plot genera +1 en la hora, es decir aparece 18 como máximo
36 ymin = '0'
37 ymax = '300'
38 dbmin = '45' #'60'#'55' #'40' #noise esf eej
38 dbmin = '45' #'60'#'55' #'40' #noise esf eej
39 dbmax = '65' #'70' #'55'
39 dbmax = '65' #'70' #'55'
40 showSPC = '1' #view plot Spectra
40 showSPC = '1' #view plot Spectra
@@ -62,10 +62,10 def main():
62
62
63 readUnitConfObj = controllerObj.addReadUnit(datatype='Spectra',
63 readUnitConfObj = controllerObj.addReadUnit(datatype='Spectra',
64 path=path,
64 path=path,
65 startDate='2021/07/11',
65 startDate='2019/12/17',#'2021/07/11',
66 endDate='2021/07/11',
66 endDate='2019/12/18',
67 startTime='07:01:30',#'07:00:00',
67 startTime='00:01:30',#'07:00:00',
68 endTime='18:00:00',#'15:00:00',
68 endTime='07:00:00',#'15:00:00',
69 walk=1,
69 walk=1,
70 timezone='lt',
70 timezone='lt',
71 online=0)
71 online=0)
@@ -94,7 +94,7 def main():
94 op3 = proc_spectra.addOperation(name='SpectraPlot', optype='external')
94 op3 = proc_spectra.addOperation(name='SpectraPlot', optype='external')
95 op3.addParameter(name='id', value='21', format='int')
95 op3.addParameter(name='id', value='21', format='int')
96 op3.addParameter(name='xaxis', value='velocity')
96 op3.addParameter(name='xaxis', value='velocity')
97 op3.addParameter(name='ymax', value=ymax, format='int')
97 #op3.addParameter(name='ymax', value=ymax, format='int')
98 op3.addParameter(name='showprofile', value='1', format='int')
98 op3.addParameter(name='showprofile', value='1', format='int')
99 op3.addParameter(name='wintitle', value='AMISR Beam 0', format='str')
99 op3.addParameter(name='wintitle', value='AMISR Beam 0', format='str')
100 op3.addParameter(name='zmin', value=dbmin, format='int')
100 op3.addParameter(name='zmin', value=dbmin, format='int')
@@ -115,7 +115,7 def main():
115 op4.addParameter(name='wintitle', value='RTI', format='str')
115 op4.addParameter(name='wintitle', value='RTI', format='str')
116 #op4.addParameter(name='xmin', value=xmin, format='int')
116 #op4.addParameter(name='xmin', value=xmin, format='int')
117 op4.addParameter(name='xmax', value=xmax, format='int') #max value =23
117 op4.addParameter(name='xmax', value=xmax, format='int') #max value =23
118 op4.addParameter(name='ymin', value=ymin, format='int')
118 #op4.addParameter(name='ymin', value=ymin, format='int')
119 #op4.addParameter(name='ymax', value=ymax, format='int')
119 #op4.addParameter(name='ymax', value=ymax, format='int')
120 op4.addParameter(name='zmin', value=dbmin, format='int')
120 op4.addParameter(name='zmin', value=dbmin, format='int')
121 op4.addParameter(name='zmax', value=dbmax, format='int')
121 op4.addParameter(name='zmax', value=dbmax, format='int')
Show file before | Show file after | Hide comments
@@ -34,11 +34,12 def main():
34 #path = '/media/soporte/AMISR_104'
34 #path = '/media/soporte/AMISR_104'
35 #outPath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/esf')
35 #outPath = os.path.join(os.environ['HOME'],'Pictures/amisr/test/proc/esf')
36 #outPath = '/media/soporte/E9F4-F053/AMISR/Data/JULIA/ESF'
36 #outPath = '/media/soporte/E9F4-F053/AMISR/Data/JULIA/ESF'
37 outPath = '/home/soporte/Data/ESF'
37 outPath = '/home/soporte/Data/OutTest/ESF'
38 remotefolder = "/home/wmaster/graficos"
38 remotefolder = "/home/wmaster/graficos"
39 #figpath = '/home/soporte/Data/OutTest/EEJ'
39
40
40 xmin = '17'
41 xmin = '18'
41 xmax = '32'
42 xmax = '07'
42 dbmin = '60' #'60'#'55' #'40' #noise esf eej
43 dbmin = '60' #'60'#'55' #'40' #noise esf eej
43 dbmax = '75' #'70' #'55'
44 dbmax = '75' #'70' #'55'
44 show = '1'
45 show = '1'
@@ -50,8 +51,8 def main():
50 nosamp = '2' # oversample
51 nosamp = '2' # oversample
51 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
52 localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc
52
53
53 startDate='2020/03/05'
54 startDate='2019/03/17'
54 endDate='2020/03/06'
55 endDate='2019/12/18'
55
56
56 # str = datetime.date.today() #ONLINE
57 # str = datetime.date.today() #ONLINE
57 # str1 = str + datetime.timedelta(days=1)
58 # str1 = str + datetime.timedelta(days=1)
@@ -286,8 +287,8 def main():
286
287
287
288
288 # print "Escribiendo el archivo XML"
289 # print "Escribiendo el archivo XML"
289 filename="/home/soporte/workspace/schain/schainpy/scripts/joab.xml"
290 #filename="/home/soporte/workspace/schain/schainpy/scripts/joab.xml"
290 controllerObj.writeXml(filename)
291 #controllerObj.writeXml(filename)
291
292
292 controllerObj.start()
293 controllerObj.start()
293 # print "Leyendo el archivo XML"
294 # print "Leyendo el archivo XML"
General Comments 0
You need to be logged in to leave comments. Login now