import os import datetime import numpy from figure import Figure, isRealtime class MomentsPlot(Figure): isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'prm' def __init__(self): self.isConfig = False self.__nsubplots = 1 self.WIDTH = 280 self.HEIGHT = 250 self.WIDTHPROF = 120 self.HEIGHTPROF = 0 self.counter_imagwr = 0 self.PLOT_CODE = 1 self.FTP_WEI = None self.EXP_CODE = None self.SUB_EXP_CODE = None self.PLOT_POS = None def getSubplots(self): ncol = int(numpy.sqrt(self.nplots)+0.9) nrow = int(self.nplots*1./ncol + 0.9) return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=True, show=True): self.__showprofile = showprofile self.nplots = nplots ncolspan = 1 colspan = 1 if showprofile: ncolspan = 3 colspan = 2 self.__nsubplots = 2 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH + self.WIDTHPROF, heightplot = self.HEIGHT + self.HEIGHTPROF, show=show) nrow, ncol = self.getSubplots() counter = 0 for y in range(nrow): for x in range(ncol): if counter >= self.nplots: break self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan, colspan, 1) if showprofile: self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan+colspan, 1, 1) counter += 1 def run(self, dataOut, id, wintitle="", channelList=None, showprofile=True, xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, save=False, figpath='', figfile=None, show=True, ftp=False, wr_period=1, server=None, folder=None, username=None, password=None, ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0, realtime=False): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, ymin : None, ymax : None, zmin : None, zmax : None """ if dataOut.flagNoData: return None if realtime: if not(isRealtime(utcdatatime = dataOut.utctime)): print 'Skipping this plot function' return if channelList == None: channelIndexList = dataOut.channelIndexList else: channelIndexList = [] for channel in channelList: if channel not in dataOut.channelList: raise ValueError, "Channel %d is not in dataOut.channelList" channelIndexList.append(dataOut.channelList.index(channel)) factor = dataOut.normFactor x = dataOut.abscissaList y = dataOut.heightList z = dataOut.data_pre[channelIndexList,:,:]/factor z = numpy.where(numpy.isfinite(z), z, numpy.NAN) avg = numpy.average(z, axis=1) noise = dataOut.noise/factor zdB = 10*numpy.log10(z) avgdB = 10*numpy.log10(avg) noisedB = 10*numpy.log10(noise) #thisDatetime = dataOut.datatime thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) title = wintitle + " Parameters" xlabel = "Velocity (m/s)" ylabel = "Range (Km)" if not self.isConfig: nplots = len(channelIndexList) self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) if xmin == None: xmin = numpy.nanmin(x) if xmax == None: xmax = numpy.nanmax(x) if ymin == None: ymin = numpy.nanmin(y) if ymax == None: ymax = numpy.nanmax(y) if zmin == None: zmin = numpy.nanmin(avgdB)*0.9 if zmax == None: zmax = numpy.nanmax(avgdB)*0.9 self.FTP_WEI = ftp_wei self.EXP_CODE = exp_code self.SUB_EXP_CODE = sub_exp_code self.PLOT_POS = plot_pos self.isConfig = True self.setWinTitle(title) for i in range(self.nplots): str_datetime = '%s %s'%(thisDatetime.strftime("%Y/%m/%d"),thisDatetime.strftime("%H:%M:%S")) title = "Channel %d: %4.2fdB: %s" %(dataOut.channelList[i]+1, noisedB[i], str_datetime) axes = self.axesList[i*self.__nsubplots] axes.pcolor(x, y, zdB[i,:,:], xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax, zmin=zmin, zmax=zmax, xlabel=xlabel, ylabel=ylabel, title=title, ticksize=9, cblabel='') #Mean Line mean = dataOut.data_param[i, 1, :] axes.addpline(mean, y, idline=0, color="black", linestyle="solid", lw=1) if self.__showprofile: axes = self.axesList[i*self.__nsubplots +1] axes.pline(avgdB[i], y, xmin=zmin, xmax=zmax, ymin=ymin, ymax=ymax, xlabel='dB', ylabel='', title='', ytick_visible=False, grid='x') noiseline = numpy.repeat(noisedB[i], len(y)) axes.addpline(noiseline, y, idline=1, color="black", linestyle="dashed", lw=2) self.draw() if figfile == None: str_datetime = thisDatetime.strftime("%Y%m%d_%H%M%S") figfile = self.getFilename(name = str_datetime) if figpath != '': self.counter_imagwr += 1 if (self.counter_imagwr>=wr_period): # store png plot to local folder self.saveFigure(figpath, figfile) # store png plot to FTP server according to RT-Web format name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) ftp_filename = os.path.join(figpath, name) self.saveFigure(figpath, ftp_filename) self.counter_imagwr = 0 class SkyMapPlot(Figure): __isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'prm' def __init__(self): self.__isConfig = False self.__nsubplots = 1 # self.WIDTH = 280 # self.HEIGHT = 250 self.WIDTH = 600 self.HEIGHT = 600 self.WIDTHPROF = 120 self.HEIGHTPROF = 0 self.counter_imagwr = 0 self.PLOT_CODE = 1 self.FTP_WEI = None self.EXP_CODE = None self.SUB_EXP_CODE = None self.PLOT_POS = None def getSubplots(self): ncol = int(numpy.sqrt(self.nplots)+0.9) nrow = int(self.nplots*1./ncol + 0.9) return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=False, show=True): self.__showprofile = showprofile self.nplots = nplots ncolspan = 1 colspan = 1 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH, #+ self.WIDTHPROF, heightplot = self.HEIGHT,# + self.HEIGHTPROF, show=show) nrow, ncol = 1,1 counter = 0 x = 0 y = 0 self.addAxes(1, 1, 0, 0, 1, 1, True) def run(self, dataOut, id, wintitle="", channelList=None, showprofile=False, xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, save=False, figpath='./', figfile=None, show=True, ftp=False, wr_period=1, server=None, folder=None, username=None, password=None, ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0, realtime=False): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, ymin : None, ymax : None, zmin : None, zmax : None """ arrayParameters = dataOut.data_param error = arrayParameters[:,-1] indValid = numpy.where(error == 0)[0] finalMeteor = arrayParameters[indValid,:] finalAzimuth = finalMeteor[:,4] finalZenith = finalMeteor[:,5] x = finalAzimuth*numpy.pi/180 y = finalZenith #thisDatetime = dataOut.datatime thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) title = wintitle + " Parameters" xlabel = "Zonal Zenith Angle (deg) " ylabel = "Meridional Zenith Angle (deg)" if not self.__isConfig: nplots = 1 self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) self.FTP_WEI = ftp_wei self.EXP_CODE = exp_code self.SUB_EXP_CODE = sub_exp_code self.PLOT_POS = plot_pos self.name = thisDatetime.strftime("%Y%m%d_%H%M%S") self.firstdate = '%s %s'%(thisDatetime.strftime("%Y/%m/%d"),thisDatetime.strftime("%H:%M:%S")) self.__isConfig = True self.setWinTitle(title) i = 0 str_datetime = '%s %s'%(thisDatetime.strftime("%Y/%m/%d"),thisDatetime.strftime("%H:%M:%S")) axes = self.axesList[i*self.__nsubplots] nevents = axes.x_buffer.shape[0] + x.shape[0] title = "Meteor Detection Sky Map\n %s - %s \n Number of events: %5.0f\n" %(self.firstdate,str_datetime,nevents) axes.polar(x, y, title=title, xlabel=xlabel, ylabel=ylabel, ticksize=9, cblabel='') self.draw() if save: self.counter_imagwr += 1 if (self.counter_imagwr==wr_period): if figfile == None: figfile = self.getFilename(name = self.name) self.saveFigure(figpath, figfile) if ftp: #provisionalmente envia archivos en el formato de la web en tiempo real name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) path = '%s%03d' %(self.PREFIX, self.id) ftp_file = os.path.join(path,'ftp','%s.png'%name) self.saveFigure(figpath, ftp_file) ftp_filename = os.path.join(figpath,ftp_file) try: self.sendByFTP(ftp_filename, server, folder, username, password) except: self.counter_imagwr = 0 raise ValueError, 'Error FTP' self.counter_imagwr = 0 class WindProfilerPlot(Figure): __isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'wind' def __init__(self): self.timerange = 2*60*60 self.__isConfig = False self.__nsubplots = 1 self.WIDTH = 800 self.HEIGHT = 150 self.WIDTHPROF = 120 self.HEIGHTPROF = 0 self.counter_imagwr = 0 self.PLOT_CODE = 0 self.FTP_WEI = None self.EXP_CODE = None self.SUB_EXP_CODE = None self.PLOT_POS = None self.tmin = None self.tmax = None self.xmin = None self.xmax = None self.figfile = None def getSubplots(self): ncol = 1 nrow = self.nplots return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=True, show=True): self.__showprofile = showprofile self.nplots = nplots ncolspan = 1 colspan = 1 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH + self.WIDTHPROF, heightplot = self.HEIGHT + self.HEIGHTPROF, show=show) nrow, ncol = self.getSubplots() counter = 0 for y in range(nrow): if counter >= self.nplots: break self.addAxes(nrow, ncol*ncolspan, y, 0, colspan, 1) counter += 1 def run(self, dataOut, id, wintitle="", channelList=None, xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, zmax_ver = None, zmin_ver = None, SNRmin = None, SNRmax = None, timerange=None, SNRthresh = None, save=False, figpath='', lastone=0,figfile=None, ftp=False, wr_period=1, show=True, server=None, folder=None, username=None, password=None, ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, ymin : None, ymax : None, zmin : None, zmax : None """ if channelList == None: channelIndexList = dataOut.channelIndexList else: channelIndexList = [] for channel in channelList: if channel not in dataOut.channelList: raise ValueError, "Channel %d is not in dataOut.channelList" channelIndexList.append(dataOut.channelList.index(channel)) if timerange != None: self.timerange = timerange tmin = None tmax = None x = dataOut.getTimeRange1() # y = dataOut.heightList y = dataOut.heightList z = dataOut.data_output.copy() nplots = z.shape[0] #Number of wind dimensions estimated nplotsw = nplots #If there is a SNR function defined if dataOut.data_SNR != None: nplots += 1 SNR = dataOut.data_SNR SNRavg = numpy.average(SNR, axis=0) SNRdB = 10*numpy.log10(SNR) SNRavgdB = 10*numpy.log10(SNRavg) if SNRthresh == None: SNRthresh = -5.0 ind = numpy.where(SNRavg < 10**(SNRthresh/10))[0] for i in range(nplotsw): z[i,ind] = numpy.nan showprofile = False # thisDatetime = dataOut.datatime thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) title = wintitle + "Wind" xlabel = "" ylabel = "Range (Km)" if not self.__isConfig: self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange) if ymin == None: ymin = numpy.nanmin(y) if ymax == None: ymax = numpy.nanmax(y) if zmax == None: zmax = numpy.nanmax(abs(z[range(2),:])) #if numpy.isnan(zmax): zmax = 50 if zmin == None: zmin = -zmax if nplotsw == 3: if zmax_ver == None: zmax_ver = numpy.nanmax(abs(z[2,:])) if zmin_ver == None: zmin_ver = -zmax_ver if dataOut.data_SNR != None: if SNRmin == None: SNRmin = numpy.nanmin(SNRavgdB) if SNRmax == None: SNRmax = numpy.nanmax(SNRavgdB) self.FTP_WEI = ftp_wei self.EXP_CODE = exp_code self.SUB_EXP_CODE = sub_exp_code self.PLOT_POS = plot_pos self.name = thisDatetime.strftime("%Y%m%d_%H%M%S") self.__isConfig = True self.setWinTitle(title) if ((self.xmax - x[1]) < (x[1]-x[0])): x[1] = self.xmax strWind = ['Zonal', 'Meridional', 'Vertical'] strCb = ['Velocity (m/s)','Velocity (m/s)','Velocity (cm/s)'] zmaxVector = [zmax, zmax, zmax_ver] zminVector = [zmin, zmin, zmin_ver] windFactor = [1,1,100] for i in range(nplotsw): title = "%s Wind: %s" %(strWind[i], thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) axes = self.axesList[i*self.__nsubplots] z1 = z[i,:].reshape((1,-1))*windFactor[i] axes.pcolorbuffer(x, y, z1, xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zminVector[i], zmax=zmaxVector[i], xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True, ticksize=9, cblabel=strCb[i], cbsize="1%", colormap="RdBu_r" ) if dataOut.data_SNR != None: i += 1 title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) axes = self.axesList[i*self.__nsubplots] SNRavgdB = SNRavgdB.reshape((1,-1)) axes.pcolorbuffer(x, y, SNRavgdB, xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax, xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True, ticksize=9, cblabel='', cbsize="1%", colormap="jet") self.draw() if self.figfile == None: str_datetime = thisDatetime.strftime("%Y%m%d_%H%M%S") self.figfile = self.getFilename(name = str_datetime) if figpath != '': self.counter_imagwr += 1 if (self.counter_imagwr>=wr_period): # store png plot to local folder self.saveFigure(figpath, self.figfile) # store png plot to FTP server according to RT-Web format name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) ftp_filename = os.path.join(figpath, name) self.saveFigure(figpath, ftp_filename) self.counter_imagwr = 0 if x[1] >= self.axesList[0].xmax: self.counter_imagwr = wr_period self.__isConfig = False self.figfile = None class ParametersPlot(Figure): __isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'prm' def __init__(self): self.timerange = 2*60*60 self.__isConfig = False self.__nsubplots = 1 self.WIDTH = 800 self.HEIGHT = 150 self.WIDTHPROF = 120 self.HEIGHTPROF = 0 self.counter_imagwr = 0 self.PLOT_CODE = 0 self.FTP_WEI = None self.EXP_CODE = None self.SUB_EXP_CODE = None self.PLOT_POS = None self.tmin = None self.tmax = None self.xmin = None self.xmax = None self.figfile = None def getSubplots(self): ncol = 1 nrow = self.nplots return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=True, show=True): self.__showprofile = showprofile self.nplots = nplots ncolspan = 1 colspan = 1 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH + self.WIDTHPROF, heightplot = self.HEIGHT + self.HEIGHTPROF, show=show) nrow, ncol = self.getSubplots() counter = 0 for y in range(nrow): for x in range(ncol): if counter >= self.nplots: break self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan, colspan, 1) if showprofile: self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan+colspan, 1, 1) counter += 1 def run(self, dataOut, id, wintitle="", channelList=None, showprofile=False, 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, 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, ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, ymin : None, ymax : None, zmin : None, zmax : None """ data_param = getattr(dataOut, parameterObject) if channelList == None: channelIndexList = numpy.arange(data_param.shape[0]) else: channelIndexList = numpy.array(channelList) nchan = len(channelIndexList) #Number of channels being plotted if timerange != None: self.timerange = timerange #tmin = None #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 # thisDatetime = dataOut.datatime if dataOut.data_SNR != None: SNRarray = dataOut.data_SNR[channelIndexList,:] SNRdB = 10*numpy.log10(SNRarray) # SNRavgdB = 10*numpy.log10(SNRavg) ind = numpy.where(SNRdB < 10**(SNRthresh/10)) z[ind] = numpy.nan thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) title = wintitle + " Parameters Plot" #: %s" %(thisDatetime.strftime("%d-%b-%Y")) xlabel = "" ylabel = "Range (Km)" if SNR: nplots = 2*nplots if onlyPositive: colormap = "jet" zmin = 0 else: colormap = "RdBu_r" if not self.__isConfig: self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange) if ymin == None: ymin = numpy.nanmin(y) if ymax == None: ymax = numpy.nanmax(y) if zmin == None: zmin = numpy.nanmin(zRange) if zmax == None: zmax = numpy.nanmax(zRange) if SNR != None: if SNRmin == None: SNRmin = numpy.nanmin(SNRdB) if SNRmax == None: SNRmax = numpy.nanmax(SNRdB) self.FTP_WEI = ftp_wei self.EXP_CODE = exp_code self.SUB_EXP_CODE = sub_exp_code self.PLOT_POS = plot_pos self.name = thisDatetime.strftime("%Y%m%d_%H%M%S") self.__isConfig = True self.figfile = figfile self.setWinTitle(title) if ((self.xmax - x[1]) < (x[1]-x[0])): 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")) if ((dataOut.azimuth!=None) and (dataOut.zenith!=None)): title = title + '_' + 'azimuth,zenith=%2.2f,%2.2f'%(dataOut.azimuth, dataOut.zenith) axes = self.axesList[j*self.__nsubplots] z1 = z[i,:].reshape((1,-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] 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, xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,colormap="jet", ticksize=9, cblabel=zlabel, cbsize="1%") self.draw() if self.figfile == None: str_datetime = thisDatetime.strftime("%Y%m%d_%H%M%S") self.figfile = self.getFilename(name = str_datetime) if figpath != '': self.counter_imagwr += 1 if (self.counter_imagwr>=wr_period): # store png plot to local folder self.saveFigure(figpath, self.figfile) # store png plot to FTP server according to RT-Web format name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) ftp_filename = os.path.join(figpath, name) self.saveFigure(figpath, ftp_filename) self.counter_imagwr = 0 if x[1] >= self.axesList[0].xmax: self.counter_imagwr = wr_period self.__isConfig = False self.figfile = None class SpectralFittingPlot(Figure): __isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'prm' N = None ippSeconds = None def __init__(self): self.__isConfig = False self.__nsubplots = 1 self.WIDTH = 450 self.HEIGHT = 250 self.WIDTHPROF = 0 self.HEIGHTPROF = 0 def getSubplots(self): ncol = int(numpy.sqrt(self.nplots)+0.9) nrow = int(self.nplots*1./ncol + 0.9) return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=False, show=True): showprofile = False self.__showprofile = showprofile self.nplots = nplots ncolspan = 5 colspan = 4 if showprofile: ncolspan = 5 colspan = 4 self.__nsubplots = 2 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH + self.WIDTHPROF, heightplot = self.HEIGHT + self.HEIGHTPROF, show=show) nrow, ncol = self.getSubplots() counter = 0 for y in range(nrow): for x in range(ncol): if counter >= self.nplots: break self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan, colspan, 1) if showprofile: self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan+colspan, 1, 1) counter += 1 def run(self, dataOut, id, cutHeight=None, fit=False, wintitle="", channelList=None, showprofile=True, xmin=None, xmax=None, ymin=None, ymax=None, save=False, figpath='./', figfile=None, show=True): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, zmin : None, zmax : None """ if cutHeight==None: h=270 heightindex = numpy.abs(cutHeight - dataOut.heightList).argmin() cutHeight = dataOut.heightList[heightindex] factor = dataOut.normFactor x = dataOut.abscissaList[:-1] #y = dataOut.getHeiRange() z = dataOut.data_pre[:,:,heightindex]/factor z = numpy.where(numpy.isfinite(z), z, numpy.NAN) avg = numpy.average(z, axis=1) listChannels = z.shape[0] #Reconstruct Function if fit==True: groupArray = dataOut.groupList listChannels = groupArray.reshape((groupArray.size)) listChannels.sort() spcFitLine = numpy.zeros(z.shape) constants = dataOut.constants nGroups = groupArray.shape[0] nChannels = groupArray.shape[1] nProfiles = z.shape[1] for f in range(nGroups): groupChann = groupArray[f,:] p = dataOut.data_param[f,:,heightindex] # p = numpy.array([ 89.343967,0.14036615,0.17086219,18.89835291,1.58388365,1.55099167]) fitLineAux = dataOut.library.modelFunction(p, constants)*nProfiles fitLineAux = fitLineAux.reshape((nChannels,nProfiles)) spcFitLine[groupChann,:] = fitLineAux # spcFitLine = spcFitLine/factor z = z[listChannels,:] spcFitLine = spcFitLine[listChannels,:] spcFitLinedB = 10*numpy.log10(spcFitLine) zdB = 10*numpy.log10(z) #thisDatetime = dataOut.datatime thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) title = wintitle + " Doppler Spectra: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S")) xlabel = "Velocity (m/s)" ylabel = "Spectrum" if not self.__isConfig: nplots = listChannels.size self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) if xmin == None: xmin = numpy.nanmin(x) if xmax == None: xmax = numpy.nanmax(x) if ymin == None: ymin = numpy.nanmin(zdB) if ymax == None: ymax = numpy.nanmax(zdB)+2 self.__isConfig = True self.setWinTitle(title) for i in range(self.nplots): # title = "Channel %d: %4.2fdB" %(dataOut.channelList[i]+1, noisedB[i]) title = "Height %4.1f km\nChannel %d:" %(cutHeight, listChannels[i]+1) axes = self.axesList[i*self.__nsubplots] if fit == False: axes.pline(x, zdB[i,:], xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax, xlabel=xlabel, ylabel=ylabel, title=title ) if fit == True: fitline=spcFitLinedB[i,:] y=numpy.vstack([zdB[i,:],fitline] ) legendlabels=['Data','Fitting'] axes.pmultilineyaxis(x, y, xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax, xlabel=xlabel, ylabel=ylabel, title=title, legendlabels=legendlabels, marker=None, linestyle='solid', grid='both') self.draw() if save: date = thisDatetime.strftime("%Y%m%d_%H%M%S") if figfile == None: figfile = self.getFilename(name = date) self.saveFigure(figpath, figfile) class EWDriftsPlot(Figure): __isConfig = None __nsubplots = None WIDTHPROF = None HEIGHTPROF = None PREFIX = 'drift' def __init__(self): self.timerange = 2*60*60 self.isConfig = False self.__nsubplots = 1 self.WIDTH = 800 self.HEIGHT = 150 self.WIDTHPROF = 120 self.HEIGHTPROF = 0 self.counter_imagwr = 0 self.PLOT_CODE = 0 self.FTP_WEI = None self.EXP_CODE = None self.SUB_EXP_CODE = None self.PLOT_POS = None self.tmin = None self.tmax = None self.xmin = None self.xmax = None self.figfile = None def getSubplots(self): ncol = 1 nrow = self.nplots return nrow, ncol def setup(self, id, nplots, wintitle, showprofile=True, show=True): self.__showprofile = showprofile self.nplots = nplots ncolspan = 1 colspan = 1 self.createFigure(id = id, wintitle = wintitle, widthplot = self.WIDTH + self.WIDTHPROF, heightplot = self.HEIGHT + self.HEIGHTPROF, show=show) nrow, ncol = self.getSubplots() counter = 0 for y in range(nrow): if counter >= self.nplots: break self.addAxes(nrow, ncol*ncolspan, y, 0, colspan, 1) counter += 1 def run(self, dataOut, id, wintitle="", channelList=None, xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None, zmaxVertical = None, zminVertical = None, zmaxZonal = None, zminZonal = None, timerange=None, SNRthresh = -numpy.inf, SNRmin = None, SNRmax = None, SNR_1 = False, save=False, figpath='', lastone=0,figfile=None, ftp=False, wr_period=1, show=True, server=None, folder=None, username=None, password=None, ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0): """ Input: dataOut : id : wintitle : channelList : showProfile : xmin : None, xmax : None, ymin : None, ymax : None, zmin : None, zmax : None """ if timerange != None: self.timerange = timerange tmin = None tmax = None x = dataOut.getTimeRange1() # y = dataOut.heightList y = dataOut.heightList z = dataOut.data_output nplots = z.shape[0] #Number of wind dimensions estimated nplotsw = nplots #If there is a SNR function defined if dataOut.data_SNR != None: nplots += 1 SNR = dataOut.data_SNR if SNR_1: SNR += 1 SNRavg = numpy.average(SNR, axis=0) SNRdB = 10*numpy.log10(SNR) SNRavgdB = 10*numpy.log10(SNRavg) ind = numpy.where(SNRavg < 10**(SNRthresh/10))[0] for i in range(nplotsw): z[i,ind] = numpy.nan showprofile = False # thisDatetime = dataOut.datatime thisDatetime = datetime.datetime.utcfromtimestamp(x[1]) title = wintitle + " EW Drifts" xlabel = "" ylabel = "Height (Km)" if not self.__isConfig: self.setup(id=id, nplots=nplots, wintitle=wintitle, showprofile=showprofile, show=show) self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange) if ymin == None: ymin = numpy.nanmin(y) if ymax == None: ymax = numpy.nanmax(y) if zmaxZonal == None: zmaxZonal = numpy.nanmax(abs(z[0,:])) if zminZonal == None: zminZonal = -zmaxZonal if zmaxVertical == None: zmaxVertical = numpy.nanmax(abs(z[1,:])) if zminVertical == None: zminVertical = -zmaxVertical if dataOut.data_SNR != None: if SNRmin == None: SNRmin = numpy.nanmin(SNRavgdB) if SNRmax == None: SNRmax = numpy.nanmax(SNRavgdB) self.FTP_WEI = ftp_wei self.EXP_CODE = exp_code self.SUB_EXP_CODE = sub_exp_code self.PLOT_POS = plot_pos self.name = thisDatetime.strftime("%Y%m%d_%H%M%S") self.__isConfig = True self.setWinTitle(title) if ((self.xmax - x[1]) < (x[1]-x[0])): x[1] = self.xmax strWind = ['Zonal','Vertical'] strCb = 'Velocity (m/s)' zmaxVector = [zmaxZonal, zmaxVertical] zminVector = [zminZonal, zminVertical] for i in range(nplotsw): title = "%s Drifts: %s" %(strWind[i], thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) axes = self.axesList[i*self.__nsubplots] z1 = z[i,:].reshape((1,-1)) axes.pcolorbuffer(x, y, z1, xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zminVector[i], zmax=zmaxVector[i], xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True, ticksize=9, cblabel=strCb, cbsize="1%", colormap="RdBu_r") if dataOut.data_SNR != None: i += 1 if SNR_1: title = "Signal Noise Ratio + 1 (SNR+1): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) else: title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) axes = self.axesList[i*self.__nsubplots] SNRavgdB = SNRavgdB.reshape((1,-1)) axes.pcolorbuffer(x, y, SNRavgdB, xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax, xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True, ticksize=9, cblabel='', cbsize="1%", colormap="jet") self.draw() if self.figfile == None: str_datetime = thisDatetime.strftime("%Y%m%d_%H%M%S") self.figfile = self.getFilename(name = str_datetime) if figpath != '': self.counter_imagwr += 1 if (self.counter_imagwr>=wr_period): # store png plot to local folder self.saveFigure(figpath, self.figfile) # store png plot to FTP server according to RT-Web format name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) ftp_filename = os.path.join(figpath, name) self.saveFigure(figpath, ftp_filename) self.counter_imagwr = 0 if x[1] >= self.axesList[0].xmax: self.counter_imagwr = wr_period self.__isConfig = False self.figfile = None