##// END OF EJS Templates
fix pulse pair spectra
avaldez -
r1283:47c7b4917f5b
parent child
Show More
@@ -324,8 +324,8 class WpowerPlot_(Figure):
324 zmin : None,
324 zmin : None,
325 zmax : None
325 zmax : None
326 """
326 """
327 print("***************PLOTEO******************")
327 #print("***************PLOTEO******************")
328 print("DATAOUT SHAPE : ",dataOut.data.shape)
328 #print("DATAOUT SHAPE : ",dataOut.data.shape)
329 if dataOut.flagNoData:
329 if dataOut.flagNoData:
330 return dataOut
330 return dataOut
331
331
@@ -344,7 +344,7 class WpowerPlot_(Figure):
344 channelIndexList.append(dataOut.channelList.index(channel))
344 channelIndexList.append(dataOut.channelList.index(channel))
345
345
346
346
347 print("channelIndexList",channelIndexList)
347 #print("channelIndexList",channelIndexList)
348 if normFactor is None:
348 if normFactor is None:
349 factor = dataOut.normFactor
349 factor = dataOut.normFactor
350 else:
350 else:
@@ -364,7 +364,7 class WpowerPlot_(Figure):
364 ylabel = "Range (km)"
364 ylabel = "Range (km)"
365
365
366 y = dataOut.getHeiRange()
366 y = dataOut.getHeiRange()
367 print("factor",factor)
367 #print("factor",factor)
368
368
369 z = dataOut.data/factor # dividido /factor
369 z = dataOut.data/factor # dividido /factor
370 z = numpy.where(numpy.isfinite(z), z, numpy.NAN)
370 z = numpy.where(numpy.isfinite(z), z, numpy.NAN)
@@ -88,6 +88,7 class SpectraProc(ProcessingUnit):
88 # calculo de self-spectra
88 # calculo de self-spectra
89 fft_volt = numpy.fft.fftshift(fft_volt, axes=(1,))
89 fft_volt = numpy.fft.fftshift(fft_volt, axes=(1,))
90 spc = fft_volt * numpy.conjugate(fft_volt)
90 spc = fft_volt * numpy.conjugate(fft_volt)
91 #print("spcch0",spc[0])
91 spc = spc.real
92 spc = spc.real
92
93
93 blocksize = 0
94 blocksize = 0
@@ -1154,6 +1155,7 class PulsePair(Operation):
1154 self.__buffer = data*numpy.conjugate(data)
1155 self.__buffer = data*numpy.conjugate(data)
1155 self.__bufferV = data[:,(self.__nProf-1):,:]*numpy.conjugate(data[:,1:,:])
1156 self.__bufferV = data[:,(self.__nProf-1):,:]*numpy.conjugate(data[:,1:,:])
1156 self.__profIndex = self.n
1157 self.__profIndex = self.n
1158 #print("spcch0",self.__buffer)
1157 return
1159 return
1158
1160
1159 def pushData(self):
1161 def pushData(self):
@@ -1162,9 +1164,13 class PulsePair(Operation):
1162 data_IV = numpy.zeros((self.__nch,self.__nHeis))
1164 data_IV = numpy.zeros((self.__nch,self.__nHeis))
1163
1165
1164 for i in range(self.__nch):
1166 for i in range(self.__nch):
1165 data_I[i,:] = numpy.sum(numpy.sum(self.__buffer[i],axis=0),axis=0)/self.n
1167 data_I[i,:] = numpy.sum(self.__buffer[i],axis=0)/self.n
1166 data_IV[i,:] = numpy.sum(numpy.sum(self.__bufferV[i],axis=0),axis=0)/(self.n-1)
1168 data_IV[i,:] = numpy.sum(self.__bufferV[i],axis=0)/(self.n-1)
1167
1169 ##print("******")
1170 #print("data_I",data_I[0])
1171 #print(self.__buffer.shape)
1172 #a=numpy.average(self.__buffer,axis=1)
1173 #print("average", a)
1168 n = self.__profIndex
1174 n = self.__profIndex
1169 ####data_intensity = numpy.sum(numpy.sum(self.__buffer,axis=0),axis=0)/self.n
1175 ####data_intensity = numpy.sum(numpy.sum(self.__buffer,axis=0),axis=0)/self.n
1170 #print("data_intensity push data",data_intensity.shape)
1176 #print("data_intensity push data",data_intensity.shape)
General Comments 0
You need to be logged in to leave comments. Login now