ispec

Return to directory of Todd's tools

Emulates AIPS' ISPEC function to produce an ASCII spectrum from a CASA image, using the region tool (rg) and the image tool (ia). Operates on a single pixel and a single polarization only. See also ispecOverlay.

Usage:

tt.ispec(image, outfile=None, region='', sep=' ', format='%e', 
        overwrite=True, blc=None, trc=None, pol=0, showplot=False, 
        source=None, velocity=False, startchan=None, 
        endchan=None, plotfile=None, plotrange=[0,0,0,0])

• outfile: the name of the ASCII file to produce
• source: append this source name to the filename
• region: the region to use (alternative to blc,trc)
• sep: the character string to use to separate columns in the output file
• blc, trc: the pixel range to use as a list (e.g. [40,50]). Currently blc must equal trc. If trc is not specified, it is assumed to equal blc. If blc not specified, then find and use the peak.
• velocity: if True, then convert Hz to km/s via restfreq keyword
• startchan: the first spectral channel to consider
• endchan: the final spectral channel to consider (if negative, then ignore this many channels at end)

Returns: three lists:

• xaxis in GHz (or km/s), flux density, blc (i.e. the pixel used)

Example:

CASA <5>: tt.ispec('/net/shadowfax/export/data_2/vext_inorth/casa_imaging/IN_VEX228.LINE_hann.image.pbcor',blc=[209,365],outfile='sma1b.228.txt',showplot=T)
  Out[5]:
(array([  2.28892724e+11,   2.28893585e+11,   2.28894446e+11, ...,
         2.30868610e+11,   2.30869471e+11,   2.30870331e+11]),
 array([ 0.06410486,  0.12293984,  0.1259677 , ...,  0.        ,
        0.        ,  0.        ]))

-- ToddHunter - 2013-12-06