sdstat – list statistics of spectral – single dish task

Description

Task sdstat computes basic statistics for each of single-dish spectrum. This task returns a Python dictionary of statistics. The return value contains the maximum and minimum intensity and their channels (‘max’, ‘max_abscissa’, ‘min’, and ‘min_abscissa’), RMS (‘rms’), mean (‘mean’), sum (‘sum’), median (‘median’), standard deviation (‘stddev’), total intensity (‘totint’), and equivalent width (‘eqw’). If you do have multiple scantable rows, then the return values will be lists.

It is possible to select channel regions to calculate spectra either non-interactively by spw parameter or interactively on a plotter by setting interactive=True.

If one of averaging parameters is set True, the spectra are averaged before calculating the statistics.

Parameters

Title

Parameter	Default	Description
infile	''
antenna	int(0)
fluxunit	''	units of the flux ['K' or 'Jy'] (''=current)
telescopeparam	''
field	''
spw	''
restfreq	''
frame	''	frequency reference frame ['LSRK', 'TOPO', 'LSRD', 'BARY', 'GALACTO', 'LGROUP', or 'CMB'] (''=current)
doppler	''	doppler convention ['RADIO', 'OPTICAL', 'Z', 'BETA', or 'GAMMA'] (''=current)
timerange	''
scan	''
pol	''
beam	''
timeaverage	False
tweight	'tintsys'	weighting for time averaging ['tintsys', 'tsys', 'tint', 'var', or 'median']
scanaverage	False
polaverage	False
pweight	'tsys'	weighting for polarization averaging ['tsys' or 'var']
interactive	False
outfile	''
format	'3.3f'
overwrite	False	overwrite the output file if already exists [True, False]
xstat	[ ]

Parameter Explanations

infile

''

name of input SD dataset

antenna

int(0)

select an antenna name or ID, e.g. 'PM03' (only effective for MS input)

fluxunit

''

units of the flux (''=current)

telescopeparam

''

parameters of telescope for flux conversion (see examples in help)

field

''

select data by field IDs and names, e.g. '3C2*' (''=all)

spw

''

select data by IF IDs (spectral windows), e.g. '3,5,7' (''=all)

restfreq

''

the rest frequency, e.g. '1.41GHz' (default unit: Hz) (see examples in help)

frame

''

frequency reference frame (''=current)

doppler

''

doppler convention (''=current). Effective only when spw selection is in velocity unit.

timerange

''

select data by time range, e.g. '09:14:0~09:54:0' (''=all) (see examples in help)

scan

''

select data by scan numbers, e.g. '21~23' (''=all)

pol

''

select data by polarization IDs, e.g. '0,1' (''=all)

beam

''

select data by beam IDs, e.g. '0,1' (''=all)

timeaverage

False

average spectra over time [True, False] (see examples in help)

tweight

'tintsys'

weighting for time averaging

scanaverage

False

average spectra within a scan number [True, False] (see examples in help)

polaverage

False

average spectra over polarizations [True, False]

pweight

'tsys'

weighting for polarization averaging

interactive

False

determines interactive masking [True, False]

outfile

''

name of output file (ASCII) to save statistics

format

'3.3f'

format string to print statistic values in file, e.g, '.7e'

overwrite

False

overwrite the output file if already exists

xstat

[ ]

RETURN ONLY: a Python dictionary of line statistics