Imaging M99

↵ CARMA Tutorials

This part of the tutorial shows an example of creating an image and a zeroth moment map of NGC4254 (M99), using the example track taken by CARMA STING.

As before, this example assumes that you have made the following global definitions in CASA:

# In CASA
msdir='./'
project='c0104I'
msfile=msdir+project+'.ms'

If your visibility data are already edited and calibrated, there's really only one step left to get you to an image: clean. Unlike in MIRIAD, which has separate tasks for inverting and cleaning an image, CASA does it in one step. If you want a dirty image, just set niter=0 so that clean doesn't actually do any...cleaning.

Here's a command-line call to clean, which demonstrates many of the interesting parameters you can tweak. This is a simple example, assuming you are making an image of the entire mosaic, and you want to clean the whole thing.

Some things of note (input parameters are italicized):

• phasecenter doesn't need to be specified, if you are happy with the phase center at the position of the first of the observed source fields. Otherwise, you can specify an RA and dec, e.g., phasecenter='J2000 12h12m12.00 12d12m12.00', or a field number, as shown in the example below.
• it is important for both imagermode and ftmachine to be set to 'mosaic' if you want your mosaicked data to be stitched together correctly (this should happen by default).
• the parameters pbcor and minpb set whether or not a primary beam correction is done, and how far out, respectively.
• always include units when specifying the threshold at which to stop cleaning.

This command produces a naturally-weighted image 512x512 arcsec in angular size, which has been cleaned down to a threshold of 54 mJy and has been primary-beam corrected (see 'minpb' and 'pbcor'):

# In CASA
clean(vis=msfile,imagename='clean_pb',spw='1~2',field='NGC4254',
      phasecenter='10',
      cell=1.,imsize=[512,512],
      mode='velocity',start='-103.00km/s',width='5.08km/s',nchan=29,
      interpolation='nearest',psfmode='clark',
      imagermode='mosaic',ftmachine='mosaic',
      scaletype='SAULT',restfreq='115.2712GHz',interactive=F,
      minpb=0.1,pbcor=True,
      niter=1000,threshold='54mJy',mosweight=False,
      weighting='natural')

To take a look at the image you produced, use the CASA viewer:

# In CASA
viewer(infile='clean_pb.image')

Desktop view of viewer display, with accompanying Data Display options GUI.

The lower part of the Viewer Display Panel provides buttons which allow you to navigate between spectral channels, and the coordinate and flux information of the pixel under the cursor is displayed at the bottom of this window. The Data Display Options window allows you to adjust the data range for the color scale as well as a number of plotting options (e.g., whether information about the third dimension of the cube is displayed in the upper right corner.)

↵ CARMA Tutorials
↵ CASAguides