5.1.2.0 sd.scantable.sinusoid_baseline

NRAO Home > CASA > CASA Toolkit Reference Manual

sd.scantable.sinusoid_baseline - Function

5.1.2 Perform a baseline subtraction using sinusoidal function
Description

Return a scan which has been baselined (all rows) with sinusoidal functions. Fit will be done with ’sigma-clipping’.

Spectral lines are detected first using linefinder and masked out to avoid them affecting the baseline solution.

If applyfft is set to True, the function performs Fourier analysis to select wave numbers for sinusoidal fitting. Currently, ’fft’ is only available to be used for the analysis. You can specify threshold for selection of wave number using fftthresh parameter. Both float and string is acceptable. Given a float value, the unit is set to sigma. For string values, allowed formats include:

any decimal number plus ’sigma’ (e.g. ’3sigma’)
’top’ plus any decimal number (e.g. ’top10’)

In addition, you can add or reject specific wave numbers from the fit using addwn and rejwn, respectively. You can specify wave numbers as an integer, string, or list of them. For string specification, syntax for those parameters are as follows:

’a-b’ (= a, a+1, a+2, ..., b-1, b)
’<a’ (= 0, 1, ..., a-2, a-1)
’>a’ (= a+1, a+2, ... up to maximum wave number corresponding to the Nyquist frequency)

You can append ’=’ after inequality sign. When both addwn and rejwn are set, rejwn will take priority of addwn.

Note: The best-fit parameter values output in logger and/or blfile are now based on specunit of ’channel’.

Arguments


Inputs
insitu	If False a new scantable is returned
	allowed:	bool
	Default:	None (use default value)
mask	An optional mask retreived from scantable
	allowed:	bool array
	Default:	None (no mask)
applyfft	Perform Fourier analysis to find appropriate sinusoidal component
	allowed:	bool
	Default:	True
fftmethod	Method to find sinusoidal component (currently only ’fft’ is available)
	allowed:	string
	Default:	’fft’
fftthresh	Threshod to select wave number in Fourier analysis
	allowed:	float, string
	Default:	3.0
addwn	Additional wave numbers to be used for fitting
	allowed:	integer, string, any array
	Default:	[]
rejwn	Waver numbers not to be used for fitting
	allowed:	integer, string, any array
	Default:	[]
clipthresh	Clipping threshold in unit of sigma
	allowed:	float
	Default:	3.0
clipniter	Maximum number of iteration of clipping
	allowed:	integer
	Default:	0
plot	Plot the fit and the residual (currently unavailable)
	allowed:	bool
	Default:	False
getresidual	If False, return best-fit value instead of residual
	allowed:	bool
	Default:	True
showprogress	Show progress status for large data
	allowed:	bool
	Default:	True
minnrow	Minimum number of spectra to show progress status
	allowed:	integer
	Default:	1000
outlog	Output the coefficients of the best-fit function to logger
	allowed:	bool
	Default:	False
blfile	Name of text file in which the best-fit parameter values to be written
	allowed:	string
	Default:	”

Returns
scantable

Example

        scan = sd.scantable(’OrionS_rawACSmod_cal’,average=False)
        # return a scan baselined by a combination of sinusoidal curves having
        # wave numbers in spectral window up to 10,
        # also with 3-sigma clipping, iteration up to 4 times
        bscan = scan.sinusoid_baseline(addwn=’<=10’,clipthresh=3.0,clipniter=4)

[next] [prev] [prev-tail] [front] [up]

More information about CASA may be found at the CASA web page

This code is available under the terms of the GNU General Public Lincense