schain-jc Commit - r1006:8363978b5396

22/08/2017

ebocanegra -

r1006:8363978b5396

parent child

Context file:

r1006:8363978b5396

Expand all files

.gitignore +3 -6

-            <<<<<<< HEAD
             # Byte-compiled / optimized / DLL files
             __pycache__/
             *.py[cod]
             # eclipse
             .project
             .pydevproject
-            =======
-            build/
-            dist/
-            schainpy.egg-info/
             .svn/
-            >>>>>>> 08c4507d6c3c48f6c52326d5dedfa1972fb26356
+            *.png
+            *.pyc
+            schainpy/scripts

schainpy/admin.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/controller.pyc binary modified 0 0

NO CONTENT: modified file, binary diff hidden

schainpy/model/data/jrodata.pyc binary modified 0 0

	1	NO CONTENT: modified file, binary diff hidden			NO CONTENT: modified file, binary diff hidden

NO CONTENT: modified file, binary diff hidden

schainpy/model/data/jroheaderIO.pyc binary modified 0 0

	1	NO CONTENT: modified file, binary diff hidden			NO CONTENT: modified file, binary diff hidden

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/io/__init__.py +1 -1

             from jroIO_bltr import *
             from jroIO_mira35c import *
-            #from io_bltr_block import *
+            from io_bltr_block import *

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/io/io_bltr_block.py +8 -6

             import warnings
             from time import gmtime
             from math import floor
             warnings.simplefilter("error")
             from numpy.lib.nanfunctions import nansum
             warnings.simplefilter('ignore', FutureWarning)
             class testBLTRReader(ProcessingUnit):
-                def __init__(self):
+                def __init__(self, **kwargs):
                     path = None
                     startDate = None
                     lon = None
                     siteFile = None
-                    ProcessingUnit.__init__(self)
+                    ProcessingUnit.__init__(self , **kwargs)
                     self.dataOut = self.createObjByDefault()
                     self.imode = 0
                     self.counter_records = 0
                     print 'Boundary Layer and Tropospheric Radar (BLTR) script, Wind velocities and SNR from *.sswma files'
                     print '{}  \n'.format(signalchainweb)
-                def run(self, path, startDate, endDate, ext, startTime, endTime):
+                def run(self, path, startDate, endDate, ext, startTime, endTime, queue=None):
                     if not(self.isConfig):
                         self.setup(path, startDate, endDate, ext)
                                       ('site', 'u1', (32,))
                                       ])
                     self.pointer = open(self.filename, 'rb')  # rb : Read Binary
+                    print self.filename
                     self.header_file = numpy.fromfile(self.pointer, self.h_file, 1)
+                    print self.header_file
                     self.nrecords = self.header_file['nrec'][0]
                     self.sizeOfFile = os.path.getsize(self.filename)
                     '''
                     self.dataOut.type = 'Parameters'
                     self.dataOut.useLocalTime = False
-                    self.dataOut.outputInterval = 157
+                    # self.dataOut.outputInterval = 157
+                    self.dataOut.paramInterval = 157
                     self.dataOut.timezone = self.timezone
                     self.dataOut.site = self.siteFile
                     self.dataOut.nrecords = self.nrecords

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/__init__.py +1 -1

             from jroproc_parameters import *
             from jroproc_spectra_lags import *
             from jroproc_spectra_acf import *
-            #from jroproc_bltr import *
+            from jroproc_bltr import *

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/jroproc_bltr.py +2 -2

                     self.dataOut.year - Experiment year
                 '''
-                def __init__(self):
+                def __init__(self, **kwargs):
                     '''
                     Inputs: None
                     '''
-                    ProcessingUnit.__init__(self)
+                    ProcessingUnit.__init__(self, **kwargs)
                     self.dataOut = Parameters()
                     # Filters

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/jroproc_parameters.py +66 -35

             import warnings
             from numpy import NaN
+            from scipy.optimize.optimize import OptimizeWarning
             warnings.filterwarnings('ignore')
             SPEED_OF_LIGHT = 299792458
                     Parameters affected:    Winds, height range, SNR
                 """
-                def run(self, dataOut, E01=None, E02=None, E12=None, N01=None, N02=None, N12=None):
+                def run(self, dataOut, E01=None, E02=None, E12=None, N01=None, N02=None, N12=None, SNRlimit=7):
                     spc = dataOut.data_pre[0].copy()
                     cspc = dataOut.data_pre[1].copy()
                     coherence=numpy.ones([nChannel,nProfiles])
                     PhaseSlope=numpy.ones(nChannel)
                     PhaseInter=numpy.ones(nChannel)
+                    dataSNR = dataOut.data_SNR
                     data = dataOut.data_pre
                     noise = dataOut.noise
                     print 'noise',noise
-                    SNRdB = 10*numpy.log10(dataOut.data_SNR)
+                    #SNRdB = 10*numpy.log10(dataOut.data_SNR)
-                    FirstMoment = []
+                    FirstMoment = numpy.average(dataOut.data_param[:,1,:],0)
-                    SNRdBMean = []
+                    #SNRdBMean = []
-                    for j in range(nHeights):
+                    #for j in range(nHeights):
-                        FirstMoment = numpy.append(FirstMoment,numpy.mean([dataOut.data_param[0,1,j],dataOut.data_param[1,1,j],dataOut.data_param[2,1,j]]))
+                    #    FirstMoment = numpy.append(FirstMoment,numpy.mean([dataOut.data_param[0,1,j],dataOut.data_param[1,1,j],dataOut.data_param[2,1,j]]))
-                        SNRdBMean = numpy.append(SNRdBMean,numpy.mean([SNRdB[0,j],SNRdB[1,j],SNRdB[2,j]]))
+                    #    SNRdBMean = numpy.append(SNRdBMean,numpy.mean([SNRdB[0,j],SNRdB[1,j],SNRdB[2,j]]))
                     data_output=numpy.ones([3,spc.shape[2]])*numpy.NaN
                     velocityY=[]
                     velocityV=[]
+                    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)
+                        [Vzon,Vmer,Vver, GaussCenter]= self.WindEstimation(spc, cspc, pairsList, ChanDist, Height, noise, VelRange, dbSNR[Height], SNRlimit)
                         if abs(Vzon)<100. and abs(Vzon)> 0.:
                             velocityX=numpy.append(velocityX, Vzon)#Vmag
                             print 'Vmer',Vmer
                             velocityY=numpy.append(velocityY, numpy.NaN)
-                        if abs(GaussCenter)<10:
+                        if dbSNR[Height] > SNRlimit:
-                            velocityV=numpy.append(velocityV, Vver)
+                            velocityV=numpy.append(velocityV, FirstMoment[Height])
                         else:
                             velocityV=numpy.append(velocityV, numpy.NaN)
                             #FirstMoment[Height]= numpy.NaN
                     data_output[0]=numpy.array(velocityX)
                     data_output[1]=numpy.array(velocityY)
-                    data_output[2]=-FirstMoment
+                    data_output[2]=-velocityV#FirstMoment
                     print ' '
                     #print 'FirstMoment'
                         if x[index]==value:
                             return index
-                def WindEstimation(self, spc, cspc, pairsList, ChanDist, Height, noise, VelRange):
+                def WindEstimation(self, spc, cspc, pairsList, ChanDist, Height, noise, VelRange, dbSNR, SNRlimit):
                     ySamples=numpy.ones([spc.shape[0],spc.shape[1]])
                     phase=numpy.ones([spc.shape[0],spc.shape[1]])
                     for i in range(spc.shape[0]):
                         '''****** Line of Data SPC ******'''
                         zline=z[i,:,Height]
                         '''****** SPC is normalized ******'''
-                        FactNorm= zline.copy() / numpy.sum(zline.copy())
+                        FactNorm= (zline.copy()-noise[i]) / numpy.sum(zline.copy())
                         FactNorm= FactNorm/numpy.sum(FactNorm)
                         SmoothSPC=self.moving_average(FactNorm,N=3)
                         xSamples = ar(range(len(SmoothSPC)))
-                        ySamples[i] = SmoothSPC-noise[i]
+                        ySamples[i] = SmoothSPC
+                    #dbSNR=10*numpy.log10(dataSNR)
                     print ' '
                     print ' '
                     print ' '
-                    print 'SmoothSPC',SmoothSPC
+                    #print 'dataSNR', dbSNR.shape, dbSNR[0,40:120]
+                    print 'SmoothSPC', SmoothSPC.shape, SmoothSPC[0:20]
                     print 'noise',noise
-                    print'zline',zline
+                    print 'zline',zline.shape, zline[0:20]
-                    print'FactNorm',FactNorm
+                    print 'FactNorm',FactNorm.shape, FactNorm[0:20]
+                    print 'FactNorm suma', numpy.sum(FactNorm)
                     for i in range(spc.shape[0]):
                         '''****** CSPC is normalized ******'''
                         chan_index0 = pairsList[i][0]
                         chan_index1 = pairsList[i][1]
-                        CSPCFactor= numpy.sum(ySamples[chan_index0]) * numpy.sum(ySamples[chan_index1])
+                        CSPCFactor= abs(numpy.sum(ySamples[chan_index0]) * numpy.sum(ySamples[chan_index1])) #
-                        CSPCNorm= cspcLine.copy() / numpy.sqrt(CSPCFactor)
+                        CSPCNorm = (cspcLine.copy() -noise[i]) / numpy.sqrt(CSPCFactor)
-                        CSPCSamples[i] = CSPCNorm-noise[i]
+                        CSPCSamples[i] = CSPCNorm
                         coherence[i] = numpy.abs(CSPCSamples[i]) / numpy.sqrt(CSPCFactor)
                         coherence[i]= self.moving_average(coherence[i],N=2)
                         phase[i] = self.moving_average( numpy.arctan2(CSPCSamples[i].imag, CSPCSamples[i].real),N=1)#*180/numpy.pi
-                    print 'CSPCSamples', CSPCSamples
+                    print 'cspcLine', cspcLine.shape, cspcLine[0:20]
+                    print 'CSPCFactor', CSPCFactor#, CSPCFactor[0:20]
+                    print numpy.sum(ySamples[chan_index0]), numpy.sum(ySamples[chan_index1]), -noise[i]
+                    print 'CSPCNorm', CSPCNorm.shape, CSPCNorm[0:20]
+                    print 'CSPCNorm suma', numpy.sum(CSPCNorm)
+                    print 'CSPCSamples', CSPCSamples.shape, CSPCSamples[0,0:20]
                     '''****** Getting fij width ******'''
                     yMean=[]
                     print '****************************'
                     print 'len(xSamples): ',len(xSamples)
-                    print 'yMean: ', yMean
+                    print 'yMean: ', yMean.shape, yMean[0:20]
-                    print 'ySamples', ySamples
+                    print 'ySamples', ySamples.shape, ySamples[0,0:20]
-                    print 'xSamples: ',xSamples
+                    print 'xSamples: ',xSamples.shape, xSamples[0:20]
                     print 'meanGauss',meanGauss
                     print 'sigma',sigma
-                    if (abs(meanGauss/sigma**2) > 0.000000001):#0.00001) :
+                    #if (abs(meanGauss/sigma**2) > 0.0001) : #0.000000001):
+                    if dbSNR > SNRlimit :
                         try:
                             popt,pcov = curve_fit(self.gaus,xSamples,yMean,p0=[1,meanGauss,sigma])
                             else:
                                 FitGauss=numpy.ones(len(xSamples))*numpy.mean(yMean)
                                 print 'Verificador:     Dentro', Height
-                        except RuntimeError:
+                        except :#RuntimeError:
                             FitGauss=numpy.ones(len(xSamples))*numpy.mean(yMean)
                         FitGauss=numpy.ones(len(xSamples))*numpy.mean(yMean)
                     Maximun=numpy.amax(yMean)
-                    eMinus1=Maximun*numpy.exp(-1)*0.8
+                    eMinus1=Maximun*numpy.exp(-1)#*0.8
                     HWpos=self.Find(FitGauss,min(FitGauss, key=lambda value:abs(value-eMinus1)))
                     HalfWidth= xFrec[HWpos]
                     if Rangpos<GCpos:
                         Range=numpy.array([Rangpos,2*GCpos-Rangpos])
-                    else:
+                    elif Rangpos< ( len(xFrec)- len(xFrec)*0.1):
                         Range=numpy.array([2*GCpos-Rangpos,Rangpos])
+                    else:
+                        Range = numpy.array([0,0])
+                    print ' '
+                    print 'GCpos',GCpos, ( len(xFrec)- len(xFrec)*0.1)
+                    print 'Rangpos',Rangpos
+                    print 'RANGE: ', Range
                     FrecRange=xFrec[Range[0]:Range[1]]
                     '''****** Getting SCPC Slope ******'''
                         if len(FrecRange)>5 and len(FrecRange)<spc.shape[1]*0.5:
                             PhaseRange=self.moving_average(phase[i,Range[0]:Range[1]],N=3)
-                            slope, intercept, r_value, p_value, std_err = stats.linregress(FrecRange,PhaseRange)
+                            print 'FrecRange', len(FrecRange) , FrecRange
-                            PhaseSlope[i]=slope
+                            print 'PhaseRange', len(PhaseRange), PhaseRange
-                            PhaseInter[i]=intercept
+                            print ' '
+                            if len(FrecRange) == len(PhaseRange):
+                                slope, intercept, r_value, p_value, std_err = stats.linregress(FrecRange,PhaseRange)
+                                PhaseSlope[i]=slope
+                                PhaseInter[i]=intercept
+                            else:
+                                PhaseSlope[i]=0
+                                PhaseInter[i]=0
                         else:
                             PhaseSlope[i]=0
                             PhaseInter[i]=0

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

NO CONTENT: modified file, binary diff hidden

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

NO CONTENT: modified file, binary diff hidden

schainpy/model/proc/jroproc_voltage.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="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/10/28" /><Parameter format="time" id="191115" name="startTime" value="09:40:00" /><Parameter format="time" id="191116" name="endTime" value="23:59: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="256" /><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="SpectraPlot" priority="4" type="external"><Parameter format="int" id="191341" name="id" value="11" /><Parameter format="str" id="191342" name="wintitle" value="SpectraPlot" /><Parameter format="str" id="191343" name="xaxis" value="velocity" /><Parameter format="int" id="191344" name="zmin" value="0" /><Parameter format="float" id="191345" name="ymin" value="1.5" /><Parameter format="int" id="191346" name="ymax" value="7" /><Parameter format="int" id="191347" name="zmax" value="40" /><Parameter format="int" id="191348" name="save" value="1" /></Operation><Operation id="19135" name="RTIPlot" priority="5" type="other"><Parameter format="int" id="191351" name="id" value="30" /><Parameter format="str" id="191352" name="wintitle" value="RTI" /><Parameter format="int" id="191353" name="zmin" value="0" /><Parameter format="int" id="191354" name="zmax" value="40" /><Parameter format="float" id="191355" name="ymin" value="1.5" /><Parameter format="int" id="191356" name="ymax" value="8" /><Parameter format="int" id="191357" name="showprofile" value="1" /><Parameter format="float" id="191358" name="xmin" value="0" /><Parameter format="float" id="191359" name="xmax" value="23.9" /><Parameter format="int" id="191360" name="save" value="1" /></Operation><Operation id="19136" name="CrossSpectraPlot" priority="6" type="other"><Parameter format="str" id="191361" name="phase_cmap" value="bwr" /><Parameter format="int" id="191362" name="id" value="2005" /><Parameter format="str" id="191363" name="wintitle" value="CrossSpectraPlot_ShortPulse" /><Parameter format="str" id="191364" name="xaxis" value="Velocity" /><Parameter format="float" id="191365" name="ymin" value="1.5" /><Parameter format="int" id="191366" name="ymax" value="7" /></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="SpectralMoments" priority="2" type="other" /><Operation id="19143" name="FullSpectralAnalysis" priority="3" type="other"><Parameter format="float" id="191431" name="E01" value="1.5" /><Parameter format="float" id="191432" name="E02" value="1.5" /><Parameter format="float" id="191433" name="E12" value="0" /><Parameter format="float" id="191434" name="N01" value="0.875" /><Parameter format="float" id="191435" name="N02" value="-0.875" /><Parameter format="float" id="191436" name="N12" value="-1.75" /></Operation><Operation id="19144" name="WindProfilerPlot" priority="4" type="other"><Parameter format="int" id="191441" name="id" value="4" /><Parameter format="str" id="191442" name="wintitle" value="Wind Profiler" /><Parameter format="bool" id="191443" name="save" value="1" /><Parameter format="float" id="191444" name="ymin" value="1.5" /><Parameter format="int" id="191445" name="ymax" value="7" /><Parameter format="float" id="191446" name="zmin" value="-0.00" /><Parameter format="float" id="191447" name="zmax" value="1.1" /></Operation></ProcUnit></Project> No newline at end of file		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

schainpy/scripts/testBLTR_block.py +8 -7

             filename = 'test1.xml'
             # path = '/home/jespinoza/workspace/data/bltr/'
-            path = '/home/erick/Documents/Data/BLTR_Data/sswma/'
+            path = '/media/erick/6F60F7113095A154/BLTR/'
             desc = "read bltr data sswma file"
-            figpath = '/home/erick/workspace'
+            figpath = '/media/erick/6F60F7113095A154/BLTR/'
             pathhdf5 = '/tmp/'
             controllerObj = Project()
             controllerObj.setup(id = '191', name='test1', description=desc)
             readUnitConfObj = controllerObj.addReadUnit(datatype='testBLTRReader',
                                                         path=path,
-                                                        startDate='2015/01/17',
+                                                        startDate='2017/01/17',
-                                                        endDate='2017/01/01',
+                                                        endDate='2018/01/01',
                                                         startTime='00:00:00',
                                                         endTime='23:59:59',
                                                         ext='sswma')
             opObj10.addParameter(name='id', value='2', format='int')
             opObj10.addParameter(name='wintitle', value='', format='str')
             #opObj10.addParameter(name='channelList', value='0', format='intlist')
-            #opObj10.addParameter(name='save', value='1', format='bool')
+            opObj10.addParameter(name='save', value='1', format='bool')
-            #opObj10.addParameter(name='figpath', value=figpath, format='str')
+            opObj10.addParameter(name='figpath', value=figpath, format='str')
             opObj10.addParameter(name='SNRmin', value='-20', format='int')
             opObj10.addParameter(name='SNRmax', value='40', format='int')
             opObj10.addParameter(name='SNRthresh', value='0', format='float')
             opObj10.addParameter(name='xmin', value='0', format='float')
             opObj10.addParameter(name='xmax', value='24', format='float')
-            #opObj10.addParameter(name='ymax', value='8', format='float')
+            opObj10.addParameter(name='ymin', value='0', format='float')
+            opObj10.addParameter(name='ymax', value='10', 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 +11 -7

             xmin = '15.5'
             xmax = '24'
             desc = "ProcBLTR Test"
             filename = "ProcBLTR.xml"
+            figpath = '/media/erick/6F60F7113095A154/BLTR'
             controllerObj = Project()
             controllerObj.setup(id='191', name='test01', description=desc)
             readUnitConfObj = controllerObj.addReadUnit(datatype='BLTRReader',
-                                                            path='/home/erick/Documents/Data/BLTR_Data/fdt/',
+                                                            path='/media/erick/6F60F7113095A154/BLTR/',
                                                             endDate='2017/10/19',
                                                             startTime='13:00:00',
             opObj11 = procUnitConfObj1.addOperation(name='IncohInt', optype='other')
-            opObj11.addParameter(name='n', value='2', format='float')
+            opObj11.addParameter(name='n', value='3', format='float')
             opObj10 = procUnitConfObj1.addOperation(name='removeDC')
             opObj22.addParameter(name='wintitle', value='Wind Profiler', format='str')
             opObj22.addParameter(name='save', value='1', format='bool')
             # opObj22.addParameter(name='figpath', value = '/home/erick/Pictures', format='str')
             opObj22.addParameter(name='zmin', value='-20', format='int')
             opObj22.addParameter(name='zmax', value='20', format='int')
             opObj22.addParameter(name='zmin_ver', value='-250', 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')
+            opObj22.addParameter(name='figpath', value=figpath, format='str')
             # opObj11.addParameter(name='pairlist', value='(1,0),(0,2),(1,2)', format='pairsList')
             #opObj10 = procUnitConfObj1.addOperation(name='selectHeights')
             #opObj10.addParameter(name='minHei', value='225', format='float')
             #opObj10.addParameter(name='maxHei', value='1000', format='float')
             # opObj11 = procUnitConfObj1.addOperation(name='CoherenceMap', optype='other')
             # opObj11.addParameter(name='id', value='102', format='int')
             # opObj11.addParameter(name='wintitle', value='Coherence', format='str')
             # opObj11.addParameter(name='ymin', value='225', format='float')
             # opObj11.addParameter(name='ymax', value='4000', format='float')
             # opObj11.addParameter(name='phase_cmap', value='jet', format='str')
             # opObj11.addParameter(name='xmin', value='8.5', format='float')
             # opObj11.addParameter(name='xmax', value='9.5', format='float')

schainpy/scripts/wind004/wind_20161108_000201.png binary modified 0 0

NO CONTENT: modified file, binary diff hidden

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