Difference between revisions of "Inspecting Data"

From CASA Guides
Jump to: navigation, search
(Checking the online flags)
Line 13: Line 13:
 
It's useful to start by checking which data have already been deleted (flagged) by the "online" flags (including time periods when the antennas were pointing off-source, or there were focus or subreflector errors), as well as data that were flagged because they were shadowed (blocked by a nearby antenna) or zero-valued (indicating a correlator error).
 
It's useful to start by checking which data have already been deleted (flagged) by the "online" flags (including time periods when the antennas were pointing off-source, or there were focus or subreflector errors), as well as data that were flagged because they were shadowed (blocked by a nearby antenna) or zero-valued (indicating a correlator error).
  
To do this, use the <tt>tflagcmd</tt> task to produce a plot:
+
To do this, use the {{tflagcmd}} task to produce a plot:
  
 
<source lang='python'>
 
<source lang='python'>
Line 27: Line 27:
  
 
== Checking the weather information ==
 
== Checking the weather information ==
 +
 +
[[Image:plotweather.png|200px|thumb|right|Weather summary plot]]
 +
 +
Another good thing to check is what the weather conditions were like during the observation.  The task that does this is {{plotweather}}:
 +
 +
<source lang='python'>
 +
# In CASA
 +
opacities = plotweather(vis='TVER0002.sb2568947.eb2579996.55518.22356400463.ms')
 +
</source>
 +
 +
The plot is automatically named vis + ".plotweather.png," so in this case, TVER0002.sb2568947.eb2579996.55518.22356400463.ms.plotweather.png.  Use your favorite graphics tool to open the file and inspect the information. 
 +
 +
The top panel shows the solar elevation -- in this case, the observation was done at night (the black disk represents the sun, below the red line, which is the horizon).  The second panel shows the wind speed and direction -- here, the wind was pretty calm (between 0.6 and 1.7 m/s).  The temperature and dewpoint are shown in the next panel.  This observation was performed in mid-November, and it was quite cold!
 +
 +
The last two panels plot the precipitable water vapor (PWV) and opacity (Tau_z) for three quantities: the value determined from the "seasonal model," the JVLA weather station, and a weighted average.  See [http://www.aoc.nrao.edu/evla/geninfo/memoseries/evlamemo143.pdf EVLA Memo 143] for more details.
 +
 +
In addition, you may notice that we set a variable <tt>opacities</tt> to capture output from {{plotweather}}.  This task outputs a list of opacities, one for each spectral window, to use for calibration.
  
 
== Initial inspection: listobs ==
 
== Initial inspection: listobs ==

Revision as of 17:06, 20 April 2012


Overview

Now that you've obtained your measurement set (MS), you will probably want to look at it carefully to see what it contains. This CASA Guide will describe a set of CASA tasks designed to help you evaluate and understand the data in an MS.

Following the example in Loading Data, we will be using the MS TVER0002.sb2568947.eb2579996.55518.22356400463.ms in this Guide. Please see Loading Data for instructions on acquiring this dataset, or use your own MS for practice.

Checking the online flags

Online, shadow, and zero flags

It's useful to start by checking which data have already been deleted (flagged) by the "online" flags (including time periods when the antennas were pointing off-source, or there were focus or subreflector errors), as well as data that were flagged because they were shadowed (blocked by a nearby antenna) or zero-valued (indicating a correlator error).

To do this, use the Template:Tflagcmd task to produce a plot:

# In CASA
tflagcmd(vis='TVER0002.sb2568947.eb2579996.55518.22356400463.ms',action='plot')

From this, we can see that antenna ea08 had some subreflector issues, and lengths of time when an antenna was not "on source" (i.e., it was slewing from one object to the next) varied a bit according to antenna. This accounts for the fact that the integration times provided by listobs are not exactly 10 s: the first integration of a given scan begins when the first antenna arrives on-source; since this does not happen simultaneously for all antennas, and the given integration time is an average across antennas, we can get values like 7.93 s (as for Scan 3).

Although the shadow and clip flags are plotted as well, the exact times affected by these flags are not known, so do not be alarmed by the fact that they appear to span the entire plot -- this is rarely actually the case.

Unfortunately, since we requested time-averaged data from the archive, the flagged data are not included in the MS. If you are concerned about the possible erroneous deletion of good data, it's best to download the complete MS, inspect it with plotms (see below), and then perform time and / or frequency averaging.

Checking the weather information

Weather summary plot

Another good thing to check is what the weather conditions were like during the observation. The task that does this is plotweather:

# In CASA
opacities = plotweather(vis='TVER0002.sb2568947.eb2579996.55518.22356400463.ms')

The plot is automatically named vis + ".plotweather.png," so in this case, TVER0002.sb2568947.eb2579996.55518.22356400463.ms.plotweather.png. Use your favorite graphics tool to open the file and inspect the information.

The top panel shows the solar elevation -- in this case, the observation was done at night (the black disk represents the sun, below the red line, which is the horizon). The second panel shows the wind speed and direction -- here, the wind was pretty calm (between 0.6 and 1.7 m/s). The temperature and dewpoint are shown in the next panel. This observation was performed in mid-November, and it was quite cold!

The last two panels plot the precipitable water vapor (PWV) and opacity (Tau_z) for three quantities: the value determined from the "seasonal model," the JVLA weather station, and a weighted average. See EVLA Memo 143 for more details.

In addition, you may notice that we set a variable opacities to capture output from plotweather. This task outputs a list of opacities, one for each spectral window, to use for calibration.

Initial inspection: listobs

Graphical display: plotms

Graphical display: viewer