schain Files · schainpy/model/graphics/jroplot_voltage.py

Probado 16 Set 2015 Ok...

Probado 16 Set 2015 Ok Modificar path de datos y graficos y de ser necesario el modo online o offline y el walk si la busqueda es por folder o por archivos

Miguel Valdez - - Load All Authors

File last commit:

r573:e93e8802031b


                r618:2f07f1815be1

Download file

             jroplot_voltage.py
        
                    180 lines
            
             | 5.7 KiB
            
                | text/x-python
            
             |
                PythonLexer
            
             / schainpy / model / graphics / jroplot_voltage.py
          
                    History
                
                 |
                  Source
                 | Raw
                 |Copy content
                 |Copy permalink

        Daniel Valdez
    
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

              r487
            
      '''

        Miguel Valdez
    
Merge with branch schain_julia_drifts from rev. 803 to 995....

              r568
            
      Created on Jul 9, 2014

        Daniel Valdez
    
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

              r487
            
        Miguel Valdez
    
Merge with branch schain_julia_drifts from rev. 803 to 995....

              r568
            
      @author: roj-idl71

      '''

      import os

        Daniel Valdez
    
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

              r487
            
      import datetime

      import numpy

      from figure import Figure

      class Scope(Figure):

          isConfig = None

          def __init__(self):

              self.isConfig = False

              self.WIDTH = 300

              self.HEIGHT = 200

              self.counter_imagwr = 0

          def getSubplots(self):

              nrow = self.nplots

              ncol = 3

              return nrow, ncol

          def setup(self, id, nplots, wintitle, show):

              self.nplots = nplots

              self.createFigure(id=id, 

                                wintitle=wintitle, 

                                show=show)

              nrow,ncol = self.getSubplots()

              colspan = 3

              rowspan = 1

              for i in range(nplots):

                  self.addAxes(nrow, ncol, i, 0, colspan, rowspan)

          def plot_iq(self, x, y, id, channelIndexList, thisDatetime, wintitle, show, xmin, xmax, ymin, ymax):

              yreal = y[channelIndexList,:].real

              yimag = y[channelIndexList,:].imag

              title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))

              xlabel = "Range (Km)"

              ylabel = "Intensity - IQ"

              if not self.isConfig:

                  nplots = len(channelIndexList)

                  self.setup(id=id,

                             nplots=nplots,

                             wintitle='',

                             show=show)

                  if xmin == None: xmin = numpy.nanmin(x)

                  if xmax == None: xmax = numpy.nanmax(x)

                  if ymin == None: ymin = min(numpy.nanmin(yreal),numpy.nanmin(yimag))

                  if ymax == None: ymax = max(numpy.nanmax(yreal),numpy.nanmax(yimag))

                  self.isConfig = True

              self.setWinTitle(title)

              for i in range(len(self.axesList)):

                  title = "Channel %d" %(i)

                  axes = self.axesList[i]

                  axes.pline(x, yreal[i,:],

                              xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,

                              xlabel=xlabel, ylabel=ylabel, title=title)

                  axes.addpline(x, yimag[i,:], idline=1, color="red", linestyle="solid", lw=2)

          def plot_power(self, x, y, id, channelIndexList, thisDatetime, wintitle, show, xmin, xmax, ymin, ymax):

              y = y[channelIndexList,:] * numpy.conjugate(y[channelIndexList,:])

              yreal = y.real

              title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))

              xlabel = "Range (Km)"

              ylabel = "Intensity"

              if not self.isConfig:

                  nplots = len(channelIndexList)

                  self.setup(id=id,

                             nplots=nplots,

                             wintitle='',

                             show=show)

                  if xmin == None: xmin = numpy.nanmin(x)

                  if xmax == None: xmax = numpy.nanmax(x)

                  if ymin == None: ymin = numpy.nanmin(yreal)

                  if ymax == None: ymax = numpy.nanmax(yreal)

                  self.isConfig = True

              self.setWinTitle(title)

              for i in range(len(self.axesList)):

                  title = "Channel %d" %(i)

                  axes = self.axesList[i]

                  ychannel = yreal[i,:]

                  axes.pline(x, ychannel,

                              xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,

                              xlabel=xlabel, ylabel=ylabel, title=title)

          def run(self, dataOut, id, wintitle="", channelList=None,

                  xmin=None, xmax=None, ymin=None, ymax=None, save=False,

                  figpath='./', figfile=None, show=True, wr_period=1,

        Miguel Valdez
    
Merge with branch schain_julia_drifts from rev. 803 to 995....

              r568
            
                  ftp=False, server=None, folder=None, username=None, password=None, type='power'):

        Daniel Valdez
    
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

              r487
            
              """

              Input:

                  dataOut         :

                  id        :

                  wintitle        :

                  channelList     :

                  xmin            :    None,

                  xmax            :    None,

                  ymin            :    None,

                  ymax            :    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))

              x = dataOut.heightList

              y = dataOut.data[channelIndexList,:] * numpy.conjugate(dataOut.data[channelIndexList,:])

              y = y.real

        Miguel Valdez
    
Merge with branch schain_julia_drifts from rev. 803 to 995....

              r568
            
              thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[0])

        Daniel Valdez
    
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

              r487
            
              if type == "power":

                  self.plot_power(dataOut.heightList, 

                                  dataOut.data,

                                  id, 

                                  channelIndexList, 

                                  thisDatetime,

                                  wintitle,

                                  show,

                                  xmin,

                                  xmax,

                                  ymin,

                                  ymax)

              if type == "iq":

                  self.plot_iq(dataOut.heightList, 

                               dataOut.data,

                               id, 

                               channelIndexList, 

                               thisDatetime,

                               wintitle,

                               show,

                               xmin,

                               xmax,

                               ymin,

                               ymax)

              self.draw()

        Miguel Valdez
    
A new SendToServer Unit has been created to upload files to a remote server....

              r573
            
              self.save(figpath=figpath,

                        figfile=figfile,

                        save=save,

                        ftp=ftp,

                        wr_period=wr_period,

                        thisDatetime=thisDatetime)

	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

Daniel Valdez This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014	r487	'''
Miguel Valdez Merge with branch schain_julia_drifts from rev. 803 to 995....	r568	Created on Jul 9, 2014
Daniel Valdez This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014	r487
Miguel Valdez Merge with branch schain_julia_drifts from rev. 803 to 995....	r568	@author: roj-idl71
		'''
		import os
Daniel Valdez This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014	r487	import datetime
		import numpy

		from figure import Figure

		class Scope(Figure):

		isConfig = None

		def __init__(self):

		self.isConfig = False
		self.WIDTH = 300
		self.HEIGHT = 200
		self.counter_imagwr = 0

		def getSubplots(self):

		nrow = self.nplots
		ncol = 3
		return nrow, ncol

		def setup(self, id, nplots, wintitle, show):

		self.nplots = nplots

		self.createFigure(id=id,
		wintitle=wintitle,
		show=show)

		nrow,ncol = self.getSubplots()
		colspan = 3
		rowspan = 1

		for i in range(nplots):
		self.addAxes(nrow, ncol, i, 0, colspan, rowspan)

		def plot_iq(self, x, y, id, channelIndexList, thisDatetime, wintitle, show, xmin, xmax, ymin, ymax):
		yreal = y[channelIndexList,:].real
		yimag = y[channelIndexList,:].imag

		title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))
		xlabel = "Range (Km)"
		ylabel = "Intensity - IQ"

		if not self.isConfig:
		nplots = len(channelIndexList)

		self.setup(id=id,
		nplots=nplots,
		wintitle='',
		show=show)

		if xmin == None: xmin = numpy.nanmin(x)
		if xmax == None: xmax = numpy.nanmax(x)
		if ymin == None: ymin = min(numpy.nanmin(yreal),numpy.nanmin(yimag))
		if ymax == None: ymax = max(numpy.nanmax(yreal),numpy.nanmax(yimag))

		self.isConfig = True

		self.setWinTitle(title)

		for i in range(len(self.axesList)):
		title = "Channel %d" %(i)
		axes = self.axesList[i]

		axes.pline(x, yreal[i,:],
		xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,
		xlabel=xlabel, ylabel=ylabel, title=title)

		axes.addpline(x, yimag[i,:], idline=1, color="red", linestyle="solid", lw=2)

		def plot_power(self, x, y, id, channelIndexList, thisDatetime, wintitle, show, xmin, xmax, ymin, ymax):
		y = y[channelIndexList,:] * numpy.conjugate(y[channelIndexList,:])
		yreal = y.real

		title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))
		xlabel = "Range (Km)"
		ylabel = "Intensity"

		if not self.isConfig:
		nplots = len(channelIndexList)

		self.setup(id=id,
		nplots=nplots,
		wintitle='',
		show=show)

		if xmin == None: xmin = numpy.nanmin(x)
		if xmax == None: xmax = numpy.nanmax(x)
		if ymin == None: ymin = numpy.nanmin(yreal)
		if ymax == None: ymax = numpy.nanmax(yreal)

		self.isConfig = True

		self.setWinTitle(title)

		for i in range(len(self.axesList)):
		title = "Channel %d" %(i)
		axes = self.axesList[i]
		ychannel = yreal[i,:]
		axes.pline(x, ychannel,
		xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax,
		xlabel=xlabel, ylabel=ylabel, title=title)


		def run(self, dataOut, id, wintitle="", channelList=None,
		xmin=None, xmax=None, ymin=None, ymax=None, save=False,
		figpath='./', figfile=None, show=True, wr_period=1,
Miguel Valdez Merge with branch schain_julia_drifts from rev. 803 to 995....	r568	ftp=False, server=None, folder=None, username=None, password=None, type='power'):
Daniel Valdez This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014	r487
		"""

		Input:
		dataOut :
		id :
		wintitle :
		channelList :
		xmin : None,
		xmax : None,
		ymin : None,
		ymax : 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))

		x = dataOut.heightList
		y = dataOut.data[channelIndexList,:] * numpy.conjugate(dataOut.data[channelIndexList,:])
		y = y.real

Miguel Valdez Merge with branch schain_julia_drifts from rev. 803 to 995....	r568	thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.getTimeRange()[0])
Daniel Valdez This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014	r487
		if type == "power":
		self.plot_power(dataOut.heightList,
		dataOut.data,
		id,
		channelIndexList,
		thisDatetime,
		wintitle,
		show,
		xmin,
		xmax,
		ymin,
		ymax)

		if type == "iq":
		self.plot_iq(dataOut.heightList,
		dataOut.data,
		id,
		channelIndexList,
		thisDatetime,
		wintitle,
		show,
		xmin,
		xmax,
		ymin,
		ymax)


		self.draw()

Miguel Valdez A new SendToServer Unit has been created to upload files to a remote server....	r573	self.save(figpath=figpath,
		figfile=figfile,
		save=save,
		ftp=ftp,
		wr_period=wr_period,
		thisDatetime=thisDatetime)