| @@ -455,7 +455,7 class WindProfilerPlot(Figure): | |||
|
|
455 | 455 | # y = dataOut.heightRange |
|
|
456 | 456 | y = dataOut.heightRange |
|
|
457 | 457 | |
|
|
458 | z = dataOut.winds | |
|
|
458 | z = dataOut.winds.copy() | |
|
|
459 | 459 | nplots = z.shape[0] #Number of wind dimensions estimated |
|
|
460 | 460 | nplotsw = nplots |
|
|
461 | 461 | |
| @@ -574,3 +574,220 class WindProfilerPlot(Figure): | |||
|
|
574 | 574 | self.counter_imagwr = wr_period |
|
|
575 | 575 | self.__isConfig = False |
|
|
576 | 576 | self.figfile = None |
|
|
577 | ||
|
|
578 | ||
|
|
579 | class RadialVelocityPlot(Figure): | |
|
|
580 | ||
|
|
581 | __isConfig = None | |
|
|
582 | __nsubplots = None | |
|
|
583 | ||
|
|
584 | WIDTHPROF = None | |
|
|
585 | HEIGHTPROF = None | |
|
|
586 | PREFIX = 'rti' | |
|
|
587 | ||
|
|
588 | def __init__(self): | |
|
|
589 | ||
|
|
590 | self.timerange = 2*60*60 | |
|
|
591 | self.__isConfig = False | |
|
|
592 | self.__nsubplots = 1 | |
|
|
593 | ||
|
|
594 | self.WIDTH = 800 | |
|
|
595 | self.HEIGHT = 150 | |
|
|
596 | self.WIDTHPROF = 120 | |
|
|
597 | self.HEIGHTPROF = 0 | |
|
|
598 | self.counter_imagwr = 0 | |
|
|
599 | ||
|
|
600 | self.PLOT_CODE = 0 | |
|
|
601 | self.FTP_WEI = None | |
|
|
602 | self.EXP_CODE = None | |
|
|
603 | self.SUB_EXP_CODE = None | |
|
|
604 | self.PLOT_POS = None | |
|
|
605 | self.tmin = None | |
|
|
606 | self.tmax = None | |
|
|
607 | ||
|
|
608 | self.xmin = None | |
|
|
609 | self.xmax = None | |
|
|
610 | ||
|
|
611 | self.figfile = None | |
|
|
612 | ||
|
|
613 | def getSubplots(self): | |
|
|
614 | ||
|
|
615 | ncol = 1 | |
|
|
616 | nrow = self.nplots | |
|
|
617 | ||
|
|
618 | return nrow, ncol | |
|
|
619 | ||
|
|
620 | def setup(self, id, nplots, wintitle, showprofile=True, show=True): | |
|
|
621 | ||
|
|
622 | self.__showprofile = showprofile | |
|
|
623 | self.nplots = nplots | |
|
|
624 | ||
|
|
625 | ncolspan = 1 | |
|
|
626 | colspan = 1 | |
|
|
627 | ||
|
|
628 | self.createFigure(id = id, | |
|
|
629 | wintitle = wintitle, | |
|
|
630 | widthplot = self.WIDTH + self.WIDTHPROF, | |
|
|
631 | heightplot = self.HEIGHT + self.HEIGHTPROF, | |
|
|
632 | show=show) | |
|
|
633 | ||
|
|
634 | nrow, ncol = self.getSubplots() | |
|
|
635 | ||
|
|
636 | counter = 0 | |
|
|
637 | for y in range(nrow): | |
|
|
638 | for x in range(ncol): | |
|
|
639 | ||
|
|
640 | if counter >= self.nplots: | |
|
|
641 | break | |
|
|
642 | ||
|
|
643 | self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan, colspan, 1) | |
|
|
644 | ||
|
|
645 | if showprofile: | |
|
|
646 | self.addAxes(nrow, ncol*ncolspan, y, x*ncolspan+colspan, 1, 1) | |
|
|
647 | ||
|
|
648 | counter += 1 | |
|
|
649 | ||
|
|
650 | def run(self, dataOut, id, wintitle="", channelList=None, showprofile=False, | |
|
|
651 | xmin=None, xmax=None, ymin=None, ymax=None, zmin=None, zmax=None,timerange=None, | |
|
|
652 | SNRmin = None, SNRmax = None, SNRthresh = None, paramIndex = None, | |
|
|
653 | save=False, figpath='', lastone=0,figfile=None, ftp=False, wr_period=1, show=True, | |
|
|
654 | server=None, folder=None, username=None, password=None, | |
|
|
655 | ftp_wei=0, exp_code=0, sub_exp_code=0, plot_pos=0): | |
|
|
656 | ||
|
|
657 | """ | |
|
|
658 | ||
|
|
659 | Input: | |
|
|
660 | dataOut : | |
|
|
661 | id : | |
|
|
662 | wintitle : | |
|
|
663 | channelList : | |
|
|
664 | showProfile : | |
|
|
665 | xmin : None, | |
|
|
666 | xmax : None, | |
|
|
667 | ymin : None, | |
|
|
668 | ymax : None, | |
|
|
669 | zmin : None, | |
|
|
670 | zmax : None | |
|
|
671 | """ | |
|
|
672 | ||
|
|
673 | if channelList == None: | |
|
|
674 | channelIndexList = dataOut.channelIndexList | |
|
|
675 | else: | |
|
|
676 | channelIndexList = [] | |
|
|
677 | for channel in channelList: | |
|
|
678 | if channel not in dataOut.channelList: | |
|
|
679 | raise ValueError, "Channel %d is not in dataOut.channelList" | |
|
|
680 | channelIndexList.append(dataOut.channelList.index(channel)) | |
|
|
681 | ||
|
|
682 | if timerange != None: | |
|
|
683 | self.timerange = timerange | |
|
|
684 | ||
|
|
685 | #tmin = None | |
|
|
686 | #tmax = None | |
|
|
687 | if paramIndex == None: | |
|
|
688 | paramIndex = 1 | |
|
|
689 | x = dataOut.getTimeRange1() | |
|
|
690 | y = dataOut.heightRange | |
|
|
691 | z = dataOut.data_param[channelIndexList,paramIndex,:].copy() | |
|
|
692 | ||
|
|
693 | zRange = dataOut.abscissaRange | |
|
|
694 | nplots = z.shape[0] #Number of wind dimensions estimated | |
|
|
695 | nplotsw = nplots | |
|
|
696 | ||
|
|
697 | if dataOut.SNR != None: | |
|
|
698 | nplots += 1 | |
|
|
699 | SNR = dataOut.SNR | |
|
|
700 | SNRavg = numpy.average(SNR, axis=0) | |
|
|
701 | ||
|
|
702 | SNRdB = 10*numpy.log10(SNR) | |
|
|
703 | SNRavgdB = 10*numpy.log10(SNRavg) | |
|
|
704 | ||
|
|
705 | if SNRthresh == None: SNRthresh = -5.0 | |
|
|
706 | ind = numpy.where(SNRavg < 10**(SNRthresh/10))[0] | |
|
|
707 | ||
|
|
708 | for i in range(nplotsw): | |
|
|
709 | z[i,ind] = numpy.nan | |
|
|
710 | # thisDatetime = dataOut.datatime | |
|
|
711 | thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[1]) | |
|
|
712 | title = wintitle + " Radial Velocity" #: %s" %(thisDatetime.strftime("%d-%b-%Y")) | |
|
|
713 | xlabel = "" | |
|
|
714 | ylabel = "Height (Km)" | |
|
|
715 | ||
|
|
716 | if not self.__isConfig: | |
|
|
717 | ||
|
|
718 | self.setup(id=id, | |
|
|
719 | nplots=nplots, | |
|
|
720 | wintitle=wintitle, | |
|
|
721 | showprofile=showprofile, | |
|
|
722 | show=show) | |
|
|
723 | ||
|
|
724 | self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange) | |
|
|
725 | ||
|
|
726 | if ymin == None: ymin = numpy.nanmin(y) | |
|
|
727 | if ymax == None: ymax = numpy.nanmax(y) | |
|
|
728 | if zmin == None: zmin = numpy.nanmin(zRange) | |
|
|
729 | if zmax == None: zmax = numpy.nanmax(zRange) | |
|
|
730 | if dataOut.SNR != None: | |
|
|
731 | if SNRmin == None: SNRmin = numpy.nanmin(SNRavgdB) | |
|
|
732 | if SNRmax == None: SNRmax = numpy.nanmax(SNRavgdB) | |
|
|
733 | ||
|
|
734 | self.FTP_WEI = ftp_wei | |
|
|
735 | self.EXP_CODE = exp_code | |
|
|
736 | self.SUB_EXP_CODE = sub_exp_code | |
|
|
737 | self.PLOT_POS = plot_pos | |
|
|
738 | ||
|
|
739 | self.name = thisDatetime.strftime("%Y%m%d_%H%M%S") | |
|
|
740 | self.__isConfig = True | |
|
|
741 | self.figfile = figfile | |
|
|
742 | ||
|
|
743 | self.setWinTitle(title) | |
|
|
744 | ||
|
|
745 | if ((self.xmax - x[1]) < (x[1]-x[0])): | |
|
|
746 | x[1] = self.xmax | |
|
|
747 | ||
|
|
748 | for i in range(nplotsw): | |
|
|
749 | title = "Channel %d: %s" %(dataOut.channelList[i]+1, thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) | |
|
|
750 | if ((dataOut.azimuth!=None) and (dataOut.zenith!=None)): | |
|
|
751 | title = title + '_' + 'azimuth,zenith=%2.2f,%2.2f'%(dataOut.azimuth, dataOut.zenith) | |
|
|
752 | axes = self.axesList[i*self.__nsubplots] | |
|
|
753 | z1 = z[i,:].reshape((1,-1)) | |
|
|
754 | axes.pcolorbuffer(x, y, z1, | |
|
|
755 | xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zmin, zmax=zmax, | |
|
|
756 | xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,colormap="RdBu_r", | |
|
|
757 | ticksize=9, cblabel='', cbsize="1%") | |
|
|
758 | ||
|
|
759 | if dataOut.SNR != None: | |
|
|
760 | i += 1 | |
|
|
761 | title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S")) | |
|
|
762 | axes = self.axesList[i*self.__nsubplots] | |
|
|
763 | ||
|
|
764 | SNRavgdB = SNRavgdB.reshape((1,-1)) | |
|
|
765 | ||
|
|
766 | axes.pcolorbuffer(x, y, SNRavgdB, | |
|
|
767 | xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax, | |
|
|
768 | xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True, | |
|
|
769 | ticksize=9, cblabel='', cbsize="1%", colormap="jet") | |
|
|
770 | ||
|
|
771 | self.draw() | |
|
|
772 | ||
|
|
773 | if self.figfile == None: | |
|
|
774 | str_datetime = thisDatetime.strftime("%Y%m%d_%H%M%S") | |
|
|
775 | self.figfile = self.getFilename(name = str_datetime) | |
|
|
776 | ||
|
|
777 | if figpath != '': | |
|
|
778 | ||
|
|
779 | self.counter_imagwr += 1 | |
|
|
780 | if (self.counter_imagwr>=wr_period): | |
|
|
781 | # store png plot to local folder | |
|
|
782 | self.saveFigure(figpath, self.figfile) | |
|
|
783 | # store png plot to FTP server according to RT-Web format | |
|
|
784 | name = self.getNameToFtp(thisDatetime, self.FTP_WEI, self.EXP_CODE, self.SUB_EXP_CODE, self.PLOT_CODE, self.PLOT_POS) | |
|
|
785 | ftp_filename = os.path.join(figpath, name) | |
|
|
786 | self.saveFigure(figpath, ftp_filename) | |
|
|
787 | ||
|
|
788 | self.counter_imagwr = 0 | |
|
|
789 | ||
|
|
790 | if x[1] >= self.axesList[0].xmax: | |
|
|
791 | self.counter_imagwr = wr_period | |
|
|
792 | self.__isConfig = False | |
|
|
793 | self.figfile = None No newline at end of file | |
| @@ -91,6 +91,7 class ParametersProc(ProcessingUnit): | |||
|
|
91 | 91 | self.buffer = None |
|
|
92 | 92 | self.firstdatatime = None |
|
|
93 | 93 | self.profIndex = 0 |
|
|
94 | return | |
|
|
94 | 95 | |
|
|
95 | 96 | #---------------------- Spectra Data --------------------------- |
|
|
96 | 97 | |
| @@ -100,10 +101,6 class ParametersProc(ProcessingUnit): | |||
|
|
100 | 101 | self.dataOut.noise = self.dataIn.getNoise() |
|
|
101 | 102 | self.dataOut.normFactor = self.dataIn.normFactor |
|
|
102 | 103 | |
|
|
103 | self.__updateObjFromInput() | |
|
|
104 | self.dataOut.flagNoData = False | |
|
|
105 | self.firstdatatime = None | |
|
|
106 | ||
|
|
107 | 104 | #---------------------- Correlation Data --------------------------- |
|
|
108 | 105 | |
|
|
109 | 106 | if self.dataIn.type == "Correlation": |
| @@ -117,9 +114,12 class ParametersProc(ProcessingUnit): | |||
|
|
117 | 114 | self.dataOut.SNR = self.dataIn.SNR |
|
|
118 | 115 | self.dataOut.pairsList = self.dataIn.pairsList |
|
|
119 | 116 | |
|
|
117 | ||
|
|
120 | 118 |
|
|
|
121 | 119 |
|
|
|
122 | 120 |
|
|
|
121 | self.dataOut.initUtcTime = self.dataIn.ltctime | |
|
|
122 | self.dataOut.windsInterval = self.dataIn.timeInterval | |
|
|
123 | 123 | |
|
|
124 | 124 | #------------------- Get Moments ---------------------------------- |
|
|
125 | 125 | def GetMoments(self, channelList = None): |
| @@ -143,12 +143,14 class ParametersProc(ProcessingUnit): | |||
|
|
143 | 143 | |
|
|
144 | 144 | data_param = numpy.zeros((data.shape[0], 4, data.shape[2])) |
|
|
145 | 145 | |
|
|
146 |
if channelList== None: |
|
|
|
146 | if channelList== None: | |
|
|
147 | channelList = self.dataIn.channelList | |
|
|
148 | self.dataOut.channelList = channelList | |
|
|
147 | 149 | |
|
|
148 | 150 | for ind in channelList: |
|
|
149 | 151 | data_param[ind,:,:] = self.__calculateMoments(data[ind,:,:], absc, noise[ind]) |
|
|
150 | 152 | |
|
|
151 | self.dataOut.data_param = data_param[:,1:] | |
|
|
153 | self.dataOut.data_param = data_param[:,1:,:] | |
|
|
152 | 154 | self.dataOut.SNR = data_param[:,0] |
|
|
153 | 155 | return |
|
|
154 | 156 | |
| @@ -1461,8 +1463,6 class WindProfiler(Operation): | |||
|
|
1461 | 1463 | |
|
|
1462 | 1464 | velRadial0 = param[:,1,:] #Radial velocity |
|
|
1463 | 1465 | dataOut.winds, dataOut.heightRange, dataOut.SNR = self.techniqueDBS(velRadial0, theta_x, theta_y, azimuth, correctFactor, horizontalOnly, heightRange, SNR) #DBS Function |
|
|
1464 | dataOut.initUtcTime = dataOut.ltctime | |
|
|
1465 | dataOut.windsInterval = dataOut.timeInterval | |
|
|
1466 | 1466 | |
|
|
1467 | 1467 | elif technique == 'SA': |
|
|
1468 | 1468 | |
| @@ -28,8 +28,8 controllerObj.setup(id = '191', name='test01', description=desc) | |||
|
|
28 | 28 | #2014051 20 Feb 2014 |
|
|
29 | 29 | path = '/home/soporte/Data/MST/DBS/d2014051' |
|
|
30 | 30 | pathFigure = '/home/soporte/workspace/Graficos/DBS/prueba1/' |
|
|
31 |
xmin = '0 |
|
|
|
32 |
xmax = ' |
|
|
|
31 | xmin = '0' | |
|
|
32 | xmax = '7.5' | |
|
|
33 | 33 | startTime = '00:00:00' |
|
|
34 | 34 | filehdf5 = "DBS_2014051.hdf5" |
|
|
35 | 35 | |
| @@ -109,6 +109,17 opObj21.addParameter(name='figpath', value=pathFigure, format='str') | |||
|
|
109 | 109 | opObj21.addParameter(name='zmin', value='5', format='int') |
|
|
110 | 110 | opObj21.addParameter(name='zmax', value='90', format='int') |
|
|
111 | 111 | |
|
|
112 | opObj21 = procUnitConfObj2.addOperation(name='RadialVelocityPlot', optype='other') | |
|
|
113 | opObj21.addParameter(name='id', value='5', format='int') | |
|
|
114 | opObj21.addParameter(name='wintitle', value='Radial Velocity Plot', format='str') | |
|
|
115 | opObj21.addParameter(name='save', value='1', format='bool') | |
|
|
116 | opObj21.addParameter(name='figpath', value=pathFigure, format='str') | |
|
|
117 | opObj21.addParameter(name='SNRmin', value='-10', format='int') | |
|
|
118 | opObj21.addParameter(name='SNRmax', value='60', format='int') | |
|
|
119 | opObj21.addParameter(name='SNRthresh', value='0', format='float') | |
|
|
120 | opObj21.addParameter(name='xmin', value=xmin, format='float') | |
|
|
121 | opObj21.addParameter(name='xmax', value=xmax, format='float') | |
|
|
122 | ||
|
|
112 | 123 | opObj22 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other') |
|
|
113 | 124 | opObj22.addParameter(name='technique', value='DBS', format='str') |
|
|
114 | 125 | opObj22.addParameter(name='correctAzimuth', value='51.06', format='float') |
| @@ -19,20 +19,20 controllerObj.setup(id = '191', name='test01', description=desc) | |||
|
|
19 | 19 | #Experimentos |
|
|
20 | 20 | |
|
|
21 | 21 | #2014050 19 Feb 2014 |
|
|
22 |
|
|
|
|
23 |
|
|
|
|
24 |
|
|
|
|
25 |
|
|
|
|
26 |
|
|
|
|
27 |
|
|
|
|
22 | path = '/home/soporte/Data/MST/SA/d2014050' | |
|
|
23 | pathFigure = '/home/soporte/workspace/Graficos/SA/new1/' | |
|
|
24 | xmin = '15.5' | |
|
|
25 | xmax = '24' | |
|
|
26 | startTime = '15:30:00' | |
|
|
27 | filehdf5 = "SA_2014050.hdf5" | |
|
|
28 | 28 | |
|
|
29 | 29 | #2014051 20 Feb 2014 |
|
|
30 | path = '/home/soporte/Data/MST/SA/d2014051' | |
|
|
31 |
pathFigure = '/home/soporte/workspace/Graficos/SA/ |
|
|
|
32 | xmin = '0.0' | |
|
|
33 | xmax = '8.0' | |
|
|
34 |
startTime = '0 |
|
|
|
35 | filehdf5 = "SA_2014051.hdf5" | |
|
|
30 | # path = '/home/soporte/Data/MST/SA/d2014051' | |
|
|
31 | # pathFigure = '/home/soporte/workspace/Graficos/SA/new/' | |
|
|
32 | # xmin = '0.0' | |
|
|
33 | # xmax = '8.0' | |
|
|
34 | # startTime = '00:00:00' | |
|
|
35 | # filehdf5 = "SA_2014051.hdf5" | |
|
|
36 | 36 | |
|
|
37 | 37 | readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader', |
|
|
38 | 38 | path=path, |
| @@ -63,7 +63,7 opObj11.addParameter(name='maxIndex', value='60', format='float') | |||
|
|
63 | 63 | #--------------------------------------------------------------------------------------------------- |
|
|
64 | 64 | procUnitConfObj1 = controllerObj.addProcUnit(datatype='CorrelationProc', inputId=procUnitConfObj0.getId()) |
|
|
65 | 65 | # procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(1,0),(2,3)', format='pairsList') |
|
|
66 | procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(0,3),(0,2),(1,3),(1,2)', format='pairsList') | |
|
|
66 | procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList') | |
|
|
67 | 67 | procUnitConfObj1.addParameter(name='fullT', value='1', format='bool') |
|
|
68 | 68 | procUnitConfObj1.addParameter(name='removeDC', value='1', format='bool') |
|
|
69 | 69 | #procUnitConfObj1.addParameter(name='lagT', value='0,1,2,3', format='intlist') |
| @@ -98,7 +98,7 opObj21.addParameter(name='technique', value='SA', format='str') | |||
|
|
98 | 98 | opObj21.addParameter(name='positionX', value='36,0,36,0', format='floatlist') |
|
|
99 | 99 | opObj21.addParameter(name='positionY', value='36,0,0,36', format='floatlist') |
|
|
100 | 100 | opObj21.addParameter(name='azimuth', value='51.06', format='float') |
|
|
101 | opObj21.addParameter(name='crosspairsList', value='(0,3),(0,2),(1,3),(1,2)', format='pairsList')#COrregir | |
|
|
101 | opObj21.addParameter(name='crosspairsList', value='(0,3),(0,2),(1,3),(1,2),(0,1),(2,3)', format='pairsList')#COrregir | |
|
|
102 | 102 | # |
|
|
103 | 103 | opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other') |
|
|
104 | 104 | opObj22.addParameter(name='id', value='4', format='int') |
General Comments 0
You need to be logged in to leave comments.
Login now