Protoplanetary Disk Simulation (CASA 3.4)

From CASA Guides
Revision as of 12:37, 4 June 2012 by Rindebet (talk | contribs)
Jump to navigationJump to search

Simulating Observations in CASA

To create a script of the Python code on this page see Extracting scripts from these tutorials.

Protoplanetary disk

  • simobserve and simanalyze version for CASA 3.4

Explanation of the script

Set sim_observe as current task and reset all parameters
# Setting everything in simobserve to original defaults
Image coordinate system can be verified
# This reports image header parameters in the Log Messages window
Image center can be identified
# These are tools found in the CASA toolkit
# They are very useful, but the interface is not as straightforward as the tasks
# You can find the tool reference manual here:
# The following command is used to open an image (this is part of the image analysis toolkit)
# Whenever data are being viewed/manipulated by tools, what is being operated on needs to be explicitly opened
# and closed (i.e. an image, a table, etc.)"input50pc_672GHz.fits")
#  Out[9]: True
# Reports the length of each axis in the opened image
#  Out[11]: [257L, 257L, 1L, 1L]
# This command converts from pixel (our source file) to world coordinates (something usable by simdata)
#  Out[12]:
#{'numeric': array([  4.71239120e+00,  -4.01423802e-01,   1.00000000e+00,
#         6.72000001e+11])}
# Formats the coordinate just converted into hms
#  Out[13]: '18:00:00.03052'
# Formats one of the other coordinates into dms
#  Out[14]: '-'
# Final housekeeping by closing the image tool
# The image tool will now be detached from the image
Brightness scale can be viewed with 'imstat' task
# Default parameters are adequate for this
# ...
#  'max': array([  6.52469971e-05]),
# ...
# that's 0.0652 mJy/pixel.
Let's call our project psim2
# This defines the root prefix for any output files from simobserve
project            =  "psim2"
We'll leave the sky model the way it is: simobserve will create psim2.skymodel CASA image since this model is a fits file, and most but not all of CASA routines can operate directly on fits
skymodel           =  "input50pc_672GHz.fits"
We need to decide where to point the telescope. The image is 2/3 arcsec in size, so we only need one pointing. We could put that in a text file ourself, or let simobserve create the ascii pointing file for us.
setpointings       =  True
direction          =  "J2000 18h00m00.031s -22d59m59.6s"
mapsize            =  "0.76arcsec"
The default pointingspacing is fine: we'll only fit one pointing in the small mapsize the default calculation maptype hexagonal is ok too since only one will fit anyway.
We do want to calculate visibilities in a measurement set: let's do a 20 min snapshot observation using the "out20" ALMA antenna configuration:
obsmode            =  "int"
totaltime          =  "1200s"
Use appropriate antenna configurations based on desired angular resolution (configuration 20 - alma.out20.cfg in this case - is the largest "compact" configuration)
antennalist        =  "alma.out20.cfg"
Deconvolve the visibilities back into an image
default ("simanalyze")
project            =  "psim2"
image              =  True
# Prior image to use in clean
modelimage         =  "input50pc_672GHz.fits"
vis                =  project+""
imsize             =  [192, 192]
Specify number of iteration of cleaning task with proper threshold and weighting
niter              =  10000
threshold          =  "1e-7Jy"
weighting          =  "natural"
We'd like to calculate a difference and fidelity image, and see some diagnostics:
analyze            =  True
And see the array but not the UV coverage:
showuv             =  False
showresidual       =  True  
showconvolved      =  True
Plot both to the screen and the png files with lots of messages:
graphics           =  "both"
verbose            =  True
overwrite = True

Run simdata

# This commands CASA to execute simanalyze

  • Output results:


Simulating Observations in CASA

Last checked on CASA Version 3.4.0.