overlayCuts
Return to directory of Todd's CASA extensions
Extract a single row or column from two CASA images and overlay them.
Usage:
au.overlayCuts(img1, img2, img3=None, row1=None, row2=None, row3=None,
column1=None, column2=None, column3=None
stokes1='XX', plotfile='', plotrange=[0, 0, 0, 0],
frequency1=None, frequency2=None, frequency3=None,
panels=1, interpolate=False, gaussian=None,
scaleByFrequency=True, plotrange2=[0,0,0,0],
scaleToArcsec=False, diameter=None, truncate=True,
obscuration=0.75)
- row1, column1: which row or column to use in the first image
- row2,column2: which row or column to use in the second image
- stokes1: which Stokes image to use for the first image
- plotfile: if set to a name, then a png file will be produced with this name
- scaleByFrequency: if True, then the xaxis of the first image (and 3rd) will be scaled to match the second by their frequency ratio.
- frequency1: the frequency of the first image (if not in header), a string with units, or floating point GHz
- frequency2: the frequency of the second image (if not in header), a string with units, or floating point GHz
- frequency3: the frequency of the second image (if not in header), a string with units, or floating point GHz
- frequencyRatio: the ratio of the frequency of img1 to img2
- panels: 1 or 2, setting this to 2 will show a second subplot with the overlay in dB units
- interpolate: this parameter is passed to extractCutFromImage
- gaussian: if set to a floating point value (in arcsec), it will also overlay a Gaussian beam (in green)
- scaleToArcsec: pass to extractCutFromImage
- diameter: the antenna diameter to use for the Airy profile to overlay (default = None = no overlay)
- truncate: if True, truncate the Airy profile at the first null
- obscuration: the central obstruction diameter for computing Airy (default = 0.75m)
Example:
--
ToddHunter - 2013-12-11