Difference between revisions of "Imaging a Mosaicked Spectral Line Dataset"

From CASA Guides
Jump to navigationJump to search
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.
 +
 +
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:
 +
* phasecenter doesn't need to be specified, if you are happy with the phase center that clean derives from the combined mosaicked 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.
 +
* always include units when specifying the threshold at which to stop cleaning.
  
 
<source lang="python">
 
<source lang="python">
 
# In CASA
 
# In CASA
clean(vis=msfile,imagename='clean3pb_07dec',spw='1,2',field='3~51',
+
clean(vis=msfile,imagename='clean_pb',spw='1,2',field='3~51',
 
       phasecenter='J2000 20h35m35.32 60d10m01.00',
 
       phasecenter='J2000 20h35m35.32 60d10m01.00',
 
       cell=2.,imsize=[400,400],
 
       cell=2.,imsize=[400,400],

Revision as of 17:11, 28 December 2009

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:

  • phasecenter doesn't need to be specified, if you are happy with the phase center that clean derives from the combined mosaicked 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.
  • 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')

CARMA Tutorials
CASAguides