##// END OF EJS Templates
schain-jc
RSS
27 de Marzo 2018
ebocanegra -
r1154:a08206fb2eae
parent child
Show More
Side by Side Unified
Context file:
r1154:a08206fb2eae Loading diff...:
Show file before | Show file after | Hide comments
@@ -10,10 +10,9 import importlib
10 import itertools
10 import itertools
11 from multiprocessing import Pool, TimeoutError
11 from multiprocessing import Pool, TimeoutError
12 from multiprocessing.pool import ThreadPool
12 from multiprocessing.pool import ThreadPool
13 import copy_reg
13
14 import cPickle
14
15 import types
15
16 from functools import partial
17 import time
16 import time
18 #from sklearn.cluster import KMeans
17 #from sklearn.cluster import KMeans
19
18
@@ -759,7 +758,7 class PrecipitationProc(Operation):
759
758
760 Numerator = ( (4*numpy.pi)**3 * aL**2 * 16 * numpy.log(2) )
759 Numerator = ( (4*numpy.pi)**3 * aL**2 * 16 * numpy.log(2) )
761 Denominator = ( Pt * Gt * Gr * Lambda**2 * SPEED_OF_LIGHT * tauW * numpy.pi * ThetaT * ThetaR)
760 Denominator = ( Pt * Gt * Gr * Lambda**2 * SPEED_OF_LIGHT * tauW * numpy.pi * ThetaT * ThetaR)
762 RadarConstant = 1e-10 * Numerator / Denominator #
761 RadarConstant = 5e-27 * Numerator / Denominator #
763 print '***'
762 print '***'
764 print '*** RadarConstant' , RadarConstant, '****'
763 print '*** RadarConstant' , RadarConstant, '****'
765 print '***'
764 print '***'
@@ -787,7 +786,7 class PrecipitationProc(Operation):
787 D_range = numpy.zeros([self.Num_Bin,self.Num_Hei])
786 D_range = numpy.zeros([self.Num_Bin,self.Num_Hei])
788 SIGMA = numpy.zeros([self.Num_Bin,self.Num_Hei])
787 SIGMA = numpy.zeros([self.Num_Bin,self.Num_Hei])
789 N_dist = numpy.zeros([self.Num_Bin,self.Num_Hei])
788 N_dist = numpy.zeros([self.Num_Bin,self.Num_Hei])
790 D_mean = numpy.zeros(self.Num_Hei)
789 V_mean = numpy.zeros(self.Num_Hei)
791 del_V = numpy.zeros(self.Num_Hei)
790 del_V = numpy.zeros(self.Num_Hei)
792 Z = numpy.zeros(self.Num_Hei)
791 Z = numpy.zeros(self.Num_Hei)
793 Ze = numpy.zeros(self.Num_Hei)
792 Ze = numpy.zeros(self.Num_Hei)
@@ -796,7 +795,7 class PrecipitationProc(Operation):
796 Range = dataOut.heightList*1000.
795 Range = dataOut.heightList*1000.
797
796
798 for R in range(self.Num_Hei):
797 for R in range(self.Num_Hei):
799
798
800 h = Range[R] + Altitude #Range from ground to radar pulse altitude
799 h = Range[R] + Altitude #Range from ground to radar pulse altitude
801 del_V[R] = 1 + 3.68 * 10**-5 * h + 1.71 * 10**-9 * h**2 #Density change correction for velocity
800 del_V[R] = 1 + 3.68 * 10**-5 * h + 1.71 * 10**-9 * h**2 #Density change correction for velocity
802
801
@@ -819,20 +818,23 class PrecipitationProc(Operation):
819
818
820 N_dist[:,R] = ETAn[:,R] / SIGMA[:,R]
819 N_dist[:,R] = ETAn[:,R] / SIGMA[:,R]
821
820
822 DMoments = self.Moments(Pr[:,R], D_range[:,R])
821 DMoments = self.Moments(Pr[:,R], Velrange[0:self.Num_Bin])
823
822
824 try:
823 try:
825 popt01,pcov = curve_fit(self.gaus, D_range[:,R] , Pr[:,R] , p0=DMoments)
824 popt01,pcov = curve_fit(self.gaus, Velrange[0:self.Num_Bin] , Pr[:,R] , p0=DMoments)
826 except:
825 except:
827 popt01=numpy.zeros(3)
826 popt01=numpy.zeros(3)
828 popt01[1]= DMoments[1]
827 popt01[1]= DMoments[1]
829 D_mean[R]=popt01[1]
828 if popt01[1]<0 or popt01[1]>20:
829 popt01[1]=numpy.NaN
830
831 V_mean[R]=popt01[1]
830
832
831 Z[R] = numpy.nansum( N_dist[:,R] * (D_range[:,R])**6 )*1e-18
833 Z[R] = numpy.nansum( N_dist[:,R] * (D_range[:,R])**6 )#*10**-18
832
834
833 RR[R] = 3.6e-6*1e-9*6*10**-4.*numpy.pi * numpy.nansum( D_range[:,R]**3 * N_dist[:,R] * Velrange[0:self.Num_Bin] ) #Rainfall rate
835 RR[R] = 0.0006*numpy.pi * numpy.nansum( D_range[:,R]**3 * N_dist[:,R] * Velrange[0:self.Num_Bin] ) #Rainfall rate
834
836
835 Ze[R] = (numpy.nansum( ETAn[:,R]) * Lambda**4) / ( numpy.pi**5 * Km)
837 Ze[R] = (numpy.nansum( ETAn[:,R]) * Lambda**4) / ( 10**-18*numpy.pi**5 * Km)
836
838
837
839
838
840
@@ -848,25 +850,26 class PrecipitationProc(Operation):
848 dataOut.channelList = [0,1,2]
850 dataOut.channelList = [0,1,2]
849
851
850 dataOut.data_param[0]=dBZ
852 dataOut.data_param[0]=dBZ
851 dataOut.data_param[1]=RR2
853 dataOut.data_param[1]=V_mean
852 dataOut.data_param[2]=RR
854 dataOut.data_param[2]=RR
853
855
854 #print 'VELRANGE', Velrange
856 #print 'VELRANGE', Velrange
855 print 'Range', len(Range)
857 #print 'Range', len(Range)
856 print 'delv',del_V
858 #print 'delv',del_V
857 #print 'DRANGE', D_range[:,50]
859 # print 'DRANGE', D_range[:,56]
858 print 'NOISE', dataOut.noise[0], 10*numpy.log10(dataOut.noise[0])
860 # #print 'NOISE', dataOut.noise[0], 10*numpy.log10(dataOut.noise[0])
859 print 'radarconstant', RadarConstant
861 # print 'radarconstant', RadarConstant
860 print 'Range', Range
862 # #print 'ETAn SHAPE', ETAn.shape
861 # print 'ETAn SHAPE', ETAn.shape
863 # # print 'ETAn ', numpy.nansum(ETAn, axis=0)
862 # print 'ETAn ', numpy.nansum(ETAn, axis=0)
864 # # print 'ETAd ', numpy.nansum(ETAd, axis=0)
863 # print 'ETAd ', numpy.nansum(ETAd, axis=0)
865 # print 'Pr ', numpy.nansum(Pr, axis=0)
864 print 'Pr ', numpy.nansum(Pr, axis=0)
866 # print 'dataOut.SPCparam[1]', numpy.nansum(dataOut.SPCparam[1][0], axis=0)
865 print 'dataOut.SPCparam[1]', numpy.nansum(dataOut.SPCparam[1][0], axis=0)
866 # print 'Ze ', dBZe
867 # print 'Ze ', dBZe
867 # print 'Z ', dBZ
868 print 'Z ', dBZ
868 print 'RR2 ', RR2
869 # print 'Ndist',N_dist[:,56]
870 # #print 'RR2 ', RR2
869 print 'RR ', RR
871 print 'RR ', RR
872 print 'Vr', V_mean
870 #print 'RR2 ', dBRR2
873 #print 'RR2 ', dBRR2
871 #print 'D_mean', D_mean
874 #print 'D_mean', D_mean
872 #print 'del_V', del_V
875 #print 'del_V', del_V
@@ -1029,7 +1032,7 class FullSpectralAnalysis(Operation):
1029 data_output[1] = numpy.array(velocityY) #self.moving_average(numpy.array(velocityY) , N=1)
1032 data_output[1] = numpy.array(velocityY) #self.moving_average(numpy.array(velocityY) , N=1)
1030 data_output[2] = -velocityV#FirstMoment
1033 data_output[2] = -velocityV#FirstMoment
1031
1034
1032 print 'FirstMoment', data_output[2]
1035 print 'data_output', data_output.shape
1033 #print FirstMoment
1036 #print FirstMoment
1034 # print 'velocityX',numpy.shape(data_output[0])
1037 # print 'velocityX',numpy.shape(data_output[0])
1035 # print 'velocityX',data_output[0]
1038 # print 'velocityX',data_output[0]
@@ -1058,7 +1061,7 class FullSpectralAnalysis(Operation):
1058 return ( Amp / ((2*numpy.pi)**0.5 * Sigma) ) * numpy.exp( -( xSamples - Mu )**2 / ( 2 * (Sigma**2) ))
1061 return ( Amp / ((2*numpy.pi)**0.5 * Sigma) ) * numpy.exp( -( xSamples - Mu )**2 / ( 2 * (Sigma**2) ))
1059
1062
1060
1063
1061
1064
1062 def Moments(self, ySamples, xSamples):
1065 def Moments(self, ySamples, xSamples):
1063 Pot = numpy.nansum( ySamples ) # Potencia, momento 0
1066 Pot = numpy.nansum( ySamples ) # Potencia, momento 0
1064 yNorm = ySamples / Pot
1067 yNorm = ySamples / Pot
Show file before | Show file after | Hide comments
@@ -1,1 +1,1
1 <Project description="Segundo Test" id="191" name="test01"><ReadUnit datatype="SpectraReader" id="1911" inputId="0" name="SpectraReader"><Operation id="19111" name="run" priority="1" type="self"><Parameter format="str" id="191111" name="datatype" value="SpectraReader" /><Parameter format="str" id="191112" name="path" value="/media/erick/6F60F7113095A154/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdataCLAIRE/Extra" /><Parameter format="date" id="191113" name="startDate" value="2018/02/01" /><Parameter format="date" id="191114" name="endDate" value="2018/02/01" /><Parameter format="time" id="191115" name="startTime" value="17:00:00" /><Parameter format="time" id="191116" name="endTime" value="20:00:00" /><Parameter format="int" id="191118" name="online" value="0" /><Parameter format="int" id="191119" name="walk" value="1" /></Operation><Operation id="19112" name="printInfo" priority="2" type="self" /><Operation id="19113" name="printNumberOfBlock" priority="3" type="self" /></ReadUnit><ProcUnit datatype="Parameters" id="1913" inputId="1912" name="ParametersProc"><Operation id="19131" name="run" priority="1" type="self" /><Operation id="19132" name="SpectralFilters" priority="2" type="other"><Parameter format="float" id="191321" name="PositiveLimit" value="1.5" /><Parameter format="float" id="191322" name="NegativeLimit" value="3.5" /></Operation><Operation id="19133" name="FullSpectralAnalysis" priority="3" type="other"><Parameter format="float" id="191331" name="SNRlimit" value="-16" /><Parameter format="float" id="191332" name="E01" value="1.500" /><Parameter format="float" id="191333" name="E02" value="1.500" /><Parameter format="float" id="191334" name="E12" value="0" /><Parameter format="float" id="191335" name="N01" value="0.875" /><Parameter format="float" id="191336" name="N02" value="-0.875" /><Parameter format="float" id="191337" name="N12" value="-1.750" /></Operation><Operation id="19134" name="WindProfilerPlot" priority="4" type="other"><Parameter format="int" id="191341" name="id" value="4" /><Parameter format="str" id="191342" name="wintitle" value="Wind Profiler" /><Parameter format="float" id="191343" name="xmin" value="17" /><Parameter format="float" id="191344" name="xmax" value="20" /><Parameter format="float" id="191345" name="ymin" value="0" /><Parameter format="int" id="191346" name="ymax" value="11" /><Parameter format="float" id="191347" name="zmin" value="-20" /><Parameter format="float" id="191348" name="zmax" value="20" /><Parameter format="float" id="191349" name="SNRmin" value="-20" /><Parameter format="float" id="191350" name="SNRmax" value="20" /><Parameter format="float" id="191351" name="zmin_ver" value="-200" /><Parameter format="float" id="191352" name="zmax_ver" value="200" /><Parameter format="float" id="191353" name="SNRthresh" value="-20" /><Parameter format="int" id="191354" name="save" value="1" /><Parameter format="str" id="191355" name="figpath" value="/media/erick/6F60F7113095A154/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdataCLAIRE/Extra" /></Operation><Operation id="19135" name="PrecipitationProc" priority="5" type="other" /><Operation id="19136" name="ParametersPlot" priority="6" type="other"><Parameter format="int" id="191361" name="id" value="10" /><Parameter format="str" id="191362" name="wintitle" value="First_gg" /><Parameter format="str" id="191363" name="colormap" value="ocean_r" /><Parameter format="int" id="191364" name="zmin" value="50" /><Parameter format="int" id="191365" name="zmax" value="100" /><Parameter format="int" id="191366" name="ymin" value="0" /><Parameter format="int" id="191367" name="ymax" value="11" /><Parameter format="int" id="191368" name="xmin" value="17" /><Parameter format="int" id="191369" name="xmax" value="20" /><Parameter format="int" id="191370" name="save" value="1" /><Parameter format="str" id="191371" name="figpath" value="/media/erick/6F60F7113095A154/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdataCLAIRE/Extra" /></Operation><Operation id="19137" name="SpcParamPlot" priority="7" type="other"><Parameter format="int" id="191371" name="id" value="21" /><Parameter format="str" id="191372" name="wintitle" value="Primer eco removido" /><Parameter format="str" id="191373" name="xaxis" value="velocity" /><Parameter format="int" id="191374" name="showprofile" value="1" /><Parameter format="int" id="191375" name="zmin" value="10" /><Parameter format="int" id="191376" name="zmax" value="40" /><Parameter format="int" id="191377" name="Selector" value="1" /></Operation><Operation id="19138" name="ParamWriter" priority="8" type="other"><Parameter format="str" id="191381" name="path" value="/media/erick/6F60F7113095A154/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdatatest/test1024" /><Parameter format="int" id="191382" name="blocksPerFile" value="100" /><Parameter format="list" id="191383" name="metadataList" value="heightList,timeZone,paramInterval" /><Parameter format="list" id="191384" name="dataList" value="data_output,data_SNR,utctime,utctimeInit" /></Operation></ProcUnit><ProcUnit datatype="SpectraProc" id="1912" inputId="1911" name="SpectraProc"><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></ProcUnit></Project> No newline at end of file
1 <Project description="Segundo Test" id="191" name="test01"><ReadUnit datatype="SpectraReader" id="1911" inputId="0" name="SpectraReader"><Operation id="19111" name="run" priority="1" type="self"><Parameter format="str" id="191111" name="datatype" value="SpectraReader" /><Parameter format="str" id="191112" name="path" value="/data/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdataCLAIRE/Extra" /><Parameter format="date" id="191113" name="startDate" value="2018/02/01" /><Parameter format="date" id="191114" name="endDate" value="2018/02/01" /><Parameter format="time" id="191115" name="startTime" value="17:00:00" /><Parameter format="time" id="191116" name="endTime" value="20:00:00" /><Parameter format="int" id="191118" name="online" value="0" /><Parameter format="int" id="191119" name="walk" value="1" /></Operation><Operation id="19112" name="printInfo" priority="2" type="self" /><Operation id="19113" name="printNumberOfBlock" priority="3" type="self" /></ReadUnit><ProcUnit datatype="Parameters" id="1913" inputId="1912" name="ParametersProc"><Operation id="19131" name="run" priority="1" type="self" /><Operation id="19132" name="SpectralFilters" priority="2" type="other"><Parameter format="float" id="191321" name="PositiveLimit" value="1.5" /><Parameter format="float" id="191322" name="NegativeLimit" value="3.5" /></Operation><Operation id="19133" name="PrecipitationProc" priority="3" type="other" /><Operation id="19134" name="ParametersPlot" priority="4" type="other"><Parameter format="int" id="191341" name="id" value="10" /><Parameter format="str" id="191342" name="wintitle" value="First_gg" /><Parameter format="str" id="191343" name="colormap" value="ocean_r" /><Parameter format="int" id="191344" name="zmin" value="00" /><Parameter format="int" id="191345" name="zmax" value="40" /><Parameter format="int" id="191346" name="ymin" value="0" /><Parameter format="int" id="191347" name="ymax" value="11" /><Parameter format="int" id="191348" name="xmin" value="17" /><Parameter format="int" id="191349" name="xmax" value="20" /><Parameter format="int" id="191350" name="save" value="1" /><Parameter format="str" id="191351" name="figpath" value="/data/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdataCLAIRE/Extra" /></Operation><Operation id="19135" name="ParamWriter" priority="5" type="other"><Parameter format="str" id="191351" name="path" value="/data/CLAIRE/CLAIRE_WINDS_2MHZ/DATA/pdatatest/test1024" /><Parameter format="int" id="191352" name="blocksPerFile" value="100" /><Parameter format="list" id="191353" name="metadataList" value="heightList,timeZone,paramInterval" /><Parameter format="list" id="191354" name="dataList" value="data_output,data_SNR,utctime,utctimeInit" /></Operation></ProcUnit><ProcUnit datatype="SpectraProc" id="1912" inputId="1911" name="SpectraProc"><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></ProcUnit></Project> No newline at end of file
General Comments 0
You need to be logged in to leave comments. Login now