##// END OF EJS Templates
scripts to processing AMISR data
scripts to processing AMISR data

File last commit:

r485:ccf1a40188a9
r485:ccf1a40188a9
Show More
amisr_proc.py
79 lines | 3.1 KiB | text/x-python | PythonLexer
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'],'Documents/amisr')
figpath = os.path.join(os.environ['HOME'],'Pictures/amisr')
readUnitConfObj = controllerObj.addReadUnit(datatype='AMISR',
path=path,
startDate='2014/08/18',
endDate='2014/08/18',
startTime='00:00:00',
endTime='23:59:59',
walk=1)
#AMISR Processing Unit
procUnitAMISRBeam0 = controllerObj.addProcUnit(datatype='AMISR', inputId=readUnitConfObj.getId())
#Beam Selector
opObj11 = procUnitAMISRBeam0.addOperation(name='BeamSelector', optype='other')
opObj11.addParameter(name='beam', value='0', format='int')
#Voltage Processing Unit
procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='Voltage', inputId=procUnitAMISRBeam0.getId())
#Coherent Integration
opObj11 = procUnitConfObjBeam0.addOperation(name='CohInt', optype='other')
opObj11.addParameter(name='n', value='128', format='int')
#Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjBeam0.getId())
procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=32, format='int')
procUnitConfObjSpectraBeam0.addParameter(name='nProfiles', value=32, format='int')
#Noise Estimation
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='getNoise')
opObj11.addParameter(name='minHei', value='100', format='float')
opObj11.addParameter(name='maxHei', value='450', 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')
#RTIPlot
title0 = 'RTI AMISR Beam 0'
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')
#Setting RTI time using xmin,xmax
opObj11.addParameter(name='xmin', value='0', format='int')
opObj11.addParameter(name='xmax', value='18', format='int')
#Setting dB range with zmin, zmax
opObj11.addParameter(name='zmin', value='45', format='int')
opObj11.addParameter(name='zmax', value='70', format='int')
#Save RTI
figfile0 = 'amisr_rti_beam0.png'
opObj11.addParameter(name='figpath', value=figpath, format='str')
opObj11.addParameter(name='figfile', value=figfile0, format='str')
print "Escribiendo el archivo XML"
controllerObj.writeXml(filename)
print "Leyendo el archivo XML"
controllerObj.readXml(filename)
controllerObj.createObjects()
controllerObj.connectObjects()
controllerObj.run()