##// END OF EJS Templates
schain
RSS
Scripts for processing of AMISR data
Julio Valdez -
r516:06661740c196
parent child
Show More
Show file before | Show file after | Hide comments
@@ -0,0 +1,123
1 import os, sys
2
3 path = os.path.split(os.getcwd())[0]
4 sys.path.append(path)
5
6 from controller import *
7
8 desc = "AMISR Experiment"
9
10 filename = "amisr_reader.xml"
11
12 controllerObj = Project()
13
14 controllerObj.setup(id = '191', name='test01', description=desc)
15
16
17 path = os.path.join(os.environ['HOME'],'amisr')
18 path = '/media/signalchain/HD-PXU2/AMISR_JULIA_MODE'
19 figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/eej')
20
21 xmin = '7'
22 xmax = '15'
23
24 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
25 path=path,
26 startDate='2014/10/07',
27 endDate='2014/10/07',
28 startTime='07:00:00',
29 endTime='15:00:00',
30 walk=0,
31 timezone='lt',
32 all=0,
33 online=0)
34
35 #AMISR Processing Unit
36 procUnitAMISRBeam0 = controllerObj.addProcUnit(datatype='AMISRProc', inputId=readUnitConfObj.getId())
37
38 opObj11 = procUnitAMISRBeam0.addOperation(name='PrintInfo', optype='other')
39
40 #Reshaper
41 opObj11 = procUnitAMISRBeam0.addOperation(name='ProfileToChannels', optype='other')
42
43 #Voltage Processing Unit
44 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=procUnitAMISRBeam0.getId())
45 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
46 opObj10.addParameter(name='frequency', value='445e6', format='float')
47
48 # opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
49 # opObj12.addParameter(name='minHei', value='0', format='float')
50
51 # 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'
52 # 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'
53 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'
54 opObj11 = procUnitConfObjBeam0.addOperation(name='Decoder', optype='other')
55 opObj11.addParameter(name='code', value=code, format='floatlist')
56 opObj11.addParameter(name='nCode', value='1', format='int')
57 opObj11.addParameter(name='nBaud', value='28', format='int')
58
59 # opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
60 # opObj12.addParameter(name='minHei', value='50', format='float')
61 # opObj12.addParameter(name='maxHei', value='150', format='float')
62 #Coherent Integration
63 #opObj11 = procUnitConfObjBeam0.addOperation(name='CohInt', optype='other')
64 #opObj11.addParameter(name='timeInterval', value='10', format='int')
65
66 #Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
67 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
68 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=64, format='int')
69 procUnitConfObjSpectraBeam0.addParameter(name='nProfiles', value=64, format='int')
70
71 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
72 # opObj11.addParameter(name='n', value='90', format='int')
73 opObj11.addParameter(name='timeInterval', value='30', format='float')
74
75 #RemoveDc
76 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
77
78 #Noise Estimation
79 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
80 opObj11.addParameter(name='minHei', value='100', format='float')
81 opObj11.addParameter(name='maxHei', value='280', format='float')
82
83 #SpectraPlot
84 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='other')
85 opObj11.addParameter(name='id', value='1', format='int')
86 opObj11.addParameter(name='wintitle', value='AMISR Beam 0', format='str')
87 opObj11.addParameter(name='zmin', value='38', format='int')
88 opObj11.addParameter(name='zmax', value='68', format='int')
89 opObj11.addParameter(name='save', value='1', format='bool')
90 opObj11.addParameter(name='figpath', value = figpath, format='str')
91 opObj11.addParameter(name='save', value='1', format='bool')
92 opObj11.addParameter(name='figpath', value = figpath, format='str')
93
94 #RTIPlot
95 #title0 = 'RTI AMISR Beam 0'
96 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='other')
97 opObj11.addParameter(name='id', value='2', format='int')
98 # opObj11.addParameter(name='wintitle', value=title0, format='str')
99 opObj11.addParameter(name='showprofile', value='0', format='int')
100 opObj11.addParameter(name='xmin', value=xmin, format='float')
101 opObj11.addParameter(name='xmax', value=xmax, format='float')
102 opObj11.addParameter(name='zmin', value='38', format='int')
103 opObj11.addParameter(name='zmax', value='68', format='int')
104 opObj11.addParameter(name='save', value='1', format='bool')
105 opObj11.addParameter(name='figpath', value = figpath, format='str')
106
107
108
109 #-----------------------------------------------------------------------------------------------
110
111
112 print "Escribiendo el archivo XML"
113 controllerObj.writeXml(filename)
114 print "Leyendo el archivo XML"
115 controllerObj.readXml(filename)
116
117 controllerObj.createObjects()
118 controllerObj.connectObjects()
119 controllerObj.run()
120
121 #21 3 pm
122
123
Show file before | Show file after | Hide comments
@@ -0,0 +1,119
1 import os, sys
2
3 path = os.path.split(os.getcwd())[0]
4 sys.path.append(path)
5
6 from controller import *
7
8 desc = "AMISR Experiment"
9
10 filename = "amisr_reader.xml"
11
12 controllerObj = Project()
13
14 controllerObj.setup(id = '191', name='test01', description=desc)
15
16
17 path = os.path.join(os.environ['HOME'],'amisr')
18 path = '/media/signalchain/HD-PXU2/AMISR_JULIA_MODE'
19 figpath = os.path.join(os.environ['HOME'],'Pictures/amisr/spreadF/filt2')
20
21 xmin = '18'
22 xmax = '31'
23
24 readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
25 path=path,
26 startDate='2014/10/06',
27 endDate='2014/10/07',
28 startTime='18:00:00',
29 endTime='07:00:00',
30 walk=0,
31 timezone='lt',
32 all=0,
33 online=0)
34
35 #AMISR Processing Unit
36 procUnitAMISRBeam0 = controllerObj.addProcUnit(datatype='AMISRProc', inputId=readUnitConfObj.getId())
37
38 opObj11 = procUnitAMISRBeam0.addOperation(name='PrintInfo', optype='other')
39
40 #Reshaper
41 opObj11 = procUnitAMISRBeam0.addOperation(name='ProfileToChannels', optype='other')
42
43 #Voltage Processing Unit
44 procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=procUnitAMISRBeam0.getId())
45 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
46 opObj10.addParameter(name='frequency', value='445e6', format='float')
47
48 opObj11 = procUnitConfObjBeam0.addOperation(name='filterByHeights')
49 opObj11.addParameter(name='window', value='2', format='int')
50
51 # 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'
52 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'
53 opObj11 = procUnitConfObjBeam0.addOperation(name='Decoder', optype='other')
54 opObj11.addParameter(name='code', value=code, format='floatlist')
55 opObj11.addParameter(name='nCode', value='1', format='int')
56 opObj11.addParameter(name='nBaud', value='28', format='int')
57
58 # opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
59 # opObj12.addParameter(name='minHei', value='0', format='float')
60 # opObj12.addParameter(name='maxHei', value='10', format='float')
61 #Coherent Integration
62 # opObj11 = procUnitConfObjBeam0.addOperation(name='CohInt', optype='other')
63 # opObj11.addParameter(name='timeInterval', value='30', format='int')
64 #Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
65 procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
66 procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=64, format='int')
67 procUnitConfObjSpectraBeam0.addParameter(name='nProfiles', value=64, format='int')
68
69 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
70 opObj11.addParameter(name='timeInterval', value='30', format='float')
71 # opObj11.addParameter(name='n', value='64', format='int')
72
73 #RemoveDc
74 #opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
75
76 #Noise Estimation
77 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
78 opObj11.addParameter(name='minHei', value='200', format='float')
79 opObj11.addParameter(name='maxHei', value='700', format='float')
80
81 #SpectraPlot
82 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='other')
83 opObj11.addParameter(name='id', value='1', format='int')
84 opObj11.addParameter(name='wintitle', value='AMISR Beam 0', format='str')
85 opObj11.addParameter(name='zmin', value='54', format='int')
86 opObj11.addParameter(name='zmax', value='70', format='int')
87 opObj11.addParameter(name='save', value='1', format='bool')
88 opObj11.addParameter(name='figpath', value = figpath, format='str')
89
90 #RTIPlot
91 #title0 = 'RTI AMISR Beam 0'
92 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='other')
93 opObj11.addParameter(name='id', value='2', format='int')
94 # opObj11.addParameter(name='wintitle', value=title0, format='str')
95 opObj11.addParameter(name='showprofile', value='0', format='int')
96 opObj11.addParameter(name='xmin', value=xmin, format='float')
97 opObj11.addParameter(name='xmax', value=xmax, format='float')
98 opObj11.addParameter(name='zmin', value='54', format='int')
99 opObj11.addParameter(name='zmax', value='70', format='int')
100 opObj11.addParameter(name='save', value='1', format='bool')
101 opObj11.addParameter(name='figpath', value = figpath, format='str')
102
103
104
105 #-----------------------------------------------------------------------------------------------
106
107
108 print "Escribiendo el archivo XML"
109 controllerObj.writeXml(filename)
110 print "Leyendo el archivo XML"
111 controllerObj.readXml(filename)
112
113 controllerObj.createObjects()
114 controllerObj.connectObjects()
115 controllerObj.run()
116
117 #21 3 pm
118
119
General Comments 0
You need to be logged in to leave comments. Login now