Plotmosaic

From CASA Guides
Jump to navigationJump to search

Return to Analysis Utilities

This page documents the plotmosaic function of Python module analysisUtils.

This function produces a plot of the pointings in a MS with primary beam FWHM circles and field names overlaid. The size of the circle plotted is 1.13*lambda/dish_diameter, which is correct for ALMA antennas. If no source ID is specified, it will use the specified intent, which defaults to OBSERVE_TARGET#ON_SOURCE. Alternatively, you can include the sourceid to limit the fields to the science target. The default coordinate system is 'relative'. You can also use 'absolute', in which case the RA and Dec are shown in decimal degrees.

Usage:

  • au.plotmosaic(vis, sourceid=2, figfile='mosaic.png', coord='relative', skipsource=-1, doplot=True, help=False, sciencespws=False, vm=", debug=False, intent='OBSERVE_TARGET#ON_SOURCE')*

Required parameters:

  • vis: the name of the measurement set

Optional parameters:

  • sourceid: may be either an integer, integer string, or the string name of the source (but it cannot be a list).
  • figfile: can be either a string or "True". In the latter case, the automatic name will be: [vis].pointings.[coord].png
  • coord: 'relative' or 'absolute'
  • skipsource: a source ID to skip
  • doplot: True or False. if False, then the central field ID is returned as an integer. Otherwise, a list is returned which is useful for determining the required image size for clean: [centralField, maxRA, minRA, minDec, maxDec] where the angles are in units of arcsec relative to the center.
  • sciencespws: if True, then run msmd (casa >= 4.1) or ValueMapping (casa < 4.1) to determine science spws.
  • vm: a ValueMapping structure (only relevant for casa < 4.1)
  • debug: if True, print extra messages
  • intent: the intent to use to select Field IDs (only used if sourceid is not specified)

Examples

# In CASA
CASA <2>: au.plotmosaic("30dor.cal.ms",sourceid='30 Doradus',doplot=False)
Field IDs with matching source ID =  [ 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23]
Got longitude = -67.755 deg
Mean frequency = 225.180484 GHz
Found 24 pointings
Field 12 is closest to the center of the area covered (3.3 arcsec away).
  Out[18]: 12
# In CASA
CASA <3>: au.plotmosaic("30dor.cal.ms",sourceid='30 Doradus',coord='abs',figfile='/tmp/absolute.png')
Field IDs with matching source ID =  [ 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23]
Got longitude = -67.755 deg
Mean frequency = 225.180484 GHz
Found 24 pointings
Field 12 is closest to the center of the area covered (3.3 arcsec away).
Wrote file = /tmp/absolute.png

Absolute.png

# In CASA
CASA <4>: au.plotmosaic("30dor.cal.ms",sourceid='30 Doradus',coord='rel',figfile='/tmp/relative.png')
Field IDs with matching source ID =  [ 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23]
Got longitude = -67.755 deg
Mean frequency = 225.180484 GHz
Found 24 pointings
Field 12 is closest to the center of the area covered (3.3 arcsec away).
Wrote file = /tmp/relative.png
  Out[4]:
[12,
 58.898348082265692,
 -56.505271415428538,
 59.801877597487355,
 -98.429477166421918]

Relative.png