| @@ -1387,7 +1387,7 class VoltageReader(JRODataReader): | |||||
| 1387 | if not( self.readNextBlock() ): |
|
1387 | if not( self.readNextBlock() ): | |
| 1388 | return 0 |
|
1388 | return 0 | |
| 1389 |
|
1389 | |||
| 1390 |
self.dataOut.dtype = |
|
1390 | self.dataOut.dtype = self.dtype | |
| 1391 |
|
1391 | |||
| 1392 | self.dataOut.nProfiles = self.processingHeaderObj.profilesPerBlock |
|
1392 | self.dataOut.nProfiles = self.processingHeaderObj.profilesPerBlock | |
| 1393 |
|
1393 | |||
| @@ -1527,7 +1527,7 class VoltageWriter(JRODataWriter): | |||||
| 1527 | self.datablock = numpy.zeros((self.systemHeaderObj.nChannels, |
|
1527 | self.datablock = numpy.zeros((self.systemHeaderObj.nChannels, | |
| 1528 | self.processingHeaderObj.profilesPerBlock, |
|
1528 | self.processingHeaderObj.profilesPerBlock, | |
| 1529 | self.processingHeaderObj.nHeights), |
|
1529 | self.processingHeaderObj.nHeights), | |
| 1530 | dtype=numpy.dtype('complex')) |
|
1530 | dtype=numpy.dtype('complex64')) | |
| 1531 |
|
1531 | |||
| 1532 |
|
1532 | |||
| 1533 | def writeBlock(self): |
|
1533 | def writeBlock(self): | |
| @@ -2045,7 +2045,7 class SpectraReader(JRODataReader): | |||||
| 2045 |
|
2045 | |||
| 2046 | self.dataOut.flagNoData = False |
|
2046 | self.dataOut.flagNoData = False | |
| 2047 |
|
2047 | |||
| 2048 |
self.dataOut.dtype = |
|
2048 | self.dataOut.dtype = self.dtype | |
| 2049 |
|
2049 | |||
| 2050 | # self.dataOut.nChannels = self.nRdChannels |
|
2050 | # self.dataOut.nChannels = self.nRdChannels | |
| 2051 |
|
2051 | |||
| @@ -2225,6 +2225,7 class SpectraWriter(JRODataWriter): | |||||
| 2225 | if not( self.processingHeaderObj.shif_fft ): |
|
2225 | if not( self.processingHeaderObj.shif_fft ): | |
| 2226 | spc = numpy.roll( spc, self.processingHeaderObj.profilesPerBlock/2, axis=2 ) #desplaza a la derecha en el eje 2 determinadas posiciones |
|
2226 | spc = numpy.roll( spc, self.processingHeaderObj.profilesPerBlock/2, axis=2 ) #desplaza a la derecha en el eje 2 determinadas posiciones | |
| 2227 | data = spc.reshape((-1)) |
|
2227 | data = spc.reshape((-1)) | |
|
|
2228 | data = data.astype(self.dtype[0]) | |||
| 2228 | data.tofile(self.fp) |
|
2229 | data.tofile(self.fp) | |
| 2229 |
|
2230 | |||
| 2230 | if self.data_cspc != None: |
|
2231 | if self.data_cspc != None: | |
| @@ -680,7 +680,7 class SpectraProc(ProcessingUnit): | |||||
| 680 | self.dataOut.systemHeaderObj = self.dataIn.systemHeaderObj.copy() |
|
680 | self.dataOut.systemHeaderObj = self.dataIn.systemHeaderObj.copy() | |
| 681 | self.dataOut.channelList = self.dataIn.channelList |
|
681 | self.dataOut.channelList = self.dataIn.channelList | |
| 682 | self.dataOut.heightList = self.dataIn.heightList |
|
682 | self.dataOut.heightList = self.dataIn.heightList | |
| 683 |
self.dataOut.dtype = |
|
683 | self.dataOut.dtype = numpy.dtype([('real','<f4'),('imag','<f4')]) | |
| 684 | # self.dataOut.nHeights = self.dataIn.nHeights |
|
684 | # self.dataOut.nHeights = self.dataIn.nHeights | |
| 685 | # self.dataOut.nChannels = self.dataIn.nChannels |
|
685 | # self.dataOut.nChannels = self.dataIn.nChannels | |
| 686 | self.dataOut.nBaud = self.dataIn.nBaud |
|
686 | self.dataOut.nBaud = self.dataIn.nBaud | |
| @@ -714,6 +714,7 class SpectraProc(ProcessingUnit): | |||||
| 714 | self.dataOut.flagNoData |
|
714 | self.dataOut.flagNoData | |
| 715 | """ |
|
715 | """ | |
| 716 | fft_volt = numpy.fft.fft(self.buffer,axis=1) |
|
716 | fft_volt = numpy.fft.fft(self.buffer,axis=1) | |
|
|
717 | fft_volt = fft_volt.astype(numpy.dtype('complex')) | |||
| 717 | dc = fft_volt[:,0,:] |
|
718 | dc = fft_volt[:,0,:] | |
| 718 |
|
719 | |||
| 719 | #calculo de self-spectra |
|
720 | #calculo de self-spectra | |
General Comments 0
You need to be logged in to leave comments.
Login now