casa  $Rev:20696$
 All Classes Namespaces Files Functions Variables
Functions
sdimaging Namespace Reference

Functions

def sdimaging

Function Documentation

def sdimaging.sdimaging (   infile = '',
  specunit = '',
  restfreq = '',
  scanlist = [],
  field = -1,
  spw = 0,
  antenna = -1,
  stokes = '',
  gridfunction = 'BOX',
  convsupport = -1,
  truncate = -1,
  gwidth = -1,
  jwidth = -1,
  outfile = '',
  overwrite = False,
  imsize = [256,
  cell = ['1.0arcmin',
  arcmin,
  dochannelmap = False,
  nchan = 1,
  start = 0,
  step = 1,
  phasecenter = '',
  ephemsrcname = '',
  pointingcolumn = 'direction' 
)
SD task: imaging for total power and spectral data
  Keyword arguments:
infile -- name of input SD (MS) dataset
specunit -- units for spectral axis
        options: (str) 'channel','km/s','GHz','MHz','kHz','Hz'
        default: 'channel'
        example: this will be the units for nchan, start, and step
restfreq -- rest frequency
        default: '' (refer input data)
        example: 1.0e11, '100GHz'
scanlist -- list of scan numbers to process
        default: [] (use all scans)
        example: [21,22,23,24]
        this selection is in addition to field and spw
field -- field id or selection string for selecting scans by name
        default: -1 (all fields)
        example: 'FLS3a', 0
        this selection is in addition to scanlist and spw
spw -- spectral window id
        default: 0
        example: 1
        this selection is in addition to scanlist and field
antenna -- select data based on antenna name(s) or id(s)
        default: -1 (all baselines, i.e. all antenna in case of auto data)
        example: 0, 'DV01'
stokes -- select data based on stokes or polarization type 
        default: '' (use all polarizations)
        example: 'XX'
gridfunction -- gridding function for imaging
        options: 'BOX' (Box-car), 'SF' (Spheroidal), 
                 'PB' (Primary-beam), 'GAUSS' (Gaussian),
                 'GJINC' (Gaussian*Jinc)
        default: 'BOX'
        example: 'SF'
    >>> gridfunction expandable parameter:
       convsupport -- convolution support for 'SF' 
           default: -1 (use default for each gridfunction)
           example: 3
       truncate -- truncattion radius of convolution kernel.
                   effective only for 'GAUSS' and 'GJINC'.
           default: '-1' (use default for each gridfunction)
           example: 3, '20arcsec', '3pixel'
       gwidth -- HWHM for gaussian. Effective only for 
                 'GAUSS' and 'GJINC'.
           default: '-1' (use default for each gridfunction)
           example: 3, '20arcsec', '3pixel'
       jwidth -- Width of jinc function. Effective only for 
                 'GJINC'.
           default: '-1' (use default for each gridfunction)
           example: 3, '20arcsec', '3pixel'
outfile -- output image name
        default: none
        example: 'mySDimage.im'
overwrite -- overwrite option for outfile
        default: False (not overwrite)
        options: True, False
        example: if True, existing file will be overwritten
imsize -- x and y image size in pixels, symmetric for single value
        default: [256,256]
        example: imsize=200 (equivalent to [200,200])
cell -- x and y cell size. default unit arcmin
        default: ['1.0arcmin', '1.0arcmin']
        example: cell=['0.2arcmin, 0.2arcmin']
                 cell='0.2arcmin' (equivalent to example above)
dochannelmap -- channel map image or total power image
        default: False (total power)
        options: True (channel map), False
    >>> dochannelmap=True expandable parameters
       nchan -- number of spectral channel for created image
           default: 1 
           example: 100
       start -- reference value of start channel (in units of specunit)
           default: 0 (0th channel if specunit='channel')
           example: 100
       step -- width of each spectral channel for created image
           default: 1 (width of 1 channel if specunit='channel')
           example: 100
phasecenter -- image phase center: direction measure or fieldid 
        default: 0
        example: 'J2000 13h44m00 -17d02m00', 'AZEL -123d48m29 15d41m41'
ephemsrcname -- ephemeris source name for moving source
        default: ''
        if the source name in the data matches one of the known 
        solar objects by the system, this task automatically set 
        the source name. 
        example: 'moon' 
pointingcolumn -- pointing data column to use
        option: 'direction', 'target', 'pointing_offset', 'source_offset', 'encoder' 
        default: 'direction'


DESCRIPTION:

Task sdimaging create image from input single-dish data.
The input can be either total power and spectral data. Currently,
this task directly accesses the Measurement Set data only because of 
the data access efficiency. So it differs from other single-dish 
tasks that mostly operate on the ASAP scantable data format.
 
Units of spectral axis can be specified via a parameter specunit.
Allowed values for specunit are 'channel', 'GHz', 'MHz', 'kHz', 'Hz', 
and 'km/s'. This parameter is also used as the units of the parameter 
start and step that specify reference value of start channel and 
width of each spectral channel for channel map, respectively. 
The parameter nchan specifies number of channels for created image. 
If you set nchan as -1, the task selects existing all channels and 
combine those data into one channel to create continuum image.

You can specify field id or name directly. By default, field is set 
to -1 that means the task selects all fields in the data.

Selection of the antennas can be made by setting antennaid(s) or 
antenna name(s) in string (e.g. '0', 'DV01',etc.) or integer 
(e.g. 0). Defalut value, -1, means that the task selects data from 
all baseline, i.e., data from all antenna when data only contains 
auto-correlation.

The parameter gridfunction sets gridding function for imaging. 
Currently, the task supports 'BOX' (Box-car), 'SF' (Prolate 
Spheroidal Wave Function), 'GAUSS' (Gaussian), 'GJINC' (Gaussian*
Jinc), where Jinc(x) = J_1(pi*x/c)/(pi*x/c) with a first order 
Bessel function J_1, and 'PB' (Primary Beam). For 'PB', correct 
antenna informations should be included in input file. 

There are four subparameters for gridfunction: convsupport, truncate, 
gwidth, and jwidth. The convsupport is an integer specifying cut-off 
radius for 'SF' in units of pixel. By default (convsupport=-1), 
the cut-off radius is set to 3 pixels. The truncate is a cut-off 
radius for 'GAUSS' or 'GJINC'. It accepts integer, float, and 
string values of numeric plus unit. Allowed units are angular 
units such as 'deg', 'arcmin', 'arcsec', and 'pixel'. Default unit 
is 'pixel' so that string without unit or numerical values (integer 
or float) will be interpreted as radius in pixel. Default value 
for truncate, which is used when negative radius is set, is 3*HWHM 
for 'GAUSS' and radius at first null for 'GJINC'. The gwidth is 
the HWHM of gaussian for 'GAUSS' and 'GJINC'. Default value is 
sqrt(log(2)) pixel for 'GAUSS' and 2.52*sqrt(log(2)) pixel for 
'GJINC'. The jwidth specifies width of the jinc function (parameter 
'c' in the definition above). Default is 1.55 pixel. Both gwidth 
jwidth allows integer, float, or string of numeric plus unit. 
Default values for gwidth and jwidth are taken from Mangum et al. 
(2007).

Definition at line 13 of file sdimaging.py.

References publish_summary.quantity, task_sdimaging.sdimaging(), and vla_uvfits_line_sf.verify.