Imaging a Mosaicked Spectral Line Dataset: Difference between revisions

From CASA Guides
Jump to navigationJump to search
No edit summary
No edit summary
Line 13: Line 13:
</source>
</source>


If your data are already edited and calibrated, there's really only one step left: '''[[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.
If your data are already edited and calibrated, there's really only one step left: '''[[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 a single field, and you want to clean the whole thing.
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 a single field, and you want to clean the whole thing.

Revision as of 10:35, 5 January 2010

CASAguides
CARMA Tutorials
Calibrating a Mosaicked Spectral Line Dataset


As before, these notes assume that you have made the following global definitions in CASA:

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

If your data are already edited and calibrated, there's really only one step left: 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 a single field, 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.
  • 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.
# In CASA
clean(vis=msfile,imagename='clean_pb',spw='1,2',field='3~51',
      phasecenter='J2000 20h35m35.32 60d10m01.00',
      cell=2.,imsize=[400,400],
      mode='velocity',start='-103.00km/s',width='5.08km/s',nchan=29,
      interpolation='linear',psfmode='clark',
      imagermode='mosaic',ftmachine='mosaic',
      scaletype='SAULT',restfreq='115.2712GHz',interactive=F,
      minpb=0.1,pbcor=True,
      niter=1000,threshold='340mJy',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