Smoothbandpass

From CASA Guides
Revision as of 16:19, 3 May 2012 by Thunter (talk | contribs)
Jump to navigationJump to search

Return to Analysis Utilities

Until casa is able to smooth the uv data prior to solving for bandpass, it may be useful to smooth the solution instead, although this is not as good as the former method. After doing this smoothing, you can use plotbandpass to overlay the smoothed solution over the original solution. Note that if the signal to noise of your original solution is too low, this method will introduce unphysical spikes in the output solution! In this case, you may have to use BPOLY option of bandpass.

au.smoothbandpass(help=True)
Usage: smoothbandpass(caltable, window_len=20, window='flat', method='ri', avoidflags=True, verbose=False, help=False, log='')
The output table will have '_smoothXXwindow_method' appended to the name,
where 'XX'=window_len and 'window' is the window type.  The window 
type options are: 'flat', 'hanning', 'hamming', 'bartlett', 'blackman'.
Method types are 'ri' = real & imaginary or 'ap' = amp & phase
avoidflags: if True (default), replace flagged channels with the value from nearest unflagged channel prior to
            passing the data to the smoothing function. The input caltable is unchanged.
verbose: if True, print each channel that is being replaced, and by which channel
log: if not blank, then write out the verbose information to the specified log file name

For more info, see the Wikipedia page on window functions.

Examples

CASA <3>: au.smoothbandpass('uid___A002_X3c7a84_X443.ms.split.bandpass_spw')
Avoided 201 flagged channels on pol=0, spw=0, antenna=0
Avoided 201 flagged channels on pol=1, spw=0, antenna=0

Smooth solution written to: uid___A002_X3c7a84_X443.ms.split.bandpass_spw_smooth20flat_ri
</code>