import os
import string
from taskinit import *
import pylab as pl
mycp=cptool()
mytp=tptool()
from parameter_check import *
def plotcal(caltable=None,xaxis=None,yaxis=None,
poln=None,field=None,antenna=None,spw=None,timerange=None,
subplot=None,overplot=None,clearpanel=None,iteration=None,
plotrange=None,showflags=None,
plotsymbol=None,plotcolor=None,markersize=None,fontsize=None,
showgui=None,figfile=None):
"""Plot calibration solutions:
An all-purpose plotter for calibration results. The values for all
calibration solutions (G, T, GSPLINE, B, BPOLY, D) can be displayed
for a variety of polarization combinations and calibrations. The
solutions may be iterated through antennas/spw/fields during one execution.
The plotter permits zooming, listing and flagging of solutions, although
the implications of flagged solutions are not yet made.
Keyword arguments:
caltable -- Name of input calibration table
default: none; example: caltable='ngc5921.gcal'
The type of calibration table is determined automatically.
xaxis -- Value to plot on the x axis
Options: 'time','chan','freq','antenna','amp','phase','real','imag','snr'
Default: cal type dependent, usually 'time'
yaxis -- Value to plot ong the y-axis
Options: 'amp','phase','real','imag','snr','antenna'
Default: cal type dependent, usually 'amp'
poln -- Polarization (or combination) to plot
default: '' (RL); all polarizations
Options: '' = ('RL'),'R','L','XY','X','Y',
'/' --> form complex poln ratio
(amp ratio and phase difference)
field -- Select field using field id(s) or field name(s).
['go listobs' to obtain the fieldt id's or names]
default: ''=all fields
If field string is a non-negative integer, it is assumed a
field index, otherwise it is assumed a field name
field='0~2'; field ids 0,1,2
field='0,4,5~7'; field ids 0,4,5,6,7
field='3C286,3C295'; field named 3C286 and 3C295
field = '3,4C*'; field id 3, all names starting with 4C
antenna -- Antenna selection (baseline syntax ignored)
default: '' (all);
example: antenna='1,3~5' means antenna
indices 1,3,4,5.
spw -- Select spectral window (channel syntax ignored, except for D)
default: ''=all spectral windows
spw='0~2,4'; spectral windows 0,1,2,4
spw='<2'; spectral windows less than 2
timerange -- Time selection
default: '' (all)
example: timerange='1995/04/13/09:15:00~1995/04/13/09:25:00'
--- Plot Options ---
subplot -- Panel number on the display screen
default: 111 (full screen display);
examples:
if iteration = 'antenna'; subplot=321 then
a plot frame will contain the first 6 antennas, in three
rows and two columns. Follow instructions on screen to
cycle through the frames
if iteration = ''; then one frame can be filled with many
plots in a piecemeal fashion; for example
antenna='0'; subplot=221; plotcal()
antenna='1'; subplot=222; plotcal()
antenna='2'; subplot=223; plotcal()
antenna='3'; subplot=224; plotcal()
overplot -- Overplot these values on current plot (if possible)
default: False;
True (overplotting) can only be done if iteration=''
clearpanel -- Clear nothing on the plot window, automatically
clear plotting area, clear the current plot area, or
clear the whole plot panel.
options: None, Auto, Current, All (None and Auto not supported)
default: Auto
example: clearpanel='Current'
iteration -- Create a sequence of plots, iterating over antenna, time,
field, and/or spw
default: '' --> create in all in one plot
example: iteration='antenna' --> create a sequence of
separate plots separated by antenna. Flagging cannot
be done in iteration mode.
plotrange -- Control the x and y ranges of the plot, as a list of
values, e.g., [xmin,xmax,ymin,ymax]
default=[] --> plot will self-scale
Note: time plotting ranges are cumbersome to use.
showflags -- If true, only flagged solutions will be plotted
default: false --> only show unflagged solutions
plotsymbol -- pylab plot symbol. See cookbook for details
default: '.': large points
',' = small points (see markersize)
'-' = connect points by line
colors are cycled automatically for multi-function plots
plotcolor -- Initial color to use on each plot
default: 'blue'
markersize -- Control the size of plot symbols
default: 5.0 --> a nice size for symbols
fontsize -- Control the font size of title (axes labels will be
80% of this size)
default: 10.0
showgui -- Whether or not to display the plotting GUI
default: True; example showgui=False
figfile -- File name to save the plotted figure to.
default: ''; example figfile=myPlot.png
"""
#Python script
#parameter_printvalues(arg_names,arg_values,arg_types)
xaxis=string.upper(xaxis)
yaxis=string.upper(yaxis)
poln=string.upper(poln)
try:
casalog.post("The plotcal task is deprecated in CASA version 5.8 and not present in CASA 6.2. All plotcal features have been made available in plotms in casa 5.8/6.2.",'WARN')
# Turn the display of the plot GUI on or off
mytp.setgui( showgui );
if ((type(caltable)==str) & (os.path.exists(caltable))):
mycp.open(caltable)
else:
raise Exception, 'Calibration table not found - please verify the name'
# If the user wants the plot window cleared then clear it.
#if ( clearpanel.lower()=='current' ) :
# mycp.clearplot( subplot=subplot );
#elif ( clearpanel.lower()=='all' ) :
# mycp.clearplot(); #RI
##if clearpanel.lower()>0: #DP
##mycp.clearplot(0) #DP
#pl.ion()
#mytp.clearplot(0,0,0)
mycp.selectcal(antenna=antenna,field=field,spw=spw,time=timerange,poln=poln)
mycp.plotoptions(subplot=subplot,
overplot=overplot,
iteration=iteration,
plotrange=plotrange,
showflags=showflags,
plotsymbol=plotsymbol,
plotcolor=plotcolor,
markersize=markersize,
fontsize=fontsize)
mycp.plot(xaxis,yaxis)
if (showgui):
pl.gcf().show()
if ( len(figfile) > 0 ) :
mycp.savefig( figfile );
except Exception, instance:
print '*** Error ***',instance
mycp.close()