amisr_eej_proc_online_v3.py
169 lines
| 8.0 KiB
| text/x-python
|
PythonLexer
r1370 | #!/usr/bin/env python | |||
import os, sys | ||||
import time | ||||
import datetime | ||||
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 EEJ Experiment" | ||||
filename = "amisr_reader.xml" | ||||
xmin = '07' | ||||
r1374 | xmax = '18' # | |||
xrange = int(xmax)-int(xmin) | ||||
r1370 | dbmin = '45' #'60'#'55' #'40' #noise esf eej | |||
dbmax = '65' #'70' #'55' | ||||
r1375 | showSPC = '0' #view plot Spectra | |||
showRTI = '0' #view plot RTI | ||||
showNOISE = '0' #view plot NOISE | ||||
r1374 | localtime='1' #para ajustar el horario en las gráficas '0' para dejar en utc | |||
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' | |||
nCode = '1' | ||||
nBaud = '28' | ||||
nosamp = '1' # oversample for EEJ | ||||
str0 = datetime.date.today() | ||||
today = str0.strftime("%Y/%m/%d") | ||||
startDate=today | ||||
endDate=today | ||||
r1375 | #startDate='2021/07/17' | |||
#endDate='2021/07/17' | ||||
#inPath= '/home/soporte/Data/OutTest/ONLINE/' | ||||
r1370 | #inPath= '/home/soporte/dataAMISR/' | |||
inPath = '/mnt/data_amisr' | ||||
outPath = '/home/soporte/Data/EEJ' | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
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+"/EEJ"+l[0]+str(DOY) | ||||
if os.path.exists(outPath): | ||||
print("outPath", outPath) | ||||
else : | ||||
os.mkdir(outPath) | ||||
print("Creating...", outPath) | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
controllerObj = Project() | ||||
controllerObj.setup(id = '10', name='eej_proc', description=desc) | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
readUnitConfObj = controllerObj.addReadUnit(datatype='AMISRReader', | ||||
path=inPath, | ||||
startDate=startDate, | ||||
endDate=endDate, | ||||
r1375 | startTime='07:01:00',#'07:00:00', | |||
r1374 | endTime='17:55:00',#'15:00:00', | |||
r1370 | walk=0, | |||
code = code, | ||||
nCode = nCode, | ||||
nBaud = nBaud, | ||||
r1374 | timezone='lt', | |||
r1370 | online=1) | |||
#AMISR Processing Unit | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId()) | ||||
r1374 | opObj10 = procUnitConfObj0.addOperation(name='setAttribute') | |||
opObj10.addParameter(name='frequency', value='445.0e6') | ||||
r1370 | ||||
opObj01 = procUnitConfObj0.addOperation(name='Decoder', optype='other') | ||||
opObj01.addParameter(name='code', value=code, format='floatlist') | ||||
opObj01.addParameter(name='nCode', value=nCode, format='int') | ||||
opObj01.addParameter(name='nBaud', value=nBaud, format='int') | ||||
opObj01.addParameter(name='osamp', value=nosamp, format='int') | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
procUnitConfObj1 = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObj0.getId()) | ||||
procUnitConfObj1.addParameter(name='nFFTPoints', value='16', format='int') | ||||
opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other') | ||||
opObj11.addParameter(name='n', value='150', format='int') #300 normal value | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
# | ||||
opObj12 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='external') | ||||
opObj12.addParameter(name='id', value='21', format='int') | ||||
opObj12.addParameter(name= 'xaxis', value='velocity') | ||||
opObj12.addParameter(name='showprofile', value='1', format='int') | ||||
opObj12.addParameter(name='wintitle', value='AMISR Beam 0', format='str') | ||||
opObj12.addParameter(name='zmin', value=dbmin, format='int') | ||||
opObj12.addParameter(name='zmax', value=dbmax, format='int') | ||||
opObj12.addParameter(name='save', value=outPath+'/plots', format='str') | ||||
opObj12.addParameter(name='colormap', value='jet', format='str') | ||||
opObj12.addParameter(name='localtime', value=localtime,format='int') | ||||
opObj12.addParameter(name='show', value = showSPC, format='int') | ||||
r1375 | opObj12.addParameter(name='exp_code', value='206', format='int') | |||
opObj12.addParameter(name='server', value='jro-realtime.igp.gob.pe:4444') | ||||
opObj12.addParameter(name='sender_period', value='120') | ||||
opObj12.addParameter(name='tag', value='AMISR') | ||||
r1370 | ||||
r1374 | opObj14 = procUnitConfObj1.addOperation(name='NoisePlot', optype='external') | |||
opObj14.addParameter(name='id', value='3', format='int') | ||||
opObj14.addParameter(name='wintitle', value='title0', format='str') | ||||
opObj14.addParameter(name='showprofile', value='0', format='int') | ||||
opObj14.addParameter(name='tmin', value=xmin, format='int') | ||||
opObj14.addParameter(name='xrange', value=xrange, format='int') | ||||
opObj14.addParameter(name='ymin', value=dbmin, format='int') | ||||
opObj14.addParameter(name='ymax', value=dbmax, format='int') | ||||
opObj14.addParameter(name='save', value=outPath+'/plots', format='str') | ||||
opObj14.addParameter(name='localtime', value=localtime,format='int') | ||||
opObj14.addParameter(name='show', value = showNOISE, format='int') | ||||
r1375 | opObj14.addParameter(name='exp_code', value='206', format='int') | |||
opObj14.addParameter(name='server', value='jro-realtime.igp.gob.pe:4444') | ||||
opObj14.addParameter(name='sender_period', value='120') | ||||
opObj14.addParameter(name='tag', value='AMISR') | ||||
r1370 | ||||
opObj15 = procUnitConfObj1.addOperation(name='RTIPlot', optype='external') | ||||
opObj15.addParameter(name='id', value='2', format='int') | ||||
opObj15.addParameter(name='localtime', value=localtime,format='int') | ||||
opObj15.addParameter(name='wintitle', value='RTI', format='str') | ||||
r1374 | opObj15.addParameter(name='tmin', value=xmin, format='int') | |||
opObj15.addParameter(name='xrange', value=xrange, format='int') | ||||
r1370 | opObj15.addParameter(name='zmin', value=dbmin, format='int') | |||
opObj15.addParameter(name='zmax', value=dbmax, format='int') | ||||
opObj15.addParameter(name='showprofile', value='1', format='int') | ||||
opObj15.addParameter(name='save', value=outPath+'/plots', format='str') | ||||
opObj15.addParameter(name='colormap', value='jet', format='str') | ||||
opObj15.addParameter(name='show', value = showRTI, format='int') | ||||
r1375 | opObj15.addParameter(name='exp_code', value='206', format='int') | |||
opObj15.addParameter(name='server', value='jro-realtime.igp.gob.pe:4444') | ||||
opObj15.addParameter(name='sender_period', value='120') | ||||
opObj15.addParameter(name='tag', value='AMISR') | ||||
r1370 | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
# | ||||
controllerObj.start() | ||||
##....................................................................................... | ||||
##....................................................................................... | ||||
if __name__ == '__main__': | ||||
import time | ||||
start_time = time.time() | ||||
main() | ||||
print("--- %s seconds ---" % (time.time() - start_time)) | ||||