schain Commit - r1027:f11c2da27e74 · Jicamarca Repository

07/09/2017

ebocanegra -

r1027:f11c2da27e74

parent child

Context file:

r1027:f11c2da27e74

Collapse all files

schainpy/model/graphics/jroplot_spectra.py +43 0

             import datetime
             import numpy
+            import matplotlib.pyplot as plt
             from figure import Figure, isRealtime, isTimeInHourRange
             from plotting_codes import *
+            from matplotlib.pyplot import savefig
             class SpectraPlot(Figure):
                     self.__xfilter_ena = False
                     self.__yfilter_ena = False
+                    self.indice=1
                 def getSubplots(self):
                     thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[0])
                     title = wintitle + " Spectra"
+                    print 'len de X',len(x), numpy.shape(x), 'len de spc line',len(dataOut.data_spc[1,:,15]), numpy.shape(dataOut.data_spc)
+                    print 'Altura:', y[0], y[1], y[13], y[14], y[10]
+                    #a=z[1,:,15]
+            #         fig = plt.figure(10+self.indice)
+            #         plt.plot( x[0:128], zdB[0,:,10] )
+            #         plt.axis([-12, 12, 15, 50])
+            #         plt.title(" %s" %(  '%s %s'%(thisDatetime.strftime("%Y/%m/%d"),thisDatetime.strftime("%H:%M:%S")))   )
+            #         plt.ylabel('Intensidad [dB]')
+            #         plt.xlabel('Velocidad [m/s]')
+            #         fig.savefig('/home/erick/Documents/Pics/to{}.png'.format(self.indice))
+            #
+            #         plt.show()
+            #
+            #         self.indice=self.indice+1
                     if ((dataOut.azimuth!=None) and (dataOut.zenith!=None)):
                         title = title + '_' + 'azimuth,zenith=%2.2f,%2.2f'%(dataOut.azimuth, dataOut.zenith)
                               ftp=ftp,
                               wr_period=wr_period,
                               thisDatetime=thisDatetime)
             class CrossSpectraPlot(Figure):
                     self.EXP_CODE = None
                     self.SUB_EXP_CODE = None
                     self.PLOT_POS = None
+                    self.indice=0
                 def getSubplots(self):
                         self.isConfig = True
                     self.setWinTitle(title)
                     for i in range(self.nplots):
                         pair = dataOut.pairsList[pairsIndexList[i]]
                         coherence = numpy.abs(coherenceComplex)
             #            phase = numpy.arctan(-1*coherenceComplex.imag/coherenceComplex.real)*180/numpy.pi
                         phase = numpy.arctan2(coherenceComplex.imag, coherenceComplex.real)*180/numpy.pi
+            #             print 'FASE', numpy.shape(phase), y[10]
+            #             fig = plt.figure(10+self.indice)
+            #             plt.plot( x[0:128],phase[:,10] )
+            #             #plt.axis([-12, 12, 15, 50])
+            #             plt.title("%s" %(  '%s %s, Channel %s'%(thisDatetime.strftime("%Y/%m/%d"),thisDatetime.strftime("%H:%M:%S") , i)))
+            #             plt.ylabel('Desfase [grados]')
+            #             plt.xlabel('Velocidad [m/s]')
+            #             fig.savefig('/home/erick/Documents/Pics/to{}.png'.format(self.indice))
+            #
+            #             plt.show()
+            #             self.indice=self.indice+1
                         title = "Coherence Ch%d * Ch%d" %(pair[0], pair[1])
                         axes0 = self.axesList[i*self.__nsubplots+2]

schainpy/model/graphics/jroplot_spectra.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/model/io/jroIO_bltr.py +1 0

                             #print ' '
                             #print ' '
                             #Height=35
                             for i in range(self.nRdPairs):
                                 chan_index0 = self.dataOut.pairsList[i][0]

schainpy/model/io/jroIO_bltr.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/jroproc_bltr.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/jroproc_parameters.py +44 -37

                     """
                     Constants:
-                    Pt:     Transmission Power               dB
+                    Pt:     Transmission Power               dB        5kW
-                    Gt:     Transmission Gain                dB
+                    Gt:     Transmission Gain                dB        24.7 dB
-                    Gr:     Reception Gain                   dB
+                    Gr:     Reception Gain                   dB        18.5 dB
-                    Lambda: Wavelenght                       m
+                    Lambda: Wavelenght                       m         0.6741 m
-                    aL:      Attenuation loses               dB
+                    aL:     Attenuation loses                dB
                     tauW:   Width of transmission pulse      s
-                    ThetaT: Transmission antenna bean angle  rad
+                    ThetaT: Transmission antenna bean angle  rad       0.1656317 rad
-                    ThetaR: Reception antenna beam angle     rad
+                    ThetaR: Reception antenna beam angle     rad       0.36774087 rad
                     """
                     Numerator = ( (4*numpy.pi)**3 * aL**2 * 16 * numpy.log(2) )
                     data = dataOut.data_pre
                     noise = dataOut.noise
-                    print 'noise',noise
+                    #print 'noise',noise
                     #SNRdB = 10*numpy.log10(dataOut.data_SNR)
                     FirstMoment = numpy.average(dataOut.data_param[:,1,:],0)
                     velocityX=[]
                     velocityY=[]
                     velocityV=[]
+                    PhaseLine=[]
                     dbSNR = 10*numpy.log10(dataSNR)
                     dbSNR = numpy.average(dbSNR,0)
                     for Height in range(nHeights):
-                        [Vzon,Vmer,Vver, GaussCenter]= self.WindEstimation(spc, cspc, pairsList, ChanDist, Height, noise, VelRange, dbSNR[Height], SNRlimit)
+                        [Vzon,Vmer,Vver, GaussCenter, PhaseSlope]= self.WindEstimation(spc, cspc, pairsList, ChanDist, Height, noise, VelRange, dbSNR[Height], SNRlimit)
+                        PhaseLine = numpy.append(PhaseLine, PhaseSlope)
                         if abs(Vzon)<100. and abs(Vzon)> 0.:
                             velocityX=numpy.append(velocityX, Vzon)#Vmag
                         else:
-                            print 'Vzon',Vzon
+                            #print 'Vzon',Vzon
                             velocityX=numpy.append(velocityX, numpy.NaN)
                         if abs(Vmer)<100. and abs(Vmer) > 0.:
                             velocityY=numpy.append(velocityY, Vmer)#Vang
                         else:
-                            print 'Vmer',Vmer
+                            #print 'Vmer',Vmer
                             velocityY=numpy.append(velocityY, numpy.NaN)
                         if dbSNR[Height] > SNRlimit:
             #                 FirstMoment[Height] = numpy.NaN
             #                 velocityX[Height] = numpy.NaN
             #                 velocityY[Height] = numpy.NaN
                     data_output[0]=numpy.array(velocityX)
                     data_output[1]=numpy.array(velocityY)
                     print 'velocityX',data_output[0]
                     print ' '
                     print 'velocityY',data_output[1]
+                    print 'PhaseLine',PhaseLine
                     #print numpy.array(velocityY)
                     print ' '
                     #print 'SNR'
                     dataOut.data_output=data_output
                     return
                     phase=numpy.ones([spc.shape[0],spc.shape[1]])
                     CSPCSamples=numpy.ones([spc.shape[0],spc.shape[1]],dtype=numpy.complex_)
                     coherence=numpy.ones([spc.shape[0],spc.shape[1]])
-                    PhaseSlope=numpy.ones(spc.shape[0])
+                    PhaseSlope=numpy.zeros(spc.shape[0])
                     PhaseInter=numpy.ones(spc.shape[0])
                     xFrec=VelRange
                     print ' '
                     #print 'dataSNR', dbSNR.shape, dbSNR[0,40:120]
-                    print 'SmoothSPC', SmoothSPC.shape, SmoothSPC[0:20]
+                    #print 'SmoothSPC', SmoothSPC.shape, SmoothSPC[0:20]
-                    print 'noise',noise
+                    #print 'noise',noise
-                    print 'zline',zline.shape, zline[0:20]
+                    #print 'zline',zline.shape, zline[0:20]
-                    print 'FactNorm',FactNorm.shape, FactNorm[0:20]
+                    #print 'FactNorm',FactNorm.shape, FactNorm[0:20]
-                    print 'FactNorm suma', numpy.sum(FactNorm)
+                    #print 'FactNorm suma', numpy.sum(FactNorm)
                     for i in range(spc.shape[0]):
                         phase[i] = self.moving_average( numpy.arctan2(CSPCSamples[i].imag, CSPCSamples[i].real),N=1)#*180/numpy.pi
-                    print 'cspcLine', cspcLine.shape, cspcLine[0:20]
+                    #print 'cspcLine', cspcLine.shape, cspcLine[0:20]
-                    print 'CSPCFactor', CSPCFactor#, CSPCFactor[0:20]
+                    #print 'CSPCFactor', CSPCFactor#, CSPCFactor[0:20]
-                    print numpy.sum(ySamples[chan_index0]), numpy.sum(ySamples[chan_index1]), -noise[i]
+                    #print numpy.sum(ySamples[chan_index0]), numpy.sum(ySamples[chan_index1]), -noise[i]
-                    print 'CSPCNorm', CSPCNorm.shape, CSPCNorm[0:20]
+                    #print 'CSPCNorm', CSPCNorm.shape, CSPCNorm[0:20]
-                    print 'CSPCNorm suma', numpy.sum(CSPCNorm)
+                    #print 'CSPCNorm suma', numpy.sum(CSPCNorm)
-                    print 'CSPCSamples', CSPCSamples.shape, CSPCSamples[0,0:20]
+                    #print 'CSPCSamples', CSPCSamples.shape, CSPCSamples[0,0:20]
                     '''****** Getting fij width ******'''
                     meanGauss=sum(xSamples*yMean) / len(xSamples)
                     sigma=sum(yMean*(xSamples-meanGauss)**2) / len(xSamples)
-                    print '****************************'
+                    #print '****************************'
-                    print 'len(xSamples): ',len(xSamples)
+                    #print 'len(xSamples): ',len(xSamples)
-                    print 'yMean: ', yMean.shape, yMean[0:20]
+                    #print 'yMean: ', yMean.shape, yMean[0:20]
-                    print 'ySamples', ySamples.shape, ySamples[0,0:20]
+                    #print 'ySamples', ySamples.shape, ySamples[0,0:20]
-                    print 'xSamples: ',xSamples.shape, xSamples[0:20]
+                    #print 'xSamples: ',xSamples.shape, xSamples[0:20]
-                    print 'meanGauss',meanGauss
+                    #print 'meanGauss',meanGauss
-                    print 'sigma',sigma
+                    #print 'sigma',sigma
                     #if (abs(meanGauss/sigma**2) > 0.0001) : #0.000000001):
-                    if dbSNR > SNRlimit :
+                    if dbSNR > SNRlimit and abs(meanGauss/sigma**2) > 0.0001:
                         try:
                             popt,pcov = curve_fit(self.gaus,xSamples,yMean,p0=[1,meanGauss,sigma])
                     else:
                         Range = numpy.array([0,0])
-                    print ' '
+                    #print ' '
-                    print 'GCpos',GCpos, ( len(xFrec)- len(xFrec)*0.1)
+                    #print 'GCpos',GCpos, ( len(xFrec)- len(xFrec)*0.1)
-                    print 'Rangpos',Rangpos
+                    #print 'Rangpos',Rangpos
                     print 'RANGE: ', Range
                     FrecRange=xFrec[Range[0]:Range[1]]
                         Vmag=numpy.sqrt(Vzon**2+Vmer**2)
                         Vang=numpy.arctan2(Vmer,Vzon)
                         Vver=xFrec[Vpos]
+                    print 'Height',Height
                     print 'vzon y vmer', Vzon, Vmer
-                    return Vzon, Vmer, Vver, GaussCenter
+                    return Vzon, Vmer, Vver, GaussCenter, PhaseSlope
             class SpectralMoments(Operation):

schainpy/model/proc/jroproc_parameters.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/scripts/schain.xml +1 -1

	@@ -1,1 +1,1
	1	<Project description="read bltr data sswma file" id="191" name="test1"><ReadUnit datatype="testBLTRReader" id="1911" inputId="0" name="testBLTRReader"><Operation id="19111" name="run" priority="1" type="self"><Parameter format="str" id="191111" name="datatype" value="testBLTRReader" /><Parameter format="str" id="191112" name="path" value="/media/erick/6F60F7113095A154/BLTR" /><Parameter format="date" id="191113" name="startDate" value="2017/01/17" /><Parameter format="date" id="191114" name="endDate" value="2018/01/01" /><Parameter format="time" id="191115" name="startTime" value="00:00:00" /><Parameter format="time" id="191116" name="endTime" value="23:59:59" /><Parameter format="str" id="191118" name="ext" value="sswma" /></Operation></ReadUnit><ProcUnit datatype="BLTRProcess" id="1912" inputId="1911" name="BLTRProcess"><Operation id="19121" name="run" priority="1" type="self" /><Operation id="19122" name="SnrFilter" priority="2" type="self"><Parameter format="float" id="191221" name="snr_val" value="-20" /><Parameter format="int" id="191222" name="modetofilter" value="2" /></Operation><Operation id="19123" name="OutliersFilter" priority="3" type="self"><Parameter format="str" id="191231" name="svalue" value="meridional" /><Parameter format="str" id="191232" name="svalue2" value="inTime" /><Parameter format="float" id="191233" name="method" value="0" /><Parameter format="float" id="191234" name="factor" value="2" /><Parameter format="float" id="191235" name="filter" value="0" /><Parameter format="float" id="191236" name="npoints" value="9" /><Parameter format="int" id="191237" name="modetofilter" value="2" /></Operation><Operation id="19124" name="OutliersFilter" priority="4" type="self"><Parameter format="str" id="191241" name="svalue" value="zonal" /><Parameter format="str" id="191242" name="svalue2" value="inTime" /><Parameter format="float" id="191243" name="method" value="0" /><Parameter format="float" id="191244" name="factor" value="2" /><Parameter format="float" id="191245" name="filter" value="0" /><Parameter format="float" id="191246" name="npoints" value="9" /><Parameter format="int" id="191247" name="modetofilter" value="2" /></Operation><Operation id="19125" name="OutliersFilter" priority="5" type="self"><Parameter format="str" id="191251" name="svalue" value="vertical" /><Parameter format="str" id="191252" name="svalue2" value="inHeight" /><Parameter format="float" id="191253" name="method" value="0" /><Parameter format="float" id="191254" name="factor" value="2" /><Parameter format="float" id="191255" name="filter" value="0" /><Parameter format="float" id="191256" name="npoints" value="9" /><Parameter format="int" id="191257" name="modetofilter" value="2" /></Operation><Operation id="19126" name="prePlot" priority="6" type="self"><Parameter format="int" id="191261" name="modeselect" value="1" /></Operation><Operation id="19127" name="WindProfilerPlot" priority="7" type="other"><Parameter format="int" id="191271" name="id" value="1" /><Parameter format="str" id="191272" name="wintitle" value="" /><Parameter format="intlist" id="191273" name="channelList" value="0" /><Parameter format="int" id="191274" name="SNRmin" value="-10" /><Parameter format="int" id="191275" name="SNRmax" value="50" /><Parameter format="float" id="191276" name="SNRthresh" value="0" /><Parameter format="float" id="191277" name="xmin" value="0" /><Parameter format="float" id="191278" name="xmax" value="24" /><Parameter format="float" id="191279" name="ymax" value="3" /><Parameter format="float" id="191280" name="zmin" value="-20" /><Parameter format="float" id="191281" name="zmax" value="20" /><Parameter format="float" id="191282" name="zmin_ver" value="-200" /><Parameter format="float" id="191283" name="zmax_ver" value="200" /></Operation><Operation id="19128" name="prePlot" priority="8" type="self"><Parameter format="int" id="191281" name="modeselect" value="2" /></Operation><Operation id="19129" name="WindProfilerPlot" priority="9" type="other"><Parameter format="int" id="191291" name="id" value="2" /><Parameter format="str" id="191292" name="wintitle" value="" /><Parameter format="bool" id="191293" name="save" value="1" /><Parameter format="str" id="191294" name="figpath" value="/media/erick/6F60F7113095A154/BLTR/" /><Parameter format="int" id="191295" name="SNRmin" value="-20" /><Parameter format="int" id="191296" name="SNRmax" value="40" /><Parameter format="float" id="191297" name="SNRthresh" value="0" /><Parameter format="float" id="191298" name="xmin" value="0" /><Parameter format="float" id="191299" name="xmax" value="24" /><Parameter format="float" id="191300" name="ymin" value="0" /><Parameter format="float" id="191301" name="ymax" value="10" /><Parameter format="float" id="191302" name="zmin" value="-4" /><Parameter format="float" id="191303" name="zmax" value="4" /><Parameter format="float" id="191304" name="zmin_ver" value="-200" /><Parameter format="float" id="191305" name="zmax_ver" value="200" /></Operation></ProcUnit></Project> No newline at end of file		1	<Project description="Segundo Test" id="191" name="test01"><ReadUnit datatype="VoltageReader" id="1911" inputId="0" name="VoltageReader"><Operation id="19111" name="run" priority="1" type="self"><Parameter format="str" id="191111" name="datatype" value="VoltageReader" /><Parameter format="str" id="191112" name="path" value="/home/erick/Documents/Data/Claire_Data/raw" /><Parameter format="date" id="191113" name="startDate" value="2017/07/26" /><Parameter format="date" id="191114" name="endDate" value="2017/07/26" /><Parameter format="time" id="191115" name="startTime" value="10:02:00" /><Parameter format="time" id="191116" name="endTime" value="10:11:00" /><Parameter format="int" id="191118" name="online" value="0" /><Parameter format="int" id="191119" name="walk" value="1" /></Operation><Operation id="19112" name="printNumberOfBlock" priority="2" type="self" /></ReadUnit><ProcUnit datatype="SpectraProc" id="1913" inputId="1912" name="SpectraProc"><Operation id="19131" name="run" priority="1" type="self"><Parameter format="int" id="191311" name="nFFTPoints" value="128" /><Parameter format="pairslist" id="191312" name="pairsList" value="(0,1),(0,2),(1,2)" /></Operation><Operation id="19132" name="removeDC" priority="2" type="self" /><Operation id="19133" name="IncohInt" priority="3" type="external"><Parameter format="float" id="191331" name="n" value="30" /></Operation><Operation id="19134" name="CrossSpectraPlot" priority="4" type="other"><Parameter format="str" id="191341" name="phase_cmap" value="bwr" /><Parameter format="int" id="191342" name="id" value="2005" /><Parameter format="str" id="191343" name="wintitle" value="CrossSpectraPlot_ShortPulse" /><Parameter format="str" id="191344" name="xaxis" value="Velocity" /><Parameter format="float" id="191345" name="ymin" value="1" /><Parameter format="int" id="191346" name="ymax" value="7" /><Parameter format="int" id="191347" name="zmin" value="15" /><Parameter format="int" id="191348" name="zmax" value="60" /><Parameter format="int" id="191349" name="save" value="2" /><Parameter format="str" id="191350" name="figpath" value="/media/erick/6F60F7113095A154/CLAIRE/CLAIRE_WINDS_2MHZ/Images" /></Operation></ProcUnit><ProcUnit datatype="VoltageProc" id="1912" inputId="1911" name="VoltageProc"><Operation id="19121" name="run" priority="1" type="self" /><Operation id="19122" name="setRadarFrequency" priority="2" type="self"><Parameter format="float" id="191221" name="frequency" value="445.09e6" /></Operation><Operation id="19123" name="selectHeights" priority="3" type="self"><Parameter format="float" id="191231" name="minHei" value="0" /><Parameter format="float" id="191232" name="maxHei" value="64" /></Operation></ProcUnit><ProcUnit datatype="Parameters" id="1914" inputId="1913" name="ParametersProc"><Operation id="19141" name="run" priority="1" type="self" /><Operation id="19142" name="GaussianFit" priority="2" type="other" /></ProcUnit></Project> No newline at end of file

schainpy/scripts/testBLTR_block.py +1 -1

             opObj10.addParameter(name='xmin', value='0', format='float')
             opObj10.addParameter(name='xmax', value='24', format='float')
             opObj10.addParameter(name='ymin', value='0', format='float')
-            opObj10.addParameter(name='ymax', value='10', format='float')
+            opObj10.addParameter(name='ymax', value='7', format='float')
             opObj10.addParameter(name='zmin', value='-4', format='float')
             opObj10.addParameter(name='zmax', value='4', format='float')
             opObj10.addParameter(name='zmin_ver', value='-200', format='float')

schainpy/scripts/testProcBLTR.py +3 -2

             procUnitConfObj2 = controllerObj.addProcUnit(datatype='Parameters', inputId=procUnitConfObj1.getId())
             opObj11 = procUnitConfObj2.addOperation(name='SpectralMoments', optype='other')
             opObj22 = procUnitConfObj2.addOperation(name='FullSpectralAnalysis', optype='other')
+            opObj22.addParameter(name='SNRlimit', value='-4', format='float')
             #
             opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other')
             opObj22.addParameter(name='id', value='4', format='int')
             opObj22.addParameter(name='SNRmin', value='-5', format='int')
             opObj22.addParameter(name='SNRmax', value='30', format='int')
             # opObj22.addParameter(name='SNRthresh', value='-3.5', format='float')
-            opObj22.addParameter(name='xmin', value=0, format='float')
+            opObj22.addParameter(name='xmin', value='0', format='float')
-            opObj22.addParameter(name='xmax', value=24, format='float')
+            opObj22.addParameter(name='xmax', value='24', format='float')
             opObj22.addParameter(name='ymin', value='225', format='float')
             #opObj22.addParameter(name='ymax', value='2000', format='float')
             opObj22.addParameter(name='save', value='1', format='int')

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