Imaging 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 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.
- 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')
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.)