schain-jc Commit - r15:731909e6d41d · Jicamarca Repository

En IO/Header.py:...

Daniel Valdez -

r15:731909e6d41d

parent child

Context file:

r15:731909e6d41d

Collapse all files

schainpy/IO/TestIO.py created 644 +49 0

@@ -0,0 +1,49
	1	'''
	2	Created on 23/01/2012
	3
	4	@author: danielangelsuarezmunoz
	5	'''
	6	import os
	7	import sys
	8	import datetime
	9	import time
	10
	11	class TestIO():
	12
	13	def __init__(self):
	14	self.setValues()
	15	self.createVoltageObjects()
	16	self.testReadVoltage()
	17	pass
	18
	19	def setValues(self):
	20
	21
	22	self.path = '/Users/danielangelsuarezmunoz/Documents/Projects'
	23	self.startDateTime = datetime.datetime(2007,5,1,17,49,0)
	24	self.endDateTime = datetime.datetime(2007,5,1,18,10,0)
	25
	26	def createVoltageObjects(self):
	27	path = os.path.split(os.getcwd())[0]
	28	sys.path.append(path)
	29
	30	from IO.Voltage import VoltageReader
	31	from Model.Voltage import Voltage
	32
	33	self.voltageModelObj = Voltage()
	34	self.voltageReaderObj = VoltageReader(self.voltageModelObj)
	35	self.voltageReaderObj.setup(self.path, self.startDateTime, self.endDateTime)
	36
	37	def testReadVoltage(self):
	38	while(not(self.voltageReaderObj.noMoreFiles)):
	39
	40	self.voltageReaderObj.getData()
	41	if self.voltageReaderObj.flagResetProcessing:
	42	print 'jump'
	43
	44	if self.voltageReaderObj.flagIsNewBlock:
	45	print 'Block No %04d, Time: %s'%(self.voltageReaderObj.nReadBlocks,
	46	datetime.datetime.fromtimestamp(self.voltageReaderObj.m_BasicHeader.utc))
	47
	48	if __name__ == '__main__':
	49	TestIO() No newline at end of file

schainpy/IO/Header.py +46 -3

             import numpy
             class PROCFLAG:
                 COHERENT_INTEGRATION = numpy.uint32(0x00000001)
                 DECODE_DATA = numpy.uint32(0x00000002)
                 SPECTRA_CALC = numpy.uint32(0x00000004)
                 INCOHERENT_INTEGRATION = numpy.uint32(0x00000008)
                 POST_COHERENT_INTEGRATION = numpy.uint32(0x00000010)
                 SHIFT_FFT_DATA = numpy.uint32(0x00000020)
                 DATATYPE_CHAR = numpy.uint32(0x00000040)
                 DATATYPE_SHORT = numpy.uint32(0x00000080)
                 DATATYPE_LONG = numpy.uint32(0x00000100)
                 DATATYPE_INT64 = numpy.uint32(0x00000200)
                 DATATYPE_FLOAT = numpy.uint32(0x00000400)
                 DATATYPE_DOUBLE = numpy.uint32(0x00000800)
                 DATAARRANGE_CONTIGUOUS_CH = numpy.uint32(0x00001000)
                 DATAARRANGE_CONTIGUOUS_H = numpy.uint32(0x00002000)
                 DATAARRANGE_CONTIGUOUS_P = numpy.uint32(0x00004000)
                 SAVE_CHANNELS_DC = numpy.uint32(0x00008000)
                 DEFLIP_DATA = numpy.uint32(0x00010000)
                 DEFINE_PROCESS_CODE = numpy.uint32(0x00020000)
                 ACQ_SYS_NATALIA = numpy.uint32(0x00040000)
                 ACQ_SYS_ECHOTEK = numpy.uint32(0x00080000)
                 ACQ_SYS_ADRXD = numpy.uint32(0x000C0000)
                 ACQ_SYS_JULIA = numpy.uint32(0x00100000)
                 ACQ_SYS_XXXXXX = numpy.uint32(0x00140000)
                 EXP_NAME_ESP = numpy.uint32(0x00200000)
                 CHANNEL_NAMES_ESP = numpy.uint32(0x00400000)
                 OPERATION_MASK = numpy.uint32(0x0000003F)
                 DATATYPE_MASK = numpy.uint32(0x00000FC0)
                 DATAARRANGE_MASK = numpy.uint32(0x00007000)
                 ACQ_SYS_MASK = numpy.uint32(0x001C0000)
             class BasicHeader:
                 size = 0
                 version = 0
                 dataBlock = 0
                 utc = 0
                 miliSecond = 0
                 timeZone = 0
                 dstFlag = 0
                 errorCount = 0
                 struct = numpy.dtype([
                                       ('nSize','<u4'),
                                       ('nVersion','<u2'),
                                       ('nDataBlockId','<u4'),
                                       ('nUtime','<u4'),
                                       ('nMilsec','<u2'),
                                       ('nTimezone','<i2'),
                                       ('nDstflag','<i2'),
                                       ('nErrorCount','<u4')
                                       ])
                 def __init__(self):
                     pass
                 def read(self, fp):
                     header = numpy.fromfile(fp, self.struct,1)
                     self.size = header['nSize'][0]
                     self.version = header['nVersion'][0]
                     self.dataBlock = header['nDataBlockId'][0]
                     self.utc = header['nUtime'][0]
                     self.miliSecond = header['nMilsec'][0]
                     self.timeZone = header['nTimezone'][0]
                     self.dstFlag = header['nDstflag'][0]
                     self.errorCount = header['nErrorCount'][0]
                     return 1
                 def copy(self):
                     obj = BasicHeader()
                     obj.size = self.size
+                    obj.version = self.version
+                    obj.dataBlock = self.dataBlock
+                    obj.utc = self.utc
+                    obj.miliSecond = self.miliSecond
+                    obj.timeZone = self.timeZone
+                    obj.dstFlag = self.dstFlag
+                    obj.errorCount = self.errorCount
                     return obj
             class SystemHeader:
                 size = 0
                 numSamples = 0
                 numProfiles = 0
                 numChannels = 0
                 adcResolution = 0
                 pciDioBusWidth = 0
                 struct = numpy.dtype([
                                     ('nSize','<u4'),
                                     ('nNumSamples','<u4'),
                                     ('nNumProfiles','<u4'),
                                     ('nNumChannels','<u4'),
                                     ('nADCResolution','<u4'),
                                     ('nPCDIOBusWidth','<u4'),
                                     ])
                 def __init__(self):
                     pass
                 def read(self, fp):
                     header = numpy.fromfile(fp,self.struct,1)
                     self.size = header['nSize'][0]
                     self.numSamples = header['nNumSamples'][0]
                     self.numProfiles = header['nNumProfiles'][0]
                     self.numChannels = header['nNumChannels'][0]
                     self.adcResolution = header['nADCResolution'][0]
                     self.pciDioBusWidth = header['nPCDIOBusWidth'][0]
                     return 1
                 def copy(self):
                     obj = SystemHeader()
                     obj.size = self.size
+                    obj.numSamples = self.numSamples
+                    obj.numProfiles = self.numProfiles
+                    obj.numChannels = self.numChannels
+                    obj.adcResolution = self.adcResolution
+                    self.pciDioBusWidth = self.pciDioBusWidth
                     return obj
             class RadarControllerHeader:
                 size = 0
                 expType = 0
                 nTx = 0
                 ipp = 0
                 txA = 0
                 txB = 0
                 numWindows = 0
                 numTaus = 0
                 codeType = 0
                 line6Function = 0
                 line5Fuction = 0
                 fClock = 0
                 prePulseBefore = 0
                 prePulserAfter = 0
                 rangeIpp = 0
                 rangeTxA = 0
                 rangeTxB = 0
                 struct = numpy.dtype([
                                      ('nSize','<u4'),
                                      ('nExpType','<u4'),
                                      ('nNTx','<u4'),
                                      ('fIpp','<f4'),
                                      ('fTxA','<f4'),
                                      ('fTxB','<f4'),
                                      ('nNumWindows','<u4'),
                                      ('nNumTaus','<u4'),
                                      ('nCodeType','<u4'),
                                      ('nLine6Function','<u4'),
                                      ('nLine5Function','<u4'),
                                      ('fClock','<f4'),
                                      ('nPrePulseBefore','<u4'),
                                      ('nPrePulseAfter','<u4'),
                                      ('sRangeIPP','<a20'),
                                      ('sRangeTxA','<a20'),
                                      ('sRangeTxB','<a20'),
                                      ])
                 def __init__(self):
                     pass
                 def read(self, fp):
                     header = numpy.fromfile(fp,self.struct,1)
                     self.size = header['nSize'][0]
                     self.expType = header['nExpType'][0]
                     self.nTx = header['nNTx'][0]
                     self.ipp = header['fIpp'][0]
                     self.txA = header['fTxA'][0]
                     self.txB = header['fTxB'][0]
                     self.numWindows = header['nNumWindows'][0]
                     self.numTaus = header['nNumTaus'][0]
                     self.codeType = header['nCodeType'][0]
                     self.line6Function = header['nLine6Function'][0]
                     self.line5Fuction = header['nLine5Function'][0]
                     self.fClock = header['fClock'][0]
                     self.prePulseBefore = header['nPrePulseBefore'][0]
                     self.prePulserAfter = header['nPrePulseAfter'][0]
                     self.rangeIpp = header['sRangeIPP'][0]
                     self.rangeTxA = header['sRangeTxA'][0]
                     self.rangeTxB = header['sRangeTxB'][0]
                     # jump Dynamic Radar Controller Header
                     jumpHeader = self.size - 116
                     fp.seek(fp.tell() + jumpHeader)
                     return 1
                 def copy(self):
                     obj = RadarControllerHeader()
                     obj.size = self.size
+                    obj.expType = self.expType
+                    obj.nTx = self.nTx
+                    obj.ipp = self.ipp
+                    obj.txA = self.txA
+                    obj.txB = self.txB
+                    obj.numWindows = self.numWindows
+                    obj.numTaus = self.numTaus
+                    obj.codeType = self.codeType
+                    obj.line6Function = self.line6Function
+                    obj.line5Fuction = self.line5Fuction
+                    obj.fClock = self.fClock
+                    obj.prePulseBefore = self.prePulseBefore
+                    obj.prePulserAfter = self.prePulserAfter
+                    obj.rangeIpp = self.rangeIpp
+                    obj.rangeTxA = self.rangeTxA
+                    obj.rangeTxB = self.rangeTxB
                     return obj
             class ProcessingHeader:
                 size = 0
                 dataType = 0
                 blockSize = 0
                 profilesPerBlock = 0
                 dataBlocksPerFile = 0
                 numWindows = 0
                 processFlags = 0
                 coherentInt = 0
                 incoherentInt = 0
                 totalSpectra = 0
                 struct = numpy.dtype([
                                        ('nSize','<u4'),
                                        ('nDataType','<u4'),
                                        ('nSizeOfDataBlock','<u4'),
                                        ('nProfilesperBlock','<u4'),
                                        ('nDataBlocksperFile','<u4'),
                                        ('nNumWindows','<u4'),
                                        ('nProcessFlags','<u4'),
                                        ('nCoherentIntegrations','<u4'),
                                        ('nIncoherentIntegrations','<u4'),
                                        ('nTotalSpectra','<u4')
                                        ])
                 samplingWindow = 0
                 structSamplingWindow = numpy.dtype([('h0','<f4'),('dh','<f4'),('nsa','<u4')])
                 numHeights = 0
                 firstHeight = 0
                 deltaHeight = 0
                 samplesWin = 0
                 spectraComb = 0
                 numCode = 0
                 codes = 0
                 numBaud = 0
                 def __init__(self):
                     pass
                 def read(self, fp):
                     header = numpy.fromfile(fp,self.struct,1)
                     self.size = header['nSize'][0]
                     self.dataType = header['nDataType'][0]
                     self.blockSize = header['nSizeOfDataBlock'][0]
                     self.profilesPerBlock = header['nProfilesperBlock'][0]
                     self.dataBlocksPerFile = header['nDataBlocksperFile'][0]
                     self.numWindows = header['nNumWindows'][0]
                     self.processFlags = header['nProcessFlags']
                     self.coherentInt = header['nCoherentIntegrations'][0]
                     self.incoherentInt = header['nIncoherentIntegrations'][0]
                     self.totalSpectra = header['nTotalSpectra'][0]
                     self.samplingWindow = numpy.fromfile(fp,self.structSamplingWindow,self.numWindows)
                     self.numHeights = numpy.sum(self.samplingWindow['nsa'])
                     self.firstHeight = self.samplingWindow['h0']
                     self.deltaHeight = self.samplingWindow['dh']
                     self.samplesWin = self.samplingWindow['nsa']
                     self.spectraComb = numpy.fromfile(fp,'u1',2*self.totalSpectra)
                     if self.processFlags & PROCFLAG.DEFINE_PROCESS_CODE == PROCFLAG.DEFINE_PROCESS_CODE:
                         self.numCode = numpy.fromfile(fp,'<u4',1)
                         self.numBaud = numpy.fromfile(fp,'<u4',1)
                         self.codes = numpy.fromfile(fp,'<f4',self.numCode*self.numBaud).reshape(self.numBaud,self.numCode)
                     return 1
                 def copy(self):
                     obj = ProcessingHeader()
                     obj.size = self.size
+                    obj.dataType = self.dataType
+                    obj.blockSize = self.blockSize
+                    obj.profilesPerBlock = self.profilesPerBlock
+                    obj.dataBlocksPerFile = self.dataBlocksPerFile
+                    obj.numWindows = self.numWindows
+                    obj.processFlags = self.processFlags
+                    obj.coherentInt = self.coherentInt
+                    obj.incoherentInt = self.incoherentInt
+                    obj.totalSpectra = self.totalSpectra
+                    obj.samplingWindow = self.samplingWindow
+                    obj.numHeights = self.numHeights
+                    obj.firstHeight = self.firstHeight
+                    obj.deltaHeight = self.deltaHeight
+                    obj.samplesWin = self.samplesWin
+                    obj.spectraComb = self.spectraComb
+                    obj.numCode = self.numCode
+                    obj.numBaud = self.numBaud
+                    obj.codes = self.codes
                     return obj

schainpy/IO/Voltage.py +17 -12

             '''
             Created on 23/01/2012
             @author: danielangelsuarezmunoz
             '''
             import os, sys
             import numpy
             import glob
             import fnmatch
             import time
             import datetime
-            from Header import *
+            from IO.Header import *
-            from Data import DataReader
+            from IO.Data import DataReader
-            from Data import DataWriter
+            from IO.Data import DataWriter
             path = os.path.split(os.getcwd())[0]
             sys.path.append(path)
             from Model.Voltage import Voltage
             class VoltageReader(DataReader):
                 __idFile = None
                 __fp = None
                 __startDateTime = None
                 __endDateTime = None
                 __dataType = None
                 __fileSizeByHeader = 0
                 __pathList = []
                 filenameList = []
                 __lastUTTime = 0
                 __maxTimeStep = 5
                 __flagIsNewFile = 0
+                __ippSeconds = 0
                 flagResetProcessing = 0
                 flagIsNewBlock = 0
                 noMoreFiles = 0
                 nReadBlocks = 0
                 online = 0
                 filename = None
                 fileSize = None
                 firstHeaderSize = 0
                 basicHeaderSize = 24
                 m_BasicHeader = BasicHeader()
                 m_SystemHeader = SystemHeader()
                 m_RadarControllerHeader = RadarControllerHeader()
                 m_ProcessingHeader = ProcessingHeader()
                 m_Voltage = None
                 __buffer = 0
                 __buffer_id = 9999
                 def __init__(self, m_Voltage = None):
                     if m_Voltage == None:
                         m_Voltage = Voltage()
                     self.m_Voltage = m_Voltage
                 def __rdSystemHeader(self,fp=None):
                     if fp == None:
                         fp = self.__fp
                     self.m_SystemHeader.read(fp)
                 def __rdRadarControllerHeader(self,fp=None):
                     if fp == None:
                         fp = self.__fp
                     self.m_RadarControllerHeader.read(fp)
                 def __rdProcessingHeader(self,fp=None):
                     if fp == None:
                         fp = self.__fp
                     self.m_ProcessingHeader.read(fp)
                 def __searchFiles(self,path, startDateTime, endDateTime, set=None, expLabel = "", ext = ".r"):
                     print "Searching files ..."
                     startUtSeconds = time.mktime(startDateTime.timetuple())
                     endUtSeconds = time.mktime(endDateTime.timetuple())
             #        startYear = startDateTime.timetuple().tm_year
             #        endYear = endDateTime.timetuple().tm_year
             #
             #        startDoy = startDateTime.timetuple().tm_yday
             #        endDoy = endDateTime.timetuple().tm_yday
             #
             #        yearRange = range(startYear,endYear+1)
             #
             #        doyDoubleList = []
             #        if startYear == endYear:
             #            doyList = range(startDoy,endDoy+1)
             #        else:
             #            for year in yearRange:
             #                if (year == startYear):
             #                    doyDoubleList.append(range(startDoy,365+1))
             #                elif (year == endYear):
             #                    doyDoubleList.append(range(1,endDoy+1))
             #                else:
             #                    doyDoubleList.append(range(1,365+1))
             #            doyList = []
             #            for list in doyDoubleList:
             #                doyList = doyList + list
             #
             #        dirList = []
             #        for thisPath in os.listdir(path):
             #            if os.path.isdir(os.path.join(path,thisPath)):
             #                #dirList.append(os.path.join(path,thisPath))
             #                dirList.append(thisPath)
             #
             #        pathList = []
             #        pathDict = {}
             #        for year in yearRange:
             #            for doy in doyList:
             #                match = fnmatch.filter(dirList, 'D' + '%4.4d%3.3d' % (year,doy))
             #                if len(match) == 0:
             #                    match = fnmatch.filter(dirList, 'd' + '%4.4d%3.3d' % (year,doy))
             #                    if len(match) == 0: continue
             #                if expLabel == '':
             #                    pathList.append(os.path.join(path,match[0]))
             #                    pathDict.setdefault(os.path.join(path,match[0]))
             #                    pathDict[os.path.join(path,match[0])] = []
             #                else:
             #                    pathList.append(os.path.join(path,os.path.join(match[0],expLabel)))
             #                    pathDict.setdefault(os.path.join(path,os.path.join(match[0],expLabel)))
             #                    pathDict[os.path.join(path,os.path.join(match[0],expLabel))] = []
                     dirList = []
                     for thisPath in os.listdir(path):
                         if os.path.isdir(os.path.join(path,thisPath)):
                             dirList.append(thisPath)
                     pathList = []
                     thisDateTime = startDateTime
                     while(thisDateTime <= endDateTime):
                         year = thisDateTime.timetuple().tm_year
                         doy = thisDateTime.timetuple().tm_yday
                         match = fnmatch.filter(dirList, '?' + '%4.4d%3.3d' % (year,doy))
                         if len(match) == 0:
                             thisDateTime += datetime.timedelta(1)
                             continue
                         pathList.append(os.path.join(path,match[0],expLabel))
                         thisDateTime += datetime.timedelta(1)
                     filenameList = []
                     for thisPath in pathList:
                         fileList = glob.glob1(thisPath, "*%s" %ext)
                         #pathDict[thisPath].append(fileList)
                         fileList.sort()
                         for file in fileList:
                             filename = os.path.join(thisPath,file)
                             if self.isThisFileinRange(filename, startUtSeconds, endUtSeconds):
                                 filenameList.append(filename)
                     self.filenameList = filenameList
                     return pathList, filenameList
                 def isThisFileinRange(self, filename, startUTSeconds=None, endUTSeconds=None):
                     try:
                         fp = open(filename,'rb')
                     except:
                         raise IOError, "The file %s can't be opened" %(filename)
                     if startUTSeconds==None:
                         startUTSeconds = self.startUTCSeconds
                     if endUTSeconds==None:
                         endUTSeconds = self.endUTCSeconds
                     m_BasicHeader = BasicHeader()
                     if not(m_BasicHeader.read(fp)):
                         return 0
                     fp.close()
                     if not ((startUTSeconds <= m_BasicHeader.utc) and (endUTSeconds >= m_BasicHeader.utc)):
                         return 0
                     return 1
                 def __readBasicHeader(self, fp=None):
                     if fp == None:
                         fp = self.__fp
                     self.m_BasicHeader.read(fp)
                 def __readFirstHeader(self):
                     self.__readBasicHeader()
                     self.__rdSystemHeader()
                     self.__rdRadarControllerHeader()
                     self.__rdProcessingHeader()
                     self.firstHeaderSize = self.m_BasicHeader.size
                     data_type=int(numpy.log2((self.m_ProcessingHeader.processFlags & PROCFLAG.DATATYPE_MASK))-numpy.log2(PROCFLAG.DATATYPE_CHAR))
                     if data_type == 0:
                         tmp=numpy.dtype([('real','<i1'),('imag','<i1')])
                     elif data_type == 1:
                         tmp=numpy.dtype([('real','<i2'),('imag','<i2')])
                     elif data_type == 2:
                         tmp=numpy.dtype([('real','<i4'),('imag','<i4')])
                     elif data_type == 3:
                         tmp=numpy.dtype([('real','<i8'),('imag','<i8')])
                     elif data_type == 4:
                         tmp=numpy.dtype([('real','<f4'),('imag','<f4')])
                     elif data_type == 5:
                         tmp=numpy.dtype([('real','<f8'),('imag','<f8')])
                     else:
                         raise ValueError, 'Data type was not defined'
                     self.__dataType = tmp
                     self.__fileSizeByHeader = self.m_ProcessingHeader.dataBlocksPerFile * self.m_ProcessingHeader.blockSize + self.firstHeaderSize + self.basicHeaderSize*(self.m_ProcessingHeader.dataBlocksPerFile - 1)
+                    c=3E8
+                    self.__ippSeconds = 2*1000*self.m_RadarControllerHeader.ipp/c
                 def __setNextFileOnline(self):
                     return 0
                 def __setNextFileOffline(self):
                     idFile = self.__idFile
                     while(True):
                         idFile += 1
                         if not(idFile < len(self.filenameList)):
                             self.noMoreFiles = 1
                             return 0
                         filename = self.filenameList[idFile]
                         fileSize = os.path.getsize(filename)
                         fp = open(filename,'rb')
                         currentSize = fileSize - fp.tell()
                         neededSize = self.m_ProcessingHeader.blockSize + self.firstHeaderSize
                         if (currentSize < neededSize):
                             continue
                         break
                     self.__flagIsNewFile = 1
                     self.__idFile = idFile
                     self.filename = filename
                     self.fileSize = fileSize
                     self.__fp = fp
                     print 'Setting the file: %s'%self.filename
                     return 1
                 def __setNextFile(self):
                     if self.online:
                         return self.__setNextFileOnline()
                     else:
                         return self.__setNextFileOffline()
                 def __setNewBlock(self):
                     if self.__fp == None:
                         return 0
                     if self.__flagIsNewFile:
                         return 1
                     currentSize = self.fileSize - self.__fp.tell()
                     neededSize = self.m_ProcessingHeader.blockSize + self.basicHeaderSize
                     # Bloque Completo
                     if (currentSize >= neededSize):
                         self.__readBasicHeader()
                         return 1
                     if not(self.__setNextFile()):
                         return 0
                     self.__readFirstHeader()
                     deltaTime = self.m_BasicHeader.utc - self.__lastUTTime # check this
                     self.flagResetProcessing = 0
                     if deltaTime > self.__maxTimeStep:
                         self.flagResetProcessing = 1
                         self.nReadBlocks = 0
                     return 1
                 def __readBlock(self):
                     """Lee el bloque de datos desde la posicion actual del puntero del archivo y
                     actualiza todos los parametros relacionados al bloque de datos (data, time,
                     etc). La data leida es almacenada en el buffer y el contador de datos leidos es
                     seteado a 0
                     """
                     pts2read = self.m_ProcessingHeader.profilesPerBlock*self.m_ProcessingHeader.numHeights*self.m_SystemHeader.numChannels
                     data = numpy.fromfile(self.__fp,self.__dataType,pts2read)
                     data = data.reshape((self.m_ProcessingHeader.profilesPerBlock, self.m_ProcessingHeader.numHeights, self.m_SystemHeader.numChannels))
                     self.__flagIsNewFile = 0
                     self.flagIsNewBlock = 1
                     self.nReadBlocks += 1
                     self.__buffer = data
                     self.__buffer_id = 0
                 def readNextBlock(self):
                     if not(self.__setNewBlock()):
                         return 0
                     self.__readBlock()
                     self.__lastUTTime = self.m_BasicHeader.utc
                     return 1
                 def __hasNotDataInBuffer(self):
                     if self.__buffer_id >= self.m_ProcessingHeader.profilesPerBlock:
                         return 1
                     return 0
                 def getData(self):
                     """Obtiene un unidad de datos del buffer de lectura y es copiada a la clase "Voltage"
                     con todos los parametros asociados a este. cuando no hay datos en el buffer de
                     lectura es necesario hacer una nueva lectura de los bloques de datos usando "readNextBlock"
                     """
                     self.flagResetProcessing = 0
                     self.flagIsNewBlock = 0
                     if self.__hasNotDataInBuffer():
                         self.readNextBlock()
                     if self.noMoreFiles == 1:
                         print 'Process finished'
                         return None
                     data = self.__buffer[self.__buffer_id,:,:]
-                    time = 111
+                    time = self.m_BasicHeader.utc + self.__buffer_id*self.__ippSeconds
-            #        self.m_Voltage.data = data
-            #        self.m_Voltage.timeProfile = time
+                    self.m_Voltage.m_BasicHeader = self.m_BasicHeader.copy()
-            #        self.m_Voltage.m_BasicHeader = self.m_BasicHeader.copy()
+                    self.m_Voltage.m_ProcessingHeader = self.m_ProcessingHeader.copy()
-            #        self.m_Voltage.m_ProcessingHeader = self.m_ProcessingHeader.copy()
+                    self.m_Voltage.m_RadarControllerHeader = self.m_RadarControllerHeader.copy()
-            #        self.m_Voltage.m_RadarControllerHeader = self.m_RadarControllerHeader.copy()
+                    self.m_Voltage.m_SystemHeader = self.m_SystemHeader.copy()
-            #        self.m_Voltage.m_SystemHeader = self.m_systemHeader.copy()
+                    self.m_Voltage.m_BasicHeader.utc = time
+                    self.m_Voltage.data = data
                     self.__buffer_id += 1
                     #call setData - to Data Object
                     return data
                 def setup(self, path, startDateTime, endDateTime, set=None, expLabel = "", ext = ".r", online = 0):
                     if online == 0:
                         pathList, filenameList = self.__searchFiles(path, startDateTime, endDateTime, set, expLabel, ext)
                         if len(filenameList) == 0:
                             self.__fp = None
                             self.noMoreFiles = 1
                             print 'Do not exist files in range: %s - %s'%(startDateTime.ctime(), endDateTime.ctime())
                             return 0
             #        for thisFile in filenameList:
             #            print thisFile
                     self.__idFile = -1
                     if not(self.__setNextFile()):
                         print "No more files"
                         return 0
                     self.__readFirstHeader()
                     self.startUTCSeconds = time.mktime(startDateTime.timetuple())
                     self.endUTCSeconds = time.mktime(endDateTime.timetuple())
                     self.startYear = startDateTime.timetuple().tm_year
                     self.endYear = endDateTime.timetuple().tm_year
                     self.startDoy = startDateTime.timetuple().tm_yday
                     self.endDoy = endDateTime.timetuple().tm_yday
                     #call fillHeaderValues() - to Data Object
                     self.__pathList = pathList
                     self.filenameList = filenameList
                     self.online = online
             class VoltageWriter(DataWriter):
                 m_BasicHeader= BasicHeader()
                 m_SystemHeader = SystemHeader()
                 m_RadarControllerHeader = RadarControllerHeader()
                 m_ProcessingHeader = ProcessingHeader()
                 m_Voltage = None
                 def __init__(self, m_Voltage):
                     self.m_Voltage = m_Voltage
   No newline at end of file

schainpy/Model/Voltage.py +11 -1

             '''
             Created on Feb 7, 2012
             @author: roj-idl71
             '''
             import os, sys
             path = os.path.split(os.getcwd())[0]
             sys.path.append(path)
             from Model.Data import Data
             from IO.Header import *
             class Voltage(Data):
                 '''
                 classdocs
                 '''
                 m_RadarControllerHeader= RadarControllerHeader()
                 m_ProcessingHeader= ProcessingHeader()
                 m_SystemHeader= SystemHeader()
                 m_BasicHeader= BasicHeader()
+                data = None
+                noData = True
                 def __init__(self):
                     '''
                     Constructor
                     '''
                     pass
                 def copy(self):
-                    pass
+                    obj = Voltage()
+                    obj.m_BasicHeader = self.m_BasicHeader.copy()
+                    obj.m_SystemHeader = self.m_SystemHeader.copy()
+                    obj.m_RadarControllerHeader = self.m_RadarControllerHeader.copy()
+                    obj.m_ProcessingHeader = self.m_ProcessingHeader.copy()
+                    return obj
   No newline at end of file

General Comments 0

Write
Preview

You need to be logged in to leave comments. Login now

No TODOs yet

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages