schain Commit - r1410:577827e9cc8b · Jicamarca Repository

fix white plot

avaldezp -

r1410:577827e9cc8b

parent child

Context file:

r1410:577827e9cc8b

Collapse all files

schainpy/model/graphics/jroplot_parameters.py +22 -13

                     data['weather'] = 10*numpy.log10(dataOut.data_360[1]/(factor))
                     ####print("weather",data['weather'])
                     data['azi']     = dataOut.data_azi
+                    data['ele']     = dataOut.data_ele
                     return data, meta
                 def get2List(self,angulos):
                         condition , value = self.search_pos(i,list1_)
                         if condition:
                             pos = tmp + c + 1
-                            for k in range(list2_[value]-1):
+                            for k in range(round(list2_[value])-1):
                                 ang_new[pos+k]  = ang_new[pos+k-1]+1
                                 ang_new2[pos+k] = numpy.nan
                             tmp = pos +k
                     return data
                 def replaceNAN(self,data_weather,data_azi,val):
-                    ####print("----------------activeNEWFUNCTION")
                     data= data_azi
                     data_T= data_weather
-                    ####print("data_azi",data_azi)
+                    if data.shape[0]> data_T.shape[0]:
-                    ####print("VAL:",val)
+                        data_N = numpy.ones( [data.shape[0],data_T.shape[1]])
-                    ####print("SHAPE",data_T.shape)
+                        c      = 0
-                    for i in range(len(data)):
+                        for i in range(len(data)):
-                        if numpy.isnan(data[i]):
+                            if numpy.isnan(data[i]):
-                           ####print("NAN")
+                                data_N[i,:]=numpy.ones(data_T.shape[1])*numpy.nan
-                           #data_T[i,:]=numpy.ones(data_T.shape[1])*val
+                            else:
-                           data_T[i,:]=numpy.ones(data_T.shape[1])*numpy.nan
+                                data_N[i,:]=data_T[c,:]
-                    return data_T
+                                sc=c+1
+                    else:
+                        for i in range(len(data)):
+                            if numpy.isnan(data[i]):
+                               data_T[i,:]=numpy.ones(data_T.shape[1])*numpy.nan
+                        return data_T
                 def const_ploteo(self,data_weather,data_azi,step,res):
                     if self.ini==0:
                         #------------------------
                         ####data_azi_new  = self.fixDATA(data_azi)
                         #ata_azi_new   = self.fixDATANEW(data_azi)
                         start = data_azi_new[-1] + res
                         end   = data_azi_new[0]  - res
                         ##### new
                         start_azi = self.res_azi[0]
                         #-----------new------------
                         data_azi ,data_azi_old= self.globalCheckPED(data_azi)
+                        print("---------------------------------------------------")
+                        print("data_azi",data_azi)
+                        print("data_azi_old",data_azi_old)
                         data_weather = self.replaceNAN(data_weather=data_weather,data_azi=data_azi_old,val=self.val_mean)
                         #--------------------------
                         ####data_azi_old = data_azi
                     self.res_weather, self.res_azi = self.const_ploteo(data_weather=data['weather'][:,r_mask],data_azi=data['azi'],step=step,res=res)
                     #numpy.set_printoptions(suppress=True)
                     #print("resultado",self.res_azi)
+                    self.res_ele   =numpy.mean(data['ele'])
                     ###########################/DATA_RM/10_tmp/ch0###############################
                     #################    PLOTEO           ###################
                     ##########################################################
                     cbar = plt.gcf().colorbar(pm, pad=0.075)
                     caax.set_xlabel('x_range [km]')
                     caax.set_ylabel('y_range [km]')
-                    plt.text(1.0, 1.05, 'azimuth '+str(thisDatetime)+"  step   "+str(self.ini), transform=caax.transAxes, va='bottom',ha='right')
+                    plt.text(1.0, 1.05, 'Azimuth '+str(thisDatetime)+"  Step   "+str(self.ini)+ " Elev: "+str(round(self.res_ele,2)), transform=caax.transAxes, va='bottom',ha='right')
                     self.ini= self.ini+1

schainpy/scripts/basic_proc00003.py +2 -2

             print("N° Profiles   : ", n)
             #---------------------------------------------------------------------------------------
             plot_rti    = 0
-            plot_ppi    = 0
+            plot_ppi    = 1
             integration = 1
-            save        = 1
+            save        = 0
             #---------------------------RANGO DE PLOTEO----------------------------------
             dBmin = '1'
             dBmax = '85'

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