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