From c9bf2c97727fa856b6e2ca628c9845c24e384774 2014-04-16 14:19:29 From: Daniel Valdez Date: 2014-04-16 14:19:29 Subject: [PATCH] se agrega metodos para almacenar en disco las fases nuevo scripts para procesamiento de jasmet 30 y 50 --- diff --git a/schainpy/model/jroplot.py b/schainpy/model/jroplot.py index e6c87c5..ec1467e 100644 --- a/schainpy/model/jroplot.py +++ b/schainpy/model/jroplot.py @@ -1078,6 +1078,8 @@ class BeaconPhase(Figure): self.SUB_EXP_CODE = None self.PLOT_POS = None + self.filename_phase = None + def getSubplots(self): ncol = 1 @@ -1103,8 +1105,27 @@ class BeaconPhase(Figure): nrow, ncol = self.getSubplots() self.addAxes(nrow, ncol*ncolspan, 0, 0, colspan, 1) - - + + def save_phase(self, filename_phase): + f = open(filename_phase,'w+') + f.write('\n\n') + f.write('JICAMARCA RADIO OBSERVATORY - Beacon Phase \n') + f.write('DD MM YYYY HH MM SS pair(2,0) pair(2,1) pair(2,3) pair(2,4)\n\n' ) + f.close() + + def save_data(self, filename_phase, data, data_datetime): + f=open(filename_phase,'a') + timetuple_data = data_datetime.timetuple() + day = str(timetuple_data.tm_mday) + month = str(timetuple_data.tm_mon) + year = str(timetuple_data.tm_year) + hour = str(timetuple_data.tm_hour) + minute = str(timetuple_data.tm_min) + second = str(timetuple_data.tm_sec) + f.write(day+' '+month+' '+year+' '+hour+' '+minute+' '+second+' '+str(data[0])+' '+str(data[1])+' '+str(data[2])+' '+str(data[3])+'\n') + f.close() + + def run(self, dataOut, id, wintitle="", pairsList=None, showprofile='True', xmin=None, xmax=None, ymin=None, ymax=None, timerange=None, @@ -1182,7 +1203,17 @@ class BeaconPhase(Figure): self.xdata = numpy.array([]) self.ydata = numpy.array([]) + + #open file beacon phase + path = '%s%03d' %(self.PREFIX, self.id) + beacon_file = os.path.join(path,'%s.txt'%self.name) + self.filename_phase = os.path.join(figpath,beacon_file) + self.save_phase(self.filename_phase) + + #store data beacon phase + self.save_data(self.filename_phase, phase_beacon, thisDatetime) + self.setWinTitle(title) diff --git a/schainpy/test/Meteor_JASMET_30mhz.py b/schainpy/test/Meteor_JASMET_30mhz.py index a3c0356..b6eb250 100644 --- a/schainpy/test/Meteor_JASMET_30mhz.py +++ b/schainpy/test/Meteor_JASMET_30mhz.py @@ -19,72 +19,28 @@ path = '/home/dsuarez/.gvfs/data on 10.10.20.13/Jasmet50' readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage', path=path, - startDate='2013/08/21', - endDate='2013/08/21', - startTime='00:00:00', + startDate='2014/04/15', + endDate='2014/04/15', + startTime='17:00:00', endTime='23:59:59', - online=1, - delay=10, + online=0, walk=1) opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') procUnitConfObj0 = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId()) -opObj11 = procUnitConfObj0.addOperation(name='setRadarFrequency') -opObj11.addParameter(name='frequency', value='30.15e6', format='float') opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') + opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') opObj11.addParameter(name='n', value='2', format='int') -procUnitConfObj1 = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) -procUnitConfObj1.addParameter(name='nProfiles', value='200', format='int') -procUnitConfObj1.addParameter(name='nFFTPoints', value='200', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other') -opObj11.addParameter(name='n', value='4', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other') -opObj11.addParameter(name='id', value='100', format='int') -opObj11.addParameter(name='wintitle', value='JASMET30MHZ', format='str') -opObj11.addParameter(name='showprofile', value='0', format='int') -opObj11.addParameter(name='timerange', value='300', format='int') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -# # opObj11.addParameter(name='xmin', value='18', format='float') -# # opObj11.addParameter(name='xmax', value='', format='float') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='1', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other') -opObj11.addParameter(name='id', value='101', format='int') -opObj11.addParameter(name='wintitle', value='JASMET30MHZ', format='str') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='1', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other') -opObj11.addParameter(name='id', value='103', format='int') -opObj11.addParameter(name='wintitle', value='JASMET30MHZ', format='str') -opObj11.addParameter(name='showprofile', value='0', format='int') -#opObj11.addParameter(name='timerange', value='300', format='int') -opObj11.addParameter(name='xmin', value='0', format='float') -opObj11.addParameter(name='xmax', value='24', format='float') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') +opObj11 = procUnitConfObj0.addOperation(name='VoltageWriter', optype='other') +opObj11.addParameter(name='path', value='/home/jasmet/jasmet30_abril') +opObj11.addParameter(name='blocksPerFile', value='100', format='int') +opObj11.addParameter(name='profilesPerBlock', value='200', format='int') """ @@ -94,7 +50,7 @@ opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril' procUnitConfObjBeacon = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) procUnitConfObjBeacon.addParameter(name='nProfiles', value='200', format='int') procUnitConfObjBeacon.addParameter(name='nFFTPoints', value='200', format='int') -procUnitConfObjBeacon.addParameter(name='pairsList', value='(0,5),(1,5),(2,5),(3,5),(4,5)', format='pairsList') +procUnitConfObjBeacon.addParameter(name='pairsList', value='(2,0),(2,1),(2,3),(2,4)', format='pairsList') opObj11 = procUnitConfObjBeacon.addOperation(name='IncohInt', optype='other') opObj11.addParameter(name='n', value='4', format='int') @@ -105,15 +61,11 @@ opObj11 = procUnitConfObjBeacon.addOperation(name='BeaconPhase', optype='other') opObj11.addParameter(name='id', value='201', format='int') opObj11.addParameter(name='wintitle', value='Beacon Phase', format='str') opObj11.addParameter(name='timerange', value='300', format='int') +opObj11.addParameter(name='xmin', value='0', format='float') +opObj11.addParameter(name='xmax', value='24', format='float') opObj11.addParameter(name='ymin', value='-180', format='float') opObj11.addParameter(name='ymax', value='180', format='float') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='0', format='int') - +opObj11.addParameter(name='figpath', value='/home/jasmet/jasmet30_phase', format='str') print "Escribiendo el archivo XML" @@ -124,4 +76,3 @@ controllerObj.readXml(filename) controllerObj.createObjects() controllerObj.connectObjects() controllerObj.run() - diff --git a/schainpy/test/Meteor_JASMET_50mhz.py b/schainpy/test/Meteor_JASMET_50mhz.py index 7692cc6..b2b59ac 100644 --- a/schainpy/test/Meteor_JASMET_50mhz.py +++ b/schainpy/test/Meteor_JASMET_50mhz.py @@ -1,9 +1,5 @@ """ - - - Se debe verficar que el disco de datos se encuentra montado en el sistema - """ import os, sys @@ -23,12 +19,11 @@ path = '/home/dsuarez/.gvfs/data on 10.10.20.13/Jasmet50' readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage', path=path, - startDate='2013/08/21', - endDate='2013/08/21', - startTime='00:00:00', + startDate='2014/04/15', + endDate='2014/04/15', + startTime='17:00:00', endTime='23:59:59', - online=1, - delay=3, + online=0, walk=1) opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') @@ -42,56 +37,10 @@ opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') opObj11.addParameter(name='n', value='2', format='int') -procUnitConfObj1 = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) -procUnitConfObj1.addParameter(name='nProfiles', value='200', format='int') -procUnitConfObj1.addParameter(name='nFFTPoints', value='200', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other') -opObj11.addParameter(name='n', value='4', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other') -opObj11.addParameter(name='id', value='100', format='int') -opObj11.addParameter(name='wintitle', value='JASMET50MHZ', format='str') -opObj11.addParameter(name='showprofile', value='0', format='int') -opObj11.addParameter(name='timerange', value='300', format='int') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -# # opObj11.addParameter(name='xmin', value='18', format='float') -# # opObj11.addParameter(name='xmax', value='', format='float') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='lastone', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='0', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='SpectraPlot', optype='other') -opObj11.addParameter(name='id', value='101', format='int') -opObj11.addParameter(name='wintitle', value='JASMET50MHZ', format='str') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -opObj11.addParameter(name='wr_period', value='10', format='int') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='0', format='int') - -opObj11 = procUnitConfObj1.addOperation(name='RTIPlot', optype='other') -opObj11.addParameter(name='id', value='103', format='int') -opObj11.addParameter(name='wintitle', value='JASMET50MHZ', format='str') -opObj11.addParameter(name='showprofile', value='0', format='int') -#opObj11.addParameter(name='timerange', value='300', format='int') -opObj11.addParameter(name='xmin', value='0', format='float') -opObj11.addParameter(name='xmax', value='24', format='float') -opObj11.addParameter(name='zmin', value='55', format='float') -opObj11.addParameter(name='zmax', value='75', format='float') -opObj11.addParameter(name='wr_period', value='10', format='int') -opObj11.addParameter(name='save', value='1', format='int') -#opObj11.addParameter(name='lastone', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') +opObj11 = procUnitConfObj0.addOperation(name='VoltageWriter', optype='other') +opObj11.addParameter(name='path', value='/media/datos/jasmet50_abril') +opObj11.addParameter(name='blocksPerFile', value='100', format='int') +opObj11.addParameter(name='profilesPerBlock', value='200', format='int') """ @@ -101,7 +50,7 @@ opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril' procUnitConfObjBeacon = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) procUnitConfObjBeacon.addParameter(name='nProfiles', value='200', format='int') procUnitConfObjBeacon.addParameter(name='nFFTPoints', value='200', format='int') -procUnitConfObjBeacon.addParameter(name='pairsList', value='(0,5),(1,5),(2,5),(3,5),(4,5)', format='pairsList') +procUnitConfObjBeacon.addParameter(name='pairsList', value='(2,0),(2,1),(2,3),(2,4)', format='pairsList') opObj11 = procUnitConfObjBeacon.addOperation(name='IncohInt', optype='other') opObj11.addParameter(name='n', value='4', format='int') @@ -112,16 +61,11 @@ opObj11 = procUnitConfObjBeacon.addOperation(name='BeaconPhase', optype='other') opObj11.addParameter(name='id', value='201', format='int') opObj11.addParameter(name='wintitle', value='Beacon Phase', format='str') opObj11.addParameter(name='timerange', value='300', format='int') +opObj11.addParameter(name='xmin', value='0', format='float') +opObj11.addParameter(name='xmax', value='24', format='float') opObj11.addParameter(name='ymin', value='-180', format='float') opObj11.addParameter(name='ymax', value='180', format='float') -opObj11.addParameter(name='wr_period', value='10', format='int') -opObj11.addParameter(name='save', value='1', format='int') -opObj11.addParameter(name='figpath', value='/home/dsuarez/Pictures/jasmet_abril', format='str') -opObj11.addParameter(name='ftp', value='1', format='int') -opObj11.addParameter(name='ftp_wei', value='1', format='int') -opObj11.addParameter(name='exp_code', value='15', format='int') -opObj11.addParameter(name='sub_exp_code', value='0', format='int') - +opObj11.addParameter(name='figpath', value='/media/datos/jasmet50_phase', format='str') print "Escribiendo el archivo XML" diff --git a/schainpy/test/jasmet30_16april.py b/schainpy/test/jasmet30_16april.py new file mode 100644 index 0000000..9ea3acd --- /dev/null +++ b/schainpy/test/jasmet30_16april.py @@ -0,0 +1,77 @@ +""" +Se debe verficar que el disco de datos se encuentra montado en el sistema +""" +import os, sys + +path = os.path.split(os.getcwd())[0] +sys.path.append(path) + +from controller import * + +desc = "Meteor Experiment Test" +filename = "meteor20130812.xml" + +controllerObj = Project() +controllerObj.setup(id = '191', name='meteor_test01', description=desc) + +path = '/home/dsuarez/.gvfs/data on 10.10.20.13/Jasmet50' + + +readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage', + path=path, + startDate='2014/04/16', + endDate='2014/04/16', + startTime='00:00:00', + endTime='23:59:59', + online=0, + walk=1) + +opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') + +procUnitConfObj0 = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId()) + + +opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') + + +opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') +opObj11.addParameter(name='n', value='2', format='int') + +opObj11 = procUnitConfObj0.addOperation(name='VoltageWriter', optype='other') +opObj11.addParameter(name='path', value='/home/jasmet/jasmet30_abril') +opObj11.addParameter(name='blocksPerFile', value='100', format='int') +opObj11.addParameter(name='profilesPerBlock', value='200', format='int') + + +""" +########################################### BEACON ########################################## +""" + +procUnitConfObjBeacon = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) +procUnitConfObjBeacon.addParameter(name='nProfiles', value='200', format='int') +procUnitConfObjBeacon.addParameter(name='nFFTPoints', value='200', format='int') +procUnitConfObjBeacon.addParameter(name='pairsList', value='(2,0),(2,1),(2,3),(2,4)', format='pairsList') + +opObj11 = procUnitConfObjBeacon.addOperation(name='IncohInt', optype='other') +opObj11.addParameter(name='n', value='4', format='int') + +opObj11 = procUnitConfObjBeacon.addOperation(name='getBeaconSignal') + +opObj11 = procUnitConfObjBeacon.addOperation(name='BeaconPhase', optype='other') +opObj11.addParameter(name='id', value='301', format='int') +opObj11.addParameter(name='wintitle', value='Beacon Phase', format='str') +opObj11.addParameter(name='xmin', value='0', format='float') +opObj11.addParameter(name='xmax', value='24', format='float') +opObj11.addParameter(name='ymin', value='-180', format='float') +opObj11.addParameter(name='ymax', value='180', format='float') +opObj11.addParameter(name='figpath', value='/home/jasmet/jasmet30_phase', format='str') + + +print "Escribiendo el archivo XML" +controllerObj.writeXml(filename) +print "Leyendo el archivo XML" +controllerObj.readXml(filename) + +controllerObj.createObjects() +controllerObj.connectObjects() +controllerObj.run() diff --git a/schainpy/test/jasmet50_16april.py b/schainpy/test/jasmet50_16april.py new file mode 100644 index 0000000..92d85c4 --- /dev/null +++ b/schainpy/test/jasmet50_16april.py @@ -0,0 +1,79 @@ +""" +Se debe verficar que el disco de datos se encuentra montado en el sistema +""" +import os, sys + +path = os.path.split(os.getcwd())[0] +sys.path.append(path) + +from controller import * + +desc = "Meteor Experiment Test" +filename = "meteor20130812.xml" + +controllerObj = Project() +controllerObj.setup(id = '191', name='meteor_test01', description=desc) + +path = '/home/dsuarez/.gvfs/data on 10.10.20.13/Jasmet50' + + +readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage', + path=path, + startDate='2014/04/16', + endDate='2014/04/16', + startTime='00:00:00', + endTime='23:59:59', + online=0, + walk=1) + +opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') + +procUnitConfObj0 = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId()) + +opObj11 = procUnitConfObj0.addOperation(name='setRadarFrequency') +opObj11.addParameter(name='frequency', value='30.15e6', format='float') + +opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') + + +opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') +opObj11.addParameter(name='n', value='2', format='int') + +opObj11 = procUnitConfObj0.addOperation(name='VoltageWriter', optype='other') +opObj11.addParameter(name='path', value='/media/datos/jasmet50_abril') +opObj11.addParameter(name='blocksPerFile', value='100', format='int') +opObj11.addParameter(name='profilesPerBlock', value='200', format='int') + + +""" +########################################### BEACON ########################################## +""" + +procUnitConfObjBeacon = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObj0.getId()) +procUnitConfObjBeacon.addParameter(name='nProfiles', value='200', format='int') +procUnitConfObjBeacon.addParameter(name='nFFTPoints', value='200', format='int') +procUnitConfObjBeacon.addParameter(name='pairsList', value='(2,0),(2,1),(2,3),(2,4)', format='pairsList') + +opObj11 = procUnitConfObjBeacon.addOperation(name='IncohInt', optype='other') +opObj11.addParameter(name='n', value='4', format='int') + +opObj11 = procUnitConfObjBeacon.addOperation(name='getBeaconSignal') + +opObj11 = procUnitConfObjBeacon.addOperation(name='BeaconPhase', optype='other') +opObj11.addParameter(name='id', value='301', format='int') +opObj11.addParameter(name='wintitle', value='Beacon Phase', format='str') +opObj11.addParameter(name='xmin', value='0', format='float') +opObj11.addParameter(name='xmax', value='24', format='float') +opObj11.addParameter(name='ymin', value='-180', format='float') +opObj11.addParameter(name='ymax', value='180', format='float') +opObj11.addParameter(name='figpath', value='/media/datos/jasmet50_phase', format='str') + + +print "Escribiendo el archivo XML" +controllerObj.writeXml(filename) +print "Leyendo el archivo XML" +controllerObj.readXml(filename) + +controllerObj.createObjects() +controllerObj.connectObjects() +controllerObj.run()