include ’imager.g’  
#  
# First make the MS from a FITS file:  
#  
m:=fitstoms(msfile=’3C273XC1.MS’, fitsfile=’3C273XC1.FITS’); m.close();  
#  
# Now make an imager tool for the MS  
#  
imgr:=imager(’3C273XC1.MS’)  
#  
# Set the imager to produce images of cellsize 0.7 and  
# 256 by 256 pixels  
#  
imgr.setimage(nx=256,ny=256, cellx=’0.7arcsec’,celly=’0.7arcsec’);  
#  
# Wait for results before proceeding to the next step  
#  
dowait:=T  
#  
# Make and display a clean image  
#  
imgr.clean(niter=1000, threshold=’30mJy’,  
model=’3C273XC1.clean.model’, image=’3C273XC1.clean.image’)  
dd.image(’3C273XC1.clean.image’)  
#  
# Fourier transform the model  
#  
imgr.ft(model=’3C273XC1.clean.model’)  
#  
# Plot the visibilities  
#  
imgr.plotvis()  
#  
# Write out the final MS and close the imager tool  
#  
imgr.close()