##// END OF EJS Templates
schain
RSS
f_
f_
joabAM - - Load All Authors

File last commit:

r1374:7243e1f8c2ee
r1374:7243e1f8c2ee
Show More
amisr_esf_proc_online_v3.py
203 lines | 9.2 KiB | text/x-python | PythonLexer
/ schainpy / scripts / amisr_esf_proc_online_v3.py
History | Source | Raw |Copy content |Copy permalink
joabAM
merged branches
 r1370 #!/usr/bin/env python
import os, sys
import time
import datetime
'''
ESF 10 BEAM ONLINE
'''
path = os.path.dirname(os.getcwd())
path = os.path.dirname(path)
sys.path.insert(0, path)
from schainpy.controller import Project
def main():
desc = "AMISR ESF 10 Beam Experiment"
filename = "amisr_esf_reader.xml"
#path = os.path.join(os.environ['HOME'],'amisr')
inPath = '/mnt/data_amisr'
#path= '/home/soporte/dataAMISR/'
#inPath= '/home/soporte/dataAMISR_test/'
outPath = '/home/soporte/Data/ESF'
remotefolder = "/home/wmaster/graficos"
joabAM
f_
 r1374 xmin = '18'
joabAM
merged branches
 r1370 xmax = '07'
dbmin = '60' #'60'#'55' #'40' #noise esf eej
dbmax = '75' #'70' #'55'
showSPC = '1' #view plot Spectra
showRTI = '1' #view plot RTI
joabAM
f_
 r1374 showNOISE = '1' #view plot NOISE
joabAM
merged branches
 r1370 # # 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'
# # 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'
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'
nCode = '1'
nBaud = '28'
nosamp = '2' # oversample
localtime='0' #para ajustar el horario en las gráficas '0' para dejar en utc
dty = datetime.date.today() #ONLINE
str1 = dty + datetime.timedelta(days=1)
str2 = dty - datetime.timedelta(days=1)
today = dty.strftime("%Y/%m/%d")
tomorrow = str1.strftime("%Y/%m/%d")
yesterday = str2.strftime("%Y/%m/%d")
startDate=today
endDate=tomorrow
#startDate='2019/12/17'
#endDate='2019/12/18'
##.......................................................................................
##.......................................................................................
l = startDate.split('/') #adding day of the year to outPath
datelist = datetime.date(int(l[0]),int(l[1]),int(l[2]))
DOY = datelist.timetuple().tm_yday
outPath= outPath+"/ESF"+l[0]+str(DOY)
if os.path.exists(outPath):
print("outPath", outPath)
else :
os.mkdir(outPath)
print("Creating...", outPath)
##.......................................................................................
##.......................................................................................
controllerObj = Project()
joabAM
f_
 r1374 controllerObj.setup(id = '21', name='esf_proc', description=desc)
joabAM
merged branches
 r1370 ##.......................................................................................
##.......................................................................................
readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader',
path=inPath,
startDate=startDate,#'2016/07/12',
endDate=endDate,#'2016/07/13',
startTime='19:01:30',#'07:00:00',
endTime='07:00:00',#'15:00:00',
walk=0,
code = code,
nCode = nCode,
nBaud = nBaud,
timezone='lt',
online=1)
#AMISR Processing Unit
##.......................................................................................
##.......................................................................................
#Voltage Processing Unit
procUnitConfObjBeam0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
opObj10 = procUnitConfObjBeam0.addOperation(name='setRadarFrequency')
#opObj10.addParameter(name='frequency', value='445e6', format='int') #changed on Dec 3, 15:40h
#Noise--> no code
opObj11 = procUnitConfObjBeam0.addOperation(name='Decoder', optype='other')
opObj11.addParameter(name='code', value=code, format='floatlist')
opObj11.addParameter(name='nCode', value=nCode, format='int')
opObj11.addParameter(name='nBaud', value=nBaud, format='int')
opObj11.addParameter(name='osamp', value=nosamp, format='int')
# opObj12 = procUnitConfObjBeam0.addOperation(name='selectHeights')
# opObj12.addParameter(name='minHei', value='50', format='float')
# opObj12.addParameter(name='maxHei', value='150', format='float')
# #Coherent Integration
# opObj11 = procUnitConfObjBeam0.addOperation(name='CohInt', optype='other')
# opObj11.addParameter(name='n', value='2', format='int')
##.......................................................................................
##.......................................................................................
#Spectra Unit Processing, getting spectras with nProfiles and nFFTPoints
procUnitConfObjSpectraBeam0 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjBeam0.getId())
procUnitConfObjSpectraBeam0.addParameter(name='nFFTPoints', value=16, format='int')
#
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='IncohInt', optype='other')
opObj11.addParameter(name='n', value='60', format='int')
#opObj11.addParameter(name='timeInterval', value='30', format='int')
# # #RemoveDc
#opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='removeDC')
#---
joabAM
f_
 r1374 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='selectChannels')
opObj11.addParameter(name='channelList', value='0,2,4,6,8', format='intlist')
joabAM
merged branches
 r1370
##.......................................................................................
##.......................................................................................
#SpectraPlot
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='SpectraPlot', optype='external')
opObj11.addParameter(name='id', value='1', format='int')
opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
opObj11.addParameter(name='showprofile', value='1', format='int')
opObj11.addParameter(name='zmin', value=dbmin, format='int')
opObj11.addParameter(name='zmax', value=dbmax, format='int')
opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
# opObj11.addParameter(name='ftp', value='1', format='int')
# opObj11.addParameter(name='wr_period', value='2', format='int')
# opObj11.addParameter(name='exp_code', value='21', format='int')
# opObj11.addParameter(name='sub_exp_code', value='4', format='int')
# opObj11.addParameter(name='ftp_wei', value='0', format='int')
# opObj11.addParameter(name='plot_pos', value='0', format='int')
opObj11.addParameter(name='localtime', value=localtime,format='int')
opObj11.addParameter(name='show', value = showSPC, format='int')
opObj11.addParameter(name='colormap', value='jet', format='str')
#RTIPlot
#title0 = 'RTI AMISR Beam 0'
opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='RTIPlot', optype='external')
opObj11.addParameter(name='id', value='2', format='int')
opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
opObj11.addParameter(name='showprofile', value='0', format='int')
joabAM
f_
 r1374 opObj11.addParameter(name='tmin', value=xmin, format='int')
opObj11.addParameter(name='tmax', value=xmax,format='int')
joabAM
merged branches
 r1370 opObj11.addParameter(name='zmin', value=dbmin, format='int')
opObj11.addParameter(name='zmax', value=dbmax, format='int')
opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
#opObj11.addParameter(name='ftp', value='1', format='int')
#opObj11.addParameter(name='wr_period', value='2', format='int')
# opObj11.addParameter(name='exp_code', value='21', format='int')
# opObj11.addParameter(name='sub_exp_code', value='4', format='int')
# opObj11.addParameter(name='ftp_wei', value='0', format='int')
# opObj11.addParameter(name='plot_pos', value='0', format='int')
opObj11.addParameter(name='localtime', value=localtime,format='int')
opObj11.addParameter(name='show', value = showRTI, format='int')
opObj11.addParameter(name='colormap', value='jet', format='str')
# # #
# #Noise
#title0 = 'RTI AMISR Beam 0'
joabAM
f_
 r1374 opObj11 = procUnitConfObjSpectraBeam0.addOperation(name='NoisePlot', optype='external')
opObj11.addParameter(name='id', value='3', format='int')
opObj11.addParameter(name='wintitle', value='ESF AMISR', format='str')
opObj11.addParameter(name='showprofile', value='0', format='int')
opObj11.addParameter(name='tmin', value=xmin, format='int')
opObj11.addParameter(name='tmax', value=xmax, format='int')
##opObj11.addParameter(name='ymin', value=dbmin, format='int')
##opObj11.addParameter(name='ymax', value=dbmax, format='int')
opObj11.addParameter(name='save', value=outPath+'/plots', format='str')
opObj11.addParameter(name='show', value = showNOISE, format='int')
opObj11.addParameter(name='localtime', value=localtime,format='int')
joabAM
merged branches
 r1370
# print "Escribiendo el archivo XML"
# controllerObj.writeXml(filename)
controllerObj.start()
# print "Leyendo el archivo XML"
# controllerObj.readXml(filename)
if __name__ == '__main__':
import time
start_time = time.time()
main()
print("--- %s seconds ---" % (time.time() - start_time))