Imaging a Mosaicked Spectral Line Dataset: Difference between revisions
No edit summary |
No edit summary |
||
Line 34: | Line 34: | ||
niter=1000,threshold='340mJy',mosweight=False, | niter=1000,threshold='340mJy',mosweight=False, | ||
weighting='natural') | weighting='natural') | ||
</source> | |||
To take a look at the image you produced, use the CASA '''[[viewer]]''': | |||
<source lang="python"> | |||
# In CASA | |||
viewer(image='clean_pb.image') | |||
</source> | </source> | ||
[[CARMA Tutorials | ↵ '''CARMA Tutorials''']] <br> | [[CARMA Tutorials | ↵ '''CARMA Tutorials''']] <br> | ||
[[Main Page | ↵ '''CASAguides''']] | [[Main Page | ↵ '''CASAguides''']] |
Revision as of 21:16, 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')
To take a look at the image you produced, use the CASA viewer:
# In CASA
viewer(image='clean_pb.image')