schain-jc Commit - r588:c1c257c97d52

parametrs plotting: SNR and DOP arguments have been added

Miguel Valdez -

r588:c1c257c97d52

parent child

Context file:

r588:c1c257c97d52

Collapse all files

schainpy/model/graphics/jroplot_parameters.py +32 -16

                         xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None,timerange=None,
                         parameterIndex = None, onlyPositive = False,
                         SNRthresh = -numpy.inf, SNR = True, SNRmin = None, SNRmax = None,
+                        DOP = True,
                         zlabel = "", parameterName = "", parameterObject = "data_param",
                         save=False, figpath='./', lastone=0,figfile=None, ftp=False, wr_period=1, show=True,
                         server=None, folder=None, username=None, password=None,
                     nchan = len(channelIndexList) #Number of channels being plotted
+                    if nchan < 1:
+                        return
+                    nGraphsByChannel = 0
+                    if SNR:
+                        nGraphsByChannel += 1
+                    if DOP:
+                        nGraphsByChannel += 1
+                    if nGraphsByChannel < 1:
+                        return
+                    nplots = nGraphsByChannel*nchan
                     if timerange != None:
                         self.timerange = timerange
                     #tmax = None
                     if parameterIndex == None:
                         parameterIndex = 1
                     x = dataOut.getTimeRange1()
                     y = dataOut.heightList
                     z = data_param[channelIndexList,parameterIndex,:].copy()
                     zRange = dataOut.abscissaList
-                    nplots = z.shape[0]    #Number of wind dimensions estimated
+            #         nChannels = z.shape[0]    #Number of wind dimensions estimated
             #        thisDatetime = dataOut.datatime
                     if dataOut.data_SNR != None:
                     xlabel = ""
                     ylabel = "Range (Km)"
-                    if SNR: nplots = 2*nplots
                     if onlyPositive:
                         colormap = "jet"
                         zmin = 0
                         if zmin == None: zmin = numpy.nanmin(zRange)
                         if zmax == None: zmax = numpy.nanmax(zRange)
-                        if SNR != None:
+                        if SNR:
                             if SNRmin == None:  SNRmin = numpy.nanmin(SNRdB)
                             if SNRmax == None:  SNRmax = numpy.nanmax(SNRdB)
                         x[1] = self.xmax
                     for i in range(nchan):
-                        if SNR: j = 2*i
-                        else: j = i
-                        title = "%s Channel %d: %s" %(parameterName, channelIndexList[i]+1, thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))
+                        j = nGraphsByChannel*i
-                        if ((dataOut.azimuth!=None) and (dataOut.zenith!=None)):
+                        if DOP:
-                            title = title + '_' + 'azimuth,zenith=%2.2f,%2.2f'%(dataOut.azimuth, dataOut.zenith)
+                            title = "%s Channel %d: %s" %(parameterName, channelIndexList[i]+1, thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))
-                        axes = self.axesList[j*self.__nsubplots]
-                        z1 = z[i,:].reshape((1,-1))
+                            if ((dataOut.azimuth!=None) and (dataOut.zenith!=None)):
-                        axes.pcolorbuffer(x, y, z1,
+                                title = title + '_' + 'azimuth,zenith=%2.2f,%2.2f'%(dataOut.azimuth, dataOut.zenith)
-                                    xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zmin, zmax=zmax,
+                            axes = self.axesList[j]
-                                    xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,colormap=colormap,
+                            z1 = z[i,:].reshape((1,-1))
-                                    ticksize=9, cblabel=zlabel, cbsize="1%")
+                            axes.pcolorbuffer(x, y, z1,
+                                        xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zmin, zmax=zmax,
+                                        xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,colormap=colormap,
+                                        ticksize=9, cblabel=zlabel, cbsize="1%")
                         if SNR:
                             title = "Channel %d Signal Noise Ratio (SNR): %s" %(channelIndexList[i]+1, thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))
-                            axes = self.axesList[(j + 1)*self.__nsubplots]
+                            axes = self.axesList[(j + nGraphsByChannel-1)]
                             z1 = SNRdB[i,:].reshape((1,-1))
                             axes.pcolorbuffer(x, y, z1,
                                     xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax,

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