Simulating Observations in CASA v3.0.1
As of release 3.0.0, CASA's tools to produce simulated interferometric data are quite powerful. As with all aspects of CASA, there is a beginner-level python task interface called simdata, which calls methods in the sm C++ tool. The task interface turns a model of the sky (2 to 4 dimensions including frequency and Stokes) into the visibilities that would be measured with ALMA, (E)VLA, CARMA, SMA, ATCA, PdB, etc. The task also can produce a cleaned image of the model visibilities, compare that image with your input convolved with the synthesized beam, and calculate a fidelity image. simdata can add thermal noise (from receiver, atmosphere, and ground) to the visibilities.
The sm tool has methods that can be used to add phase delay variations, gain fluctuations and drift, cross-polarization, and (coming soon) bandpass and pointing errors to your simulated data. sm also has more flexibility in adding thermal noise than simdata, for example for new observatories that are unknown to simdata.
CASA simulation uses the aatm atmospheric model, a thin wrapper of Juan Pardo's ATM library, to accurately calculate all atmospheric corruption terms (noise, phase delay) accurately as a function of frequency and site characteristics.
Part of CASA's simulation routines are generic ephemeris and geodesy calculations available in python - see simutil.py.
simdata is under active development, documentation may lag reality, please email rindebet at nrao.edu with any questions - It's my job to help you use this software
Steps to simulation
1. Install CASA
the links below describe the various sections of inputs
2. Input Model - Preparing a patch of sky for simdata to pseudoobserve.
- Sample Model Images
- Simulation Recipes - a nearby galaxy with ALMA, a protoplanetary disk with SMA, etc
3. Antenna List - how to specify the positions and diameters of your antennas or stations.
4. Specifying Observation - how to set up what/when/how you want to observe, and the output image details
5. Corrupting Observation - (Optional) For added realism, corrupt your visibilities with thermal and phase noise.
6. Deconvolve Image (Optional) Go back from the calculated visibilities to a synthesis image