##// END OF EJS Templates
schain
RSS
test
test
jespinoza - - Load All Authors

File last commit:

r1343:95b6dc6e7058
r1376:3ce340933772
Show More
jroplot_voltage.py
302 lines | 9.6 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
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 from schainpy.model.graphics.jroplot_base import Plot, plt
George Yong
Multiprocessing for voltage (all operations) working
 r1173
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285
class ScopePlot(Plot):
'''
Plot for Scope
avaldez
update for weather radar options
 r1296 '''
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285
CODE = 'scope'
plot_type = 'scatter'
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 def setup(self):
self.xaxis = 'Range (Km)'
self.ncols = 1
self.nrows = 1
self.nplots = 1
self.ylabel = 'Intensity [dB]'
self.titles = ['Scope']
self.colorbar = False
self.width = 6
self.height = 4
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 def update(self, dataOut):
data = {}
meta = {
'nProfiles': dataOut.nProfiles,
'flagDataAsBlock': dataOut.flagDataAsBlock,
'profileIndex': dataOut.profileIndex,
}
if self.CODE == 'scope':
data[self.CODE] = dataOut.data
elif self.CODE == 'pp_power':
data[self.CODE] = dataOut.dataPP_POWER
elif self.CODE == 'pp_signal':
data[self.CODE] = dataOut.dataPP_POW
elif self.CODE == 'pp_velocity':
data[self.CODE] = dataOut.dataPP_DOP
elif self.CODE == 'pp_specwidth':
data[self.CODE] = dataOut.dataPP_WIDTH
return data, meta
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 def plot_iq(self, x, y, channelIndexList, thisDatetime, wintitle):
avaldez
update for weather radar options
 r1296
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 yreal = y[channelIndexList,:].real
yimag = y[channelIndexList,:].imag
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y"))
self.xlabel = "Range (Km)"
self.ylabel = "Intensity - IQ"
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 self.y = yreal
self.x = x
avaldez
update for weather radar options
 r1296
self.titles[0] = title
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
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 for i,ax in enumerate(self.axes):
title = "Channel %d" %(i)
if ax.firsttime:
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.xmin = min(x)
self.xmax = max(x)
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 ax.plt_r = ax.plot(x, yreal[i,:], color='b')[0]
ax.plt_i = ax.plot(x, yimag[i,:], color='r')[0]
else:
ax.plt_r.set_data(x, yreal[i,:])
ax.plt_i.set_data(x, yimag[i,:])
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 def plot_power(self, x, y, channelIndexList, thisDatetime, wintitle):
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 y = y[channelIndexList,:] * numpy.conjugate(y[channelIndexList,:])
yreal = y.real
avaldez
update for weather radar options
 r1296 yreal = 10*numpy.log10(yreal)
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 self.y = yreal
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 title = wintitle + " Power: %s" %(thisDatetime.strftime("%d-%b-%Y"))
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 self.xlabel = "Range (Km)"
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.ylabel = "Intensity [dB]"
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 self.titles[0] = title
for i,ax in enumerate(self.axes):
title = "Channel %d" %(i)
ychannel = yreal[i,:]
avaldez
update for weather radar options
 r1296
if ax.firsttime:
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.xmin = min(x)
self.xmax = max(x)
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 ax.plt_r = ax.plot(x, ychannel)[0]
else:
ax.plt_r.set_data(x, ychannel)
avaldez
update for weather radar options
 r1296
def plot_weatherpower(self, x, y, channelIndexList, thisDatetime, wintitle):
y = y[channelIndexList,:]
yreal = y.real
yreal = 10*numpy.log10(yreal)
self.y = yreal
title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))
self.xlabel = "Range (Km)"
self.ylabel = "Intensity"
self.xmin = min(x)
self.xmax = max(x)
self.titles[0] =title
for i,ax in enumerate(self.axes):
title = "Channel %d" %(i)
ychannel = yreal[i,:]
if ax.firsttime:
ax.plt_r = ax.plot(x, ychannel)[0]
else:
#pass
ax.plt_r.set_data(x, ychannel)
def plot_weathervelocity(self, x, y, channelIndexList, thisDatetime, wintitle):
x = x[channelIndexList,:]
yreal = y
self.y = yreal
title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))
self.xlabel = "Velocity (m/s)"
self.ylabel = "Range (Km)"
self.xmin = numpy.min(x)
self.xmax = numpy.max(x)
self.titles[0] =title
for i,ax in enumerate(self.axes):
title = "Channel %d" %(i)
xchannel = x[i,:]
if ax.firsttime:
ax.plt_r = ax.plot(xchannel, yreal)[0]
else:
#pass
ax.plt_r.set_data(xchannel, yreal)
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
avaldez
graficos de ploteo pulsepair
 r1305 def plot_weatherspecwidth(self, x, y, channelIndexList, thisDatetime, wintitle):
x = x[channelIndexList,:]
yreal = y
self.y = yreal
title = wintitle + " Scope: %s" %(thisDatetime.strftime("%d-%b-%Y %H:%M:%S"))
self.xlabel = "width "
self.ylabel = "Range (Km)"
self.xmin = numpy.min(x)
self.xmax = numpy.max(x)
self.titles[0] =title
for i,ax in enumerate(self.axes):
title = "Channel %d" %(i)
xchannel = x[i,:]
if ax.firsttime:
ax.plt_r = ax.plot(xchannel, yreal)[0]
else:
#pass
ax.plt_r.set_data(xchannel, yreal)
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 def plot(self):
if self.channels:
channels = self.channels
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 else:
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 channels = self.data.channels
thisDatetime = datetime.datetime.utcfromtimestamp(self.data.times[-1])
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343
scope = self.data[-1][self.CODE]
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 if self.data.flagDataAsBlock:
avaldez
update for weather radar options
 r1296
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 for i in range(self.data.nProfiles):
wintitle1 = " [Profile = %d] " %i
avaldez
update for weather radar options
 r1296 if self.CODE =="scope":
if self.type == "power":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_power(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope[:,i,:],
channels,
thisDatetime,
wintitle1
)
if self.type == "iq":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_iq(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope[:,i,:],
channels,
thisDatetime,
wintitle1
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 )
avaldez
update for weather radar options
 r1296 if self.CODE=="pp_power":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_weatherpower(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope[:,i,:],
channels,
thisDatetime,
wintitle
)
avaldez
UPDATE del jroplot_voltage.py
 r1312 if self.CODE=="pp_signal":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_weatherpower(self.data.yrange,
avaldez
UPDATE del jroplot_voltage.py
 r1312 scope[:,i,:],
channels,
thisDatetime,
wintitle
)
avaldez
update for weather radar options
 r1296 if self.CODE=="pp_velocity":
self.plot_weathervelocity(scope[:,i,:],
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.data.yrange,
avaldez
update for weather radar options
 r1296 channels,
thisDatetime,
wintitle
)
avaldez
graficos de ploteo pulsepair
 r1305 if self.CODE=="pp_spcwidth":
self.plot_weatherspecwidth(scope[:,i,:],
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.data.yrange,
avaldez
graficos de ploteo pulsepair
 r1305 channels,
thisDatetime,
wintitle
)
avaldez
update for weather radar options
 r1296 else:
wintitle = " [Profile = %d] " %self.data.profileIndex
if self.CODE== "scope":
if self.type == "power":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_power(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope,
channels,
thisDatetime,
wintitle
)
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285
if self.type == "iq":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_iq(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope,
channels,
thisDatetime,
wintitle
Juan C. Espinoza
Update and fix plot modules #TODO: correlation & spectraheis
 r1285 )
avaldez
update for weather radar options
 r1296 if self.CODE=="pp_power":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_weatherpower(self.data.yrange,
avaldez
update for weather radar options
 r1296 scope,
channels,
thisDatetime,
wintitle
)
avaldez
UPDATE del jroplot_voltage.py
 r1312 if self.CODE=="pp_signal":
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.plot_weatherpower(self.data.yrange,
avaldez
UPDATE del jroplot_voltage.py
 r1312 scope,
channels,
thisDatetime,
wintitle
)
avaldez
update for weather radar options
 r1296 if self.CODE=="pp_velocity":
self.plot_weathervelocity(scope,
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.data.yrange,
avaldez
update for weather radar options
 r1296 channels,
thisDatetime,
wintitle
)
avaldez
graficos de ploteo pulsepair
 r1305 if self.CODE=="pp_specwidth":
self.plot_weatherspecwidth(scope,
jespinoza
Add update method to plots to pass data (no more changes in jrodata)
 r1343 self.data.yrange,
avaldez
graficos de ploteo pulsepair
 r1305 channels,
thisDatetime,
wintitle
)
avaldez
update for weather radar options
 r1296
class PulsepairPowerPlot(ScopePlot):
'''
avaldez
UPDATE del jroplot_voltage.py
 r1312 Plot for P= S+N
avaldez
update for weather radar options
 r1296 '''
CODE = 'pp_power'
plot_type = 'scatter'
class PulsepairVelocityPlot(ScopePlot):
'''
avaldez
UPDATE del jroplot_voltage.py
 r1312 Plot for VELOCITY
avaldez
update for weather radar options
 r1296 '''
CODE = 'pp_velocity'
plot_type = 'scatter'
avaldez
graficos de ploteo pulsepair
 r1305
class PulsepairSpecwidthPlot(ScopePlot):
'''
avaldez
UPDATE del jroplot_voltage.py
 r1312 Plot for WIDTH
avaldez
graficos de ploteo pulsepair
 r1305 '''
CODE = 'pp_specwidth'
plot_type = 'scatter'
avaldez
UPDATE del jroplot_voltage.py
 r1312
class PulsepairSignalPlot(ScopePlot):
'''
Plot for S
'''
CODE = 'pp_signal'
plot_type = 'scatter'