##// END OF EJS Templates
schain
RSS
First Spectral Fitting and EW Drifts operative module inside Signal Chain TRUNK
First Spectral Fitting and EW Drifts operative module inside Signal Chain TRUNK
Julio Valdez - - Load All Authors

File last commit:

r507:5d24178b6c75
r513:19f921674eb5
Show More
amisr_windEstimation.py
147 lines | 6.6 KiB | text/x-python | PythonLexer
/ schainpy / test / amisr_windEstimation.py
History | Source | Raw |Copy content |Copy permalink
Julio Valdez
Processing Modules added:...
 r502 import os, sys
path = os.path.split(os.getcwd())[0]
sys.path.append(path)
from controller import *
desc = "AMISR Experiment"
filename = "amisr_reader.xml"
controllerObj = Project()
controllerObj.setup(id = '191', name='test01', description=desc)
path = os.path.join(os.environ['HOME'],'Development/amisr/data')
path = '/home/soporte/Data/AMISR'
figpath = os.path.join(os.environ['HOME'],'Pictures/amisr')
Julio Valdez
AMISR wind estimation script updated
 r507 pathFigure = '/home/operaciones/Documents/AMISR_windprofiler/20141023'
xmin = '8.0'
xmax = '12.0'
Julio Valdez
Processing Modules added:...
 r502
readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
path=path,
Julio Valdez
AMISR wind estimation script updated
 r507 startDate='2014/10/23',
endDate='2014/10/23',
Julio Valdez
Processing Modules added:...
 r502 startTime='00:00:00',
endTime='23:59:59',
walk=1,
timezone='lt')
#AMISR Processing Unit
procUnitAMISRBeam0 = controllerObj.addProcUnit(datatype='AMISRProc', inputId=readUnitConfObj.getId())
opObj11 = procUnitAMISRBeam0.addOperation(name='PrintInfo', optype='other')
#Reshaper
opObj11 = procUnitAMISRBeam0.addOperation(name='ProfileToChannels', optype='other')
#Beam Selector
#opObj11 = procUnitAMISRBeam0.addOperation(name='BeamSelector', optype='other')
#opObj11.addParameter(name='beam', value='0', format='int')
#Voltage Processing Unit
procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=procUnitAMISRBeam0.getId())
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
opObj10.addParameter(name='frequency', value='445e6', format='float')
opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
opObj12.addParameter(name='minHei', value='0', format='float')
opObj12.addParameter(name='maxHei', value='10', format='float')
Julio Valdez
Processing Modules added:...
 r502 #Coherent Integration
opObj11 = procUnitConfObjBeam0.addOperation(name='CohInt', optype='other')
opObj11.addParameter(name='n', value='8', format='int')
#Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=32, format='int')
procUnitConfObjSpectraBeam0.addParameter(name='nProfiles', value=32, format='int')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
opObj11.addParameter(name='n', value='16', format='int')
Julio Valdez
Processing Modules added:...
 r502 #RemoveDc
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
#Noise Estimation
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
opObj11.addParameter(name='minHei', value='5', format='float')
opObj11.addParameter(name='maxHei', value='20', format='float')
#SpectraPlot
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='other')
opObj11.addParameter(name='id', value='100', format='int')
opObj11.addParameter(name='wintitle', value='AMISR Beam 0', format='str')
opObj11.addParameter(name='zmin', value='30', format='int')
opObj11.addParameter(name='zmax', value='80', format='int')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj11.addParameter(name='save', value='1', format='bool')
opObj11.addParameter(name='figpath', value = pathFigure, format='str')
Julio Valdez
Processing Modules added:...
 r502 #RTIPlot
#title0 = 'RTI AMISR Beam 0'
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='other')
opObj11.addParameter(name='id', value='200', format='int')
# opObj11.addParameter(name='wintitle', value=title0, format='str')
opObj11.addParameter(name='showprofile', value='0', format='int')
opObj11.addParameter(name='xmin', value=xmin, format='float')
opObj11.addParameter(name='xmax', value=xmax, format='float')
Julio Valdez
AMISR wind estimation script updated
 r507 opObj11.addParameter(name='zmin', value='30', format='int')
opObj11.addParameter(name='zmax', value='80', format='int')
opObj23.addParameter(name='save', value='1', format='bool')
opObj23.addParameter(name='figpath', value = pathFigure, format='str')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 # Save RTI
# figfile0 = 'amisr_rti_beam0.png'
# opObj11.addParameter(name='figpath', value=figpath, format='str')
# opObj11.addParameter(name='figfile', value=figfile0, format='str')
Julio Valdez
Processing Modules added:...
 r502
#-----------------------------------------------------------------------------------------------
procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObjSpectraBeam0 .getId())
opObj20 = procUnitConfObj2.addOperation(name='GetMoments')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj21 = procUnitConfObj2.addOperation(name='MomentsPlot', optype='other')
opObj21.addParameter(name='id', value='3', format='int')
opObj21.addParameter(name='wintitle', value='Moments Plot', format='str')
opObj21.addParameter(name='save', value='1', format='bool')
opObj21.addParameter(name='figpath', value=pathFigure, format='str')
opObj21.addParameter(name='zmin', value='30', format='int')
opObj21.addParameter(name='zmax', value='80', format='int')
Julio Valdez
Processing Modules added:...
 r502
opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other')
opObj22.addParameter(name='technique', value='DBS', format='str')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj22.addParameter(name='correctAzimuth', value='51.06', format='float')
Julio Valdez
Processing Modules added:...
 r502 opObj22.addParameter(name='correctFactor', value='-1', format='float')
Julio Valdez
AMISR wind estimation script updated
 r507 opObj22.addParameter(name='azimuth', value='0,-90,0,90,180', format='floatlist')
opObj22.addParameter(name='elevation', value='74.53,75.90.0,75.60,75.60', format='floatlist')
Julio Valdez
Processing Modules added:...
 r502 # opObj22.addParameter(name='horizontalOnly', value='1', format='bool')
# opObj22.addParameter(name='channelList', value='1,2', format='intlist')
opObj23 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
opObj23.addParameter(name='id', value='4', format='int')
opObj23.addParameter(name='wintitle', value='Wind Profiler', format='str')
opObj23.addParameter(name='save', value='1', format='bool')
opObj23.addParameter(name='figpath', value = pathFigure, format='str')
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503 opObj23.addParameter(name='zmin', value='-20', format='int')
opObj23.addParameter(name='zmax', value='20', format='int')
opObj23.addParameter(name='zmin_ver', value='-100', format='float')
opObj23.addParameter(name='zmax_ver', value='100', format='float')
Julio Valdez
AMISR wind estimation script updated
 r507 opObj23.addParameter(name='SNRmin', value='-20', format='int')
opObj23.addParameter(name='SNRmax', value='30', format='int')
opObj23.addParameter(name='SNRthresh', value='-50', format='float')
Julio Valdez
Processing Modules added:...
 r502 opObj23.addParameter(name='xmin', value=xmin, format='float')
opObj23.addParameter(name='xmax', value=xmax, format='float')
print "Escribiendo el archivo XML"
controllerObj.writeXml(filename)
print "Leyendo el archivo XML"
controllerObj.readXml(filename)
controllerObj.createObjects()
controllerObj.connectObjects()
controllerObj.run()
Julio Valdez
-Modifications to DBS module to accept azimuth and elevation angles as inputs....
 r503
#21 3 pm